A new integral representation for quasi-periodic scattering problems in two dimensions

Alex Barnett, Leslie Greengard

Research output: Contribution to journalArticle

Abstract

Boundary integral equations are an important class of methods for acoustic and electromagnetic scattering from periodic arrays of obstacles. For piecewise homogeneous materials, they discretize the interface alone and can achieve high order accuracy in complicated geometries. They also satisfy the radiation condition for the scattered field, avoiding the need for artificial boundary conditions on a truncated computational domain. By using the quasi-periodic Green's function, appropriate boundary conditions are automatically satisfied on the boundary of the unit cell. There are two drawbacks to this approach: (i) the quasi-periodic Green's function diverges for parameter families known as Wood's anomalies, even though the scattering problem remains well-posed, and (ii) the lattice sum representation of the quasi-periodic Green's function converges in a disc, becoming unwieldy when obstacles have high aspect ratio. In this paper, we bypass both problems by means of a new integral representation that relies on the free-space Green's function alone, adding auxiliary layer potentials on the boundary of the unit cell strip while expanding the linear system to enforce quasi-periodicity. Summing nearby images directly leaves auxiliary densities that are smooth, hence easily represented in the Fourier domain using Sommerfeld integrals. Wood's anomalies are handled analytically by deformation of the Sommerfeld contour. The resulting integral equation is of the second kind and achieves spectral accuracy. Because of our image structure, inclusions which intersect the unit cell walls are handled easily and automatically. We include an implementation and simple code example with a freely-available MATLAB toolbox.

Original languageEnglish (US)
Pages (from-to)67-90
Number of pages24
JournalBIT Numerical Mathematics
Volume51
Issue number1
DOIs
StatePublished - Mar 2011

Fingerprint

Periodic Problem
Scattering Problems
Green's function
Integral Representation
Two Dimensions
Periodic Functions
Scattering
Unit
Anomaly
Wood
Boundary conditions
Layer Potentials
Artificial Boundary Conditions
Quasiperiodicity
Radiation Condition
High Order Accuracy
Acoustic Scattering
Boundary integral equations
Electromagnetic Scattering
Cell

ASJC Scopus subject areas

  • Computer Networks and Communications
  • Software
  • Applied Mathematics
  • Computational Mathematics

Cite this

A new integral representation for quasi-periodic scattering problems in two dimensions. / Barnett, Alex; Greengard, Leslie.

In: BIT Numerical Mathematics, Vol. 51, No. 1, 03.2011, p. 67-90.

Research output: Contribution to journalArticle

@article{c8f9e68d6f1d4a598e38522bfb88683d,
title = "A new integral representation for quasi-periodic scattering problems in two dimensions",
abstract = "Boundary integral equations are an important class of methods for acoustic and electromagnetic scattering from periodic arrays of obstacles. For piecewise homogeneous materials, they discretize the interface alone and can achieve high order accuracy in complicated geometries. They also satisfy the radiation condition for the scattered field, avoiding the need for artificial boundary conditions on a truncated computational domain. By using the quasi-periodic Green's function, appropriate boundary conditions are automatically satisfied on the boundary of the unit cell. There are two drawbacks to this approach: (i) the quasi-periodic Green's function diverges for parameter families known as Wood's anomalies, even though the scattering problem remains well-posed, and (ii) the lattice sum representation of the quasi-periodic Green's function converges in a disc, becoming unwieldy when obstacles have high aspect ratio. In this paper, we bypass both problems by means of a new integral representation that relies on the free-space Green's function alone, adding auxiliary layer potentials on the boundary of the unit cell strip while expanding the linear system to enforce quasi-periodicity. Summing nearby images directly leaves auxiliary densities that are smooth, hence easily represented in the Fourier domain using Sommerfeld integrals. Wood's anomalies are handled analytically by deformation of the Sommerfeld contour. The resulting integral equation is of the second kind and achieves spectral accuracy. Because of our image structure, inclusions which intersect the unit cell walls are handled easily and automatically. We include an implementation and simple code example with a freely-available MATLAB toolbox.",
author = "Alex Barnett and Leslie Greengard",
year = "2011",
month = "3",
doi = "10.1007/s10543-010-0297-x",
language = "English (US)",
volume = "51",
pages = "67--90",
journal = "BIT Numerical Mathematics",
issn = "0006-3835",
publisher = "Springer Netherlands",
number = "1",

}

TY - JOUR

T1 - A new integral representation for quasi-periodic scattering problems in two dimensions

AU - Barnett, Alex

AU - Greengard, Leslie

PY - 2011/3

Y1 - 2011/3

N2 - Boundary integral equations are an important class of methods for acoustic and electromagnetic scattering from periodic arrays of obstacles. For piecewise homogeneous materials, they discretize the interface alone and can achieve high order accuracy in complicated geometries. They also satisfy the radiation condition for the scattered field, avoiding the need for artificial boundary conditions on a truncated computational domain. By using the quasi-periodic Green's function, appropriate boundary conditions are automatically satisfied on the boundary of the unit cell. There are two drawbacks to this approach: (i) the quasi-periodic Green's function diverges for parameter families known as Wood's anomalies, even though the scattering problem remains well-posed, and (ii) the lattice sum representation of the quasi-periodic Green's function converges in a disc, becoming unwieldy when obstacles have high aspect ratio. In this paper, we bypass both problems by means of a new integral representation that relies on the free-space Green's function alone, adding auxiliary layer potentials on the boundary of the unit cell strip while expanding the linear system to enforce quasi-periodicity. Summing nearby images directly leaves auxiliary densities that are smooth, hence easily represented in the Fourier domain using Sommerfeld integrals. Wood's anomalies are handled analytically by deformation of the Sommerfeld contour. The resulting integral equation is of the second kind and achieves spectral accuracy. Because of our image structure, inclusions which intersect the unit cell walls are handled easily and automatically. We include an implementation and simple code example with a freely-available MATLAB toolbox.

AB - Boundary integral equations are an important class of methods for acoustic and electromagnetic scattering from periodic arrays of obstacles. For piecewise homogeneous materials, they discretize the interface alone and can achieve high order accuracy in complicated geometries. They also satisfy the radiation condition for the scattered field, avoiding the need for artificial boundary conditions on a truncated computational domain. By using the quasi-periodic Green's function, appropriate boundary conditions are automatically satisfied on the boundary of the unit cell. There are two drawbacks to this approach: (i) the quasi-periodic Green's function diverges for parameter families known as Wood's anomalies, even though the scattering problem remains well-posed, and (ii) the lattice sum representation of the quasi-periodic Green's function converges in a disc, becoming unwieldy when obstacles have high aspect ratio. In this paper, we bypass both problems by means of a new integral representation that relies on the free-space Green's function alone, adding auxiliary layer potentials on the boundary of the unit cell strip while expanding the linear system to enforce quasi-periodicity. Summing nearby images directly leaves auxiliary densities that are smooth, hence easily represented in the Fourier domain using Sommerfeld integrals. Wood's anomalies are handled analytically by deformation of the Sommerfeld contour. The resulting integral equation is of the second kind and achieves spectral accuracy. Because of our image structure, inclusions which intersect the unit cell walls are handled easily and automatically. We include an implementation and simple code example with a freely-available MATLAB toolbox.

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

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

U2 - 10.1007/s10543-010-0297-x

DO - 10.1007/s10543-010-0297-x

M3 - Article

AN - SCOPUS:79953054189

VL - 51

SP - 67

EP - 90

JO - BIT Numerical Mathematics

JF - BIT Numerical Mathematics

SN - 0006-3835

IS - 1

ER -