Abstract:In order to conquer the spurious reflections from the truncated edges and maintain the stability in the long-time simulation of elastic wave propagation, several perfectly matched layer (PML) methods have been proposed in the first-order (e.g., velocity–stress equations) and the second-order (e.g., energy equation with displacement unknown only) formulations. The multiaxial perfectly matched layer (M-PML) holds the excellent stability for the long-time simulation of wave propagation, even though it is not perfectly matched in the discretized M-PML equation system. This absorbing boundary approach can offer an alternative way to solve the problem of the late-time instability, especially for anisotropic media, which is also suffered by the convolutional perfectly matched layer (C-PML) that is supposed to be competent to handle most stable problems. The M-PML termination implementation in the first-order formulations is well proposed. The common drawback of the implementation of the first-order M-PML formulations is that it necessitates fundamental reconstruction of the existing codes of the second-order spectral element method (SEM) or finite element method (FEM). Therefore, we propose a nonconvolutional second-order M-PML absorbing boundary condition approach for the wave propagation simulation in elastic media that has not yet been developed before. Two-dimensional numerical simulation validations demonstrate that the proposed second-order M-PML has good performances: 1) superior efficiency and stability of absorbing the spurious elastic wavefields, both the surface waves and body waves, reflected on the boundaries; 2) superior stability in the long-time simulation even in the isotropic medium with a high Poisson's ratio; 3) superior efficiency and stability in the long-time simulation for anisotropic media. This method hence makes the SEM and FEM in the second-order wave equation formulation more efficient and stable for the long-time simulation.

Efficient TDFEM Schemes with Second-Order PML Based on Different Temporal Basis Functions

A temporal stable smoothed particle finite element method for large deformation problems in geomechanics

Implementation of Perfectly Matched Layer (PML) for the WCS-FDTD Method Using DSP Techniques

A Modified PML Implementation in FDTD for Handling Dispersive Media

Efficiency of Perfectly Matched Layers for Seismic Wave Modeling in Second-Order Viscoelastic Equations

Time‐domain Integral Equation Solvers Using Quadratic B‐spline Temporal Basis Functions

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

A Simple Implementation of PML for Second-Order Elastic Wave Equations

A Circular-Shaped Perfectly Matched Layer Strategy for Rectangular Fdtd Grids

Efficient Implementation of CFS-CPML in FDTD Solutions of Second-Order Seismic Wave Equations

Numerical implementations of PML boundary conditions in the TLM-based SCN FDTD grid

Cubic Lattice-Based Spherical Uniaxial Perfectly Matched Layer For The Fdtd Method

Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method

Efficient Interpolation Scheme for Multilevel Fast Multipole Algorithm

Three-Dimensional Spherical-Shaped UPML for FDTD with Cubic Lattices

A Memory-efficient Implementation of Perfectly Matched Layer with Smoothly-varying Coefficients in Discontinuous Galerkin Time-Domain Method

The Second-Order Numerical Approximation for a Modified Ericksen–Leslie Model

Variable-order optimal implicit finite-difference schemes for explicit time-marching solutions to wave equations

Temporal Second-Order Fast Finite Difference/Compact Difference Schemes for Time-Fractional Generalized Burgers' Equations

Modified symplectic scheme with finite element method for seismic wavefield modeling

Compact Matrix-Exponential-Based FDTD with Second-Order PML and Direct Z-Transform for Modeling Complex Subsurface Sensing and Imaging Problems