feat: Enable aperture modifications without a full mechanics solver

host-configs/ubuntu.cmake

@@ -0,0 +1,32 @@
+set( CONFIG_NAME "ubuntu" ) 
+
+# Set compilers path
+set(CMAKE_C_COMPILER "/usr/bin/gcc" CACHE PATH "")   # This is typically something like /usr/bin/gcc ... or clang
+set(CMAKE_CXX_COMPILER "/usr/bin/g++" CACHE PATH "") # This is typically something like /usr/bin/g++ ... or clang++
+set(ENABLE_FORTRAN OFF CACHE BOOL "" FORCE)
+
+# Set paths to mpi
+set(ENABLE_MPI ON CACHE PATH "")
+set(MPI_C_COMPILER "/usr/bin/mpicc" CACHE PATH "")    # This is typically something like /usr/bin/mpicc
+set(MPI_CXX_COMPILER "/usr/bin/mpicxx" CACHE PATH "") # This is typically something like /usr/bin/mpicxx
+set(MPIEXEC "/usr/bin/mpirun" CACHE PATH "")          # This is typically something like /usr/bin/mpirun
+
+# Set paths to blas and lapack
+set( BLAS_LIBRARIES "/usr/lib/x86_64-linux-gnu/libblas.so" CACHE PATH "" FORCE )     # This is typically something like /usr/lib64/libblas.so 
+set( LAPACK_LIBRARIES "/usr/lib/x86_64-linux-gnu/liblapack.so" CACHE PATH "" FORCE ) # This is typically something like /usr/lib64/liblapack.so
+
+# Cuda and openMP
+set( ENABLE_CUDA OFF CACHE PATH "" FORCE )
+set( ENABLE_OPENMP OFF CACHE PATH "" FORCE )
+
+# TPLs
+set( ENABLE_TRILINOS OFF CACHE PATH "" FORCE )
+set( ENABLE_CALIPER OFF CACHE PATH "" FORCE )
+set( ENABLE_DOXYGEN OFF CACHE BOOL "" FORCE)
+set( ENABLE_MATHPRESSO OFF CACHE BOOL "" FORCE )
+
+if(NOT ( EXISTS "${GEOS_TPL_DIR}" AND IS_DIRECTORY "${GEOS_TPL_DIR}" ) )
+   set(GEOS_TPL_DIR "${CMAKE_SOURCE_DIR}/../../../thirdPartyLibs/install-${CONFIG_NAME}-release" CACHE PATH "" FORCE )
+endif()
+
+include(${CMAKE_CURRENT_LIST_DIR}/tpls.cmake)

inputFiles/stressPathDrivenGeomechanics/fractureMatrixFlow_edfm_BBP_base.xml

@@ -0,0 +1,158 @@
+<?xml version="1.0" ?> 
+
+<Problem>
+  <!-- SPHINX_TUT_STRESS_PATH_DRIVEN -->
+  <Solvers
+    gravityVector="{ 0.0, 0.0, -9.81 }">
+    <SinglePhaseFVM
+      name="SinglePhaseFlow"
+      logLevel="1"
+      discretization="singlePhaseTPFA"
+      targetRegions="{ RockMatrix, Fracture }">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-6"
+        newtonMaxIter="30"/>
+      <LinearSolverParameters
+        directParallel="0"/>
+    </SinglePhaseFVM>
+	<!-- SPHINX_TUT_STRESS_PATH_DRIVEN_END -->
+
+    <EmbeddedSurfaceGenerator
+      name="SurfaceGenerator"
+      logLevel="1"
+      discretization="FE1"
+      targetRegions="{ RockMatrix, Fracture }"
+      targetObjects="{ FracturePlane }"
+      fractureRegion="Fracture"
+      mpiCommOrder="1"/>
+  </Solvers>
+
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"/>
+    </FiniteVolume>
+    <FiniteElements>
+    <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="RockMatrix"
+      cellBlocks="{ * }" 
+      materialList="{ water, rock }"/>
+
+    <SurfaceElementRegion
+      name="Fracture"
+      faceBlock="embeddedSurfaceSubRegion"
+      defaultAperture="1.0e-4"
+      materialList="{ water, fractureFilling }"
+      subRegionType="embeddedElement"/>
+  </ElementRegions>
+
+  <Constitutive>
+    <CompressibleSinglePhaseFluid
+      name="water"
+      defaultDensity="2000"
+      defaultViscosity="0.001"
+      referencePressure="0.0"
+      compressibility="5e-12"
+      viscosibility="0.0"/>
+
+    <CompressibleSolidBartonBandisPermeability
+      name="fractureFilling"
+      solidModelName="nullSolid"
+      porosityModelName="fracturePorosity"
+      permeabilityModelName="fracturePerm"/>
+
+    <CompressibleSolidConstantPermeability
+      name="rock"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"/>
+
+    <NullModel
+      name="nullSolid"/>
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.1"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-15, 1.0e-15, 1.0e-15 }"/>
+
+    <PressurePorosity
+      name="fracturePorosity"
+      defaultReferencePorosity="1.00"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+	<!-- SPHINX_TUT_STRESS_PATH_DRIVEN_LAW -->
+    <BartonBandisPermeability
+      name="fracturePerm"
+      biotCoefficient="1.0" 
+      poissonRatio="0.3"	
+      normalStiffness="10.0e9" 
+      referenceAperture="1.0e-4"
+      referencePressure="1e7" 
+      referenceTotalStress="{ 10.0e6, 10.0e6, 30.0e6 }"/>
+	<!-- SPHINX_TUT_STRESS_PATH_DRIVEN_LAW_END -->
+
+  </Constitutive>
+
+
+  <FieldSpecifications>
+
+    <FieldSpecification
+      name="Porosity"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/RockMatrix/cb1"
+      fieldName="rockPorosity_referencePorosity"
+      scale="0.1"/>
+
+    <FieldSpecification
+      name="fracPorosity"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fracture"
+      fieldName="fracturePorosity_porosity"
+      scale="1.00"/>
+
+    <FieldSpecification
+      name="initialFracturePressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="pressure"
+      scale="1.0e7"/>
+
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/RockMatrix/cb1"
+      fieldName="pressure"
+      scale="0.0"/>
+
+    <FieldSpecification
+      name="sourceTerm"
+      objectPath="ElementRegions/Fracture"
+      fieldName="pressure"
+      scale="1.0e7"
+      setNames="{ source }"/>
+
+    <FieldSpecification
+      name="sinkTerm"
+      objectPath="ElementRegions/Fracture"
+      fieldName="pressure"
+      scale="1.904100e7"
+      setNames="{ sink }"/>
+
+  </FieldSpecifications>
+</Problem>

