refactor: Preparatory for well to frac connection

inputFiles/poromechanicsFractures/singlePhasePoromechanics_FaultModel_well_fim_new_smoke.xml

@@ -0,0 +1,116 @@
+<Problem>
+  <Included>
+    <File name="./singlePhasePoromechanics_FaultModel_well_new_base.xml"/>
+  </Included>
+
+  <Solvers
+    gravityVector="{ 0.0, 0.0, 0.0 }">
+
+    <SinglePhasePoromechanicsConformingFracturesReservoir
+      name="reservoirSolver"
+      targetRegions="{ Region, Fault, wellRegion1, wellRegion2 }" 
+      initialDt="1e-3"
+      poromechanicsConformingFracturesSolverName="poroFractureSolver"
+      wellSolverName="singlePhaseWellSolver"
+      logLevel="1">
+      <NonlinearSolverParameters
+        lineSearchAction="None"
+        maxTimeStepCuts="10"
+        newtonTol="1e-3"
+        newtonMaxIter="20"
+        maxAllowedResidualNorm="1e+25"/>
+      <LinearSolverParameters
+        solverType="direct"/> 
+    </SinglePhasePoromechanicsConformingFracturesReservoir>
+
+    <SinglePhasePoromechanicsConformingFractures
+      name="poroFractureSolver"
+      targetRegions="{ Region, Fault }" 
+      flowSolverName="singlePhaseFlowSolver"
+      solidSolverName="fractureMechSolver"
+      logLevel="1">
+    </SinglePhasePoromechanicsConformingFractures>
+
+    <SolidMechanicsLagrangeContact
+      name="fractureMechSolver"
+      stabilizationName="singlePhaseTPFA_stab"
+      logLevel="1"
+      targetRegions="{ Region, Fault }" 
+      timeIntegrationOption="QuasiStatic"
+      discretization="FE1"/>  
+
+    <SinglePhaseFVM 
+      name="singlePhaseFlowSolver"
+      logLevel="1"
+      allowNegativePressure="1"
+      discretization="singlePhaseTPFA"
+      targetRegions="{ Region, Fault}" />
+
+    <SinglePhaseWell
+      name="singlePhaseWellSolver"
+      logLevel="1"
+      targetRegions="{ wellRegion1, wellRegion2 }">
+      <WellControls
+        name="wellControls1"
+        logLevel="2"
+        type="producer"
+        control="BHP"
+        referenceElevation="-2500"
+        targetBHP="1e7"
+        targetTotalRate="1e10" />
+      <WellControls
+        name="wellControls2"
+        logLevel="2"
+        type="injector"
+        control="BHP"
+        referenceElevation="-2500"
+        targetBHP="10e7"
+        targetTotalRate="1e10" />
+    </SinglePhaseWell>
+  </Solvers>
+
+  <Tasks>
+      <SinglePhasePoromechanicsConformingFracturesInitialization 
+        logLevel="1" 
+        name="SinglePhasePoroelasticityPreEquilibrationStep"  
+        poromechanicsSolverName="poroFractureSolver"/>
+  </Tasks>
+
+  <Events
+    minTime="-1e11"
+    maxTime="3.0e7">
+
+    <!-- Initialization -->
+    <SoloEvent 
+      name="SinglePhasePoroelasticityPreEquilibrationStep" 
+      targetTime="-1e11" 
+      beginTime="-1e11" 
+      target="/Tasks/SinglePhasePoroelasticityPreEquilibrationStep"
+    />
+    <!-- Injection -->
+    <SoloEvent
+      name="outputPostEquilibrationStep1"
+      targetTime="0"
+      targetExactTimestep="1"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="outputs"
+      beginTime="0.0"
+      endTime="3.0e7"
+      forceDt="3.0e7"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="solverApplication"
+      endTime="3.0e7"
+      forceDt="3.0e7"
+      target="/Solvers/reservoirSolver"/>
+
+    <PeriodicEvent
+      name="restarts"
+      beginTime="0.0"
+      endTime="3.0e7"
+      target="/Outputs/restartOutput"/>
+  </Events>
+</Problem>

inputFiles/poromechanicsFractures/singlePhasePoromechanics_FaultModel_well_new_base.xml

