liu.seSearch for publications in DiVA
Change search
ReferencesLink to record
Permanent link

Direct link
Simulation of Earthquake Rupture Dynamics in Complex Geometries Using Coupled Finite Difference and Finite Volume Methods
Department of Geophysics, Stanford University, CA 94305-2215, USA.
Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.
Department of Applied Mathematics, Naval Postgraduate School, Monterey, CA 93943-5216, USA.
Institute for Computational and Mathematical Engineering, Stanford University, CA 94305-4042, USA.
2013 (English)Report (Other academic)
Abstract [en]

A numerical method suitable for wave propagation problems in complex geometries is developed for simulating dynamic earthquake ruptures with realistic friction laws. The numerical method couples an unstructured, node-centered finite volume method to a structured, high order finite difference method. In this work we our focus attention on 2-D antiplane shear problems. The finite volume method is used on unstructured triangular meshes to resolve earthquake ruptures propagating along a nonplanar fault. Outside the small region containing the geometrically complex fault, a high order finite difference method, having superior numerical accuracy, is used on a structured grid.

The finite difference method is coupled weakly to the finite volume method along interfaces of collocated grid points. Both methods are on summation-by-parts form. The simultaneous approximation term method is used to weakly enforce the interface conditions. At fault interfaces, fault strength is expressed as a nonlinear function of sliding velocity (the jump in particle velocity across the fault) and a state variable capturing the history dependence of frictional resistance. Energy estimates are used to prove that both types of interface conditions are imposed in a stable manner.

Stability and accuracy of the numerical implementation are verified through numerical experiments, and efficiency of the hybrid approach is confirmed through grid coarsening tests. Finally, the method is used to study earthquake rupture propagation along the margins of a volcanic plug.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2013. , 36 p.
LiTH-MAT-R, ISSN 0348-2960 ; 11
Keyword [en]
elastic waves, earthquake, high order finite difference finite volume, summation-by-parts, simultaneous approximation term, weak nonlinear boundary conditions
National Category
Computational Mathematics
URN: urn:nbn:se:liu:diva-98812ISRN: LiTH-MAT-R--2013/11--SEOAI: diva2:656225
Available from: 2013-10-15 Created: 2013-10-14 Last updated: 2013-10-15

Open Access in DiVA

Simulation of Earthquake Rupture Dynamics in Complex Geometries Using Coupled Finite Difference and Finite Volume Methods(2579 kB)141 downloads
File information
File name FULLTEXT01.pdfFile size 2579 kBChecksum SHA-512
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Nordström, Jan
By organisation
Computational MathematicsThe Institute of Technology
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 141 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 118 hits
ReferencesLink to record
Permanent link

Direct link