Surface scattering in three dimensions: an accelerated high–order solver

Oscar P. Bruno,Leonid A. Kunyansky
DOI: https://doi.org/10.1098/rspa.2001.0882
2001-12-08
Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences
Abstract:We present a new algorithm for the numerical solution of problems of acoustic scattering by surfaces in threedimensional space. This algorithm evaluates scattered fields through fast, highorder, accurate solution of the corresponding boundary integral equation. The highorder accuracy of our solver is achieved through use of partitions of unityI together with analytical resolution of kernel singularities. The acceleration, in turn, results from use of highorder equivalent source approximations, which allow for fast evaluation of nonadjacent interactions by means of the threedimensional fast Fourier transform (FFT). Our acceleration scheme has dramatically lower memory requirements and yields much higher accuracy than existing FFTaccelerated techniques. The present algorithm computes one matrixvector multiply in O(N6/5logN) to O(N4/3logN) operations (depending on the geometric characteristics of the scattering surface), it exhibits superalgebraic convergence, and it does not suffer from accuracy breakdowns of any kind. We demonstrate the efficiency of our method through a variety of examples. In particular, we show that the present algorithm can evaluate accurately, on a personal computer, scattering from bodies of acoustical sizes (ka) of several hundreds.
What problem does this paper attempt to address?