@@ -0,0 +1,235 @@
+<Problem>
+  <Constitutive>
+    <PorousElasticIsotropic
+      name="porousRock"
+      permeabilityModelName="rockPerm"
+      porosityModelName="rockPorosity"
+      solidModelName="rock" />
+
+    <ElasticIsotropic
+      name="rock"
+      defaultDensity="2927"
+      defaultPoissonRatio="0.15"
+      defaultYoungModulus="14950000000" />
+
+    <CompressibleSinglePhaseFluid
+      name="water"
+      compressibility="2.09028227021e-10"
+      defaultDensity="1325"
+      defaultViscosity="0.001"
+      referenceDensity="1325"
+      referencePressure="0e6"
+      referenceViscosity="0.001"
+      viscosibility="0.0" />
+
+    <BiotPorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.3"
+      defaultGrainBulkModulus="7.12e30" />
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-18, 1.0e-18, 1.0e-18 }" />
+      <!-- Material inside the fault -->
+    <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="40.0e6"
+      frictionCoefficient="0.01"/>
+
+    <HydraulicApertureTable
+    name="hApertureModel"
+    apertureTableName="apertureTable"/>   
+
+    <NullModel
+      name="nullSolid"/>
+
+  </Constitutive>
+
+  <FieldSpecifications>
+
+    <FieldSpecification
+      name="fracPorosity"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Fault"
+      fieldName="fracturePorosity_porosity"
+      scale="1.00"/>
+
+    <FieldSpecification
+      name="xconstraint"
+      component="0"
+      fieldName="totalDisplacement"
+      objectPath="nodeManager"
+      scale="0.0"
+      setNames="{ xnegFace, xposFace }" />
+
+    <FieldSpecification
+      name="yconstraint"
+      component="1"
+      fieldName="totalDisplacement"
+      objectPath="nodeManager"
+      scale="0.0"
+      setNames="{ ynegFace, yposFace }" />
+
+    <FieldSpecification
+      name="zconstraint"
+      component="2"
+      fieldName="totalDisplacement"
+      objectPath="nodeManager"
+      scale="0.0"
+      setNames="{ znegFace }" /> 
+
+    <Traction
+      name="tractionTop"
+      objectPath="faceManager"
+      tractionType="normal"
+      scale="-4e7"      
+      setNames="{ zposFace }"/>
+
+    <FieldSpecification
+      name="init_pressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="pressure"
+      scale="4e7"/>
+  </FieldSpecifications>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Region"
+      cellBlocks="{ 2_hexahedra, 1_hexahedra, 0_hexahedra }"
+      materialList="{ water, porousRock }"/>  
+
+    <SurfaceElementRegion
+      name="Fault"
+      faceBlock="fracture"
+      materialList="{ water, faultFilling, fractureContact, hApertureModel }"
+      defaultAperture="1e-3"/>
+
+    <WellElementRegion
+      name="wellRegion1"
+      materialList="{ water }"/>
+
+    <WellElementRegion
+      name="wellRegion2"
+      materialList="{ water }"/>
+  </ElementRegions>
+
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1" />
+    </FiniteElements>
+
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"/>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA_stab"/>  
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput" />
+
+    <Restart
+      name="restartOutput"/>
+  </Outputs>
+
+  <Functions>
+    <TableFunction
+      name="apertureTable"
+      coordinates="{ -1.0e-3, 0.0 }"
+      values="{ 1.0e-3, 1.0e-3 }"
+    />
+  </Functions>
+
+  <Mesh>
+    <VTKMesh
+      name="mesh1"
+      useGlobalIds="1"
+      faceBlocks="{ fracture }"
+      file="verticalFault_ExternalMesh.vtm">
+      <InternalWell
+        name="well_producer1"
+        logLevel="1"
+        wellRegionName="wellRegion1"
+        wellControlsName="wellControls1"
+        polylineNodeCoords="{ { 5400, 5400, 0 },
+                              { 5400, 5400, -2500 } }"
+        polylineSegmentConn="{ { 0, 1 } }"
+        radius="0.1"
+        numElementsPerSegment="1">
+        <Perforation
+          name="producer1_perf1"
+          targetRegion="Region"
+          distanceFromHead="2500"/>
+      </InternalWell>
+      <InternalWell
+        name="well_injector1"
+        logLevel="1"
+        wellRegionName="wellRegion2"
+        wellControlsName="wellControls2"
+        polylineNodeCoords="{ { -5400, -5400, 0 },
+                              { -5400, -5400, -2500 } }"
+        polylineSegmentConn="{ { 0, 1 } }"
+        radius="0.1"
+        numElementsPerSegment="1">
+        <Perforation
+          name="injector1_perf1"
+          distanceFromHead="2500"/>
+      </InternalWell>
+    </VTKMesh>
+  </Mesh>
+
+  <Geometry>
+    <Box
+      name="xposFace"
+      xMin="{ 5999.99,  -6000.01,  -4000.01 }"
+      xMax="{ 6000.01,   6000.01,      0.01 }"/>
+
+    <Box
+      name="xnegFace"
+      xMin="{-6000.01,  -6000.01,  -4000.01 }"
+      xMax="{-5999.99,   6000.01,      0.01 }"/>
+
+    <Box
+      name="yposFace"
+      xMin="{-6000.01,   5999.99,  -4000.01 }"
+      xMax="{ 6000.01,   6000.01,      0.01 }"/>
+
+    <Box
+      name="ynegFace"
+      xMin="{-6000.01,  -6000.01,  -4000.01 }"
+      xMax="{ 6000.01,  -5999.99,      0.01 }"/>
+
+    <Box
+      name="zposFace"
+      xMin="{-6000.01,  -6000.01,     -0.01 }"
+      xMax="{ 6000.01,   6000.01,      0.01 }"/>
+
+    <Box
+      name="znegFace"
+      xMin="{-6000.01,  -6000.01,  -4000.01 }"
+      xMax="{ 6000.01,   6000.01,  -3999.99 }"/>
+  </Geometry>
+
+</Problem>