A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form (Q2217869)

This article discusses a a stable discontinuous Galerkin method with a perfectly matched layer (PML) for 2D and 3D space dimensional linear elastodynamics, in velocity-stress formulation, subject to well-posed linear boundary conditions. First, the elastodynamics equation in a cuboidal domain is discusses and an unsplit PML is derived by truncating the domain using complex coordinate stretching. Second, the authors develop a discontinuous Galerkin numerical method for the linear elastodynamics equation using physically motivated numerical flux and penalty parameters, which are compatible with all well-posed, internal and external, boundary conditions. Third, to ensure numerical stability of the discretization in the presence of PML damping, the numerical DG fluxes and the numerical inter-element and boundary procedures are extended to the PML auxiliary differential equations. This step turns out to be crucial for deriving discrete energy estimates analogous to the continuous energy case. Numerical experiments included in the paper support the theoretical findings.