GEOS-DEV · jafranc · Nov 4, 2025 · Nov 5, 2025 · Nov 5, 2025 · Nov 6, 2025
@@ -0,0 +1,32 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <!-- SPHINX_TASK -->
+  <Tasks>
+    <FrictionDriver
+      name="frictionDriver"
+      friction="fractureContact0"
+      dispControl="jFunction"
+      dJumpControl="djFunction"
+      stressControlsL="sLFunction"
+      stressControlsR="sRFunction"
+      yTiltAngle="0"
+      zTiltAngle="0"
+      steps="1"
+      simultaneous="1"
+      tolJumpN="1e-1"
+      tolJumpT="1e-5"
+      tolNormalTrac="1e-1"
+      iterPenNFac="10"
+      iterPenTFac="1"
+      faceArea="1"
+      neighborsVolume="{1,1}"
+      neighborsShear="{1e9,1e9}"
+      neighborsBulk="{1e7,1e7}"
+      output="frictionDriver_Coulomb.txt" />
+  </Tasks>
+
+  <Included>
+    <File name="./frictionDriver_base.xml"/>
+  </Included>
+</Problem>
@@ -0,0 +1,63 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Outputs>
+    <Restart
+      name="restartOutput"/>
+  </Outputs>
+
+  <!-- SPHINX_EVENTS -->
+  <Events
+    maxTime="1">
+    <SoloEvent
+      name="frictionDriver"
+      targetTime="0"
+      target="/Tasks/frictionDriver"/>
+
+    <PeriodicEvent
+      name="restarts"
+      timeFrequency="1"
+      targetExactTimestep="0"
+      target="/Outputs/restartOutput"/>          
+  </Events>
+  <!-- SPHINX_EVENTS_END -->
+
+  <Constitutive>  
+    <Coulomb
+      name="fractureContact0"
+      defaultCohesion="0.0e6"
+      defaultFrictionCoefficient="0.01"/>
+    <Coulomb
+      name="fractureContact1"
+      defaultCohesion="1.0e6"
+      defaultFrictionCoefficient="0.01"/>
+
+  </Constitutive>
+
+  <!-- SPHINX_FUNCTION -->
+  <Functions>
+    <TableFunction
+      name="jFunction"
+      inputVarNames="{ time }"
+      coordinateFiles="{ tables/time.geos }"
+      voxelFile="tables/jumps.geos"/>
+
+    <TableFunction
+      name="djFunction"
+      inputVarNames="{ time }"
+      coordinateFiles="{ tables/time.geos }"
+      voxelFile="tables/djumps.geos"/>
+
+    <TableFunction
+      name="sLFunction"
+      inputVarNames="{ time }"
+      coordinateFiles="{ tables/time.geos }"
+      voxelFile="tables/tractions.geos"/>
+    <TableFunction
+      name="sRFunction"
+      inputVarNames="{ time }"
+      coordinateFiles="{ tables/time.geos }"
+      voxelFile="tables/constant.geos"/>
+  </Functions>
+  <!-- SPHINX_FUNCTION_END -->
+</Problem>
@@ -0,0 +1,2 @@
+1e2
+1e2
@@ -0,0 +1,2 @@
+0
+0
@@ -0,0 +1,2 @@
+0
+-5e-3
@@ -0,0 +1,2 @@
+0
+5
@@ -0,0 +1,2 @@
+0
+-5e2
@@ -189,6 +189,15 @@ class CoulombFriction : public FrictionBase
    */
   KernelWrapper createKernelUpdates() const;
 
+  // TODO (jafranc) make sure it is only helpers for FrictionDriver
+  /// getting cohesion value
+  real64 getCohesion() const
+  { return m_cohesion[0]; }
+
+  /// getting friction coeff
+  real64 getFrictionCoeff() const
+  { return m_frictionCoefficient[0]; }
+
   /**
    * @struct Set of "char const *" and keys for data specified in this class.
    */

@@ -29,12 +29,16 @@ set( constitutiveDrivers_headers
      relativePermeability/RelpermDriver.hpp
      relativePermeability/RelpermDriverRunTest.hpp
      solid/TriaxialDriver.hpp
+     contact/FrictionDriver.hpp
+     contact/FrictionDriverRunTest.hpp
    )
 #
 # Specify all sources
 #
 set( constitutiveDrivers_sources
      ConstitutiveDriver.cpp
+     contact/FrictionDriver.cpp
+     contact/FrictionDriverRunTest.cpp
      fluid/multiFluid/PVTDriver.cpp
      fluid/multiFluid/constant/PVTDriverRunTestInvariantImmiscibleFluid.cpp
      fluid/multiFluid/blackOil/PVTDriverRunTestDeadOilFluid.cpp

@@ -145,9 +145,9 @@ void ConstitutiveDriver::outputToFile() const
 
   integer const precision = LvArray::math::max( LvArray::math::min( m_precision, maxPrecision ), minPrecision );
   file << std::scientific << std::setprecision( precision );
-  integer const width =  precision+7;
+  integer const width =  precision+7+1;//spacing for minus sign
 
-  for( integer step = 0; step <= m_numSteps; ++step )
+  for( integer step = 0; step < m_table.size(0); ++step )
   {
     file << std::setw( width ) << m_table( step, 0 );
     for( integer col = 1; col < numColumns; ++col )
@@ -213,7 +213,7 @@ void ConstitutiveDriver::outputToConsole() const
   TableData tableData;
 
   stdVector< TableData::CellData > tableRow( numColumns );
-  for( integer step = 0; step <= m_numSteps; ++step )
+  for( integer step = 0; step < m_table.size(0); ++step )
   {
     for( integer col = 0; col < numColumns; ++col )
     {