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.
- Publication:
-
arXiv Mathematics e-prints
- Pub Date:
- November 2003
- DOI:
- 10.48550/arXiv.math/0311515
- arXiv:
- arXiv:math/0311515
- Bibcode:
- 2003math.....11515H
- Keywords:
-
- Numerical Analysis;
- 65Rxx;
- 78-xx
- E-Print:
- 33 pages, 4 figures