Abstract:The absorbing boundary conditions play a key role in numerical simulations. The complex frequency-shifted perfectly matched layer (CFS-PML), one of the absorbing boundary conditions, has attracted widespread attention because it can efficiently attenuate near-grazing incident waves and evanescent waves. In practical application, the distinct boundary reflections still exist using the CFS-PML in a thin PML region. To further enhance the absorbing performance, we have developed a novel unified high-order unsplit CFS-PML technique. The high-order CFS-PML combined with the discontinuous Galerkin (DG) method is applied to solve the second-order seismic wave equation on an unstructured mesh. In this work, we establish a unified framework for high-order unsplit CFS-PML formulations. The classical CFS-PML is the first-order form of our proposed unsplit high-order approach. By using auxiliary variables and auxiliary ordinary differential equations (AODEs), a higher-order CFS-PML formulation can be reduced to a first-order CFS-PML formula. For an Nth-order AODE CFS-PML, it will lead to (2N-1) additional variables and (2N-1) auxiliary differential equations. Taking second- and third-order as examples, we have derived the high-order unsplit AODE CFS-PML expressions in detail. The DG method can effectively deal with complex geological structures, such as caves, cracks, faults, etc., which yields more accurate numerical solutions. We incorporate the high-order unsplit AODE CFS-PML technique into the DG method to approximate seismic wave equations. Their strengths are combined to calculate the propagating wavefield in more complex finite models and obtain better results. In addition, the stability conditions of the second- and third-order CFS-PMLs are provided. We perform the high-order unsplit CFS-PML implementation for seismic wave modeling in several models. The isotropic and anisotropic examples prove that the high-order unsplit CFS-PML enjoys better absorbing performance than the classical CFS-PML. The numerical result for the graben model certifies the feasibility and flexibility of our proposed high-order CFS-PML combined with the DG method.

Low- and High-Order Unsplit ADE CFS-PML Boundary Conditions with Discontinuous Galerkin Method for Wavefield Simulation in Multiporosity Media.

A Unified Higher-Order Unsplit CFS-PML Technique for Solving Second-Order Seismic Equations Using Discontinuous Galerkin Method

Double-pole Unsplit Complex-Frequency-shifted Multiaxial Perfectly Matched Layer Combined with Strong-Stability-preserved Runge-Kutta Time Discretization for Seismic Wave Equation Based on the Discontinuous Galerkin Method

Oscillatory Responses of a Spatially Random Seabed in the Vicinity of Structures: Validation and Application of the Integrated CFD-SFEM

A CIP-based Numerical Simulation of Wave Interaction with a Fluid-Filled Membrane Submerged Breakwater

Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling

An Unsplit CFS-PML Scheme for the Second-Order Wave Equation with Its Application in Fractional Viscoacoustic Simulation

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

Wavefield Simulation with the Discontinuous Galerkin Method for a Poroelastic Wave Equation in Triple-Porosity Media

Implementation of ADE-CFS-PML for the Single Field WCS-FDTD Method

Nonsplit Complex-Frequency-shifted Perfectly Matched Layer Combined with Symplectic Methods for Solving Second-Order Seismic Wave Equations - Part 2: Wavefield Simulations

A hybrid PML formulation for the 2D three-field dynamic poroelastic equations

Coupled Finite Particle Method for Simulations of Wave and Structure Interaction

Nonsplit Complex-Frequency Shifted Perfectly Matched Layer Combined with Symplectic Methods for Solving Second-Order Seismic Wave Equations — Part 1: Method

A Novel P/S Decoupling Scheme with an Exact Riemann Solver on Coupling Fluid–Solid Media

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

Unified Riemann Solution for Multi-Physics Coupling: Anisotropic Poroelastic/elastic/fluid Interfaces

Full-anisotropic Poroelastic Wave Modeling: A Discontinuous Galerkin Algorithm with a Generalized Wave Impedance

Scaled Boundary Perfectly Matched Layer for Wave Propagation in a Three-Dimensional Poroelastic Medium

Amplitude-preserving P-/S-wavefield Separation with the Discontinuous Galerkin Method on Unstructured Meshes

A Compact Upwind Flux with More Physical Insight for Wave Propagation in 3-D Poroelastic Media