A High-Order, Fast Algorithm for Scattering Calculation in Two Dimensions

J. C. Aguilar, Yu Chen

Research output: Contribution to journalArticle

Abstract

We present a high-order, fast, iterative solver for the direct scattering calculation for the Helmholtz equation in two dimensions. Our algorithm solves the scattering problem formulated as the Lippmann-Schwinger integral equation for compactly supported, smoothly vanishing scatterers. There are two main components to this algorithm. First, the integral equation is discretized with quadratures based on high-order corrected trapezoidal rules for the logarithmic singularity present in the kernel of the integral equation. Second, on the uniform mesh required for the trapezoidal rule we rewrite the discretized integral operator as a composition of two linear operators: a discrete convolution followed by a diagonal multiplication; therefore, the application of these operators to an arbitrary vector, required by an iterative method for the solution of the discretized linear system, will cost N 2 log(N) for a N-by-N mesh, with the help of FFT. We will demonstrate the performance of the algorithm for scatterers of complex structures and at large wave numbers. For numerical implementations, GMRES iterations will be used, and corrected trapezoidal rules up to order 20 will be tested.

Original languageEnglish (US)
Pages (from-to)1-11
Number of pages11
JournalComputers and Mathematics with Applications
Volume47
Issue number1
DOIs
StatePublished - Jan 2004

Fingerprint

Trapezoidal Rule
Fast Algorithm
Integral equations
Integral Equations
Two Dimensions
Scattering
Higher Order
Mesh
Iteration
Iterative Solver
Helmholtz equation
GMRES
Scattering Problems
Helmholtz Equation
Iterative methods
Convolution
Complex Structure
Integral Operator
Quadrature
Fast Fourier transforms

Keywords

  • Correction coefficients
  • Integral equation
  • Large wave numbers
  • Logarithmic singularity
  • Quadrature rules

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics
  • Modeling and Simulation

Cite this

A High-Order, Fast Algorithm for Scattering Calculation in Two Dimensions. / Aguilar, J. C.; Chen, Yu.

In: Computers and Mathematics with Applications, Vol. 47, No. 1, 01.2004, p. 1-11.

Research output: Contribution to journalArticle

@article{72ade98d828346858cb27ab215858d8e,
title = "A High-Order, Fast Algorithm for Scattering Calculation in Two Dimensions",
abstract = "We present a high-order, fast, iterative solver for the direct scattering calculation for the Helmholtz equation in two dimensions. Our algorithm solves the scattering problem formulated as the Lippmann-Schwinger integral equation for compactly supported, smoothly vanishing scatterers. There are two main components to this algorithm. First, the integral equation is discretized with quadratures based on high-order corrected trapezoidal rules for the logarithmic singularity present in the kernel of the integral equation. Second, on the uniform mesh required for the trapezoidal rule we rewrite the discretized integral operator as a composition of two linear operators: a discrete convolution followed by a diagonal multiplication; therefore, the application of these operators to an arbitrary vector, required by an iterative method for the solution of the discretized linear system, will cost N 2 log(N) for a N-by-N mesh, with the help of FFT. We will demonstrate the performance of the algorithm for scatterers of complex structures and at large wave numbers. For numerical implementations, GMRES iterations will be used, and corrected trapezoidal rules up to order 20 will be tested.",
keywords = "Correction coefficients, Integral equation, Large wave numbers, Logarithmic singularity, Quadrature rules",
author = "Aguilar, {J. C.} and Yu Chen",
year = "2004",
month = "1",
doi = "10.1016/S0898-1221(04)90001-6",
language = "English (US)",
volume = "47",
pages = "1--11",
journal = "Computers and Mathematics with Applications",
issn = "0898-1221",
publisher = "Elsevier Limited",
number = "1",

}

TY - JOUR

T1 - A High-Order, Fast Algorithm for Scattering Calculation in Two Dimensions

AU - Aguilar, J. C.

AU - Chen, Yu

PY - 2004/1

Y1 - 2004/1

N2 - We present a high-order, fast, iterative solver for the direct scattering calculation for the Helmholtz equation in two dimensions. Our algorithm solves the scattering problem formulated as the Lippmann-Schwinger integral equation for compactly supported, smoothly vanishing scatterers. There are two main components to this algorithm. First, the integral equation is discretized with quadratures based on high-order corrected trapezoidal rules for the logarithmic singularity present in the kernel of the integral equation. Second, on the uniform mesh required for the trapezoidal rule we rewrite the discretized integral operator as a composition of two linear operators: a discrete convolution followed by a diagonal multiplication; therefore, the application of these operators to an arbitrary vector, required by an iterative method for the solution of the discretized linear system, will cost N 2 log(N) for a N-by-N mesh, with the help of FFT. We will demonstrate the performance of the algorithm for scatterers of complex structures and at large wave numbers. For numerical implementations, GMRES iterations will be used, and corrected trapezoidal rules up to order 20 will be tested.

AB - We present a high-order, fast, iterative solver for the direct scattering calculation for the Helmholtz equation in two dimensions. Our algorithm solves the scattering problem formulated as the Lippmann-Schwinger integral equation for compactly supported, smoothly vanishing scatterers. There are two main components to this algorithm. First, the integral equation is discretized with quadratures based on high-order corrected trapezoidal rules for the logarithmic singularity present in the kernel of the integral equation. Second, on the uniform mesh required for the trapezoidal rule we rewrite the discretized integral operator as a composition of two linear operators: a discrete convolution followed by a diagonal multiplication; therefore, the application of these operators to an arbitrary vector, required by an iterative method for the solution of the discretized linear system, will cost N 2 log(N) for a N-by-N mesh, with the help of FFT. We will demonstrate the performance of the algorithm for scatterers of complex structures and at large wave numbers. For numerical implementations, GMRES iterations will be used, and corrected trapezoidal rules up to order 20 will be tested.

KW - Correction coefficients

KW - Integral equation

KW - Large wave numbers

KW - Logarithmic singularity

KW - Quadrature rules

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

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

U2 - 10.1016/S0898-1221(04)90001-6

DO - 10.1016/S0898-1221(04)90001-6

M3 - Article

AN - SCOPUS:0942300720

VL - 47

SP - 1

EP - 11

JO - Computers and Mathematics with Applications

JF - Computers and Mathematics with Applications

SN - 0898-1221

IS - 1

ER -