inputFiles/stressPathDrivenGeomechanics/fractureMatrixFlow_edfm_horizontalFrac_BBP_smoke.xml

@@ -0,0 +1,122 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Included>
+    <File name="./fractureMatrixFlow_edfm_BBP_base.xml" />
+  </Included>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ -0.5, 0.5 }"
+      yCoords="{ 0, 11 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 1 }"
+      ny="{ 11 }"
+      nz="{ 1 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>
+    <Box
+      name="source"
+      xMin="{ -1.01, -0.01, -0.01 }"
+      xMax="{ 1.01, 1.1, 1.01 }"/>
+
+    <Box
+      name="sink"
+      xMin="{ -1.01, 9.9, -0.01 }"
+      xMax="{ 1.01, 11.01, 1.01 }"/>
+
+    <Rectangle
+      name="FracturePlane"
+      normal="{ 0, 0, 1 }"
+      origin="{ 0.0, 0.5, 0.5 }"
+      lengthVector="{ 0, 1, 0 }"
+      widthVector="{ 1, 0, 0 }"
+      dimensions="{ 30, 1 }"/>
+  </Geometry>
+
+  <Tasks>
+    <PackCollection
+      name="pressureCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="pressure"/> 
+    <PackCollection
+      name="apertureCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="hydraulicAperture"/> 
+    <PackCollection
+      name="permeabilityCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="fracturePerm_permeability"/>           
+  </Tasks>
+
+  <Events
+    maxTime="50"
+    logLevel="3">
+    <SoloEvent
+      name="preFracture"
+      target="/Solvers/SurfaceGenerator"/>
+
+    <PeriodicEvent
+      name="solverApplications"
+      forceDt="1"
+      target="/Solvers/SinglePhaseFlow"/>
+
+    <PeriodicEvent
+      name="outputs"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="pressureHistoryCollection"
+      target="/Tasks/pressureCollection" />
+    <PeriodicEvent
+      name="pressureHistoryOutput"
+      targetExactTimestep="0"
+      target="/Outputs/pressureHistoryOutput"/> 
+
+    <PeriodicEvent
+      name="apertureHistoryCollection"
+      target="/Tasks/apertureCollection" />
+    <PeriodicEvent
+      name="apertureHistoryOutput"
+      targetExactTimestep="0"
+      target="/Outputs/apertureHistoryOutput"/> 
+
+    <PeriodicEvent
+      name="permeabilityHistoryCollection"
+      target="/Tasks/permeabilityCollection" />
+    <PeriodicEvent
+      name="permeabilityHistoryOutput"
+      targetExactTimestep="0"
+      target="/Outputs/permeabilityHistoryOutput"/> 
+  </Events>
+
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="stressPathDriven"
+      plotLevel="2"
+      format="ascii"/>
+
+    <Restart
+      name="sidreRestart"/>
+	<!-- SPHINX_TUT_STRESS_PATH_DRIVEN_TIMEHISTORY -->
+    <TimeHistory
+      name="pressureHistoryOutput"
+      sources="{/Tasks/pressureCollection}"
+      filename="fracPressure_history" />
+    <TimeHistory
+      name="apertureHistoryOutput"
+      sources="{/Tasks/apertureCollection}"
+      filename="fracAperture_history" />
+    <TimeHistory
+      name="permeabilityHistoryOutput"
+      sources="{/Tasks/permeabilityCollection}"
+      filename="fracPermeability_history" />
+	<!-- SPHINX_TUT_STRESS_PATH_DRIVEN_TIMEHISTORY_END -->
+  </Outputs>
+</Problem>

