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.

A Finite Volume Method for Elastic Waves in Heterogeneous, Anisotropic and Fractured Media

A finite-volume discretization for deformation of fractured media

A high order discontinuous Galerkin method for the symmetric form of the anisotropic viscoelastic wave equation

An efficient multi-field dynamic model for 3D wave propagation in saturated anisotropic porous media

Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media

A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex media using physics based numerical fluxes

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

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

Hybrid-Dimensional Finite Volume Discretizations for Fractured Porous Media

An Exact Riemann Solver for Wave Propagation in Arbitrary Anisotropic Elastic Media with Fluid Coupling

A Nearly Analytical Discrete Method For Wave-Field Simulations In 2d Porous Media

Optimally Accurate Nearly Analytic Discrete Scheme for Wave-Field Simulation in 3D Anisotropic Media

3-D Simulation of Wave Propagation in the Fractured Porous Media for Subsurface Sensing: A Rotated Staggered Finite-Difference Algorithm Based on Equivalent Medium Theory

A hybrid numerical method for elastic wave propagation in discontinuous media with complex geometry

A modified time-space optimized symplectic method for solving 3D elastic wave equations

Modeling 3-D Elastic Wave Propagation in TI Media Using Discontinuous Galerkin Method on Tetrahedral Meshes

A Simple Multi-Directional Absorbing Layer Method to Simulate Elastic Wave Propagation in Unbounded Domains

Symplectic Stereomodelling Method for Solving Elastic Wave Equations in Porous Media

A discontinuous Galerkin based multiscale method for heterogeneous elastic wave equations

A meshless method for stress-wave propagation in anisotropic and cracked media

Efficient scattering approach to seismic full-waveform inversion in anisotropic elastic media with variable density