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.

Improvement of accuracy of the spectral element method for elastic wave computation using modified numerical integration operators

Accuracy and Stability Enhancements in the Incompressible Finite-Volume-particle Method for Multiphase Flow Simulations.

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

Reducing the Dispersion Error in Acoustic Wave Modeling with Modified Integration Rules

Implicit Spectral Element Method for Wave Equation and Some Factors Influencing Numerical Accuracy

K-space dispersion error compensators for the fractional spatial derivatives based constant-Q viscoelastic wave equation modeling

Application of modified integration rule to time-domain finite-element acoustic simulation of rooms.

A Legendre spectral element method with optimal mass matrix for seismic wave modeling

A low-dispersive method using the high-order stereo-modelling operator for solving 2-D wave equations

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

A 2.5d Finite-Difference Algorithm for Elastic Wave Modeling Using Near-Optimal Quadratures

High Accuracy Wave Simulation – Revised Derivation, Numerical Analysis and Testing of a Nearly Analytic Integration Discrete Method for Solving Acoustic Wave Equation

An novel finite difference dispersion error elimination mechanism in the Lax–Wendroff high‐order time discretization

An Effective Method to Suppress Numerical Dispersion in 2D Acoustic and Elastic Modelling Using a High-Order Padé Approximation

Elastic Wave Modelling Method Based on the Displacement–velocity Fields: an Improving Nearly Analytic Discrete Approximation

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

A 2.5D Elastic Wave Modeling Algorithm Using Near-optimal Quadratures

A 2.5 D finite-difference algorithm for elastic wave modeling using near-optimal quadratures2. 5D elastic wave modeling

Spectral-element Method with an Optimal Mass Matrix for Seismic Wave Modelling

A new discontinuous Galerkin spectral element method for elastic waves with physically motivated numerical fluxes

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