TY - JOUR

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

AU - Barnett, Alex

AU - Greengard, Leslie

N1 - Funding Information:
Acknowledgements We thank Stephen Shipman and Simon Chandler-Wilde for invaluable discussions. The work of AHB was supported by NSF grant DMS-0811005, and by the Class of 1962 Fellowship at Dartmouth College. The work of LG was supported by the Department of Energy under contract DEFG0288ER25053 and by AFOSR under MURI grant FA9550-06-1-0337.

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 -