A Lax-Wendroff lumped mass finite element method for seismic simulations
Liu Shaolin1, Li Xiaofan1, Wang Wenshuai1,2, Liu Youshan1
1. Key Laboratory of Earth and Planetary Physics, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
2. School of Mathematics and Computer Science, Ningxia University, Yinchuan, Ningxia 750021, China
Abstract：Using lumped mass finite element method with triangle mesh to solve elastic wave equations and Lax-Wendroff method to obtain fourth-order temporal accuracy, we develop Lax-Wendroff lumped mass finite element method (LWFEM) for elastic wave simulations. To avoid artificial boundary reflections, the second-order PML absorbing boundary condition (PML ABC) in terms of displacement is constructed. The stability criterion on a period grid is analyzed based on constructing generalized eigenvalue problems of the finite element method. In numerical experiments, the accuracy and efficiency are discussed by comparing with conventional methods such as central difference finite element method (CDFEM), Runge-Kutta finite element method (RKFEM) and Newmark spectral element method (NSEM). Numerical results demonstrate the validity of LWFEM in complex models.
刘少林, 李小凡, 汪文帅, 刘有山. Lax-Wendroff集中质量有限元法地震波场模拟[J]. 石油地球物理勘探, 2015, 50(5): 905-911,924.
Liu Shaolin, Li Xiaofan, Wang Wenshuai, Liu Youshan. A Lax-Wendroff lumped mass finite element method for seismic simulations. OGP, 2015, 50(5): 905-911,924.
Bielak J,Ghattas O,Kim E J. Parallel octree-based finite element method for large-scale earthquake ground motion simulation. Computer Modeling in Engineering and Sciences,2005,10(2):99-112.
Li Weihua,Liu Qinghua,Zhao Chenggang. Three-dimensional viscous-spring boundaries in time domain and dynamic analysis using explicit finite element method of saturated porous medium. Chinese Journal of Geophysics,2010,53(10):2460-2469.
[7]
Marfurt K J. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics,1984,49(5):533-549.
[8]
Pasquetti R,Rapetti F. Spectral element methods on triangles and quadrilaterals:comparisons and applications. Journal of Computational Physics,2004,198(1):349-362.
[9]
Pasquetti R,Rapetti F. Spectral element methods on unstructured meshes:comparisons and recent advances. Journal of Scientific Computing,2006,27(1-3):377-387.
Liu Youshan,Teng Jiwen,Liu Shaolin et al. Explicit finite element method with triangle meshes stored by spare format and its perfectly match layer absorbing boundary condition. Chinese Journal of Geophysics,2013,56(9):3085-3099.
[11]
Komatitsch D,Erlebacher G,Göddeke D et al. High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster. Journal of Computational Physics,2010,229(20):7692-7714.
[12]
Wu S R. Lumped mass matrix in explicit finite element method for transient dynamics of elasticity. Computer Methods in Applied Mechanics and Engineering,2006,195(44-47):5983-5994
[13]
Etienne V,Chaljub E,Virieux J et al. An hp-adaptive discontinuous Galerkin finite-element method for 3-D elastic wave modelling. Geophysical Journal International,2010,183(2):941-962.
[14]
Zhang J F,Gao H W. Irregular perfectly matched layers for 3D elastic wave modeling. Geophysics,2011,76(2):T27-T36.
[15]
Liu H,Dai N,Niu F et al. An explicit time evolution method for acoustic wave propagation. Geophysics,2014,79(3):T117-T124.
[16]
Lax P D,Wendroff B. Difference schemes for hyperbolic equations with high order of accuracy. Communications on Pure and Applied Mathematics,1964,17(3):381-398.
[17]
Dablain M. The application of high-order differencing to the scalar wave equation. Geophysics,1986,51(1):54-66.
[18]
Liu Y,Sen M K. Time-space domain dispersion-relation-based finite-difference method with arbitrary even-order for the 2D acoustic wave equation. Journal of Computational Physics,2013,232(1):327-345.
[19]
Tan S R,Huang L J. An efficient finite-difference method with high-order accuracy in both time and space domains for modelling scalar-wave propagation. Geophysical Journal International,2014,197(1):1250-1267.
[20]
Long G H,Zhao Y B,Zhou J. A temporal fourth-order scheme for the first-order acoustic wave equations. Geophysical Journal International,2013,194(3):1473-1485.
[21]
Tong P,Yang D H,Hua B L et al. A high-order stereo-modeling method for solving wave equations. Bulletin of Seismological Society of America,2013,103(2A):811-833.
[22]
De Basabe J D,Sen M K. Stability of the high-order finite elements for acoustic or elastic wave propagation with high-order time stepping. Geophysical Journal International,2010,181(1):577-590.
Chen Shan,Yang Dinghui,Deng Xiaoying. An improved algorithm of the fourth-order Runge-Kutta method and seismic wave-field simulation. Chinese Journal of Geophysics,2010,53(5):1196-1206.
[24]
Liu Y S,Teng J W,Lan H Q et al. A comparative study of finite element and spectral element methods in seismic wavefield modeling. Geophysics, 2014, 79(2):T91-T104.
Zhang Meigen,Wang Miaoyue,Li Xiaofan et al. Full wavefield inversion of anisotropic elastic parameters in the time domain. Chinese Journal of Geophysics,2003,46(1):94-100.
[26]
Komatitsch D,Tromp J. A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation. Geophysical Journal International,2003,154(1):146-153.
[27]
Chew W C and Liu Q H. Perfectly matched layers for elastodynamics:A new absorbing boundary condition. Journal of Computational Acoustics,1996,4(4):341-359.
[28]
Komatitsch D,Barnes C,Tromp J. Wave propagation near a fluid-solid interface:A spectral-element approach. Geophysics,2000,65(2):623-631.
[29]
Komatitsch D,Barnes C,Tromp J. Simulation of anisotropic wave propagation based upon a spectral element method. Geophysics,2000,65(4):1251-1260.
[30]
De Basabe J D,Sen M K. Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations. Geophysics,2007,72(6):T81-T95.
[31]
de Hoop A T. A modification of Cagniard's method for solving seismic pulse problems. Applied Scientific Research,1960,8(1):349-356.
