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

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

Year:    2010

Communications in Computational Physics, Vol. 7 (2010), Iss. 5 : pp. 1027–1048

Abstract

In this paper, we propose a weighted Runge-Kutta (WRK) method to solve the 2D acoustic and elastic wave equations. This method successfully suppresses the numerical dispersion resulted from discretizing the wave equations. In this method, the partial differential wave equation is first transformed into a system of ordinary differential equations (ODEs), then a third-order Runge-Kutta method is proposed to solve the ODEs. Like the conventional third-order RK scheme, this new method includes three stages. By introducing a weight to estimate the displacement and its gradients in every stage, we obtain a weighted RK (WRK) method. In this paper, we investigate the theoretical properties of the WRK method, including the stability criteria, numerical error, and the numerical dispersion in solving the 1D and 2D scalar wave equations. We also compare it against other methods such as the high-order compact or so-called Lax-Wendroff correction (LWC) and the staggered-grid schemes. To validate the efficiency and accuracy of the method, we simulate wave fields in the 2D homogeneous transversely isotropic and heterogeneous isotropic media. We conclude that the WRK method can effectively suppress numerical dispersions and source noises caused in using coarse grids and can further improve the original RK method in terms of the numerical dispersion and stability condition.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.2009.09.088

Communications in Computational Physics, Vol. 7 (2010), Iss. 5 : pp. 1027–1048

Published online:    2010-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    22

Keywords:   

  1. Dispersion-relation preserving stereo-modeling method beyond Nyquist frequency for acoustic wave equation

    Jing, Hao | Chen, Yushu | Yang, Dinghui | Wu, Rushan

    GEOPHYSICS, Vol. 82 (2017), Iss. 1 P.T1

    https://doi.org/10.1190/geo2016-0104.1 [Citations: 16]
  2. A nearly analytic symplectically partitioned Runge-Kutta method for 2-D seismic wave equations

    Ma, Xiao | Yang, Dinghui | Liu, Faqi

    Geophysical Journal International, Vol. 187 (2011), Iss. 1 P.480

    https://doi.org/10.1111/j.1365-246X.2011.05158.x [Citations: 44]
  3. Using a stereo-modeling operator to reduce aliasing during migration

    Li, Jingshuang | Yang, Dinghui | Fehler, Michael | Huang, Xueyuan

    GEOPHYSICS, Vol. 83 (2018), Iss. 5 P.S437

    https://doi.org/10.1190/geo2017-0514.1 [Citations: 1]
  4. 3D weak-dispersion reverse time migration using a stereo-modeling operator

    Li, Jingshuang | Fehler, Michael | Yang, Dinghui | Huang, Xueyuan

    GEOPHYSICS, Vol. 80 (2015), Iss. 1 P.S19

    https://doi.org/10.1190/geo2013-0472.1 [Citations: 13]
  5. Dynamic inversion method based on the time-staggered stereo-modeling scheme and its acceleration

    Jing, Hao | Yang, Dinghui | Wu, Hao

    Geophysical Journal International, Vol. 207 (2016), Iss. 3 P.1675

    https://doi.org/10.1093/gji/ggw367 [Citations: 3]
  6. An anti‐dispersion reverse‐time migration method with local nearly‐analytic operators and its application

    Li, Jingshuang | Yang, Dinghui | Liu, Faqi | Hua, Biaolong

    SEG Technical Program Expanded Abstracts 2011, (2011), P.3340

    https://doi.org/10.1190/1.3627891 [Citations: 1]
  7. An efficient symplectic reverse time migration method using a local nearly analytic discrete operator in acoustic transversely isotropic media with a vertical symmetry axis

    Feng, Lanlan | Yang, Dinghui | Xie, Wei

    GEOPHYSICS, Vol. 80 (2015), Iss. 4 P.S103

    https://doi.org/10.1190/geo2014-0248.1 [Citations: 2]
  8. A nearly analytic exponential time difference method for solving 2D seismic wave equations

    Zhang, Xiao | Yang, Dinghui | Song, Guojie

    Earthquake Science, Vol. 27 (2014), Iss. 1 P.57

    https://doi.org/10.1007/s11589-013-0056-6 [Citations: 6]
  9. A new high-order scheme based on numerical dispersion analysis of the wave phase velocity for semidiscrete wave equations

    Ma, Xiao | Yang, Dinghui | He, Xijun | Li, Jingshuang | Zheng, Yongchang

    GEOPHYSICS, Vol. 83 (2018), Iss. 3 P.T123

    https://doi.org/10.1190/geo2017-0441.1 [Citations: 10]
  10. An efficient reverse time migration method using local nearly analytic discrete operator

    Li, Jingshuang | Yang, Dinghui | Liu, Faqi

    GEOPHYSICS, Vol. 78 (2013), Iss. 1 P.S15

    https://doi.org/10.1190/geo2012-0247.1 [Citations: 17]
  11. An efficient fourth-order low dispersive finite difference scheme for a 2-D acoustic wave equation

    Das, Sambit | Liao, Wenyuan | Gupta, Anirudh

    Journal of Computational and Applied Mathematics, Vol. 258 (2014), Iss. P.151

    https://doi.org/10.1016/j.cam.2013.09.006 [Citations: 22]
  12. A Nearly‐Analytic Symplectic PRK method for Solving 2D Acoustic and Elastic Equations

    Ma, Xiao | Yang, Dinghui

    SEG Technical Program Expanded Abstracts 2010, (2010), P.2993

    https://doi.org/10.1190/1.3513468 [Citations: 4]
  13. 3D weak-dispersion reverse time migration with a stereomodeling method

    Jingshuang, Li | Yang, Dinghui | Huang, Xueyuan | Fehler, Mike

    SEG Technical Program Expanded Abstracts 2013, (2013), P.3995

    https://doi.org/10.1190/segam2013-0996.1 [Citations: 0]
  14. Full waveform inversion based on deep learning and optimal nearly analytic discrete method

    Fan, Lu | Yan-Jie, Zhou | Xi-Jun, He | Xiao, Ma | Xue-Yuan, Huang

    Applied Geophysics, Vol. 18 (2021), Iss. 4 P.483

    https://doi.org/10.1007/s11770-021-0912-4 [Citations: 3]
  15. 3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

    Zou, Qiang | Huang, Jian-Ping | Yong, Peng | Li, Zhen-Chun

    Petroleum Science, Vol. 17 (2020), Iss. 4 P.967

    https://doi.org/10.1007/s12182-020-00477-3 [Citations: 10]
  16. A Strong Stability-Preserving Predictor-Corrector Method for the Simulation of Elastic Wave Propagation in Anisotropic Media

    Yang, D. H. | Wang, N. | Liu, E.

    Communications in Computational Physics, Vol. 12 (2012), Iss. 4 P.1006

    https://doi.org/10.4208/cicp.010111.230911a [Citations: 13]
  17. A low-dispersive method using the high-order stereo-modelling operator for solving 2-D wave equations

    Li, Jingshuang | Yang, Dinghui | Wu, Hao | Ma, Xiao

    Geophysical Journal International, Vol. 210 (2017), Iss. 3 P.1938

    https://doi.org/10.1093/gji/ggx277 [Citations: 13]