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.

MRFD Method for Numerical Solution of Wave Propagation in Layered Media with General Boundary Condition

SIMULATION OF ACOUSTIC WAVE PROPAGATION IN COMPLEX MEDIA USING MRFD METHOD

2D Frequency-Domain Finite-Difference Acoustic Wave Modeling Using Optimized Perfectly Matched Layers

Numerical Investigation of 3D multichannel analysis of surface wave method

Finite element modeling of wave propagation problems in multilayered soils resting on a rigid base

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

Wavelet-Galerkin Method for Computations of Electromagnetic Fields

Improved HIE-FDTD Method for Simulating Thin Layer Structure in Linear Dispersive Media

Multilevel adaptive mesh modeling for wave propagation in layer media

Efficient 2D Acoustic Wave Finite-Difference Numerical Simulation in Strongly Heterogeneous Media Using the Adaptive Mesh Refinement Technique

On Free Wave Propagation in Anisotropic Layered Media

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

A novel CFS-PML boundary condition for transient electromagnetic simulation using a fictitious wave domain method

A semianalytical solution for one‐dimensional transient wave propagation in a saturated single‐layer porous medium with a fluid surface layer

A finite-element algorithm with a perfectly matched layer boundary condition for seismic modelling in a diffusive-viscous medium

FDwave3D: a MATLAB solver for the 3D anisotropic wave equation using the finite-difference method

Fast Multipole Method For 3-D Helmholtz Equation In Layered Media

2-D wave equation modeling and migration by a new finite difference scheme based on the Galerkin method

Numerical Simulation of Elastic Wave Propagation in Azimuthally Anisotropic Media Using Multi-Level Finite Difference Method

NEW NUMERICAL METHOD FOR SURFACE WAVE PROPAGATION IN MULTILAYER ANISOTROPY MATERIALS

Two Dimensional Linear Irregular Wave Propagation by Fmbem