feat: SCEC benchmarks BP6 and BP7

inputFiles/inducedSeismicity/SCEC_BP6_QD_S_base.xml

@@ -18,9 +18,9 @@
       name="mesh1"
       elementTypes="{ C3D8 }"
       xCoords="{ -10.0e3, 10.0e3 }"
-      yCoords="{ -0.5e3, 0.5e3 }"
-      zCoords="{ -20.0e3, -10.0e3, 10.0e3, 20.0e3 }"
-      nx="{ 100 }"
+      yCoords="{ -0.1e3, 0.1e3 }"
+      zCoords="{ -40.0e3, -10.0e3, 10.0e3, 40.0e3 }"
+      nx="{ 50 }"
       ny="{ 1 }"
       nz="{ 20, 201, 20 }"
       cellBlockNames="{ cb1 }"/>
@@ -29,13 +29,58 @@
   <Geometry>    
     <Box
       name="faultPlane"
-      xMin="{ -0.01, -1.01e3, -10.0e3 }"
-      xMax="{ 0.01,   1.01e3,  10.0e3 }"/>
+      xMin="{ -0.01, -1.01e3, -20.0e3 }"
+      xMax="{ 0.01,   1.01e3,  20.0e3 }"/>
 
+    <!-- Boxes for single-cell sources/receivers on fault -->
+
+    <!-- Source at z = 0 m -->
     <Box
       name="source"
-      xMin="{ -0.01, -1.01e3, -50.0 }"
-      xMax="{ 0.01,   1.01e3,  50.0 }"/>
+      xMin="{ -0.01, -1.01e3, -0.1e3 }"
+      xMax="{ 0.01,   1.01e3,  0.1e3 }"/>
+
+    <!-- Receiver at  z = 500 m-->
+    <Box 
+      name="receiver1"
+      xMin="{ -0.01, -1.01e3,  0.4e3 }"
+      xMax="{ 0.01,   1.01e3,  0.6e3 }"/>
+
+    <!-- Receiver at  z = 1500 m-->
+    <Box
+      name="receiver2"
+      xMin="{ -0.01, -1.01e3,  1.4e3 }"
+      xMax="{ 0.01,   1.01e3,  1.6e3 }"/>
+
+    <!-- Receiver at  z = 2500 m-->
+    <Box
+      name="receiver3"
+      xMin="{ -0.01, -1.01e3,  2.4e3 }"
+      xMax="{ 0.01,   1.01e3,  2.6e3 }"/>
+
+    <!-- Receiver at  z = 3500 m-->
+    <Box
+      name="receiver4"
+      xMin="{ -0.01, -1.01e3,  3.4e3 }"
+      xMax="{ 0.01,   1.01e3,  3.6e3 }"/>
+
+    <!-- Receiver at  z = 5000 m-->
+    <Box
+      name="receiver5"
+      xMin="{ -0.01, -1.01e3,  4.9e3 }"
+      xMax="{ 0.01,   1.01e3,  5.1e3 }"/> 
+
+    <!-- Receiver at  z = 7500 m-->
+    <Box
+      name="receiver6"
+      xMin="{ -0.01, -1.01e3,  7.4e3 }"
+      xMax="{ 0.01,   1.01e3,  7.6e3 }"/>
+
+    <!-- Receiver at  z = -1500 m-->
+    <Box
+      name="receiver7"
+      xMin="{ -0.01, -1.01e3,  -1.6e3 }"
+      xMax="{ 0.01,   1.01e3,  -1.4e3 }"/>     
   </Geometry>
 
   <ElementRegions>
@@ -47,7 +92,7 @@
     <SurfaceElementRegion
       name="Fault"
       materialList="{frictionLaw, faultMaterial, water}"
-      defaultAperture="1e-3"/>  
+      defaultAperture="5e-3"/>  
     </ElementRegions>
 
     <Constitutive>
@@ -83,7 +128,7 @@
      name="rock"
      defaultDensity="2670"
      defaultBulkModulus="50.0e9"
