Phase Field Fracture (#683)

* updated versions of Reaction-Diffusion solver

* ReactionDiffusion is now PhaseField. Also, created a copy of PoroElastic and reamed it PhaseField

* enable split operator algorithm in PhaseFieldFractureSolver. Add test problem

* add damage and stess updates to solid mechanics kernels. Still not quite right but moving in the right direction

* moved strain energy density calculation into constitutive relation

* enable mapping of damage to quadrature from nodal values

* enforce strain energy density is always increasing

* add solidMaterialName in the PhaseFieldDamageFem.hpp

* Adjusted parameters in PF Coupled Solver xml file

* Added load controlled xml file for PF fracture

* made call to map nodal damage to quadrature points

* add displacement controlled sample problem

* testing coupled phase-field solver after BCs were fixed, looks like there is some issue with the damage passed to the solid solver

* Rewrote Damage Solver in Residual Form

* Implemented computation of L2 norm of the residual

* Moved phasefield damage and strain energy density state variables into a Damage constitutive type

Co-authored-by: Randolph Settgast <[email protected]>


Co-authored-by: Andre Costa <[email protected]>
Co-authored-by: Andre Macieira Braga Costa <[email protected]>
Co-authored-by: Andre Macieira Braga Costa <[email protected]>
Co-authored-by: Tao Jin <[email protected]>

Author: 4 people

geosx_xml_tools_package/geosx_xml_tools/tests/generate_test_xml.py

@@ -57,7 +57,7 @@ def generate_test_xml_files(root_dir):
                      fluidName="water"
                      solidName="rock"
                      targetRegions="{Region1}">
-      <SystemSolverParameters krylovTol="1.0e-10"
+      <LinearSolverParameters krylovTol="1.0e-10"
                               newtonTol="1.0e-6"
                               maxIterNewton="8"/>
     </SinglePhaseFlow>