Superalgebraically convergent smoothly windowed lattice sums for doubly periodic Green functions in three-dimensional space

This work, part I in a two-part series, presents: (i) a simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) an associated boundary-integral equation method for the numerical solution of problems of scattering of waves by doubly periodic arrays of scatterers in three-dimensional space. Except for certain ‘Wood frequencies’ at which the quasi-periodic Green function ceases to exist, the proposed approach, which is based on smooth windowing functions, gives rise to tapered lattice sums which converge superalgebraically fast to the Green function—that is, faster than any power of the number of terms used. This is in sharp contrast to the extremely slow convergence exhibited by the lattice sums in the absence of smooth windowing. (The Wood-frequency problem is treated in part II.) This paper establishes rigorously the superalgebraic convergence of the windowed lattice sums. A variety of numerical results demonstrate the practical efficiency of the proposed approach.

[1]  Gérard Tayeb,et al.  Combined Method for the Computation of the Doubly Periodic Green's Function , 2001 .

[2]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[3]  Fernando Reitich,et al.  Solution of a boundary value problem for the Helmholtz equation via variation of the boundary into the complex domain , 1992, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[4]  L. Rayleigh.  III. Note on the remarkable case of diffraction spectra described by Prof. Wood , 1907 .

[5]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[6]  O. Bruno,et al.  A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications , 2001 .

[7]  R. Wood.  XLII. On a remarkable case of uneven distribution of light in a diffraction grating spectrum , 1902 .

[8]  O. Bruno,et al.  Efficient Solution of Acoustic and Electromagnetic Scattering Problems in Three-Dimensional Periodic Media , 2011 .

[9]  Jonathan M. Borwein,et al.  Lattice Sums Then and Now , 2013 .

[10]  Randy C. Paffenroth,et al.  Electromagnetic integral equations requiring small numbers of Krylov-subspace iterations , 2009, J. Comput. Phys..

[11]  Jin Au Kong,et al.  Polarimetric Passive Remote Sensing of Periodic Surfaces , 1991 .

[12]  Vassilis G. Papanicolaou,et al.  Ewald's Method Revisited: Rapidly Convergent Series Representations of Certain Green's Functions , 1999 .

[13]  Chris M. Linton,et al.  Lattice Sums for the Helmholtz Equation , 2010, SIAM Rev..

[14]  Leslie Greengard,et al.  A new integral representation for quasi-periodic scattering problems in two dimensions , 2011 .

[15]  R. Kress,et al.  Integral equation methods in scattering theory , 1983 .

[16]  Oscar P. Bruno,et al.  Regularized integral equations and fast high‐order solvers for sound‐hard acoustic scattering problems , 2012 .

[17]  Donald R. Wilton,et al.  Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method , 2007, J. Comput. Phys..

[18]  Oscar P. Bruno,et al.  Rapidly convergent two-dimensional quasi-periodic Green function throughout the spectrum - including Wood anomalies , 2014, J. Comput. Phys..

[19]  R. Wood,et al.  On a Remarkable Case of Uneven Distribution of Light in a Diffraction Grating Spectrum , 1902 .

[20]  Avner Friedman,et al.  Maxwell’s equations in a periodic structure , 1991 .

[21]  Oscar P. Bruno,et al.  Surface scattering in three dimensions: an accelerated high–order solver , 2001, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[22]  John Anderson Monro,et al.  A Super-Algebraically Convergent, Windowing-Based Approach to the Evaluation of Scattering from Periodic Rough Surfaces , 2008 .