feat: QD rate and state friction coupled with Mechanics and Poromechanics solvers.

inputFiles/inducedSeismicity/SCEC_BP6_QD_S_base.xml

@@ -0,0 +1,217 @@
+<?xml version="1.0" ?>
+<Problem>
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <Mesh>
+    <InternalMesh
+      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 }"
+      ny="{ 1 }"
+      nz="{ 20, 201, 20 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>    
+    <Box
+      name="faultPlane"
+      xMin="{ -0.01, -1.01e3, -10.0e3 }"
+      xMax="{ 0.01,   1.01e3,  10.0e3 }"/>
+
+    <Box
+      name="source"
+      xMin="{ -0.01, -1.01e3, -50.0 }"
+      xMax="{ 0.01,   1.01e3,  50.0 }"/>
+  </Geometry>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ cb1 }"
+      materialList="{rock}"/>
+
+    <SurfaceElementRegion
+      name="Fault"
+      materialList="{frictionLaw, faultMaterial, water}"
+      defaultAperture="1e-3"/>  
+    </ElementRegions>
+
+    <Constitutive>
+
+    <CompressibleSinglePhaseFluid
+      name="water"
+      defaultDensity="1000"
+      defaultViscosity="1.0e-3"
+      referencePressure="0.0"
+      compressibility="1.0e-8"
+      viscosibility="0.0"/>
+
+    <CompressibleSolidConstantPermeability
+      name="faultMaterial"
+      solidModelName="nullSolid"
+      porosityModelName="faultPorosity"
+      permeabilityModelName="faultPermeability"/>  
+
+    <PressurePorosity
+      name="faultPorosity"
+      defaultReferencePorosity="0.100"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+
+    <NullModel
+      name="nullSolid"/>  
+
+    <ConstantPermeability
+      name="faultPermeability"
+      permeabilityComponents="{ 1.0e-13, 1.0e-13, 1.0e-13 }"/>
+
+    <ElasticIsotropic
+     name="rock"
+     defaultDensity="2670"
+     defaultBulkModulus="50.0e9"
+     defaultShearModulus="32.04e9"/>
+
+    <RateAndStateFriction
+      name="frictionLaw"
+      defaultA="0.007"
+      defaultB="0.005"
+      defaultDc="4.0e-3"  
+      defaultReferenceVelocity="1.0e-6"
+      defaultReferenceFrictionCoefficient="0.6"/>  
+  </Constitutive>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="BP6_QD_S"/>
+
+    <TimeHistory
+      name="timeHistoryOutput"
+      sources="{/Tasks/pressureCollection, /Tasks/slipRateCollection}"
+      filename="BP6_DQ_S" />  
+  </Outputs>
+
+  <Tasks>
+    <PackCollection
+      name="slipRateCollection"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      fieldName="slipRate"
+      setNames="{source}"/>
+
+    <PackCollection
+      name="pressureCollection"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      fieldName="pressure"
+      setNames="{source}"/>    
+   </Tasks>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ cb1 }"
+      materialList="{rock}"/>
+
+    <SurfaceElementRegion
+      name="Fault"
+      materialList="{frictionLaw}"
+      defaultAperture="1e-3"/>  
+    </ElementRegions>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="fault"
+      fieldName="ruptureState"
+      initialCondition="1"
+      objectPath="faceManager"
+      scale="1"
+      setNames="{ faultPlane }"/>
+
+    <FieldSpecification
+      name="xconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ xneg, xpos }"/>
+
+    <FieldSpecification
+      name="yconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ ypos, yneg }"/>
+
+    <FieldSpecification
+      name="zconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zpos, zneg }"/>
+
+    <FieldSpecification
+      name="initialStateVariable"
+      objectPath="ElementRegions/Fault"
+      fieldName="stateVariable"
+      scale="0.68070857"
+      setNames="{ all }"
+      initialCondition="1"/>  
+
+    <FieldSpecification
+      name="slipVelocity1"
+      fieldName="slipVelocity"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="0"
+      scale="1.0e-12"
+      setNames="{all}"/>
+
+    <FieldSpecification
+      name="slipVelocity2"
+      fieldName="slipVelocity"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="1"
+      scale="0."
+      setNames="{all}"/>
+
+    <FieldSpecification
+      name="backgroundShearStress"
+      fieldName="backgroundShearStress"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="0"
+      scale="29.2e6"
+      setNames="{all}"/> 
+
+    <FieldSpecification
+      name="backgroundNormalStress"
+      fieldName="backgroundNormalStress"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      scale="50.0e6"
+      setNames="{all}"/>                      
+
+    <SourceFlux
+      name="sourceTerm"
+      objectPath="ElementRegions/Fault"
+      scale="-1.25e-3"
+      setNames="{ source }"
+      endTime="8.64e6"/>
+  </FieldSpecifications>
+</Problem>

