feat: add rate- and state-dependent friction model & quasi-dynamic solver.

inputFiles/inducedSeismicity/SpringSlider_base.xml

@@ -0,0 +1,165 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Solvers>
+    <QuasiDynamicEQ
+      name="SpringSlider"
+      targetRegions="{ Fault }"
+      shearImpedance="4.41"
+      initialDt="1"
+      logLevel="1"
+      discretization="FE1">
+    </QuasiDynamicEQ>
+
+    <SurfaceGenerator
+      name="SurfaceGen"
+      targetRegions="{ Domain }"
+      rockToughness="1.0"
+      mpiCommOrder="1"
+      fractureRegion="Fault"/>  
+  </Solvers>
+
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+  </NumericalMethods>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0, 1 }"
+      yCoords="{ 0, 2 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 100 }"
+      ny="{ 2 }"
+      nz="{ 1 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>    
+    <ThickPlane
+      name="faultPlane"
+      normal="{ 0, 1, 0 }"
+      origin="{ 0, 1, 0 }"
+      thickness="0.1"/>
+  </Geometry>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ cb1 }"
+      materialList="{}"/>
+
+    <SurfaceElementRegion
+      name="Fault"
+      materialList="{frictionLaw}"
+      defaultAperture="1e-3"/>  
+    </ElementRegions>
+
+  <Constitutive>
+    <RateAndStateFriction
+      name="frictionLaw"
+      defaultA="0.01"
+      defaultB="0.015"
+      defaultDc="1.0e-5"  
+      defaultReferenceVelocity="1.0e-6"
+      defaultReferenceFrictionCoefficient="0.6"/>  
+  </Constitutive>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="fault"
+      fieldName="ruptureState"
+      initialCondition="1"
+      objectPath="faceManager"
+      scale="1"
+      setNames="{ faultPlane }"/>
+
+    <FieldSpecification
+      name="normalTraction"
+      fieldName="traction"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="0"
+      scale="50"
+      setNames="{all}"/>
+
+    <FieldSpecification
+      name="shearTraction1"
+      fieldName="traction"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="1"
+      scale="21.2132034356"
+      setNames="{all}"/>
+    <FieldSpecification
+      name="shearTraction2"
+      fieldName="traction"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="2"
+      scale="21.2132034356"
+      setNames="{all}"/>  
+    <FieldSpecification
+      name="stateVariable"
+      fieldName="stateVariable"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      scale="0.6"
+      setNames="{all}"/>
+    <FieldSpecification
+      name="slipRate"
+      fieldName="slipRate"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      scale="1.0e-6"
+      setNames="{all}"/>
+    <FieldSpecification
+      name="slipVelocity1"
+      fieldName="slipVelocity"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="0"
+      scale="0."
+      setNames="{all}"/>               
+    <FieldSpecification
+      name="slipVelocity2"
+      fieldName="slipVelocity"
+      initialCondition="1"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      component="1"
+      scale="0."
+      setNames="{all}"/>               
+  </FieldSpecifications>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="springSlider"/>
+
+    <TimeHistory
+      name="timeHistoryOutput"
+      sources="{/Tasks/slipCollection,/Tasks/slipRateCollection,/Tasks/stateVariableCollection}"
+      filename="springSlider"/>
+  </Outputs>
+
+  <Tasks>
+    <PackCollection
+      name="slipCollection"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      fieldName="displacementJump"/>
+
+    <PackCollection
+      name="slipRateCollection"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      fieldName="slipRate"/>
+
+    <PackCollection
+      name="stateVariableCollection"
+      objectPath="ElementRegions/Fault/FractureSubRegion"
+      fieldName="stateVariable"/>  
+  </Tasks>
+</Problem>

inputFiles/inducedSeismicity/SpringSlider_smoke.xml