-     defaultShearModulus="32.04e9"/>
+     defaultShearModulus="24.03e9"/> <!-- Note: Adjusted to translate antiplane shear to plane strain  -->
 
     <RateAndStateFriction
       name="frictionLaw"
@@ -110,13 +155,13 @@
       name="slipRateCollection"
       objectPath="ElementRegions/Fault/FractureSubRegion"
       fieldName="slipRate"
-      setNames="{source}"/>
+      setNames="{ receiver1, receiver2, receiver3, receiver4, receiver5, receiver6, receiver7, source}"/>
 
     <PackCollection
       name="pressureCollection"
       objectPath="ElementRegions/Fault/FractureSubRegion"
       fieldName="pressure"
-      setNames="{source}"/>    
+      setNames="{ receiver1, receiver2, receiver3, receiver4, receiver5, receiver6, receiver7, source}"/> 
    </Tasks>
 
   <ElementRegions>
@@ -169,7 +214,7 @@
       objectPath="ElementRegions/Fault"
       fieldName="stateVariable"
       scale="0.68070857"
-      setNames="{ all }"
+      setNames="{ faultPlane }"
       initialCondition="1"/>  
 
     <FieldSpecification
@@ -179,7 +224,7 @@
       objectPath="ElementRegions/Fault/FractureSubRegion"
       component="0"
       scale="1.0e-12"
-      setNames="{all}"/>
+      setNames="{ faultPlane }"/>
 
     <FieldSpecification
       name="slipVelocity2"
@@ -188,7 +233,7 @@
       objectPath="ElementRegions/Fault/FractureSubRegion"
       component="1"
       scale="0."
-      setNames="{all}"/>
+      setNames="{ faultPlane }"/>
 
     <FieldSpecification
       name="backgroundShearStress"
@@ -197,21 +242,21 @@
       objectPath="ElementRegions/Fault/FractureSubRegion"
       component="0"
       scale="29.2e6"
-      setNames="{all}"/> 
+      setNames="{ faultPlane }"/> 
 
     <FieldSpecification
       name="backgroundNormalStress"
       fieldName="backgroundNormalStress"
       initialCondition="1"
       objectPath="ElementRegions/Fault/FractureSubRegion"
       scale="50.0e6"
-      setNames="{all}"/>                      
+      setNames="{ faultPlane }"/>                      
 
     <SourceFlux
       name="sourceTerm"
       objectPath="ElementRegions/Fault"
       scale="-1.25e-3"
-      setNames="{ source }"
+      setNames="{source}"
       endTime="8.64e6"/>
   </FieldSpecifications>
 </Problem>

inputFiles/inducedSeismicity/SCEC_BP6_QD_S_explicit.xml

@@ -56,8 +56,8 @@
         lineSearchMaxCuts="2"
         maxTimeStepCuts="1"/>
       <LinearSolverParameters
-        solverType="direct"
-        directParallel="0"
+        solverType="gmres"
+        preconditionerType="mgr"
         logLevel="0"/>
     </SolidMechanicsLagrangeContactBubbleStab>
 
@@ -77,14 +77,16 @@
         newtonTol = "1.0e-6"
         newtonMaxIter ="8"/>
       <LinearSolverParameters
-        solverType="direct"
-        directParallel="0"/>   
+        solverType="gmres"
+        preconditionerType="amg"
+        logLevel="0"/>
     </SinglePhaseFVM>
   </Solvers>
 
-
+  <!-- <Events
+    maxTime="6.312e7"> -->
   <Events
-    maxTime="6.312e7">
+    maxTime="1e6">
 
    <!-- fault generation  -->
    <SoloEvent 
@@ -100,21 +102,9 @@
 
   <!-- solver application  -->
    <PeriodicEvent
-     name="solverApplications1"
+     name="solverApplications"
      maxEventDt="1e4"
-     endTime="8.64e6"
      target="/Solvers/QDSolver"/> 
-   <PeriodicEvent
-     name="solverApplications2"
-     maxEventDt="10"
-     beginTime="8.64e6"
-     endTime="8.6499e6"
-     target="/Solvers/QDSolver"/>
-   <PeriodicEvent
-     name="solverApplications3"
-     maxEventDt="1e3"
-     beginTime="8.6499e6"
-     target="/Solvers/QDSolver"/>
 
   <!-- Time history collection -->
     <PeriodicEvent

