fix: update TFrac tutorial example (#3595)

Author: jhuang2601

inputFiles/lagrangianContactMechanics/ContactMechanics_TFrac_benchmark.xml

@@ -16,16 +16,17 @@
    discretization="FE1"
    targetRegions="{ Region, Fracture }">
    <NonlinearSolverParameters
-     newtonTol="1.0e-8"
+     newtonTol="1.0e-5"
      logLevel="2"
      maxNumConfigurationAttempts="10"
      newtonMaxIter="10"
      lineSearchAction="Require"
      lineSearchMaxCuts="2"
      maxTimeStepCuts="2"/>
    <LinearSolverParameters
-     solverType="direct"
-     directParallel="0"
+     solverType="gmres"
+     preconditionerType="mgr"
+     krylovTol="1e-8"
      logLevel="0"/>
   </SolidMechanicsLagrangeContact>
   <!-- SPHINX_SOLVER_END -->
@@ -39,15 +40,15 @@
   </NumericalMethods>
 
   <Events
-    maxTime="0.2">
+    maxTime="1.0">
     <SoloEvent
       name="preFracture"
       target="/Solvers/SurfaceGen"/>
 
     <PeriodicEvent
       name="solverApplications"
       beginTime="0.0"
-      endTime="0.2"
+      endTime="1.0"
       forceDt="0.2"
       target="/Solvers/lagrangiancontact"/>  
 
@@ -77,12 +78,6 @@
       name="displacementHistoryOutput"
       timeFrequency="0.2"
       targetExactTimestep="0"
-      target="/Outputs/displacementOutput"/>
-      
-    <PeriodicEvent
-      name="restarts"
-      timeFrequency="0.2"
-      targetExactTimestep="0"
-      target="/Outputs/restartOutput"/>     
+      target="/Outputs/displacementOutput"/>    
   </Events>
 </Problem>

src/docs/sphinx/advancedExamples/validationStudies/faultMechanics/intersectFrac/Example.rst

@@ -170,7 +170,7 @@ The next step is to specify fields, including:
   - The boundary conditions (traction loaded on the vertical fracture and the constraints of the outer boundaries have to be set).
 
 In this tutorial, we specify an uniaxial vertical stress ``SigmaY`` (:math:`\sigma_y` = -1.0e8 Pa). 
-A compressive traction ``NormalTraction`` (:math:`P_in` = -1.0e8 Pa) is applied at the surface of vertical fracture.
+A compressive traction ``NormalTraction`` (:math:`P_{in}` = -1.0e8 Pa) is applied at the surface of vertical fracture.
 The remaining parts of the outer boundaries are subjected to roller constraints.  
 These boundary conditions are set up through the ``FieldSpecifications`` section.
 

src/docs/sphinx/advancedExamples/validationStudies/faultMechanics/intersectFrac/intersectFracFigure.py

@@ -175,7 +175,7 @@ def main():
     cmap = plt.get_cmap("tab10")
 
     ax[0, 0].plot(r1 - length2, tn_literature, color=cmap(-1), label='Phan et al.(2003)', lw=lw)
-    ax[0, 0].plot(xlist, tnlist, 'o', alpha=0.6, color=cmap(2), mec='k', label='GEOSX Results', markersize=msize)
+    ax[0, 0].plot(xlist, tnlist, 'o', alpha=0.6, color=cmap(2), mec='k', label='GEOS Results', markersize=msize)
     ax[0, 0].grid()
     ax[0, 0].set_xlabel('Horizontal Frac Length [m]', size=fsize, weight="bold")
     ax[0, 0].set_ylabel('Normal Traction [MPa]', size=fsize, weight="bold")
@@ -184,7 +184,7 @@ def main():
     ax[0, 0].yaxis.set_tick_params(labelsize=fsize)
 
     ax[0, 1].plot(r2 - length2, gt_literature, color=cmap(-1), label='Phan et al.(2003)', lw=lw)
-    ax[0, 1].plot(xlist, gtlist, 'o', alpha=0.6, color=cmap(2), mec='k', label='GEOSX Results', markersize=msize)
+    ax[0, 1].plot(xlist, gtlist, 'o', alpha=0.6, color=cmap(2), mec='k', label='GEOS Results', markersize=msize)
     ax[0, 1].grid()
     ax[0, 1].set_xlabel('Horizontal Frac Length [m]', size=fsize, weight="bold")
     ax[0, 1].set_ylabel('Slip [mm]', size=fsize, weight="bold")
@@ -200,9 +200,8 @@ def main():
                   alpha=0.6,
                   color=cmap(2),
                   mec='k',
-                  label='GEOSX Results',
+                  label='GEOS Results',
                   markersize=msize)
-    ax[1, 0].plot(x_analytical, aperture_analytical * 1.0e3, '--', color=cmap(1), label='Sneddon Solution', lw=lw)
     ax[1, 0].grid()
     ax[1, 0].set_xlabel('Vertical Frac Length [m]', size=fsize, weight="bold")
     ax[1, 0].set_ylabel('Aperture [mm]', size=fsize, weight="bold")