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.Communicated by Per Lötstedt.

Propagation of guided waves through weak penetrable scatterers

“…This latter problem has its own interest and it has been extensively studied, notably to handle for several mathematical difficulties and to develop accurate and efficient numerical methods. [11][12][13] Also, a recent increasing interest for these periodic structures is inspired by the success of studies of photonic crystals 14 and possible applications in other contexts of waves. The literature is vast, see Refs.…”

On the use of perfectly matched layers in the presence of long or backward propagating guided elastic waves

“…It therefore bears some strong similarities with the method proposed by Skelton et al in [8] to overcome the very same issue, which uses the biorthogonality relations satisfied by the Rayleigh-Lamb modes to separate the forward propagating waves from the backward ones in order to treat them appropriately within the PML. It also shares a bond with the approach proposed by Barnett and Greengard in [11] for an integral representation for quasi-periodic scattering problems, in the sense that it involves the computation of a finite number of "corrections", which measure in some way the failure of the approximate solution to satisfy a radiation condition.…”

A method to compute periodic sums

“…Even though it is dedicated to the Helmholtz equation in two dimensions, a method for "periodization" of free-space solutions similar in spirit to that presented here was proposed and tested. Moreover, the cited paper contains an additional "periodization" method based on boundary integrals, which can be tried for different kernels and space dimensionality (also see [28], where periodization along one dimension is performed).…”