Defmod - Parallel multiphysics finite element code for modeling crustal deformation during the earthquake/rifting cycle

In this article, we present Defmod, an open source, fully unstructured, two or three dimensional, parallel finite element code for modeling crustal deformation over time scales ranging from milliseconds to thousands of years. Unlike existing public domain numerical codes, Defmod can simulate deformation due to all major processes that make up the earthquake/rifting cycle, in non-homogeneous media. Specifically, it can be used to model deformation due to dynamic and quasistatic processes such as co-seismic slip or dike intrusion(s), poroelastic rebound due to fluid flow and post-seismic or post-rifting viscoelastic relaxation. It can also be used to model deformation due to processes such as post-glacial rebound, hydrological (un)loading, injection and/or withdrawal of fluids from subsurface reservoirs etc. Defmod is written in Fortran 95 and uses PETSc's parallel sparse data structures and implicit solvers. Problems can be solved using (stabilized) linear triangular, quadrilateral, tetrahedral or hexahedral elements on shared or distributed memory machines with hundreds or even thousands of processor cores. In the current version of the code, prescribed loading is supported. Results are written in ASCII VTK format for easy visualization. The source code is released under the terms of GNU General Public License (v3.0) and is freely available from this https URL

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

[2]  William Roshan Quadros,et al.  CUBIT geometry and mesh generation toolkit 15.1 user documentation , 2016 .

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

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

[5]  R. Taylor.  The Finite Element Method, the Basis , 2000 .

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

[7]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[8]  Message Passing Interface Forum.  MPI: A message - passing interface standard , 1994 .

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

[10]  Ronaldo I. Borja,et al.  Stabilized low-order finite elements for coupled solid-deformation/fluid-diffusion and their application to fault zone transients , 2008 .

[11]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[12]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[13]  J. Z. Zhu,et al.  The finite element method , 1977 .

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