@@ -0,0 +1,46 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Included>
+    <File
+      name="./SpringSlider_base.xml"/>
+  </Included>
+
+  <Events
+    maxTime="4.5e4">
+
+    <SoloEvent 
+      name="generateFault"
+      target="/Solvers/SurfaceGen"/>
+
+    <PeriodicEvent
+      name="vtkOutput"
+      cycleFrequency="1"
+      targetExactTimestep="0"
+      target="/Outputs/vtkOutput"/>  
+
+    <PeriodicEvent
+      name="solverApplications"
+      maxEventDt="1e4"
+      target="/Solvers/SpringSlider"/>
+
+    <PeriodicEvent
+      name="seismicRateCollection"
+      cycleFrequency="1"
+      target="/Tasks/slipCollection"/>
+
+    <PeriodicEvent
+      name="slipRateCollection"
+      cycleFrequency="1"
+      target="/Tasks/slipRateCollection"/>
+
+    <PeriodicEvent
+      name="stateVariableCollection"
+      cycleFrequency="1"
+      target="/Tasks/stateVariableCollection"/>  
+
+    <PeriodicEvent
+      name="timeHistoryOutput"
+      cycleFrequency="1"
+      target="/Outputs/timeHistoryOutput" />    
+  </Events>
+</Problem>

inputFiles/inducedSeismicity/inducedSeismicity.ats

@@ -15,12 +15,19 @@ decks = [
         restart_step=20,
         check_step=30,
         restartcheck_params=RestartcheckParameters(atol=1e-4, rtol=1e-3)),
-     TestDeck(
+    TestDeck(
         name="SeismicityRate_analytical_verification_smoke",
         description="Prescribed logarithmic stressing history",
         partitions=((1, 1, 1), ),
         restart_step=0,
         check_step=100,
         curvecheck_params=CurveCheckParameters(**curvecheck_params))
+    TestDeck(
+        name="SpringSlider_smoke",
+        description="Spring slider 0D system",
+        partitions=((1, 1, 1), ),
+        restart_step=0,
+        check_step=3265,
+        restartcheck_params=RestartcheckParameters(atol=1e-4, rtol=1e-3) )
         ]
 generate_geos_tests(decks)

src/coreComponents/constitutive/CMakeLists.txt

@@ -25,6 +25,7 @@ set( constitutive_headers
      contact/HydraulicApertureBase.hpp
      contact/HydraulicApertureRelationSelector.hpp
      contact/HydraulicApertureTable.hpp
+     contact/RateAndStateFriction.hpp
      diffusion/ConstantDiffusion.hpp
      diffusion/DiffusionBase.hpp
      diffusion/DiffusionFields.hpp
@@ -194,6 +195,7 @@ set( constitutive_sources
      contact/FrictionlessContact.cpp
      contact/HydraulicApertureBase.cpp
      contact/HydraulicApertureTable.cpp
+     contact/RateAndStateFriction.cpp
      diffusion/ConstantDiffusion.cpp
      diffusion/DiffusionBase.cpp
      dispersion/DispersionBase.cpp

src/coreComponents/constitutive/contact/CoulombFriction.cpp

@@ -50,7 +50,6 @@ CoulombFriction::CoulombFriction( string const & name, Group * const parent ):
   registerWrapper( viewKeyStruct::elasticSlipString(), &m_elasticSlip ).
     setApplyDefaultValue( 0.0 ).
     setDescription( "Elastic Slip" );
-
 }
 
 CoulombFriction::~CoulombFriction()

src/coreComponents/constitutive/contact/CoulombFriction.hpp

@@ -169,20 +169,6 @@ class CoulombFriction : public FrictionBase
   virtual void allocateConstitutiveData( dataRepository::Group & parent,
                                          localIndex const numConstitutivePointsPerParentIndex ) override final;
 
-  /**
-   * @brief Const accessor for cohesion
-   * @return A const reference to arrayView1d<real64 const> containing the
-   *         cohesions (at every element).
-   */
-  real64 const & cohesion() const { return m_cohesion; }
-
-  /**
-   * @brief Const accessor for friction angle
-   * @return A const reference to arrayView1d<real64 const> containing the
-   *         friction coefficient (at every element).
-   */
-  real64 const & frictionCoefficient() const { return m_frictionCoefficient; }
-
   /// Type of kernel wrapper for in-kernel update
   using KernelWrapper = CoulombFrictionUpdates;