Perfectly matched layer for second-order time-domain elastic wave equation: formulation and stability

A time domain system of equations is proposed to model elastic wave propaga-tion in an unbounded two-dimensional anisotropic solid using perfectly matched layer (PML). Starting from a system of first-order frequency domain stress-velocity equations and using complex coordinate stretching approach with a two-parameter stretch function, a second-order formulation is obtained. The final system, which consists of just two second order equations along with four aux-iliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. The discrete stability of the so-lutions for a given mesh size is examined with the help of a plane-wave analysis of the corresponding continuous problem. It is shown that increasing the scaling parameter of the stretch function leads to significant stability improvements for certain anisotropic media that have known issues. Numerical computations for different isotropic and anisotropic media are used to illustrate the results