Development Plan#
Future implementation of features is guided by several target applications, including
Earthquake cycle modeling with quasi-static simulation of interseismic deformation and dynamic simulation of coseismic deformation.
Inversion of geodetic data for slow slip events, fault creep, and long-term fault slip rates.
Quasistatic and dynamic modeling of fluids and faulting.
Version 3.0.4 (October 2022)#
Updates for examples and documentation along with bugfixes.
Version 3.1.0 (November 2022)#
Parallel mesh loading [15%]
Diagonal Jacobian for explicit part of IMEX formulation [0%]
Improve robustness of HDF5 output by opening/closing at each time step [0%]
Additional minor cleanup of code internals to improve maintainability.
Version 3.2 (February 2022)#
Spontaneous rupture for quasistatic and dynamic simulations [20%]
Reimplementation of small strain formulation for elasticity [20%]
Reimplementation of Drucker-Prager elastoplastic bulk rheology [0%]
Add support for GeoModelGrids implementation of spatial databases for 3D seismic velocity models.
Line/point fluid sources in poroelasticity [20%]
Integration with libCEED for fast high order residual evaluation
Contribution led by Jed Brown.Add ability to output residual field during nonlinear solve for debugging [0%]
Version 3.3 (June 2023)#
Coupling of problems with compatible meshes [10%]
Implement “injectors” for solution and state variables.Migrate examples to Jupyter notebooks
Update to current version of Pyre
Moment tensor point sources [5%]
Moment tensor point sources provide a mesh independent deformation source that is better suited for Green’s function calculations than slip on a fault surface via cohesive cells.Adaptive mesh refinement
Features for Future Releases#
Consolidate HDF5 output into a single file
Elasticity with self-gravitation
Drucker-Prager bulk rheology with relaxation to yield surface
Drucker-Prager bulk rheology with strain hardening/softening
Adjoint for data assimilation
Fault with both prescribed slip and spontaneous rupture
Use fault constitutive model to control slip on fault except during episodes of prescribed slip. Need some way to describe when to turn on/off prescribed slip.