An efficient solver for problems of scattering by bodies of revolution

Abstract

This proposal relates to the design, analysis and application of a novel numerical scheme for the solution of axisymmetric scattering problems. To this end, a procedure is introduced to iteratively evaluate the solution of the Lippmann-Schwinger integral equation in $O(N\log^2 N)$ operations, where $N$ is the number of the discretization points. The method achieves its efficiency through the use of the addition theorem and Fast Legendre Transforms (FLT). For globally smooth sound velocities/refractive indexes the method is spectrally accurate. More generally the order of convergence is tied to and in fact, limited by, the smoothness of the solution.