inputFiles/inducedSeismicity/SCEC_BP6_QD_S_implicit.xml

@@ -57,8 +57,8 @@
         lineSearchMaxCuts="2"
         maxTimeStepCuts="1"/>
       <LinearSolverParameters
-        solverType="direct"
-        directParallel="0"
+        solverType="gmres"
+        preconditionerType="mgr"
         logLevel="0"/>
     </SolidMechanicsLagrangeContactBubbleStab>
 
@@ -78,14 +78,15 @@
         newtonTol = "1.0e-6"
         newtonMaxIter ="8"/>
       <LinearSolverParameters
-        solverType="direct"
-        directParallel="0"/>   
+        solverType="gmres"
+        preconditionerType="amg"
+        logLevel="0"/>  
     </SinglePhaseFVM>
   </Solvers>
 
 
   <Events
-    maxTime="6.312e7">
+    maxTime="1e6">
 
    <!-- fault generation  -->
    <SoloEvent 

inputFiles/inducedSeismicity/SCEC_BP7_QD_base.xml

@@ -1,54 +1,5 @@
 <?xml version="1.0" ?>
 <Problem>
-  <Solvers
-    gravityVector="{0.0, 0.0, 0.0}">
-    <!-- <ImplicitQuasiDynamicEQ
-      name="QDSolver"
-      targetRegions="{ Domain, Fault }"
-      shearImpedance="4.6247113164e6"
-      initialDt="1000"
-      logLevel="1"
-      discretization="FE1"
-      stressSolverName="stressSolver"
-      targetSlipIncrement="1.0e-6">
-      <NonlinearSolverParameters
-        newtonTol="1.0e-8"
-        logLevel="2"
-        newtonMaxIter="20"
-        maxNumConfigurationAttempts="1"
-        lineSearchAction="Require"
-        lineSearchMaxCuts="2"
-        maxTimeStepCuts="1"/>
-     </QuasiDynamicEQ> -->
-
-    <SolidMechanicsLagrangeContactBubbleStab
-      name="stressSolver"
-      timeIntegrationOption="QuasiStatic"
-      logLevel="2"
-      writeLinearSystem="0"
-      discretization="FE1"
-      targetRegions="{ Domain, Fault }">
-      <NonlinearSolverParameters
-        newtonTol="1.0e-8"
-        logLevel="2"
-        newtonMaxIter="10"
-        maxNumConfigurationAttempts="1"
-        lineSearchAction="Require"
-        lineSearchMaxCuts="2"
-        maxTimeStepCuts="1"/>
-      <LinearSolverParameters
-        solverType="direct"
-        directParallel="0"
-        logLevel="0"/>
-    </SolidMechanicsLagrangeContactBubbleStab>
-
-    <SurfaceGenerator
-      name="SurfaceGen"
-      targetRegions="{ Domain }"
-      initialRockToughness="1.0"
-      mpiCommOrder="1"
-      fractureRegion="Fault"/>  
-  </Solvers>
 
   <NumericalMethods>
     <FiniteElements>
@@ -222,30 +173,4 @@
       fieldName="pressure"
       setNames="{source}"/>    
    </Tasks> -->
-
-  <Events
-    maxTime="4.5e4"
-    maxCycle="1">
-
-    <SoloEvent 
-      name="generateFault"
-      target="/Solvers/SurfaceGen"/>
-
-    <PeriodicEvent
-      name="vtkOutput"
-      cycleFrequency="1"
-      targetExactTimestep="0"
-      target="/Outputs/vtkOutput"/>  
-
-    <PeriodicEvent
-      name="solverApplications"
-      maxEventDt="1e4"
-      target="/Solvers/stressSolver"/>
-
-    <!-- <PeriodicEvent
-      name="resarts"
-      timeFrequency="2e4"
-      targetExactTimestep="0"
-      target="/Outputs/restart"/>   -->
-  </Events> 
 </Problem>