Abstract
The advantage of the finite element method (FEM) lies in its flexibility in addressing rugged interfaces in complex geological models. However, the efficiency of the FEM is relatively low for large-scale seismic wave modelling. Here, we introduce an order-corrected symplectic FEM (OCSFEM) with structure-preserving properties and parsimonious memory requirements for the elastic wave equation. In this method, the storage of the large sparse stiffness matrix is changed to the storage of the element Jacobian matrix. An efficient order-corrected symplectic method with third-order temporal accuracy is combined with a triangle-based FEM to construct the OCSFEM. The structure-preserving characteristics and high efficiency of the OCSFEM facilitate the high-fidelity modelling of large-scale and long-term wave phenomena. Complex and large-scale numerical examples show that the OCSFEM exhibits low numerical dispersion and high stability compared with conventional methods, such as the second-order symplectic FEM.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
