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.
Note
Because we strictly follow the semantic versioning guidelines, a minor release may get promoted to a major releases if we make changes to the public API (parameters). This can happen if realize that we should modify the parameters to improve maintainability or prepare for future changes.
Version 4.0.0 (December 2023)#
Analytical function spatial database [100%]
Use
pod
initial guess to improve convergence [100%]Output of fault tractions [100%]
Output of boundary and fault orientation [100%]
Convert from CppUnit to Catch2 [100%]
Version 4.1.0 (March 2024)#
Output of fault rupture auxiliary subfields [0%]
Dynamic prescribed slip with diagonal Jacobian for explicit part of IMEX formulation [75%]
Parallel mesh loading [15%]
Add 2D and 3D examples for crustal faults with complex fault geometry [50%]
Better preconditioners [25%]
elasticity with fault
incompressible elasticity
poroelasticity
Update VTK output to use
vtu
files rather than legacyvtk
files [0%]Finish updating
examples/subduction-3d
[40%]
Version 5.0.0 (June 2024)#
Spontaneous rupture for quasistatic and dynamic simulations [20%]
Convert from SWIG to pybind11 [0%]
Reimplementation of Drucker-Prager elastoplastic bulk rheology [0%]
Add support for GeoModelGrids implementation of spatial databases for 3D seismic velocity models. [0%]
Improve robustness of HDF5 output by opening/closing at each time step [0%]
Add
examples/barwaves-2d
[15%]Update coordinates with solution [50%]
Version 6.0.0 (TBD)#
Migrate examples to Jupyter notebooks
Update to current version of Pyre
More flexible specification of time-dependent boundary conditions. [0%]
Dirichlet boundary conditions with constraints on normal and tangential components. [0%]
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%]
Elasticity with self-gravitation [0%]
Features for Future Releases#
Coupling of problems with compatible meshes [10%]
Implement “injectors” for solution and state variables.Reimplementation of small strain formulation for elasticity [20%]
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
Line/point fluid sources in poroelasticity [20%]
Consolidate HDF5 output into a single file
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.