Abstract:This paper studies the perfectly matched layer (PML) method for wave scattering in a half space of a homogeneous medium bounded by a two-dimensional, perfectly conducting, and locally defected periodic surface, and develops a high-accuracy boundary integral equation (BIE) solver. Along the vertical direction, we place a PML to truncate the unbounded domain onto a strip and prove that the PML solution converges to the true solution in the physical subregion of the strip with an error bounded by the reciprocal PML thickness. Laterally, we divide the unbounded strip into three regions: one region containing the defect and two semi-waveguide regions, separated by two vertical line segments. In both semi-waveguides, we prove the well-posedness of an associated scattering problem so as to well define a Neumann-to-Dirichlet (NtD) operator on the associated vertical segment. The two NtD operators, serving as exact lateral boundary conditions, reformulate the unbounded strip problem as a boundary value problem over the defected region. Due to the periodicity of the semi-waveguides, the two NtD operators are related to two different Neumann-marching operators, and either Neumann-marching operator is governed by a nonlinear Riccati equation. It is proved that the Neumann-marching operators are contracting, so that the PML solution decays exponentially fast along both lateral directions. The consequences culminate in two opposite aspects. Negatively, the PML solution cannot converge exponentially to the true solution in the whole physical region of the strip. Positively, from a numerical perspective, the Riccati equations can now be efficiently solved by a recursive doubling procedure and a high-accuracy PML-based BIE method so that the boundary value problem on the defected region can be solved efficiently and accurately. Numerical experiments demonstrate that the PML solution converges exponentially fast to the true solution in any compact subdomain of the strip.

A Bivariate Spline Solution to the Exterior Helmholtz Equation and Its Applications

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

Boundary Knot Method for Some Inverse Problems Associated with the Helmholtz Equation

A Multiaxial Perfectly Matched Layer (M-PML) for the Long-Time Simulation of Elastic Wave Propagation in the Second-Order Equations

An iterative method for Helmholtz boundary value problems arising in wave propagation

Coordinate complexification for the Helmholtz equation with Dirichlet boundary conditions in a perturbed half-space

Numerical solution to the PML problem of the biharmonic wave scattering in periodic structures

PML and High-Accuracy Boundary Integral Equation Solver for Wave Scattering by a Locally Defected Periodic Surface

SOLVING HELMHOLTZ EQUATION BY MESHLESS RADIAL BASIS FUNCTIONS METHOD

Sharp preasymptotic error bounds for the Helmholtz $h$-FEM

Convergence of overlapping domain decomposition methods with PML transmission conditions applied to nontrapping Helmholtz problems

Mixed boundary value problems for the Helmholtz equation in a model 2D angular domain

The Helmholtz Dirichlet and Neumann problems on piecewise smooth open curves

The modified method of fundamental solutions for exterior problems of the Helmholtz equation; spurious eigenvalues and their removals

Schwarz methods with PMLs for Helmholtz problems: fast convergence at high frequency

The $hp$-FEM applied to the Helmholtz equation with PML truncation does not suffer from the pollution effect

Quadratic Spline Collocation Method and Efficient Preconditioner for the Helmholtz Equation with the Sommerfeld Boundary Conditions

An immersed boundary hierarchical B-spline method for flexoelectricity

Multivariate Splines and Their Applications

WAVE SCATTERING IN LAYERED ORTHOTROPIC MEDIA I: A STABLE PML AND A HIGH-ACCURACY BOUNDARY INTEGRAL EQUATION METHOD

Method of virtual sources using on-surface radiation conditions for the Helmholtz equation