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.

Numerical analysis of the diffusive-viscous wave equation

NUMERICAL WAVE OF FINITE ELEMENT SOLUTION IN FLUID MECHANICS USING DIFFERENT INTERPOLATION FUNCTIONS

Numerical Solution of a Two-Dimensional Nonlocal Wave Equation on Unbounded Domains

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

Numerical analysis for high-order methods for variable-exponent fractional diffusion-wave equation

A finite volume scheme employing the multipoint flux approximation with diamond stencil for the diffusive‐viscous wave equation on general polyhedral meshes

A high order discontinuous Galerkin method for the symmetric form of the anisotropic viscoelastic wave equation

A finite-element algorithm with a perfectly matched layer boundary condition for seismic modelling in a diffusive-viscous medium

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

Numerical approximation of discontinuous solutions of the semilinear wave equation

Discontinuous Galerkin method for solving wave equations in two-phase and viscoelastic media

A high order finite difference method for the elastic wave equation in bounded domains with nonconforming interfaces

A Modified Numerical-Flux-based Discontinuous Galerkin Method for 2D Wave Propagations in Isotropic and Anisotropic Media

Efficient Ordinary Differential Equation-Based Discontinuous Galerkin Method for Viscoelastic Wave Modeling

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

Numerical Studies on the Classical Long Wave System Describing Shallow Water Waves With Dispersion

Dispersion-dissipation Analysis of Triangular Numerical-Flux-based Discontinuous Galerkin Method for Elastic Wave Equations

The Numerical Solution of Parabolic Volterra Integro-Differential Equations on Unbounded Spatial Domains

Numerical methods for the forward and backward problems of a time-space fractional diffusion equation

Numerical Simulation of Thermoelastic Nonlinear Waves in Fluid Saturated Porous Media with Non-local Darcy Law

Artificial Boundary Conditions and Finite Difference Approximations for a Time-Fractional Diffusion-Wave Equation on a Two-Dimensional Unbounded Spatial Domain