You are here: Home / Groups / Short-Term Crustal Dynamics / Wiki / Work Plans / 2016 PyLith Dev Plans
  • Discoverability Visible
  • Join Policy Invite Only
  • Created 05 Jan 2021

Work Plans /

2016 PyLith Dev Plans

PyLith Development Plans, Jun 2016

Priorities for PyLith software development, such as new features and enhancements. This a draft for community comment (Jun 24, 2016).

This plan attempts to balance meeting short-term objectives of delivering high priority, new features and meeting long-term objectives of extending the code to solve a broader range of scientific problems.

Current development is slower than in the past because Brad Aagaard (who is not funded by CIG) is unable to devote as much time to PyLith development due to a substantial increase of other obligations.

Version 3.0 (Spring 2017)

  1. Multiphysics expert.png 15%
    • Setup modular approach for specifying governing equations and computing residuals and Jacobians. expert.png
    • Incompressible elasticity via a pressure field intermediate.png
    • Elasticity + heat flow difficult.png
    • Poroelasticity difficult.png
  2. Higher order basis functions difficult.png 20%
    • Allow user to select order of basis functions independent of the mesh (which defines the geometry). This permits higher resolution for a given mesh.
  3. Switch to using PETSc time-stepping (TS) algorithms. intermediate.png 25%
    • Replace simple Python-based time-stepping implementations with PETSc time-stepping algorithms that provide support for higher order discretization in time and real adaptive time stepping.
  4. Improve fault formulation for spontaneous rupture intermediate.png 10%
    • Removes inner solve associated with updating Lagrange multipliers. This will significantly accelerate the nonlinear solve.
  5. Add Drucker-Prager with relaxation to yield surface bulk rheology intermediate.png 0%

Version 3.1 (late 2017)

  1. Reorganize top-level code to allow different startup cases intermediate.png 0%
    • Elastic prestep
    • User-specified initial solution
    • Checkpoint via special spatial database?
    • Reorganization for time-dependent Green’s functions and adjoints
  2. Multilevel nonlinear solve
  3. Radial basis functions for spatial databases intermediate.png 0%

Version 4.0 (TBD)

  1. Earthquake cycle modeling difficult.png
    • Same mesh for dynamic and quasi-static parts (dynamic -> quasi-static, quasi-static -> dynamic, complete cycle)
  2. Create strain hardening/softening 2-D and 3-D Drucker-Prager elastoplastic models.intermediate.png
  3. Moment tensor point sources via equivalent body forces difficult.png 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.

Features for Future Releases

Major features

  1. Earthquake Cycle Modeling
    • Different meshes for dynamic and quasi-static parts expert.png
      • Requires interpolation of fields between different meshes/discretizations and may require extrapolation of solutions when quasi-static problems span a larger domain than the dynamic problems.
  2. Data assimilation
    • Use flexibility of multiphysics implementation to support inclusion of data assimilation expert.png

Minor features

  1. GUI interface for specifying parameters difficult.png 25%
  2. Begin implementation of data assimilation capabilities via adjoint equation.
  3. Combined prescribed slip / spontaneous rupture fault condition difficult.png
    • 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.
  4. Use threading to accelerate integrations on multi-core machines. difficult.png

difficulty rating system

intermediate.png intermediate difficult.png difficult expert.png expert

based on the ski trail rating system

Created on , Last modified on