inputFiles/inducedSeismicity/SCEC_BP6_QD_S_explicit.xml

@@ -0,0 +1,139 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Included>
+    <File
+      name="./SCEC_BP6_QD_S_base.xml"/>
+  </Included>
+
+  <Solvers
+    gravityVector="{0.0, 0.0, 0.0}">
+    <ExplicitQuasiDynamicEQ
+      name="QDSolver"
+      targetRegions="{ Domain, Fault }"
+      shearImpedance="4.6247113164e6"
+      initialDt="1"
+      logLevel="1"
+      discretization="FE1"
+      stressSolverName="stressSolver">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-8"
+        logLevel="2"
+        newtonMaxIter="20"
+        maxNumConfigurationAttempts="1"
+        lineSearchMaxCuts="2"
+        maxTimeStepCuts="1"/>
+     </ExplicitQuasiDynamicEQ>
+
+    <OneWayCoupledFractureFlowContactMechanics
+      name="stressSolver"
+      logLevel="1"
+      solidSolverName="contactSolver"
+      flowSolverName="singlePhaseFlowSolver"
+      targetRegions="{ Domain, Fault }">
+       <NonlinearSolverParameters
+        newtonMaxIter="2"
+        newtonTol="1.0e-2"
+	      couplingType="Sequential"  
+        lineSearchAction="None"
+	      subcycling="0"
+        maxTimeStepCuts="1"
+        lineSearchMaxCuts="0"/>
+    </OneWayCoupledFractureFlowContactMechanics> 
+
+    <SolidMechanicsLagrangeContactBubbleStab
+      name="contactSolver"
+      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 }"
+      rockToughness="1.0"
+      mpiCommOrder="1"
+      fractureRegion="Fault"/>  
+
+    <SinglePhaseFVM 
+      name="singlePhaseFlowSolver"
+      logLevel="1"
+      discretization="singlePhaseTPFA"
+      targetRegions="{ Fault }"> 
+      <NonlinearSolverParameters
+        newtonTol = "1.0e-6"
+        newtonMaxIter ="8"/>
+      <LinearSolverParameters
+        solverType="direct"
+        directParallel="0"/>   
+    </SinglePhaseFVM>
+  </Solvers>
+
+
+  <Events
+    maxTime="6.312e7">
+
+   <!-- fault generation  -->
+   <SoloEvent 
+     name="generateFault"
+     target="/Solvers/SurfaceGen"/>
+
+  <!-- vtk output -->
+   <PeriodicEvent
+     name="vtkOutput"
+     cycleFrequency="20"
+     targetExactTimestep="0"
+     target="/Outputs/vtkOutput"/>
+
+  <!-- solver application  -->
+   <PeriodicEvent
+     name="solverApplications1"
+     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
+      name="slipRateCollection"
+      cycleFrequency="1"
+      targetExactTimestep="0"
+      target="/Tasks/slipRateCollection" />
+
+    <PeriodicEvent
+      name="pressureCollection"
+      cycleFrequency="1"
+      targetExactTimestep="0"
+      target="/Tasks/pressureCollection" />  
+
+  <!-- Time history output -->
+   <PeriodicEvent
+      name="timeHistoryOutput"
+      cycleFrequency="10"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput"/>
+  </Events>
+</Problem>