Simulation of Earthquake Rupture Dynamics in Complex Geometries Using Coupled Finite Difference and Finite Volume Methods
2015 (English)In: Communications in Computational Physics, ISSN 1815-2406, E-ISSN 1991-7120, Vol. 17, no 2, 337-370 p.Article in journal (Refereed) Published
We couple a node-centered finite volume method to a high order finite difference method to simulate dynamic earthquake ruptures along nonplanar faults in two dimensions. The finite volume method is implemented on an unstructured mesh, providing the ability to handle complex geometries. The geometric complexities are limited to a small portion of the overall domain and elsewhere the high order finite difference method is used, enhancing efficiency. Both the finite volume and finite difference methods are in summation-by-parts form. Interface conditions coupling the numerical solution across physical interfaces like faults, and computational ones between structured and unstructured meshes, are enforced weakly using the simultaneousapproximation-term technique. The fault interface condition, or friction law, provides a nonlinear relation between fields on the two sides of the fault, and allows for the particle velocity field to be discontinuous across it. Stability is proved by deriving energy estimates; stability, accuracy, and efficiency of the hybrid method are confirmed with several computational experiments. The capabilities of the method are demonstrated by simulating an earthquake rupture propagating along the margins of a volcanic plug.
Place, publisher, year, edition, pages
Global Science Press, 2015. Vol. 17, no 2, 337-370 p.
elastic waves, earthquake, high order finite difference finite volume, summation-by-parts, simultaneous approximation term, nonlinear boundary conditions
Computational Mathematics Mathematics
IdentifiersURN: urn:nbn:se:liu:diva-113745DOI: 10.4208/cicp.111013.120914aISI: 000353693400002OAI: oai:DiVA.org:liu-113745DiVA: diva2:784471