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.

Long-time Momentum and Actions Behaviour of Energy-Preserving Methods for Semi-Linear Wave Equations Via Spatial Spectral Semi-Discretisations

An Energy-Preserving Algorithm for Nonlinear Hamiltonian Wave Equations with Neumann Boundary Conditions

NUMERICAL ANALYSIS OF AN ENERGY-CONSERVATION SCHEME FOR TWO-DIMENSIONAL HAMILTONIAN WAVE EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

Numerical Conservation of Energy, Momentum and Actions for Extended RKN Methods when Applied to Nonlinear Wave Equations Via Spatial Spectral Semi-Discretizations

Numerical solutions of the nonlinear wave equations with energy-preserving sixth-order finite difference schemes

Long-time Oscillatory Energy Conservation of Total Energy-Preserving Methods for Highly Oscillatory Hamiltonian Systems

An Energy-Preserving and Symmetric Scheme for Nonlinear Hamiltonian Wave Equations

A variant of the discrete gradient method for the solution of the semilinear wave equation under different boundary conditions

A numerical method with properties of consistency in the energy domain for a class of dissipative nonlinear wave equations with applications to a Dirichlet boundary-value problem

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

Energy conserving methods for Hamiltonian PDEs based on spectral space decomposition

Long-term Analysis of Symplectic or Symmetric Extended RKN Methods for Nonlinear Wave Equations

The Conservative Time High-Order AVF Compact Finite Difference Schemes for Two-Dimensional Variable Coefficient Acoustic Wave Equations

Two linear energy-preserving compact finite difference schemes for coupled nonlinear wave equations

A novel explicit fully-discrete momentum-preserving scheme of damped nonlinear stochastic wave equation influenced by multiplicative space-time noise

Energy-preserving splitting finite element method for nonlinear stochastic space-fractional wave equations with multiplicative noise

Structure-Preserving Numerical Methods for Two Nonlinear Systems of Dispersive Wave Equations

A conservative and energy stable discontinuous spectral element method for the shifted wave equation in second order form

A Linearly-Fitted Conservative (dissipative) Scheme for Efficiently Solving Conservative (dissipative) Nonlinear Wave PDEs

Linear High-Order Energy-Preserving Schemes for the Nonlinear Schrödinger Equation with Wave Operator Using the Scalar Auxiliary Variable Approach

Modified Symplectic Schemes with Nearly-Analytic Discrete Operators for Acoustic Wave Simulations