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

O. P. Bruno, Leonid Kunyansky

Research output: Contribution to journalArticle

45 Citations (Scopus)

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 of London, Series A: Mathematical and Physical Sciences
Volume457
Issue number2016
DOIs
StatePublished - Dec 8 2001
Externally publishedYes

Fingerprint

Surface scattering
Three-dimension
Scattering
Higher Order
Fast Fourier transform
scattering
Fast Fourier transforms
High Order Accuracy
Three-dimensional
Acoustic Scattering
acoustic scattering
Partition of Unity
Boundary integral equations
Evaluate
personal computers
Personal Computer
Boundary Integral Equations
Personal computers
Breakdown
integral equations

Keywords

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

ASJC Scopus subject areas

  • General

Cite this

@article{a13e7687bbdc4f5babb695e0d030fbf4,
title = "Surface scattering in three dimensions: An accelerated high-order solver",
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.",
keywords = "Equivalent sources, Fast algorithm, Fast Fourier transform, Integral equation, Spherical wave expansion, Wave scattering",
author = "Bruno, {O. P.} and Leonid Kunyansky",
year = "2001",
month = "12",
day = "8",
doi = "10.1098/rspa.2001.0882",
language = "English (US)",
volume = "457",
pages = "2921--2934",
journal = "Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences",
issn = "0962-8444",
publisher = "Royal Society of London",
number = "2016",

}

TY - JOUR

T1 - Surface scattering in three dimensions

T2 - An accelerated high-order solver

AU - Bruno, O. P.

AU - Kunyansky, Leonid

PY - 2001/12/8

Y1 - 2001/12/8

N2 - 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.

AB - 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.

KW - Equivalent sources

KW - Fast algorithm

KW - Fast Fourier transform

KW - Integral equation

KW - Spherical wave expansion

KW - Wave scattering

UR - http://www.scopus.com/inward/record.url?scp=3042796907&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=3042796907&partnerID=8YFLogxK

U2 - 10.1098/rspa.2001.0882

DO - 10.1098/rspa.2001.0882

M3 - Article

AN - SCOPUS:3042796907

VL - 457

SP - 2921

EP - 2934

JO - Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

JF - Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

SN - 0962-8444

IS - 2016

ER -