A Domain Decomposition Approach to Implementing Fault Slip in Finite-Element Models of Quasi-static and Dynamic Crustal Deformation

Abstract

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.

14 pages, 15 figures

Cited by in corpus (6)

• Three-Dimensional Numerical Modeling of Shear Stimulation of Naturally Fractured Reservoirs
• Deformation of a Half-Space from Anelastic Strain Confined in a Tetrahedral Volume
• Experiences with efficient methodologies for teaching computer programming to geoscientists
• A Reverse Augmented Constraint preconditioner for Lagrange multiplier methods in contact mechanics
• A scalable preconditioning framework for stabilized contact mechanics with hydraulically active fractures
• Gamra: Simple Meshes for Complex Earthquakes