Journal
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH
Volume 118, Issue 6, Pages 3059-3079Publisher
AMER GEOPHYSICAL UNION
DOI: 10.1002/jgrb.50217
Keywords
fault slip; finite-element modeling; crustal deformation; earthquake physics
Categories
Funding
- Earthquake Hazards Program of the U.S. Geological Survey
- Computational Infrastructure for Geodynamics (NSF) [EAR-0949446]
- NSF [EAR-0529922, 07HQAG0008, EAR/ITR-0313238, EAR-0745391]
Ask authors/readers for more resources
We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available