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

Direct link
Simulation of Dynamic Earthquake Ruptures in Complex Geometries Using High-Order Finite Difference Methods
Department of Geophysics, Stanford University, Stanford, CA, USA.
Department of Geophysics and Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA, USA.
Linköping University, Department of Mathematics, Scientific Computing. Linköping University, The Institute of Technology.ORCID iD: 0000-0002-7972-6183
2012 (English)Report (Other academic)
Abstract [en]

We develop a stable and high-order accurate finite difference method for problems in earthquake rupture dynamics capable of handling complex geometries and multiple faults. The bulk material is an isotropic elastic solid cut by preexisting fault interfaces. The fields across the interfaces are related through friction laws which depend on the sliding velocity, tractions acting on the interface, and state variables which evolve according to ordinary differential equations involving local fields.

The method is based on summation-by-parts finite difference operators with irregular geometries handled through coordinate transforms and multi-block meshes. Boundary conditions as well as block interface conditions (whether frictional or otherwise) are enforced weakly through the simultaneous approximation term method, resulting in a provably stable discretization.

The theoretical accuracy and stability results are confirmed with the method of manufactured solutions. The practical benefits of the new methodology are illustrated in a simulation of a subduction zone megathrust earthquake, a challenging application problem involving complex free-surface topography, nonplanar faults, and varying material properties.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2012. , 34 p.
LiTH-MAT-R, ISSN 0348-2960 ; 2012:2
Keyword [en]
High-order finite difference, nonlinear boundary conditions, simultaneous approximation term method, elastodynamics, summation-by-parts, friction, wave propagation, multi-block, coordinate transforms, weak, multi-block, coordinate transforms, weak boundary conditions
National Category
URN: urn:nbn:se:liu:diva-74594ISRN: LiTH-MAT-R--2012/02--SEOAI: diva2:488731
Available from: 2012-02-02 Created: 2012-02-02 Last updated: 2013-08-30Bibliographically approved

Open Access in DiVA

fulltext(2372 kB)207 downloads
File information
File name FULLTEXT01.pdfFile size 2372 kBChecksum SHA-512
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Nordström, Jan
By organisation
Scientific ComputingThe Institute of Technology

Search outside of DiVA

GoogleGoogle Scholar
Total: 207 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: 206 hits
ReferencesLink to record
Permanent link

Direct link