A new second-order formulation is obtained for elastic wave propagation in 2D media bounded by a perfectly matched layer (PML). The formulation uses a complex coordinate stretching approach with a two-parameter stretch function. The final system, consisting of just two second-order displacement equations along with four auxiliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. With the help of a plane-wave analysis, the stability of the continuous formulation is examined. It is shown that by increasing the scaling parameter in the stretch function, any existing instability is moved to higher spatial frequencies. Since discrete models cannot resolve frequencies beyond a certain limit, this can lead to significant computational stability improvements. Numerical results are shown to validate our formulation and to illustrate the improved stability that can be achieved with certain anisotropic media that have known issues.
BérengerJP. A perfectly matched layer for the absorption of electromagnetic waves. J Comput Phys1994; 114(2): 185–200.
2.
SacksZSKingslandDMLeeR. A perfectly matched anisotropic absorber for use as an absorbing boundary condition. IEEE Trans Antennas Propag1995; 43(12): 1460–1463.
3.
RodenJAGedneySD. Efficient implementation of the uniaxial-based PML media in three-dimensional nonorthogonal coordinates with the use of the FDTD technique. Microwave Opt Technol Lett1997; 14(2): 71–75.
4.
GedneySD. An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices. IEEE Trans Antennas Propag1996; 44(12): 1630–1639.
5.
ChewWCWeedonWH. A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates. Microwave Opt Technol Lett1994; 7(13): 599–604.
6.
KuzuogluMMittraR. Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers. IEEE Microwave Guided Wave Lett1996; 6(12): 447–449.
7.
TeixeiraFLChewWC. On causality and dynamic stability of perfectly matched layers for FDTD simulations. IEEE Trans Microwave Theory Tech1999; 47(6): 775–785.
8.
Meza-FajardoKCPapageorgiouAS. A nonconvolutional, split-field, perfectly matched layer for wave propagation in isotropic and anisotropic elastic media: stability analysis. Bull Seismol Soc Am2008; 98(4): 1811–1836.
9.
BérengerJP. Application of the CFS PML to the absorption of evanescent waves in waveguides. IEEE J Sel Areas Commun2002; 12(6): 218–220.
10.
BécacheEPetropoulosPGGedneySD. On the long-time behavior of unsplit perfectly matched layers. IEEE Trans Antennas Propag2004; 52(5): 1335–1342.
11.
BécacheEFauqueuxSJolyP. Stability of perfectly matched layers, group velocities and anisotropic waves. J Comput Phys2003; 188(2): 399–433.
12.
AppelöDKreissG. A new absorbing layer for elastic waves. J Comput Phys2006; 215(2): 642–660.
13.
DuruKKreissG. A well-posed and discretely stable perfectly matched layer for elastic wave equations in second order formulation. Commun Comput Phys2012; 11(5): 1643–1672.
14.
CollinoFTsogkaC. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics2001; 66(1): 294–307.
15.
HastingsFDSchneiderJBBroschatSL. Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation. J Acoust Soc Am1996; 100(5): 3061–3069.
16.
DrossaertFHGiannopoulosA. Complex frequency shifted convolution PML for FDTD modelling of elastic waves. Wave Motion2007; 44(7–8): 593–604.
17.
ChewWCLiuQH. Perfectly matched layers for elastodynamics: A new absorbing boundary condition. J Comput Acoust1996; 4(4): 341–359.
18.
KucukcobanSKallivokasLF. Mixed perfectly-matched-layers for direct transient analysis in 2D elastic heterogeneous media. Comput Meth Appl Mech Eng2011; 200(1–4): 57–76.
19.
AppelöDHagstromTKreissG. Perfectly matched layers for hyperbolic systems: general formulation, well-posedness, and stability. SIAM J Appl Math2006; 67(1): 1–23.
20.
RodenJAGedneySD. Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media. Microwave Opt Technol Lett2000; 27(5): 334–339.
21.
LohPROskooiAFIbanescuMSkorobogatiyMJohnsonSG. Fundamental relation between phase and group velocity, and application to the failure of perfectly matched layers in backward-wave structures. Phys Rev E2009; 79: 065601.
22.
KreissGDuruK. Discrete stability of perfectly matched layers for anisotropic wave equations in first and second order formulation. BIT Numer Math2013; 53(3): 641–663.
23.
KomatitschDMartinR. An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation. Geophysics2007; 72(5): SM155.
24.
BeltzerAI. Acoustics of solids. New York: Springer-Verlag, 1988.
25.
CollinoFMonkP. The perfectly matched layer in curvilinear coordinates. SIAM J Sci Comput1998; 19(6): 2061–2090.
26.
ChewWCJinJMMichielssenE. Complex coordinate stretching as a generalized absorbing boundary condition. Microwave Opt Technol Lett1997; 15(6): 363–369.
27.
TeixeiraFLChewWC. Complex space approach to perfectly matched layers: a review and some new developments. Int J Numer Modell Electron Networks Devices Fields2000; 13(5): 441–455.
28.
JohnsonSG. Notes on perfectly matched layers (PMLs). Technical Report, Massachusetts Institute of Technology, Cambridge, MA.
29.
ZhangWShenY. Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling. Geophysics2010; 75(4): T141–T154.
30.
DrossaertFHGiannopoulosA. A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves. Geophysics2007; 72(2): T9–T17.
31.
PetropoulosPG. Reflectionless sponge layers as absorbing boundary conditions for the numerical solution of Maxwell equations in rectangular, cylindrical, and spherical coordinates. SIAM J Appl Math2000; 60(3): 1037–1058.
32.
HarrisGMartinC. The roots of a polynomial vary continuously as a function of the coefficients. Proc Am Math Soc1987; 100(2): 390–392.
33.
ChungJHulbertGM. A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized-α method. J Appl Mech1993; 60(2): 371–375.
