Surface scattering in three dimensions: An accelerated high-order solver

O. P. Bruno, L. A. Kunyansky

Research output: Contribution to journalArticle

46 Scopus citations

Abstract

We present a new algorithm for the numerical solution of problems of acoustic scattering by surfaces in three-dimensional space. This algorithm evaluates scattered fields through fast, high-order, accurate solution of the corresponding boundary integral equation. The high-order accuracy of our solver is achieved through use of partitions of unity together with analytical resolution of kernel singularities. The acceleration, in turn, results from use of high-order equivalent source approximations, which allow for fast evaluation of non-adjacent interactions by means of the three-dimensional fast Fourier transform (FFT). Our acceleration scheme has dramatically lower memory requirements and yields much higher accuracy than existing FFT-accelerated techniques. The present algorithm computes one matrix-vector multiply in O(N6/5 log N) to O(N4/3 log N) operations (depending on the geometric characteristics of the scattering surface), it exhibits super-algebraic 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.

Original languageEnglish (US)
Pages (from-to)2921-2934
Number of pages14
JournalProceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume457
Issue number2016
DOIs
StatePublished - Dec 8 2001
Externally publishedYes

Keywords

  • Equivalent sources
  • Fast Fourier transform
  • Fast algorithm
  • Integral equation
  • Spherical wave expansion
  • Wave scattering

ASJC Scopus subject areas

  • Mathematics(all)
  • Engineering(all)
  • Physics and Astronomy(all)

Fingerprint Dive into the research topics of 'Surface scattering in three dimensions: An accelerated high-order solver'. Together they form a unique fingerprint.

  • Cite this