Abstract:The complex-frequency-shifted convolutional perfectly matched layer (CFS-CPML) has been widely used in numerical simulations of the first-order seismic wave equations to avoid artificial reflections caused by truncated boundaries. However, numerically solving the second-order seismic wave equation is preferred for some large-scale geophysical algorithms, such as reverse-time migration (RTM). The extension of CFS-CPML to the second-order seismic wave equations has been realized by a strict derivation approach in the existing literature. However, the derived CFS-CPML formulations are complex, due to the introduction of a large number of intermediate variables, which decreases the overall computational efficiency, visibly. We present an efficient strategy to implement CFS-CPML in finite-difference time-domain (FDTD) modeling of the second-order seismic wave equations. We do not derive any continuous CFS-CPML formulations. Instead, we propose to split the high-order centered-grid finite-difference (CGFD) stencil that is applied to approximate the spatial second-order derivatives into two levels of CGFD stencils. The inner stencils are viewed as the second-order CGFD approximations of the first-order derivatives at symmetric grid points that are centered at the grid point where the second-order derivative is specified. With this viewpoint, the convolution term related to CFS-CPML can be introduced to each inner CGFD approximation, naturally. The outer CGFD stencil is also viewed as a CGFD approximation of a first-order derivative and it is augmented by a convolution term as well. In this way, CFS-CPML is incorporated into the high-order CGFD scheme, efficiently. We take three-dimensional (3-D) acoustic, viscoacoustic, and elastic wavefields FDTD modeling, for example, to verify the feasibility of the new implementation strategy of CFS-CPML.

A higher-order spatial FDTD scheme with CFS PML for 3D numerical simulation of wave propagation in cold plasma

Development of a Novel FDTD (2, 4)-Compatible Conformal Scheme for Electromagnetic Computations of Complex Curved PEC Objects

A Modified PML Implementation in FDTD for Handling Dispersive Media

A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

A modified high-order FDTD conformal technique for modeling perfectly conducting objects

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Finite-difference computation of transient electromagnetic waves for cylindrical geometries in complex media

Finite-difference Time-Domain Studies of Plasma Resonating Structures

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

High-order accurate schemes for Maxwell's equations with nonlinear active media and material interfaces

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

Perfectly Matched Layer for Accurate FDTD for Anisotropic Magnetized Plasma

Analysis of electromagnetic wave propagation and scattering characteristics of plasma shealth via high order ADE-ADI FDTD

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

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

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

ANSYS HFSS as a new numerical tool to study wave propagation inside anisotropic magnetized plasmas in the Ion Cylotron Range of Frequencies

Tridiagonalized WLP-FDTD Implementation for Wave Propagation in Magnetized Plasma

High-order FDTD schemes for Maxwell's interface problems with discontinuous coefficients and complex interfaces based on the Correction Function Method

Research on the moving plasma photonic crystals based on the novel symplectic finite-difference time-domain method