A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation
Reads0
Chats0
TLDR
In this article, a domain decomposition approach with Lagrange multipliers is employed to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications.Abstract:
[1] 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.read more
Citations
More filters
Journal ArticleDOI
Coupled multiphase flow and poromechanics: A computational model of pore pressure effects on fault slip and earthquake triggering
Birendra Jha,Ruben Juanes +1 more
TL;DR: In this article, the authors presented a new computational approach to model coupled multiphase flow and geomechanics of faulted reservoirs, where faults were modeled as surfaces embedded in a three-dimensional medium by using zero-thickness interface elements to accurately model fault slip under dynamically evolving fluid pressure and fault strength.
Journal ArticleDOI
A Suite of Exercises for Verifying Dynamic Earthquake Rupture Codes
Ruth A. Harris,Michael Barall,Brad T. Aagaard,Shuo Ma,D. Roten,Kim B. Olsen,Benchun Duan,Dunyu Liu,Bin Luo,Kangchen Bai,Jean-Paul Ampuero,Yoshihiro Kaneko,Alice-Agnes Gabriel,Kenneth Duru,Thomas Ulrich,Stephanie Wollherr,Zheqiang Shi,Eric M. Dunham,S. A. Bydlon,Zhenguo Zhang,Xiaofei Chen,Surendra Nadh Somala,Christian Pelties,Josué Tago,Víctor M. Cruz-Atienza,Jeremy E. Kozdon,Eric G. Daub,Khurram S. Aslam,Yuko Kase,K. Withers,Luis A. Dalguer +30 more
TL;DR: In this paper, the authors describe a set of benchmark exercises that are designed to test if computer codes that simulate dynamic earthquake rupture are working as intended, and they produce simulation results that include earthquake size, amounts of fault slip, and the patterns of ground shaking and crustal deformation.
Journal ArticleDOI
Green’s functions for geophysics: a review
TL;DR: This article provides a relatively complete review on GFs for geophysics, which can be scaled and applied to large- scale problems such as those involved in Earth sciences as well as to nano-scale problems associated with quantum nanostructures.
Journal ArticleDOI
Revisiting viscoelastic effects on interseismic deformation and locking degree: A case study of the Peru‐North Chile subduction zone
TL;DR: In this paper, the authors use finite element method (FEM) models to investigate the control of viscoelasticity on interseismic deformation and to highlight the pitfalls of interpreting the data with purely elastic models for both the forward and inverse problems.
Journal ArticleDOI
Earthquakes and Tremor Linked to Seamount Subduction During Shallow Slow Slip at the Hikurangi Margin, New Zealand
Erin Todd,Erin Todd,Susan Y. Schwartz,Kimihiro Mochizuki,Laura M. Wallace,Laura M. Wallace,Anne F. Sheehan,Spahr C. Webb,Charles Williams,J. S. Nakai,J. Yarce,Bill Fry,Stuart Henrys,Yoshihiro Ito +13 more
References
More filters
Journal ArticleDOI
Matplotlib: A 2D Graphics Environment
TL;DR: Matplotlib is a 2D graphics package used for Python for application development, interactive scripting, and publication-quality image generation across user interfaces and operating systems.
Book
The finite element method
TL;DR: In this article, the methodes are numeriques and the fonction de forme reference record created on 2005-11-18, modified on 2016-08-08.
Book
Iterative Methods for Sparse Linear Systems
TL;DR: This chapter discusses methods related to the normal equations of linear algebra, and some of the techniques used in this chapter were derived from previous chapters of this book.
Journal ArticleDOI
Tectonic stress and the spectra of seismic shear waves from earthquakes
TL;DR: In this paper, an earthquake model is derived by considering the effective stress available to accelerate the sides of the fault, and the model describes near and far-field displacement-time functions and spectra and includes the effect of fractional stress drop.
Journal ArticleDOI
A Method of Computation for Structural Dynamics
TL;DR: In this paper, the method is capable of application to structures of any degree of complication, with any relationship between force and displacement, from linear elastic behavior through various degrees of inelastic behavior or plastic response, up to failure; any type of dynamic loading, due to shock or impact, vibration, earthquake, or nuclear blast can be considered.