A High-Order NAD-SPRK Algorithm Based on the Two Dimensions Elastic Wave Equation

Abstract

Abstract:

Based on the two dimensions elastic wave equation with the Hamiltonian system， an eighth-order NAD-SPRK algorithm is obtained by combining the NAD operator for higher order partial derivatives of discrete space and the symplectic partitioned Runge-Kutta algorithm of discrete time derivative. The stability conditions， numerical dispersion and computational efficiency of the algorithm are studied theoretically and numerically. The results show that compared with the fourth-order NSPRK algorithm， the eighth-order Lax-Wendroff algorithm and the eighth-order staggered grid algorithm， the eighth-order NAD-SPRK algorithm is significantly superior to the traditional numerical algorithm in suppressing the numerical dispersion， and has the minimum numerical error and the highest computational efficiency.

Key words: elastic wave equation, NAD operator, stability analysis, numerical dispersion, computational efficiency

CLC Number:

P631
O241

Chen Li, Zhu Xingwen, Zhang Chaoyuan. A High-Order NAD-SPRK Algorithm Based on the Two Dimensions Elastic Wave Equation[J]. Journal of Dali University, 2022, 7(12): 1-7.

References

〔1〕 DABLAIN M A. The Application of High-Order Differen-cing to Scalar Wave Equation〔J〕. Geophysics， 1986， 51（1）： 54-66.

〔2〕 YANG D H， LIU E， ZHANG Z J， et al. Finite-Difference Modelling in Two-Dimensional Anisotropic Media Using a Flux-Corrected Transport Technique〔J〕. Geophysical Jour-nal International， 2002， 148（2）： 320-328.

〔3〕杨顶辉. 双相各向异性介质中弹性波方程的有限元解法及波场模拟〔J〕. 地球物理学报， 2002， 45（4）： 575-583. 〔4〕孙耿. 波动方程的一类显式辛格式〔J〕. 计算数学， 1997， 19（1）： 1-10.

〔5〕董良国，马在田，曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究〔J〕. 地球物理学报，2000，43（6）：856-864.

〔6〕 YANG D H， TENG J W， ZHANG Z J， et al. A Nearly-Analytic Discrete Algorithm for Acoustic and Elastic Wave Equations in Anisotropic Media〔J〕. Bulletin of the Seismo-logical Society of America，2003，93（2）：882-890.

〔7〕 CHEN S， YANG D H， DENG X Y. A Weighted Runge-Kutta Method with Weak Numerical Dispersion for Solving Wave Equations〔J〕.Communications in Computational Phy-sics， 2010， 7（5）： 1027-1048.

〔8〕王磊，杨顶辉，邓小英. 非均匀介质中地震波应力场的WNAD方法及其数值模拟〔J〕. 地球物理学报， 2009， 52（6）： 1526-1535.

〔9〕 YANG D H， SONG G J， CHEN S， et al. An Improved Nearly Analytical Discrete Algorithm： An Efficient Tool to Simulate the Seismic Response of 2-D Porous Structures〔J〕. Journal of Geophysics and Engineering，2007，4（1）：40-52.

〔10〕 YANG D H， WANG N， CHEN S， et al. An Explicit Algorithm Based on the Implicit Runge-Kutta Algorithm for Solving the Wave Equations〔J〕. Bulletin of the Seis-mological Society of America， 2009， 99： 3340-3354.

〔11〕马啸，杨顶辉，张锦华.求解声波方程的辛可分Runge-Kutta方法〔J〕.地球物理学报，2010，53（8）：1993-2003.

〔12〕 ZHANG C Y， LI X， MA X， et al. A Runge-Kutta Method with Using Eighth-Order Nearly-Analytic Spatial Discretization Operator for Solving a 2D Acoustic Wave Equation〔J〕. Journal of Seismic Exploration， 2014，23： 279-302.

〔13〕 TONG P， YANG D H， HUA B L， et al. A High-Order Stereo-Modeling Method for Solving Wave Equations〔J〕. Bulletin of the Seismological Society of America，2013， 103（2A）：811-833.