Abstract:The construction of low dispersive and strong stable numerical schemes is critical for the investigation of wave propagation. Solving seismic wave equations in the time domain, the dispersion and stability are related to the temporal and spatial discretizations. Here we focus on the temporal discretization and develop a high-efficient time integration scheme for the elastic wave equation. Following the transformation of the elastic seismic wave equation into a Hamiltonian system, an explicit second-order symplectic Runge-Kutta-Nystrom (RKN) scheme is proposed for elastic wave modeling. The order conditions are obtained by the rooted trees theory. To obtain symmetric solutions, we assign a constraint condition to the symplectic coefficients. In order to minimize the error in temporal discretization, we propose a method that compares numerical angular frequency with real frequency, and obtain a constraint equation based on Taylor expansion. A symplectic RKN scheme with minimized numerical dispersion is developed by solving the order conditions associating with constraint equations. Through analyzing the time advance equation, we construct another constraint equation to allow stability limit to attain its maximal value, and then another symplectic RKN scheme with the same symplectic coefficients as the first scheme is developed. The optimization process for determining the coefficients of the RKN scheme is independent of spatial discretizations. For the ease of comparison, we simply choose finite difference method to approximate spatial derivatives. Generally, the numerical dispersion amount is lower than the conventional second-order schemes, and is slightly larger than the high-order schemes. The maximal Courant number for the new scheme is the largest compared to those for the conventional schemes. The largest time increment for the new scheme is nearly two times as large as those for the conventional second-order schemes, and is about 1.5 times as large as those for the high-order ones. In a practical simulation, all these schemes discussed in the paper consume approximately the same computer memory, but the numerical accuracy and computational speed are quite different. The numerical accuracy and complexity of our scheme indicate that it may be a good candidate for the balance of numerical accuracy and complexity. Both theoretical analyses and numerical experiments show that our scheme is superior to conventional symplectic schemes in some aspects including suppressing dispersion and increasing stability range. The optimization procedure proposed in the paper may be used to define the coefficients in other symplectic or non-symplectic schemes. Though striking numerical results are obtained by our scheme with a finite difference operator, a low-dispersive spatial operator is still needed to be combined with our scheme to form a more powerful tool for wave simulations.

A symplectic RKN scheme for solving elastic wave equations

A modified symplectic scheme for seismic wave modeling

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

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

A New Kind of Optimal Second-Order Symplectic Scheme for Seismic Wave Simulations

A Modified Symplectic PRK Scheme for Seismic Wave Modeling

Symplectic partitioned Runge-Kutta method for solving the acoustic wave equation

Modified Symplectic Schemes with Nearly-Analytic Discrete Operators for Acoustic Wave Simulations

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

A Runge-Kutta method using high-order interpolation approximation for solving 2D acoustic and elastic wave equations

Modified symplectic scheme with finite element method for seismic wavefield modeling

3D numerical simulation of frequency-domain elastic wave equation with a compact second-order scheme

A modified time-space optimized symplectic method for solving 3D elastic wave equations

An affine generalized optimal scheme with improved free-surface expression using adaptive strategy for frequency-domain elastic wave equation

High-order Temporal and Implicit Spatial Staggered-Grid Finite-Difference Operators for Modelling Seismic Wave Propagation

A nearly analytic symplectically partitioned Runge-Kutta method for 2-D seismic wave equations

An Improved Algorithm of the Fourth-Order Runge-Kutta Method and Seismic Wave-Field Simulation

Multi-symplectic exponentially fitted RKN methods for wave equations

THE SEM BASED ON SYMPLECTICAL SCHEMES AND ITS APPLICATION IN MODELING THE WAVE PROPAGATION IN TRANSVERSELY ISOTROPIC MEDIA

Modeling of the Acoustic Wave Equation by Staggered‐Grid Finite‐Difference Schemes with High‐Order Temporal and Spatial Accuracy

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