fractureMatrixFlow_conforming_2d.xml

<?xml version="1.0" ?>

<Problem>
  <Solvers
    gravityVector="{ 0.0, 0.0, -9.81 }">
    <SinglePhaseFVM
      name="SinglePhaseFlow"
      logLevel="1"
      discretization="singlePhaseTPFA"
      targetRegions="{ Region2, Fracture }">
      <NonlinearSolverParameters
        newtonTol="1.0e-2"
        newtonMaxIter="8"/>
      <LinearSolverParameters
        directParallel="0"/>
    </SinglePhaseFVM>

    <SurfaceGenerator
      name="SurfaceGen"
      targetRegions="{ Fracture }"
      initialRockToughness="1e6"
      mpiCommOrder="1"/>
  </Solvers>

  <Mesh>
    <InternalMesh
      name="mesh1"
      elementTypes="{ C3D8 }"
      xCoords="{ -5, 5 }"
      yCoords="{ 0, 10 }"
      zCoords="{ 0, 1 }"
      nx="{ 10 }"
      ny="{ 10 }"
      nz="{ 1 }"
      cellBlockNames="{ cb1 }"/>
  </Mesh>

  <Geometry>
    <Box
      name="left"
      xMin="{ -0.01, -0.01, -0.01 }"
      xMax="{ 1.01, 1.01, 1.01 }"/>

    <Box
      name="right"
      xMin="{ 8.99, -0.01, -0.01 }"
      xMax="{ 10.01, 1.01, 1.01 }"/>

    <Box
      name="bottom"
      xMin="{ -0.01, -0.01, -0.01 }"
      xMax="{ 10.01, 1.11, 1.01 }"/>

    <Box
      name="top"
      xMin="{ -0.01, 8.99, -0.01 }"
      xMax="{ 10.01, 10.01, 1.01 }"/>

    <Box
      name="fracture"
      xMin="{ -0.01, -0.01, -0.01 }"
      xMax="{ 0.01, 10.01, 1.01 }"/>

    <Box
      name="source"
      xMin="{ -0.01, -0.01, -0.01 }"
      xMax="{ 0.01, 1.01, 1.01 }"/>
  </Geometry>

  <Events
    maxTime="50000">
    <!--SoloEvent
            name="output0"
            target="/Outputs/siloOutput" / -->
    <SoloEvent
      name="preFracture"
      target="/Solvers/SurfaceGen"/>

    <!--SoloEvent
            name="output1"
            target="/Outputs/siloOutput" / -->
    <PeriodicEvent
      name="solverApplications"
      forceDt="1e3"
      target="/Solvers/SinglePhaseFlow"/>

    <PeriodicEvent
      name="outputs"
      target="/Outputs/siloOutput"/>

    <PeriodicEvent
      name="restarts"
      timeFrequency="25000"
      targetExactTimestep="0"
      target="/Outputs/restartOutput"/>
  </Events>

  <NumericalMethods>
    <FiniteVolume>
      <TwoPointFluxApproximation
        name="singlePhaseTPFA"/>
    </FiniteVolume>
  </NumericalMethods>

  <ElementRegions>
    <CellElementRegion
      name="Region2"
      cellBlocks="{ * }"
      materialList="{ water, rock }"/>

    <SurfaceElementRegion
      name="Fracture"
      defaultAperture="1.0e-5"
      materialList="{ water, fractureFilling }"/>
  </ElementRegions>

  <Constitutive>
    <CompressibleSinglePhaseFluid
      name="water"
      defaultDensity="1000"
      defaultViscosity="0.001"
      referencePressure="0.0"
      compressibility="5e-10"
      viscosibility="0.0"/>

    <CompressibleSolidParallelPlatesPermeability
      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.01"
      referencePressure="0.0"
      compressibility="1.0e-9"/>

    <ConstantPermeability
      name="rockPerm"
      permeabilityComponents="{ 5.0e-19, 5.0e-19, 5.0e-19 }"/>

    <PressurePorosity
      name="fracturePorosity"
      defaultReferencePorosity="1.00"
      referencePressure="0.0"
      compressibility="0.0"/>

    <ParallelPlatesPermeability
      name="fracturePerm"/>
  </Constitutive>

  <FieldSpecifications>
    <FieldSpecification
      name="frac"
      initialCondition="1"
      setNames="{ fracture }"
      objectPath="faceManager"
      fieldName="ruptureState"
      scale="1"/>

    <FieldSpecification
      name="Porosity"
      initialCondition="1"
      setNames="{ all }"
      objectPath="ElementRegions/Region2/cb1"
      fieldName="rockPorosity_referencePorosity"
      scale="0.01"/>

    <FieldSpecification
      name="fracPorosity"
      initialCondition="1"
      setNames="{ all }"
      objectPath="ElementRegions/Fracture"
      fieldName="fracturePorosity_porosity"
      scale="1.00"/>

    <FieldSpecification
      name="initialPressure"
      initialCondition="1"
      setNames="{ all }"
      objectPath="ElementRegions/Region2/cb1"
      fieldName="pressure"
      scale="0.0"/>

    <FieldSpecification
      name="source"
      objectPath="ElementRegions/Fracture"
      fieldName="pressure"
      scale="2.0e6"
      setNames="{ source }"/>
  </FieldSpecifications>

  <Outputs>
    <VTK
      name="siloOutput"
      plotFileRoot="fractureMatrixFlow"/>

    <Restart
      name="restartOutput"/>
  </Outputs>
</Problem>