Three-dimensional quasi-periodic shifted Green function throughout the spectrum, including Wood anomalies.

This work, part II in a series, presents an efficient method for evaluation of wave scattering by doubly periodic diffraction gratings at or near what are commonly called 'Wood anomaly frequencies'. At these frequencies, there is a grazing Rayleigh wave, and the quasi-periodic Green function ceases to exist. We present a modification of the Green function by adding two types of terms to its lattice sum. The first type are transversely shifted Green functions with coefficients that annihilate the growth in the original lattice sum and yield algebraic convergence. The second type are quasi-periodic plane wave solutions of the Helmholtz equation which reinstate certain necessary grazing modes without leading to blow-up at Wood anomalies. Using the new quasi-periodic Green function, we establish, for the first time, that the Dirichlet problem of scattering by a smooth doubly periodic scattering surface at a Wood frequency is uniquely solvable. We also present an efficient high-order numerical method based on this new Green function for scattering by doubly periodic surfaces at and around Wood frequencies. We believe this is the first solver able to handle Wood frequencies for doubly periodic scattering problems in three dimensions. We demonstrate the method by applying it to acoustic scattering.

Rapidly convergent quasi-periodic Green functions for scattering by arrays of cylinders-including Wood anomalies.

This paper presents a full-spectrum Green-function methodology (which is valid, in particular, at and around Wood-anomaly frequencies) for evaluation of scattering by periodic arrays of cylinders of arbitrary cross section-with application to wire gratings, particle arrays and reflectarrays and, indeed, general arrays of conducting or dielectric bounded obstacles under both transverse electric and transverse magnetic polarized illumination. The proposed method, which, for definiteness, is demonstrated here for arrays of perfectly conducting particles under transverse electric polarization, is based on the use of the shifted Green-function method introduced in a recent contribution (Bruno & Delourme 2014 J. Computat. Phys.262, 262-290 (doi:10.1016/j.jcp.2013.12.047)). A certain infinite term arises at Wood anomalies for the cylinder-array problems considered here that is not present in the previous rough-surface case. As shown in this paper, these infinite terms can be treated via an application of ideas related to the Woodbury-Sherman-Morrison formulae. The resulting approach, which is applicable to general arrays of obstacles even at and around Wood-anomaly frequencies, exhibits fast convergence and high accuracies. For example, a few hundreds of milliseconds suffice for the proposed approach to evaluate solutions throughout the resonance region (wavelengths comparable to the period and cylinder sizes) with full single-precision accuracy.

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.

High-order integral equation methods for problems of scattering by bumps and cavities on half-planes.

This paper presents high-order integral equation methods for the evaluation of electromagnetic wave scattering by dielectric bumps and dielectric cavities on perfectly conducting or dielectric half-planes. In detail, the algorithms introduced in this paper apply to eight classical scattering problems, namely, scattering by a dielectric bump on a perfectly conducting or a dielectric half-plane, and scattering by a filled, overfilled, or void dielectric cavity on a perfectly conducting or a dielectric half-plane. In all cases field representations based on single-layer potentials for appropriately chosen Green functions are used. The numerical far fields and near fields exhibit excellent convergence as discretizations are refined-even at and around points where singular fields and infinite currents exist.

Fourier continuation methods for high-fidelity simulation of nonlinear acoustic beams.

On the basis of recently developed Fourier continuation (FC) methods and associated efficient parallelization techniques, this text introduces numerical algorithms that, due to very low dispersive errors, can accurately and efficiently solve the types of nonlinear partial differential equation (PDE) models of nonlinear acoustics in hundred-wavelength domains as arise in the simulation of focused medical ultrasound. As demonstrated in the examples presented in this text, the FC approach can be used to produce solutions to nonlinear acoustics PDEs models with significantly reduced discretization requirements over those associated with finite-difference, finite-element and finite-volume methods, especially in cases involving waves that travel distances that are orders of magnitude longer than their respective wavelengths. In these examples, the FC methodology is shown to lead to improvements in computing times by factors of hundreds and even thousands over those required by the standard approaches. A variety of one-and two-dimensional examples presented in this text demonstrate the power and capabilities of the proposed methodology, including an example containing a number of scattering centers and nonlinear multiple-scattering events.

Efficient high-order evaluation of scattering by periodic surfaces: deep gratings, high frequencies, and glancing incidences.

We present a superalgebraically convergent integral equation algorithm for evaluation of TE and TM electromagnetic scattering by smooth perfectly conducting periodic surfaces z=f(x). For grating-diffraction problems in the resonance regime (heights and periods up to a few wavelengths) the proposed algorithm produces solutions with full double-precision accuracy in single-processor computing times of the order of a few seconds. The algorithm can also produce, in reasonable computing times, highly accurate solutions for very challenging problems, such as (a) a problem of diffraction by a grating for which the peak-to-trough distance equals 40 times its period that, in turn, equals 20 times the wavelength; and (b) a high-frequency problem with very small incidence, up to 0.01 degrees from glancing. The algorithm is based on the concurrent use of Floquet and Chebyshev expansions together with certain integration weights that are computed accurately by means of an asymptotic expansion as the number of integration points tends to infinity.

Prescribed error tolerances within fixed computational times for scattering problems of arbitrarily high frequency: the convex case.

We present a new algorithm for the numerical solution of problems of electromagnetic or acoustic scattering by large, convex obstacles. This algorithm combines the use of an ansatz for the unknown density in a boundary-integral formulation of the scattering problem with an extension of the ideas of the method of stationary phase. We include numerical results illustrating the high-order convergence of our algorithm as well as its asymptotically bounded computational cost as the frequency increases.

Inverse scattering problem for optical coherence tomography.

We deal with the imaging problem of determining the internal structure of a body from backscattered laser light and low-coherence interferometry. Specifically, using the interference fringes that result when the backscattering of low-coherence light is made to interfere with the reference beam, we obtain maps detailing the values of the refractive index within the sample. Our approach accounts fully for the statistical nature of the coherence phenomenon; the numerical experiments that we present, which show image reconstructions of high quality, were obtained under noise floors exceeding those present for various experimental setups reported in the literature.