feat: add MGR for single phase poromechanics with conforming fractures and wells

inputFiles/thermoPoromechanicsFractures/THM_conforming_wellbore/displacedFault.vtm

@@ -0,0 +1,7 @@
+<?xml version="1.0"?>
+<VTKFile type="vtkMultiBlockDataSet" version="1.0">
+    <vtkMultiBlockDataSet>
+        <DataSet name="main" file="Domain_Mesh.vtu"/>
+        <DataSet name="fracture" file="Fault_Mesh.vtu"/>
+    </vtkMultiBlockDataSet>
+</VTKFile>

inputFiles/thermoPoromechanicsFractures/THM_conforming_wellbore/domain_boundaries.xml

@@ -0,0 +1,30 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Geometry>
+     <Box
+      name="xmin"
+      xMin="{ -4.0000000e-03, -2.0000000e-03, -1.0000010e+03 }"
+      xMax="{ 4.0000000e-03, 2.0000020e+03, 1.0000000e-03 }"/>
+     <Box
+      name="xmax"
+      xMin="{ 3.9999960e+03, -2.0000000e-03, -1.0000010e+03 }"
+      xMax="{ 4.0000040e+03, 2.0000020e+03, 1.0000000e-03 }"/>
+     <Box
+      name="ymin"
+      xMin="{ -4.0000000e-03, -2.0000000e-03, -1.0000010e+03 }"
+      xMax="{ 4.0000040e+03, 2.0000000e-03, 1.0000000e-03 }"/>
+     <Box
+      name="ymax"
+      xMin="{ -4.0000000e-03, 1.9999980e+03, -1.0000010e+03 }"
+      xMax="{ 4.0000040e+03, 2.0000020e+03, 1.0000000e-03 }"/>
+     <Box
+      name="zmin"
+      xMin="{ -4.0000000e-03, -2.0000000e-03, -1.0000010e+03 }"
+      xMax="{ 4.0000040e+03, 2.0000020e+03, -9.9999900e+02 }"/>
+     <Box
+      name="zmax"
+      xMin="{ -4.0000000e-03, -2.0000000e-03, -1.0000000e-03 }"
+      xMax="{ 4.0000040e+03, 2.0000020e+03, 1.0000000e-03 }"/>
+  </Geometry>
+</Problem>

inputFiles/thermoPoromechanicsFractures/THM_conforming_wellbore/faultPoroelastic_isothermal_base.xml

@@ -0,0 +1,197 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <!-- SPHINX_ELEMENTREGIONS -->
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ 0_hexahedra, 1_hexahedra, 2_hexahedra, 3_hexahedra, 4_hexahedra, 5_hexahedra }"
+      materialList="{ water, porousRock }"/>  
+
+    <SurfaceElementRegion
+      name="Fault"
+      faceBlock="fracture"
+      materialList="{ water, faultFilling, fractureContact, hApertureModel }"
+      defaultAperture="1e-3"/>
+
+    <WellElementRegion
+      name="WELL.W001"
+      materialList="{ water }" />
+  </ElementRegions> 
+  <!-- SPHINX_ELEMENTREGIONS_END -->
+
+  <!-- SPHINX_NUMERICAL -->
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"/>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA_stab"/>  
+    </FiniteVolume>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+  </NumericalMethods>
+  <!-- SPHINX_NUMERICAL_END -->
+
+  <!-- SPHINX_MATERIAL -->
+  <Constitutive>
+    <PorousElasticIsotropic
+      name="porousRock"
+      solidModelName="rock"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"
+    />
+
+    <ElasticIsotropic
+      name="rock"
+      defaultDensity="2700"
+      defaultYoungModulus="14.95e9"
+      defaultPoissonRatio="0.15"
+    />
+
+    <CompressibleSinglePhaseFluid 
+      name="water"
+      defaultDensity="1000"
+      defaultViscosity="0.001"
+      referencePressure="0e6"
+      referenceDensity="1000"
+      compressibility="2.09028227021e-10"
+      referenceViscosity="0.001"
+      viscosibility="0.0"
+    />
+
+    <BiotPorosity
+      name="rockPorosity"
+      defaultGrainBulkModulus="7.12e10"
+      defaultReferencePorosity="0.3"
+    />  
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{1.0e-12, 1.0e-12, 1.0e-12}"
+    /> 
+
+    <CompressibleSolidParallelPlatesPermeability
+      name="faultFilling"
+      solidModelName="nullSolid"
+      porosityModelName="fracturePorosity"
+      permeabilityModelName="fracturePerm"/>
+
+    <PressurePorosity
+      name="fracturePorosity"
+      defaultReferencePorosity="1.00"
+      referencePressure="0e6"
+      compressibility="0.0"/>
+
+    <ParallelPlatesPermeability
+      name="fracturePerm"/>
+
+    <Coulomb
+      name="fractureContact"
+      cohesion="1.0e7"
+      frictionCoefficient="0.5"/>
+
+    <HydraulicApertureTable
+      name="hApertureModel"
+      apertureTableName="apertureTable"/>   
+
+    <NullModel
+      name="nullSolid"/>
+  </Constitutive>
+  <!-- SPHINX_MATERIAL_END -->
+
+  <Functions>
+    <TableFunction
+      name="apertureTable"
+      coordinates="{ -1.0e-3, 0.0 }"
+      values="{ 1.0e-3, 1.0e-3 }"/>
+  </Functions>
+
+  <!-- SPHINX_BC --> 
+  <FieldSpecifications>
+    <FieldSpecification
+      name="fracPorosity"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fault"
+      fieldName="fracturePorosity_porosity"
+      scale="1.00"/>
+
+    <FieldSpecification
+      name="xconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ xmin, xmax }"/>
+
+    <FieldSpecification
+      name="yconstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ ymin, ymax }"/>
+
+    <FieldSpecification
+      name="zconstraintFront"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zmin }"/>
+
+    <Traction
+      name="NormalTraction"
+      objectPath="faceManager"
+      tractionType="normal"
+      scale="-3.0e6"      
+      setNames="{ zmax }"/>  
+
+    <HydrostaticEquilibrium
+      name="equil"
+      objectPath="ElementRegions"
+      maxNumberOfEquilibrationIterations="100"
+      datumElevation="0"
+      datumPressure="0"/>	  
+  </FieldSpecifications>
+  <!-- SPHINX_BC_END -->
+
+  <Tasks>
+    <SinglePhasePoromechanicsConformingFracturesInitialization
+      logLevel="1"      
+      name="ELASTICITY.PRE.INIT.STEP"
+      poromechanicsSolverName="poroFractureSolver"/>
+
+    <PackCollection
+      name="tractionCollection"
+      objectPath="ElementRegions/Fault/fracture"
+      fieldName="traction"/>
+
+    <PackCollection
+      name="displacementJumpCollection"
+      objectPath="ElementRegions/Fault/fracture"
+      fieldName="displacementJump"/>   
+  </Tasks> 
+
+  <!-- SPHINX_OUTPUT -->  
+  <Outputs>
+    <VTK 
+      name="vtkOutput"
+      plotFileRoot="faultModel"/>
+
+    <TimeHistory
+      name="tractionOutput"
+      sources="{/Tasks/tractionCollection}"
+      filename="traction_history" />
+
+    <TimeHistory
+      name="displacementJumpOutput"
+      sources="{/Tasks/displacementJumpCollection}"
+      filename="displacementJump_history" /> 
+  </Outputs>
+  <!-- SPHINX_OUTPUT_END -->
+</Problem>