src/coreComponents/constitutive/CMakeLists.txt

@@ -151,6 +151,7 @@ set( constitutive_headers
      permeability/DamagePermeability.hpp
      permeability/ExponentialDecayPermeability.hpp
      permeability/ParallelPlatesPermeability.hpp
+     permeability/BartonBandisPermeability.hpp
      permeability/PermeabilityBase.hpp
      permeability/PermeabilityFields.hpp
      permeability/PressurePermeability.hpp
@@ -303,6 +304,7 @@ set( constitutive_sources
      permeability/DamagePermeability.cpp
      permeability/ExponentialDecayPermeability.cpp
      permeability/ParallelPlatesPermeability.cpp
+     permeability/BartonBandisPermeability.cpp
      permeability/PermeabilityBase.cpp
      permeability/PressurePermeability.cpp
      permeability/ProppantPermeability.cpp

src/coreComponents/constitutive/ConstitutivePassThru.hpp

@@ -49,6 +49,7 @@
 #include "permeability/CarmanKozenyPermeability.hpp"
 #include "permeability/ExponentialDecayPermeability.hpp"
 #include "permeability/ParallelPlatesPermeability.hpp"
+#include "permeability/BartonBandisPermeability.hpp"
 #include "permeability/PressurePermeability.hpp"
 #include "permeability/ProppantPermeability.hpp"
 #include "permeability/SlipDependentPermeability.hpp"
@@ -375,6 +376,7 @@ struct ConstitutivePassThru< CompressibleSolidBase >
                                  CompressibleSolid< PressurePorosity, CarmanKozenyPermeability >,
                                  CompressibleSolid< PressurePorosity, ExponentialDecayPermeability >,
                                  CompressibleSolid< PressurePorosity, ParallelPlatesPermeability >,
+                                 CompressibleSolid< PressurePorosity, BartonBandisPermeability >,
                                  CompressibleSolid< PressurePorosity, PressurePermeability >,
                                  CompressibleSolid< PressurePorosity, SlipDependentPermeability >,
                                  CompressibleSolid< PressurePorosity, WillisRichardsPermeability >
@@ -389,6 +391,7 @@ struct ConstitutivePassThru< CompressibleSolidBase >
                                  CompressibleSolid< PressurePorosity, CarmanKozenyPermeability >,
                                  CompressibleSolid< PressurePorosity, ExponentialDecayPermeability >,
                                  CompressibleSolid< PressurePorosity, ParallelPlatesPermeability >,
+                                 CompressibleSolid< PressurePorosity, BartonBandisPermeability >,
                                  CompressibleSolid< PressurePorosity, PressurePermeability >,
                                  CompressibleSolid< PressurePorosity, SlipDependentPermeability >,
                                  CompressibleSolid< PressurePorosity, WillisRichardsPermeability >
@@ -462,6 +465,7 @@ struct ConstitutivePassThru< CoupledSolidBase >
                                  CompressibleSolid< PressurePorosity, CarmanKozenyPermeability >,
                                  CompressibleSolid< PressurePorosity, ExponentialDecayPermeability >,
                                  CompressibleSolid< PressurePorosity, ParallelPlatesPermeability >,
+                                 CompressibleSolid< PressurePorosity, BartonBandisPermeability >,
                                  CompressibleSolid< PressurePorosity, PressurePermeability >,
                                  CompressibleSolid< PressurePorosity, SlipDependentPermeability >,
                                  CompressibleSolid< PressurePorosity, WillisRichardsPermeability >,
@@ -503,6 +507,7 @@ struct ConstitutivePassThru< CoupledSolidBase >
                                  CompressibleSolid< PressurePorosity, CarmanKozenyPermeability >,
                                  CompressibleSolid< PressurePorosity, ExponentialDecayPermeability >,
                                  CompressibleSolid< PressurePorosity, ParallelPlatesPermeability >,
+                                 CompressibleSolid< PressurePorosity, BartonBandisPermeability >,
                                  CompressibleSolid< PressurePorosity, PressurePermeability >,
                                  CompressibleSolid< PressurePorosity, SlipDependentPermeability >,
                                  CompressibleSolid< PressurePorosity, WillisRichardsPermeability >,