Resolve dynamic + fault IMEX issues - Githubissues

geodynamics / pylith

PyLith is a finite element code for the solution of dynamic and quasi-static tectonic deformation problems.

Other

152 stars 96 forks source link

Resolve dynamic + fault IMEX issues #203

Open baagaard-usgs opened 4 years ago

baagaard-usgs commented 4 years ago

Resolve implementation of dynamic + fault IMEX issues.

Need consistent Jacobian and residuals for IMEX formulation.
Implement matrix-free Jacobian.

baagaard-usgs commented 3 years ago

Current branch: baagaard/feature-prescribed-slip-dynamic

Implemented and tested

Slip, slip rate, and slip acceleration time history functions.
FaultCohesive::_createIntegrationPatches() Create integration patches for pair of material ids on negative and positive sides of the fault. Used in dynamic case.

Implemented, not tested

Pointwise function for slip acceleration.
Pointwise functions for negative and positive sides of the fault.
Pointwise functions for negative and positive sides of the fault in IsotropicLinearElasticity

Not implemented

Compute LHS Jacobian for displacement and velocity fields (J_ii=stshift)
Set pointwise functions for negative and positive sides of the fault in other elastic rheologies materials