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.

Stability And Accuracy Of The Cubic Interpolated Propagation Scheme

An Optimized Dispersion-Relation-Preserving Combined Compact Difference Scheme to Solve Advection Equations

Error estimates of the cubic interpolated pseudo-particle scheme for one-dimensional advection equations

A Numerical Method For Solving The One-Dimensional Vlasov-Poisson Equation In Phase Space

Stability and numerical dispersion of symplectic fourth-order time-domain schemes for optical field simulation

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

A 2DH fully dispersive and weakly nonlinear Boussinesq-type model based on a finite-volume and finite-difference TVD-type scheme

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

An Efficient and high accuracy finite-difference scheme for the acoustic wave equation in 3D heterogeneous media

A conservative semi-Lagrangian method for oscillation-free computation of advection processes

Efficient and Stable Finite Difference Modelling of Acoustic Wave Propagation in Variable-density Media

An Optimized Time-Space-domain Finite Difference Method with Piecewise Constant Interpolation Coefficients for Scalar Wave Propagation

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

An Efficient Fourth-Order Three-Point Scheme for Solving Some Nonlinear Dispersive Wave Equations

Generalized Transfer Matrix Method for Propagation of Surface Waves in Layered Azimuthally Anisotropic Half-Space

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

A time-space domain stereo finite difference method for 3D scalar wave propagation

A space-time generalized finite difference scheme for long wave propagation based on high-order Korteweg-de Vries type equations

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

Rusanov-type Schemes for Hyperbolic Equations: Wave-Speed Estimates, Monotonicity and Stability

A Finite Difference Scheme based on Cubic Trigonometric B-splines for Time Fractional Diffusion-wave Equation