Abstract:In the numerical computation of wave equations, numerical dispersion is a persistent problem arising from inadequate discretization of the continuous wave equation. To thoroughly understand the mechanism of numerical dispersion, we separately analyze the numerical dispersion relations of time-stepping and spatial discretization schemes by Fourier analysis. The relevant results show that the numerical dispersion errors of time-marching schemes depend on the time step length or the Courant number, whereas the numerical dispersion errors of spatial discretization schemes are determined by the error between the eigenvalues of the numerical spatial differential operator and the continuous spatial differential operator. We also find that the much better numerical dispersion accuracy of the stereo-modeling discrete (SMD)-type operator can be attributed to the inclusion of diversified basis functions for its eigenvalue. Based on these findings, we combine the optimal four-stage symplectic partitioned Runge-Kutta and eighth-order SMD as a new fully discrete scheme. The subsequent analysis of its normalized phase velocity is consistent with the numerical dispersion analysis in semidiscrete forms. This is followed by an acoustic wave simulation in a homogeneous model and corresponding computational efficiency comparison. The results show that the new scheme is much more accurate and antidispersive on a coarse grid. In the final two numerical experiments, we use the new scheme to model the acoustic wave propagation in a three-layer model and the Marmousi model. The convolutional perfectly matched layer is applied to eliminate artificial boundary reflections. Our semidiscrete numerical dispersion analysis provides an efficient tool to quantitatively evaluate the time-stepping and spatial discretization schemes. It can facilitate the development of more accurate fully discrete numerical schemes for solving seismic wave equations.

Space-time Spectral Element Method Solution for the Acoustic Wave Equation and Its Dispersion Analysis

The research of space-time coupled spectral element method for acoustic wave equation

Time-Space Coupled Spectral Element Method for Acoustic Wave Propagation on a Rigid Cylinder

An Accurate Triangular Spectral Element Method-Based Numerical Simulation for Acoustic Problems in Complex Geometries

Spectral Element Method for Wave Propagation on Irregular Domains

Solutions of Surface-Wave Dispersion and Attenuation in Stratified Viscoelastic Media Using a Spectral-Element Approach

SPECTRAL ELEMENT METHOD FOR ACOUSTIC PROPAGATION PROBLEMS BASED ON LINEARIZED EULER EQUATIONS

Implicit Chebyshev spectral element method for acoustics wave equations

Spectral element method for acoustic wave simulation in heterogeneous media

A Wideband Fast Multipole Boundary Element Method for Half-Space/plane-symmetric Acoustic Wave Problems

Seismic Wave Modeling in Viscoelastic VTI Media Using Spectral Element Method

Application of Spectral Element Method Based on Convolution Filtering in Long-Term Wavefield Modeling

Chebyshev Spectral Elements Method for 2-Dimensional Acoustic Propagation Problem in a Uniform Mean Flow with Absorbing Boundary Condition

Chebychev Spectral Elements Method for Wave Equation with Absorbing Boundary Conditions

A Spectral Method for Depth-Separated Solution of a Wavenumber Integration Model in Horizontally Stratified Fluid Acoustic Waveguides

Dispersion-Relation Preserving Stereo-Modeling Method Beyond Nyquist Frequency for Acoustic Wave Equation

A cell‐based smoothed radial point interpolation method with virtual nodes for three‐dimensional mid‐frequency acoustic problems

Simulation of Acoustic Wavefields in Heterogeneous Media: A Robust Method for Automatic Suppression of Numerical Dispersion

A New High-Order Scheme Based on Numerical Dispersion Analysis of the Wave Phase Velocity for Semidiscrete Wave Equations

Symplectic Interior Penalty Discontinuous Galerkin Method for Solving the Seismic Scalar Wave Equation

Fast Multipole Boundary Element Method for the Acoustic Analysis of Finite Periodic Structures