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.

The posteriori analysis of the spectral element discretization of the wave equation

Time and space adaptivity of the wave equation discretized in time by a second order scheme

A posteriori error estimates for the wave equation with mesh change in the leapfrog method

Semi and fully-discrete analysis of lowest-order nonstandard finite element methods for the biharmonic wave problem

A posteriori error analysis of a space-time hybridizable discontinuous Galerkin method for the advection-diffusion problem

Adaptive mesh strategies for the spectral element method

A posteriori error analysis for a continuous space-time finite element method for a hyperbolic integro-differential equation

A posteriori error estimates and space--time adaptivity for parabolic partial differential equations on stationary surfaces

Spectral Difference method with a posteriori limiting: Application to the Euler equations in one and two space dimensions

An easily computable error estimator in space and time for the wave equation

Robust a Posteriori Error Estimates of Time-Dependent Poisson–Nernst–Planck Equations

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

A posteriori error analysis of semilinear parabolic interface problems using elliptic reconstruction

Frequency-Explicit A Posteriori Error Estimates for Discontinuous Galerkin Discretizations of Maxwell's Equations

Error analysis of an implicit-explicit time discretization scheme for semilinear wave equations with application to multiscale problems

A finite difference - discontinuous Galerkin method for the wave equation in second order form

Error analysis for a spectral element method for solving two-parameter singularly perturbed diffusion equation

An alternating direction implicit finite element Galerkin method for the linear Schrödinger equation

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

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

A-posteriori diffusion analysis of higher-order numerical schemes for application to propagating linear waves