Fault integration and projection cases

[baagaard-usgs]

Related issues: #578, #636, #646

Fault integration cases

• Lagrange multiplier in elasticity equation (prescribed slip) [Implemented]
  • FIeld: displacement
  • Inputs: Solution (Lagrange multiplier)
  • Output: Traction on positive and negative fault faces
• Lagrange multiplier in slip constraint equation (quasistatic; implcit) [Implemented]
  • FIeld: Lagrange multiplier
  • Inputs: Solution (displacements on positive and negative sides of the fault), fault auxiliary field with current prescribed slip
  • Output: Difference between current slip and prescribed slip (slip constraint)
• Lagrange multiplier in slip acceleration equation (dynamic; IMEX) [In progress]
  • Field: Lagrange multiplier
  • Inputs: Solution (displacement on negative or positive side of the fault), bulk auxiliary field with material properties
  • Output: Bulk traction and Lagrange multiplier on positive and negative fault faces
• Fault traction in elasticity equation (spontaneous rupture) [Not implemented]
  • Field: displacement
  • Inputs: Solution (displacements and velocities on positive and negative sides of the fault; Lagrange multiplier) and fault auxiliary field with fault friction parameters
  • Output: Traction on positive and negative fault faces

Fault projection cases

Output mesh is always the fault mesh.

• Output of prescribed slip [Implemented]
  • Inputs: fault auxiliary field (fault mesh)
  • Output of auxiliary subfield of prescribed slip
• Output of fault orientation [Not implemented]
  • Inputs: fault normal and reference directions
  • Output of strike direction, dip direction, and fault normal (separate kernels for each one)
• Output of fault traction [Not implemented]
  • Inputs: solution (Lagrange multiplier) (domain mesh)
  • Output: Lagrange multiplier
• Output of slip for spontaneous rupture [Not implemented]
  • Inputs: solution (displacements on positive and negative sides of the fault)
  • Output: Slip

[baagaard-usgs]

Fault integration case "Fault traction in elasticity equation (spontaneous rupture)" is now implemented in PETSc branch baagaard/update-plex-ex5-fault-projection.

[baagaard-usgs]

Getting the normal works in projection (computing the derived field) in PETSc branch baagaard/update-plex-ex5-fault-projection. See src/dm/impls/plex/tests/ex5.c.