Benchmarking Defmod, an open source FEM code for modeling episodic fault rupture

We present Defmod, an open source (linear) finite element code that enables us to efficiently model the crustal deformation due to (quasi-)static and dynamic loadings, poroelastic flow, viscoelastic flow and frictional fault slip. Ali (2015) provides the original code introducing an implicit solver for (quasi-)static problem, and an explicit solver for dynamic problem. The fault constraint is implemented via Lagrange Multiplier. Meng (2015) combines these two solvers into a hybrid solver that uses failure criteria and friction laws to adaptively switch between the (quasi-)static state and dynamic state. The code is capable of modeling episodic fault rupture driven by quasi-static loadings, e.g. due to reservoir fluid withdraw or injection. Here, we focus on benchmarking the Defmod results against some establish results. The Defmod code is validated to model the static and/or dynamic crustal deformations.The code adaptively switches between the (quasi-)static and dynamic modes.The code enables one to model earthquakes induced by fluid withdraw and/or injection.

[1]  Ruben Juanes,et al.  Coupled multiphase flow and poromechanics: A computational model of pore pressure effects on fault slip and earthquake triggering , 2014 .

[2]  K. Satō,et al.  Mandel and Cryer problems of fluid-saturated foams with negative Poisson’s ratio , 2005 .

[3]  L. John,et al.  Finite dynamic model for infinite media , 1969 .

[4]  Annalisa Fregolent,et al.  Dynamic rupture at a frictional interface between dissimilar materials with asperities , 2010 .

[5]  J. Ampuero,et al.  Earthquake nucleation on rate and state faults - Aging and slip laws , 2008 .

[6]  Dhabaleswar K. Panda,et al.  Scalable Earthquake Simulation on Petascale Supercomputers , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[7]  J. C. Savage,et al.  Asthenosphere readjustment and the earthquake cycle , 1978 .

[8]  R. Taylor,et al.  Lagrange constraints for transient finite element surface contact , 1991 .

[9]  Ladislav Frýba.  History of Winkler Foundation , 1995 .

[10]  D. V. Griffiths,et al.  Programming the finite element method , 1982 .

[11]  I. Gladwell,et al.  PCG methods in transient FE analysis. Part II: Second order problems , 1989 .

[12]  Clark R. Dohrmann,et al.  A Computational Study of Stabilized, Low-order C0 Finite Element Approximations of Darcy Equations , 2006 .

[13]  Stephen P. Timoshenko,et al.  Vibration Problems in Engineering - fourth edition , 1974 .

[14]  Sylvain Barbot,et al.  A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow , 2010 .

[15]  S. Day,et al.  Comparison of Fault Representation Methods in Finite Difference Simulations of Dynamic Rupture , 2006 .

[16]  A. Pitarka,et al.  The SCEC/USGS Dynamic Earthquake Rupture Code Verification Exercise , 2012 .

[17]  Stephen P. Timoshenko,et al.  Vibration problems in engineering , 1928 .

[18]  H. J. Melosh,et al.  The dynamical origin of subduction zone topography , 1980 .

[19]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

[20]  S. Tabrez Ali,et al.  Defmod - Parallel multiphysics finite element code for modeling crustal deformation during the earthquake/rifting cycle , 2014, ArXiv.

[21]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[22]  Luis A. Dalguer,et al.  Staggered-grid split-node method for spontaneous rupture simulation , 2007 .