Skip to content

Commit

Permalink
Modified parallel plates permeability and updated discretization (#2623)
Browse files Browse the repository at this point in the history
---------

Co-authored-by: Matteo Cusini <[email protected]>
  • Loading branch information
karimifard and CusiniM authored Sep 14, 2023
1 parent 3cd5d63 commit 7f0065a
Show file tree
Hide file tree
Showing 28 changed files with 655 additions and 101 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -14,9 +14,9 @@
targetPhaseVolFractionChangeInTimeStep="1"
maxCompFractionChange="1"
maxRelativePressureChange="1"
initialDt="1e4">
initialDt="1e2">
<NonlinearSolverParameters
newtonTol="1.0e-6"
newtonTol="1.0e-5"
newtonMaxIter="15"
maxTimeStepCuts="2"
lineSearchMaxCuts="2"/>
Expand Down Expand Up @@ -62,20 +62,20 @@
</NumericalMethods>

<Events
maxTime="5e5">
maxTime="1e5">
<SoloEvent
name="preFracture"
target="/Solvers/SurfaceGenerator"/>

<PeriodicEvent
name="outputs"
timeFrequency="1e5"
timeFrequency="1e4"
targetExactTimestep="1"
target="/Outputs/vtkOutput"/>

<PeriodicEvent
name="restarts"
timeFrequency="2e5"
timeFrequency="0.5e5"
targetExactTimestep="0"
target="/Outputs/sidreRestart"/>

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -108,7 +108,7 @@

<ConstantPermeability
name="rockPerm"
permeabilityComponents="{ 1.0e-15, 1.0e-15, 1.0e-15 }"/>
permeabilityComponents="{ 1.0e-12, 1.0e-12, 1.0e-12 }"/>

<CompressibleSinglePhaseFluid
name="water"
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -13,16 +13,21 @@
xCoords="{ -1, 1 }"
yCoords="{ -1, 1 }"
zCoords="{ 0, 1 }"
nx="{ 2 }"
ny="{ 4 }"
nx="{ 8 }"
ny="{ 10 }"
nz="{ 1 }"
cellBlockNames="{ cb1 }"/>
</Mesh>

<Geometry>
<Box
name="fracture"
name="fractureTip"
xMin="{ -0.01, -0.8, -0.01 }"
xMax="{ 0.01, -0.59, 1.01 }"/>

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

<Box name="bottomFace"
Expand All @@ -32,7 +37,6 @@
<Box name="topFace"
xMax="{1.01, 1.01, 1.01 }"
xMin="{-1.01, -1.01, 0.99 }"/>

</Geometry>

<Events
Expand Down Expand Up @@ -65,15 +69,15 @@
<FieldSpecification
name="frac"
initialCondition="1"
setNames="{ fracture }"
setNames="{ fracture, fractureTip }"
objectPath="faceManager"
fieldName="ruptureState"
scale="1"/>

<FieldSpecification
name="separableFace"
initialCondition="1"
setNames="{ fracture }"
setNames="{ fracture, fractureTip}"
objectPath="faceManager"
fieldName="isFaceSeparable"
scale="1"/>
Expand All @@ -88,11 +92,10 @@

<FieldSpecification
name="prescribedFracturePressure"
setNames="{ all }"
setNames="{ fractureTip }"
objectPath="ElementRegions/Fracture"
fieldName="pressure"
scale="1.0e5"
functionName="injectionSchedule"/>
scale="1.0e5"/>

<FieldSpecification
name="initialStress"
Expand All @@ -101,7 +104,7 @@
objectPath="ElementRegions/Domain"
fieldName="rock_stress"
component="0"
scale="-1.0e5"/>
scale="-1.0e0"/>

<FieldSpecification
name="xconstraint"
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -60,6 +60,10 @@
</Mesh>

<Geometry>
<Box
name="fractureTip"
xMin="{ 1., -0.01, -0.1 }"
xMax="{ 2., .3, 1.1 }"/>
<Rectangle
name="FracturePlane"
normal="{ 1, 0, 0 }"
Expand Down Expand Up @@ -192,7 +196,7 @@

<FieldSpecification
name="initialPressureFracture"
setNames="{ all }"
setNames="{ fractureTip }"
objectPath="ElementRegions/Fracture"
fieldName="pressure"
scale="1.0e5"/>
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,136 @@
<?xml version="1.0" ?>

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

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

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

<ElementRegions>
<CellElementRegion
name="RockMatrix"
cellBlocks="{ cb1 }"
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"/>

<CompressibleSolidConstantPermeability
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-12, 5.0e-12, 5.0e-12 }"/>

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

<ConstantPermeability
name="fracturePerm"
permeabilityComponents="{ 1.0e-18, 1.0e-18, 1.0e-18 }"/>
</Constitutive>

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

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

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

<FieldSpecification
name="leftPressure"
objectPath="faceManager"
fieldName="pressure"
scale="1.0e6"
setNames="{ xneg }"/>

<FieldSpecification
name="rightPressure"
objectPath="faceManager"
fieldName="pressure"
scale="0.0e6"
setNames="{ xpos }"/>
</FieldSpecifications>

<Outputs>
<VTK
name="vtkOutput"
plotFileRoot="fractureMatrixFlow_impermeableFault"/>
<Restart name="restartOutput"/>
</Outputs>
</Problem>
Original file line number Diff line number Diff line change
@@ -0,0 +1,43 @@
<?xml version="1.0" ?>

<Problem>
<Included>
<File name="./impermeableFault_conforming_base.xml" />
</Included>

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

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

<Events
maxTime="5000">
<SoloEvent
name="preFracture"
target="/Solvers/SurfaceGen"/>

<PeriodicEvent
name="solverApplications"
forceDt="1e3"
target="/Solvers/SinglePhaseFlow"/>

<PeriodicEvent
name="outputs"
target="/Outputs/vtkOutput"/>
</Events>
</Problem>
Original file line number Diff line number Diff line change
@@ -0,0 +1,48 @@
<?xml version="1.0" ?>

<Problem>
<Included>
<File name="./impermeableFault_conforming_base.xml" />
</Included>

<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="fracture"
xMin="{ -0.01, 2.5, -0.01 }"
xMax="{ 0.01, 7.5, 1.01 }"/>
</Geometry>

<Events
maxTime="5000">
<SoloEvent
name="preFracture"
target="/Solvers/SurfaceGen"/>

<PeriodicEvent
name="solverApplications"
forceDt="1e3"
target="/Solvers/SinglePhaseFlow"/>

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

<PeriodicEvent
name="restarts"
timeFrequency="2500"
target="/Outputs/restartOutput"/>
</Events>
</Problem>
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@
isThermal="1">
<NonlinearSolverParameters
newtonTol="1.0e-5"
newtonMaxIter="10"
newtonMaxIter="20"
maxTimeStepCuts="4"/>
<LinearSolverParameters
directParallel="0"
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,7 @@

<PeriodicEvent
name="restarts"
timeFrequency="1.0e6"
timeFrequency="5.0e5"
target="/Outputs/restartOutput"/>
</Events>
</Problem>
Original file line number Diff line number Diff line change
Expand Up @@ -47,7 +47,7 @@

<PeriodicEvent
name="restarts"
timeFrequency="1e6"
timeFrequency="5e5"
target="/Outputs/restartOutput"/>
</Events>

Expand Down
2 changes: 1 addition & 1 deletion integratedTests
Submodule integratedTests updated 378 files
Loading

0 comments on commit 7f0065a

Please sign in to comment.