[Feature] implement fault permeability models

inputFiles/poromechanics/SlipPermeability_embeddedFrac.xml

@@ -0,0 +1,377 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Solvers
+    gravityVector="{0.0, 0.0, 0.0}">
+    <SinglePhasePoromechanicsEmbeddedFractures
+      name="poroSolver"
+      targetRegions="{ Domain, Fracture }"
+      initialDt="1e-3"
+      fracturesSolverName="fractureMechSolver"
+      fluidSolverName="flowSolver"
+      solidSolverName="matrixMechSolver"
+      porousMaterialNames="{ porousRock, fractureFilling }"
+      logLevel="1">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-5"
+        newtonMaxIter="10"
+        maxTimeStepCuts="1"/>
+      <LinearSolverParameters
+        directParallel="0"/>
+    </SinglePhasePoromechanicsEmbeddedFractures>
+
+    <SolidMechanicsEmbeddedFractures
+      name="fractureMechSolver"
+      targetRegions="{ Domain }"
+      solidSolverName="matrixMechSolver"
+      fractureRegionName="Fracture"
+      contactRelationName="fractureContact"/>
+
+    <SolidMechanicsLagrangianSSLE
+      name="matrixMechSolver"
+      timeIntegrationOption="QuasiStatic"
+      discretization="FE1"
+      targetRegions="{ Domain }"
+      solidMaterialNames="{ rock }"/>
+
+    <SinglePhaseFVM
+      name="flowSolver"
+      discretization="singlePhaseTPFA"
+      targetRegions="{ Domain, Fracture }"
+      fluidNames="{ water }"
+      solidNames="{ porousRock, fractureFilling }"
+      permeabilityNames="{ rockPerm, fracturePerm }"/>
+
+    <EmbeddedSurfaceGenerator
+      name="SurfaceGenerator"
+      solidMaterialNames="{ rock }"
+      targetRegions="{ Domain, Fracture }"
+      fractureRegion="Fracture"
+      logLevel="2"
+      mpiCommOrder="1"/>
+  </Solvers>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ -40, -3.5, 3.5, 40 }"
+      yCoords="{ -40, -3.5, 3.5, 40 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 10, 21, 10 }"
+      ny="{ 10, 21, 10 }"
+      nz="{ 1 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+
+  <Geometry>
+    <BoundedPlane
+      name="FracturePlane"
+      normal="{-0.342020143325669, 0.939692620785908, 0.0}"
+      origin="{0.0, 0.0, 0.5}"
+      lengthVector="{0.939692620785908, 0.342020143325669, 0.0}"
+      widthVector="{0, 0, 1}"
+      dimensions="{ 2,  10 }"/>
+
+    <Box
+      name="rightPoint"
+      xMin="{ 39.9, -0.5, -0.001}"
+      xMax="{ 40.1,  0.5,  1.001}"/>
+
+    <Box
+      name="leftPoint"
+      xMin="{-40.1, -0.5, -0.001}"
+      xMax="{-39.9,  0.5,  1.001}"/>
+
+    <Box
+      name="topPoint"
+      xMin="{-0.5, 39.9, -0.001}"
+      xMax="{ 0.5, 40.1,  1.001}"/>
+
+    <Box
+      name="bottomPoint"
+      xMin="{-0.5, -40.1, -0.001}"
+      xMax="{ 0.5, -39.9,  1.001}"/>
+  </Geometry>
+
+  <Events
+    maxTime="1.0e-2">
+    <SoloEvent
+      name="preFracture"
+      target="/Solvers/SurfaceGenerator"/>  
+
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="5.0e-03"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="solverApplication0"
+      endTime="1.e-10"
+      forceDt="1.e-10"
+      target="/Solvers/poroSolver"/>
+
+    <PeriodicEvent
+      name="solverApplication1"
+      beginTime="1.e-10"
+      endTime="5.0000000e-03"
+      forceDt="4.9999999e-03"
+      target="/Solvers/poroSolver"/>
+
+    <PeriodicEvent
+      name="solverApplication2"
+      beginTime="5.0e-03"
+      endTime="0.05"
+      forceDt="5.0e-03"
+      target="/Solvers/poroSolver"/>
+
+    <PeriodicEvent
+      name="solverApplication3"
+      beginTime="0.05"
+      endTime="10.0"
+      forceDt="0.01"
+      target="/Solvers/poroSolver"/>
+
+    <PeriodicEvent
+      name="timeHistoryCollection"      
+      target="/Tasks/tractionCollection" />
+
+    <PeriodicEvent
+      name="timeHistoryOutput"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput"/>
+
+    <PeriodicEvent
+      name="timeHistoryCollection1"
+      target="/Tasks/displacementJumpCollection" />
+
+    <PeriodicEvent
+      name="timeHistoryOutput1"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput1"/> 
+
+    <PeriodicEvent
+      name="timeHistoryCollection2"
+      target="/Tasks/fracPermeabilityCollection" />
+
+    <PeriodicEvent
+      name="timeHistoryOutput2"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput2"/>     
+  </Events>
+
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"
+        fieldName="pressure"
+        coefficientName="permeability"
+        coefficientModelNames="{ rockPerm, fracturePerm, fractureContact }"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ cb1 }"
+      materialList="{ porousRock, water }"/>
+
+    <SurfaceElementRegion
+      name="Fracture"
+      subRegionType="embeddedElement"
+      materialList="{ water, fractureFilling, fractureContact }"
+      defaultAperture="1e-3"/>
+  </ElementRegions>
+
+  <Constitutive>
+    <PorousElasticIsotropic
+      name="porousRock"
+      solidModelName="rock"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"/>
+
+    <ElasticIsotropic
+      name="rock"
+      defaultDensity="2700"
+      defaultBulkModulus="16.66666666666666e9"
+      defaultShearModulus="1.0e10"/>
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-12, 1.0e-12, 1.0e-12 }"/>
+
+    <CompressibleSinglePhaseFluid
+      name="water"
+      defaultDensity="1000"
+      defaultViscosity="1.0e-3"
+      referencePressure="1e6"
+      compressibility="0.0e0"
+      referenceViscosity="1e-3"
+      viscosibility="0.0"/>
+
+    <BiotPorosity
+      name="rockPorosity"
+      grainBulkModulus="1.0e27"
+      defaultReferencePorosity="0.2"/>
+
+    <CompressibleSolidSlipDependentPermeability
+      name="fractureFilling"
+      solidModelName="nullSolid"
+      porosityModelName="fracturePorosity"
+      permeabilityModelName="fracturePerm"/>
+
+    <!--CompressibleSolidParallelPlatesPermeability
+      name="fractureFilling"
+      solidModelName="nullSolid"
+      porosityModelName="fracturePorosity"
+      permeabilityModelName="fracturePerm2"/-->  
+
+    <NullModel
+      name="nullSolid"/>
+
+    <PressurePorosity
+      name="fracturePorosity"
+      defaultReferencePorosity="1.00"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+
+    <SlipDependentPermeability
+      name="fracturePerm"
+      shearDispThreshold="0.05"
+      maxPermMultiplier="100"
+      initialPermeability="{1e-15, 1e-15, 1e-15}"/>
+
+    <ParallelPlatesPermeability
+      name="fracturePerm2"/>  
+
+    <FrictionlessContact
+      name="fractureContact"
+      penaltyStiffness="2.0e11"
+      apertureTableName="aperTable"/>
+  </Constitutive>
+
+  <Functions>
+    <TableFunction
+        name="aperTable"
+        coordinates="{ -1.0e-3, 0.0 }"
+        values="{ 1.0e-6, 1.0e-4 }"/>    
+  </Functions>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="fracPorosity"
+      initialCondition="0"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fracture"
+      fieldName="fracturePorosity_porosity"
+      scale="1.00"/>
+
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Domain/cb1"
+      fieldName="pressure"
+      scale="0.0"/>
+
+    <FieldSpecification
+      name="initialPressureFracture"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="pressure"
+      scale="0.0"/>
+
+    <FieldSpecification
+      name="xnegconstraint"
+      objectPath="nodeManager"
+      fieldName="TotalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ bottomPoint, topPoint }"/>
+
+    <FieldSpecification
+      name="yposconstraint"
+      objectPath="nodeManager"
+      fieldName="TotalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ leftPoint, rightPoint }"/>
+
+    <FieldSpecification
+      name="zconstraint"
+      objectPath="nodeManager"
+      fieldName="TotalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zneg, zpos }"/>
+
+    <Traction
+      name="xload1"
+      objectPath="faceManager"
+      tractionType="vector"
+      direction="{1,0,0}"
+      scale="-1.0e7"
+      setNames="{ xpos }"/>
+
+    <Traction
+      name="xload2"
+      objectPath="faceManager"
+      tractionType="vector"
+      direction="{1,0,0}"
+      scale="+1.0e7"
+      setNames="{ xneg }"/> 
+
+   <FieldSpecification
+      name="boundaryPressure"
+      objectPath="faceManager"
+      fieldName="pressure"
+      scale="0.0"
+      setNames="{ ypos, yneg }"/>    
+  </FieldSpecifications>
+
+  <Tasks>
+    <PackCollection
+      name="tractionCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="fractureTraction"/>
+
+    <PackCollection
+      name="displacementJumpCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="displacementJump"/>  
+
+    <PackCollection
+      name="fracPermeabilityCollection"
+      objectPath="ElementRegions/Fracture/embeddedSurfaceSubRegion"
+      fieldName="fracturePerm_permeability"/>             
+  </Tasks>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotLevel="3"
+      writeBinaryData="1"/>
+
+   <TimeHistory
+      name="timeHistoryOutput"
+      sources="{/Tasks/tractionCollection}"
+      filename="traction_history" />
+
+    <TimeHistory
+      name="timeHistoryOutput1"
+      sources="{/Tasks/displacementJumpCollection}"
+      filename="displacementJump_history" />
+
+    <TimeHistory
+      name="timeHistoryOutput2"
+      sources="{/Tasks/fracPermeabilityCollection}"
+      filename="fracPermeability_history" />
+  </Outputs>
+</Problem>