Abstract:Symplectic integrators are well known for their excellent performance in solving partial differential equation of dynamical systems because they are capable of preserving some conservative properties of dynamic equations. However, there are not enough high-order, for example third-order symplectic schemes, which are suitable for seismic wave equations. Here, we propose a strategy to construct a symplectic scheme that is based on a so-called high-order operator modification method. We first employ a conventional two-stage Runge–Kutta–Nyström (RKN) method to solve the ordinary differential equations, which are derived from the spatial discretization of the seismic wave equations. We then add a high-order term to the RKN method. Finally, we obtain a new third-order symplectic scheme with all positive symplectic coefficients, and it is defined based on the order condition, the symplectic condition, the stability condition and the dispersion relation. It is worth noting that the new scheme is independent of the spatial discretization type used, and we simply apply some finite difference operators to approximate the spatial derivatives of the isotropic elastic equations for a straightforward discussion. For the theoretical analysis, we obtain the semi-analytic stability conditions of our scheme with various orders of spatial approximation. The stability and dispersion properties of our scheme are also compared with conventional schemes to illustrate the favorable numerical behaviors of our scheme in terms of precision, stability and dispersion characteristics. Finally, three numerical experiments are employed to further demonstrate the validity of our method. The modified strategy that is proposed in this paper can be used to construct other explicit symplectic schemes.

A General Symplectic Scheme with Three Free Parameters and Its Applications.

A General Symplectic Integrator for Canonical Hamiltonian Systems

A Symplectic Approximation with Nonlinear Stability and Convergence Analysis for Efficiently Solving Semi-Linear Klein–Gordon Equations

Convergence of Time-Splitting Energy-Conserved Symplectic Schemes for 3d Maxwell'S Equations

Applications of the Multi-Symplectic Euler-Box Scheme

Symplectic Properties Research for Finite Element Methods of Hamiltonian System

High order symplectic schemes for the sine-Gordon equation

Symplectic method of nonlinear Schr(o)dinger equation Mathematics Teaching Group

Local Energy-Preserving Scalar Auxiliary Variable Approaches for General Multi-Symplectic Hamiltonian PDEs

Two Novel Classes of Arbitrary High-Order Structure-Preserving Algorithms for Canonical Hamiltonian Systems

The symplectic schemes of Schr(o)dinger equation

A Novel Energy-conserving Scheme for Eight-dimensional Hamiltonian Problems

The Exponential Invariant Energy Quadratization Approach for General Multi-Symplectic Hamiltonian PDEs

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

A Note on the Generating Function Method

A modified symplectic scheme for seismic wave modeling

Novel semi-explicit symplectic schemes for nonseparable stochastic Hamiltonian systems

A New Class of Symplectic Methods for Stochastic Hamiltonian Systems

Symplectic Method Based on Dual Variable Principle

Symmetric-Adjoint and Symplectic-Adjoint Runge-Kutta Methods and Their Applications

Generalized flow-composed symplectic methods for post-Newtonian Hamiltonian systems