refactor: Phase-field and damage solid model cleanup, refactor (#4096)

* PhaseFieldFractureSolver, PhaseFieldPoromechanicsSolver, and PhaseFieldDamageFEM cleanup
* added a check for nucleation model to see if the dissipation option is linear
* move localDissipationoption to constitutive model and added FractureModelType option
* updated xmls
* renamed sed to crackDrivingForce
* apply Damage default field values after array allocation
* refactored and added spectral/vol-dev split test decks for all fracture-model types
* fix save history of crack driving force
* enabled multi-region run with certain defined regions have no damage dof
* register nodal damage as a field trait instead of an XML fieldName
* Add phase-field benchmark
* add viscous regularization
* removed fieldName in benchmarks, updated schemas
* use internal meshes for single notch examples and stop tracking lfs
* removed useless coeffField and call setConstitutiveNamesCallSuper
* Error out on unimplemented ExplicitTransient time integration
* moved fine mesh files to GEOSDATA, created coarse meshes for integrated tests
* fix a bug to save the crack driving force history also on the host
* rebaseline
* make the fracture model type and dissipation option clear in examples
* merged damage interpolation kernels

Author: frankfeifan

.integrated_tests.yaml

@@ -1,6 +1,6 @@
 baselines:
   bucket: geosx
-  baseline: integratedTests/baseline_integratedTests-pr4127-17399-63445db
+  baseline: integratedTests/baseline_integratedTests-pr4096-17439-18734f3
 
 allow_fail:
   all: ''

BASELINE_NOTES.md

@@ -5,6 +5,9 @@ This file is designed to track changes to the integrated test baselines.
 Any developer who updates the baseline ID in the .integrated_tests.yaml file is expected to create an entry in this file with the pull request number, date, and their justification for rebaselining.
 These notes should be in reverse-chronological order, and use the following time format: (YYYY-MM-DD).
 
+PR #4096 (2026-08-27) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4096-17439-18734f3.tar.gz>
+Refactor the phase-field module and add new benchmarks for baseline
+
 PR #4127 (2026-08-23) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4127-17399-63445db.tar.gz>
 =====================
 Rebaseline five restart checks after the TPL update changed VTK/Scotch mesh partitioning. Global mesh topology and fields are unchanged.

inputFiles/phaseField/PhaseFieldFracture.ats

@@ -28,21 +28,45 @@ decks = [
         check_step=10,
         restartcheck_params=RestartcheckParameters(**restartcheck_params)),
     TestDeck(
-        name="PhaseFieldFracture_SpectralSplit",
+        name="PhaseFieldFracture_SpectralSplit_Cohesive",
         description=
         'Testing the spectral split of the stress tensor in PF Fracture',
         partitions=((1, 1, 1), (2, 2, 1)),
         restart_step=15,
         check_step=30,
         restartcheck_params=RestartcheckParameters(**restartcheck_params)),
     TestDeck(
-        name="PhaseFieldFracture_VolDevSplit",
+        name="PhaseFieldFracture_SpectralSplit_Brittle",
+        description=
+        'Testing the spectral split with the brittle degradation function in PF Fracture',
+        partitions=((1, 1, 1), (2, 2, 1)),
+        restart_step=15,
+        check_step=30,
+        restartcheck_params=RestartcheckParameters(**restartcheck_params)),
+    TestDeck(
+        name="PhaseFieldFracture_VolDevSplit_Brittle",
         description=
         'Testing the volumetric-deviatoric split of the stress tensor in PF Fracture',
         partitions=((1, 1, 1), (2, 2, 1)),
         restart_step=15,
         check_step=30,
         restartcheck_params=RestartcheckParameters(**restartcheck_params)),
+    TestDeck(
+        name="PhaseFieldFracture_VolDevSplit_Cohesive",
+        description=
+        'Testing the volumetric-deviatoric split with the cohesive degradation function in PF Fracture',
+        partitions=((1, 1, 1), (2, 2, 1)),
+        restart_step=15,
+        check_step=30,
+        restartcheck_params=RestartcheckParameters(**restartcheck_params)),
+    TestDeck(
+        name="PhaseFieldFracture_MultiRegion",
+        description=
+        'Testing PF Fracture on one region of a multi-region problem, with mechanics only elsewhere',
+        partitions=((1, 1, 1), (2, 2, 1)),
+        restart_step=10,
+        check_step=20,
+        restartcheck_params=RestartcheckParameters(**restartcheck_params)),
     TestDeck(name="PhaseFieldFracture_Nucleation_smoke",
              description='Testing the nucleation model for PF Fracture',
              partitions=((1, 1, 1), (2, 2, 1)),

inputFiles/phaseField/PhaseFieldFracture_CohesiveModel.xml

@@ -31,8 +31,6 @@
       name="DamageSolve"
       discretization="FE1"
       timeIntegrationOption="SteadyState"
-      fieldName="Damage"
-      localDissipation="Linear"
       logLevel="1"
       targetRegions="{ Region1 }">
       <NonlinearSolverParameters
@@ -110,7 +108,9 @@
       defaultShearModulus="8.3333e9"
       lengthScale="0.2"
       defaultCriticalFractureEnergy="2.7"
-      criticalStrainEnergy="1.0"/>
+      criticalStrainEnergy="1.0"
+      fractureModelType="Cohesive"
+      localDissipationOption="Linear"/>
   </Constitutive>
 
   <FieldSpecifications>
@@ -119,14 +119,6 @@
                         objectPath="nodeManager"
                         scale="0"
                         setNames="{xneg, xpos, yneg, ypos, zneg, zpos}"/> -->
-    <FieldSpecification
-       name="Gc"
-       initialCondition="1"
-       setNames="{ all }"
-       objectPath="ElementRegions/Region1"
-       fieldName="shale_criticalFractureEnergy"
-       scale="2.7"/>
-
     <FieldSpecification
       name="xpos"
       objectPath="nodeManager"
@@ -167,25 +159,17 @@
       setNames="{ all }"/>
 
     <!-- <FieldSpecification name="pointSource"
-                        fieldName="Damage"
+                        fieldName="damage"
                         initialCondition="0"
                         objectPath="nodeManager"
                         scale="1"
                         setNames="{centerLine}"/> -->
     <!-- FieldSpecification name="pointSource2"
-                        fieldName="Damage"
+                        fieldName="damage"
                         initialCondition="0"
                         objectPath="nodeManager"
                         scale="0"
                         setNames="{xneg,xpos}"/> -->
-    <!-- <FieldSpecification name="diffusionCoeff"
-                        initialCondition="1"
-                        setNames="{all}"
-                        objectPath="ElementRegions"
-                        fieldName="coeffField"
-                        scale="2e3"
-                        functionName="f_b"
-                        /> -->
   </FieldSpecifications>
 
   <Functions>

inputFiles/phaseField/PhaseFieldFracture_DamageAndLoad.xml

@@ -34,9 +34,6 @@
       name="DamageSolve"
       discretization="FE1"
       timeIntegrationOption="SteadyState"
-      fieldName="Damage"
-      
-      localDissipation="Quadratic"
       logLevel="1"
       targetRegions="{ Region1 }">
       <NonlinearSolverParameters
@@ -114,18 +111,12 @@
       defaultShearModulus="3.00e5"
       lengthScale="0.2"
       defaultCriticalFractureEnergy="2.7"
-      criticalStrainEnergy="0"/>
+      criticalStrainEnergy="0"
+      fractureModelType="Brittle"
+      localDissipationOption="Quadratic"/>
   </Constitutive>
 
   <FieldSpecifications>
-    <FieldSpecification
-       name="Gc"
-       initialCondition="1"
-       setNames="{ all }"
-       objectPath="ElementRegions/Region1"
-       fieldName="shale_criticalFractureEnergy"
-       scale="2.7"/>
-
     <FieldSpecification
       name="xtraction"
       fieldName="totalDisplacement"
@@ -158,7 +149,7 @@
 
     <FieldSpecification
       name="pointSource"
-      fieldName="Damage"
+      fieldName="damage"
       objectPath="nodeManager"
       scale="1"
       setNames="{ middle }"/>

inputFiles/phaseField/PhaseFieldFracture_DamageOnly.xml

@@ -34,9 +34,6 @@
       name="DamageSolve"
       discretization="FE1"
       timeIntegrationOption="SteadyState"
-      fieldName="Damage"
-      
-      localDissipation="Linear"
       logLevel="1"
       targetRegions="{ Region1 }">
       <NonlinearSolverParameters
@@ -116,18 +113,12 @@
       defaultShearModulus="8.0769e4"
       lengthScale="0.2"
       defaultCriticalFractureEnergy="2.7"
-      criticalStrainEnergy="1.5"/>
+      criticalStrainEnergy="1.5"
+      fractureModelType="Brittle"
+      localDissipationOption="Linear"/>
   </Constitutive>
 
   <FieldSpecifications>
-    <FieldSpecification
-       name="Gc"
-       initialCondition="1"
-       setNames="{ all }"
-       objectPath="ElementRegions/Region1"
-       fieldName="shale_criticalFractureEnergy"
-       scale="2.7"/>
-
     <FieldSpecification
       name="xboundary"
       fieldName="totalDisplacement"
@@ -151,7 +142,7 @@
 
     <FieldSpecification
       name="pointSource"
-      fieldName="Damage"
+      fieldName="damage"
       objectPath="nodeManager"
       scale="1"
       setNames="{ centerPoint }"/>

…Field/PhaseFieldFracture_VolDevSplit.xml …ieldFracture_DecompositionSplit_base.xmlinputFiles/phaseField/PhaseFieldFracture_VolDevSplit.xml renamed to inputFiles/phaseField/PhaseFieldFracture_DecompositionSplit_base.xml inputFiles/phaseField/PhaseFieldFracture_VolDevSplit.xml renamed to inputFiles/phaseField/PhaseFieldFracture_DecompositionSplit_base.xml

@@ -7,11 +7,11 @@
       name="PhaseFieldSolve"
       solidSolverName="SolidMechSolve"
       damageSolverName="DamageSolve"
-      logLevel="2"
+      logLevel="1"
       targetRegions="{ Region1 }">
       <NonlinearSolverParameters
         couplingType="Sequential"
-	lineSearchAction="None"  	  
+        lineSearchAction="None"
         newtonMaxIter="100"/>
     </PhaseFieldFracture>
 
@@ -32,9 +32,6 @@
       name="DamageSolve"
       discretization="FE1"
       timeIntegrationOption="SteadyState"
-      fieldName="Damage"
-      
-      localDissipation="Quadratic"
       logLevel="1"
       targetRegions="{ Region1 }">
       <NonlinearSolverParameters
@@ -73,37 +70,6 @@
       xMax="{  1.01, 1.01,  0.11 }"/>
   </Geometry>
 
-  <Events
-    maxTime="30.0">
-    <!-- <SoloEvent name="siloOutput0"
-               target="/Outputs/SiloOutput"/> -->
-    <!-- This event is applied every cycle, and overrides the
-    solver time-step request -->
-    <PeriodicEvent
-      name="solverApplications"
-      forceDt="1.0"
-      target="/Solvers/PhaseFieldSolve"/>
-
-    <!-- This event is applied every 1.0s.  The targetExactTimestep
-    flag allows this event to request a dt modification to match an
-    integer multiple of the timeFrequency. -->
-    <PeriodicEvent
-      name="vtkOutput"
-      timeFrequency="1.0"
-      target="/Outputs/Mode2VolDev_VTKOutput"/>
-
-    <PeriodicEvent
-      name="siloOutputs"
-      timeFrequency="1.0"
-      target="/Outputs/Mode2VolDev_SiloOutput"/>
-
-    <PeriodicEvent
-      name="restarts"
-      timeFrequency="15.0"
-      targetExactTimestep="0"
-      target="/Outputs/restartOutput"/>
-  </Events>
-
   <NumericalMethods>
     <FiniteElements>
       <FiniteElementSpace
@@ -119,26 +85,7 @@
       materialList="{ shale }"/>
   </ElementRegions>
 
-  <Constitutive>
-    <DamageVolDevElasticIsotropic
-      name="shale"
-      defaultDensity="2700"
-      defaultBulkModulus="1.750e5"
-      defaultShearModulus="8.077e4"
-      lengthScale="0.25"
-      defaultCriticalFractureEnergy="54.0"
-      criticalStrainEnergy="14.88"/>
-  </Constitutive>
-
   <FieldSpecifications>
-    <FieldSpecification
-       name="Gc"
-       initialCondition="1"
-       setNames="{ all }"
-       objectPath="ElementRegions/Region1"
-       fieldName="shale_criticalFractureEnergy"
-       scale="54.0"/>
-    
     <FieldSpecification
       name="shearLoad"
       fieldName="totalDisplacement"
@@ -178,22 +125,4 @@
       values="{ -1, 100 }"/>
   </Functions>
 
-  <Outputs>
-    <Silo
-      name="Mode2VolDev_SiloOutput"/>
-
-    <VTK
-      name="Mode2VolDev_VTKOutput"
-      plotFileRoot="plot"/>
-
-    <Restart
-      name="restartOutput"/>
-  </Outputs>
-
-  <!-- <Geometry> -->
-  <!--   <Box name="bottom" xMin="{-0.51, -0.51, -0.01" xMax="0.51, -0.49, +0.11"/> -->
-  <!--   <Box name="top" xMin="{-0.51, 0.49, -0.01" xMax="0.51, 0.51, +0.11"/> -->
-  <!--   <Box name="left" xMin="{-0.51, -0.51, -0.01" xMax="-0.49, 0.51, +0.11"/> -->
-  <!--   <Box name="right" xMin="0.49, -0.51, -0.01" xMax="0.51, 0.51, +0.11"/> -->
-  <!-- </Geometry> -->
-</Problem>
+</Problem>