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.

Compact ETDRK Scheme for Nonlinear Dispersive Wave Equations

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

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

A fourth-order exponential time differencing scheme with dimensional splitting for non-linear reaction-diffusion systems

A fast compact finite difference scheme for the fourth-order diffusion-wave equation

An Efficient High-Order Explicit Scheme for Solving Hamiltonian Nonlinear Wave Equations

A Compact Exponential Scheme for Solving 1D Unsteady Convection-Diffusion Equation with Neumann Boundary Conditions

Fourth-Order Accurate Compact Scheme for First-Order Maxwell's Equations

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

Linearly-fitted energy-mass-preserving schemes for Korteweg–de Vries equations

An Exponential Time Differencing Runge–Kutta Method ETDRK32 for Phase Field Models

A COMPACT DIFFERENCE SCHEME FOR FOURTH-ORDER FRACTIONAL SUB-DIFFUSION EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

High-order exponential time differencing multi-resolution alternative finite difference WENO methods for nonlinear degenerate parabolic equations

A Weighted Runge-Kutta Method with Weak Numerical Dispersion for Solving Wave Equations

High-Order Dissipation-Preserving Methods for Nonlinear Fractional Generalized Wave Equations

A fourth-order compact solver for fractional-in-time fourth-order diffusion equations

Compact difference scheme for two‐dimensional fourth‐order nonlinear hyperbolic equation

High-Order Compact Implicit Difference Method For the Wave Equation

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

Efficient second-order semi-implicit finite element method for fourth-order nonlinear diffusion equations