adding integrated test

GEOS-DEV · Aug 12, 2024 · f8df31b · f8df31b
1 parent 798ee87
commit f8df31b
Show file tree

Hide file tree

Showing 2 changed files with 280 additions and 0 deletions.
diff --git a/inputFiles/compositionalMultiphaseFlow/4comp_2ph_1d_DD.xml b/inputFiles/compositionalMultiphaseFlow/4comp_2ph_1d_DD.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Solvers>
+    <CompositionalMultiphaseFVM
+      name="compflow"
+      logLevel="1"
+      discretization="fluidTPFA"
+      targetRegions="{ Region1 }"
+      initialDt="1e3"
+      targetFlowCFL="2"
+      temperature="297.15">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-6"
+        newtonMaxIter="15"/>
+      <LinearSolverParameters
+        directParallel="0"/>
+    </CompositionalMultiphaseFVM>
+  </Solvers>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0, 10 }"
+      yCoords="{ 0, 1 }"
+      zCoords="{ 0, 1 }"
+      nx="{ 100 }"
+      ny="{ 1 }"
+      nz="{ 1 }"
+      cellBlockNames="{ block1 }"/>
+  </Mesh>
+
+  <Geometry>
+    <Box
+      name="source"
+      xMin="{ -0.01, -0.01, -0.01 }"
+      xMax="{ 0.11, 1.01, 1.01 }"/>
+
+    <Box
+      name="sink"
+      xMin="{ 9.9, -0.01, -0.01 }"
+      xMax="{ 10.1, 1.01, 1.01 }"/>
+  </Geometry>
+
+  <Events
+    maxTime="2e6">
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="2e5"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="solverApplications1"
+      endTime="2e6"
+      target="/Solvers/compflow"/>
+
+    <PeriodicEvent
+      name="restarts"
+      timeFrequency="1e6"
+      targetExactTimestep="0"
+      target="/Outputs/restartOutput"/>
+  </Events>
+
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="fluidTPFA"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Region1"
+      cellBlocks="{ block1 }"
+      materialList="{ fluid1, rock, relperm, diffusion, dispersion}"/>
+  </ElementRegions>
+
+  <Constitutive>
+    <CompositionalMultiphaseFluid
+      name="fluid1"
+      phaseNames="{ oil, gas }"
+      equationsOfState="{ PR, PR }"
+      componentNames="{ N2, C10, C20, H2O }"
+      componentCriticalPressure="{ 34e5, 25.3e5, 14.6e5, 220.5e5 }"
+      componentCriticalTemperature="{ 126.2, 622.0, 782.0, 647.0 }"
+      componentAcentricFactor="{ 0.04, 0.443, 0.816, 0.344 }"
+      componentMolarWeight="{ 28e-3, 134e-3, 275e-3, 18e-3 }"
+      componentVolumeShift="{ 0, 0, 0, 0 }"
+      componentBinaryCoeff="{ { 0, 0, 0, 0 },
+                              { 0, 0, 0, 0 },
+                              { 0, 0, 0, 0 },
+                              { 0, 0, 0, 0 } }"/>
+
+    <CompressibleSolidConstantPermeability
+      name="rock"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"/>
+
+    <NullModel
+      name="nullSolid"/>
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.2"
+      referencePressure="0.0"
+      compressibility="1.0e-9"/>
+
+    <BrooksCoreyRelativePermeability
+      name="relperm"
+      phaseNames="{ oil, gas }"
+      phaseMinVolumeFraction="{ 0.1, 0.15 }"
+      phaseRelPermExponent="{ 2.0, 2.0 }"
+      phaseRelPermMaxValue="{ 0.8, 0.9 }"/>
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-16, 1.0e-16, 1.0e-16 }"/>
+
+    <ConstantDiffusion
+      name="diffusion"
+      phaseNames="{ oil, gas }"
+      defaultPhaseDiffusivityMultipliers="{ 0.1, 0.3 }"
+      diffusivityComponents="{ 1e-5, 1e-5, 1e-5 }"/>
+
+    <LinearIsotropicDispersion
+      name="dispersion"
+      phaseNames="{ oil, gas }"
+      longitudinalDispersivity="1e-1" /> 
+
+  </Constitutive>
+
+  <FieldSpecifications>
+    <!-- Initial pressure: ~5 bar -->
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="pressure"
+      scale="5e6"/>
+
+    <!-- Initial composition: no water, only heavy hydrocarbon components and N2 -->
+    <FieldSpecification
+      name="initialComposition_N2"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="0"
+      scale="0.099"/>
+
+    <FieldSpecification
+      name="initialComposition_C10"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="1"
+      scale="0.001"/>
+
+    <FieldSpecification
+      name="initialComposition_C20"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="2"
+      scale="0.899"/>
+
+    <FieldSpecification
+      name="initialComposition_H20"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="3"
+      scale="0.001"/>
+
+    <!-- Injection pressure: ~10 bar -->
+    <FieldSpecification
+      name="sourceTermPressure"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="pressure"
+      scale="1e7"
+      setNames="{ source }"/>
+
+    <!-- Injection stream: mostly water -->
+    <FieldSpecification
+      name="sourceTermComposition_N2"
+      setNames="{ source }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="0"
+      scale="0.001"/>
+
+    <FieldSpecification
+      name="sourceTermComposition_C10"
+      setNames="{ source }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="1"
+      scale="0.997"/>
+
+    <FieldSpecification
+      name="sourceTermComposition_C20"
+      setNames="{ source }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="2"
+      scale="0.001"/>
+
+    <FieldSpecification
+      name="sourceTermComposition_H20"
+      setNames="{ source }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="3"
+      scale="0.001"/>
+
+    <!-- Production pressure: ~1 bar, close to stock tank conditions -->
+    <FieldSpecification
+      name="sinkTerm"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="pressure"
+      scale="1e5"
+      setNames="{ sink }"/>
+
+    <!-- Production stream: same as initial (should not matter due to upwinding) -->
+    <FieldSpecification
+      name="sinkTermComposition_N2"
+      setNames="{ sink }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="0"
+      scale="0.099"/>
+
+    <FieldSpecification
+      name="sinkTermComposition_C10"
+      setNames="{ sink }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="1"
+      scale="0.001"/>
+
+    <FieldSpecification
+      name="sinkTermComposition_C20"
+      setNames="{ sink }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="2"
+      scale="0.899"/>
+
+    <FieldSpecification
+      name="sinkTermComposition_H20"
+      setNames="{ sink }"
+      objectPath="ElementRegions/Region1/block1"
+      fieldName="globalCompFraction"
+      component="3"
+      scale="0.001"/>
+  </FieldSpecifications>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      plotFileRoot="4comp_2ph_1d"/>
+
+    <Restart
+      name="restartOutput"/>
+  </Outputs>
+</Problem>
diff --git a/inputFiles/compositionalMultiphaseFlow/compositionalMultiphaseFlow.ats b/inputFiles/compositionalMultiphaseFlow/compositionalMultiphaseFlow.ats
@@ -21,6 +21,14 @@ decks = [
         restart_step=118,
         check_step=218,
         restartcheck_params=RestartcheckParameters(**restartcheck_params)),
+    TestDeck(
+        name="4comp_2ph_cap_1d_DD",
+        description=
+            "Compositional multiphase flow test (1D displacement, 2-phase 4-component, capillary pressure) with dispersion and diffusion",
+        partitions=((1, 1, 1), (2, 1, 1), (3, 1, 1)),
+        restart_step=118,
+        check_step=218,
+        restartcheck_params=RestartcheckParameters(**restartcheck_params)),
     TestDeck(
         name="deadoil_3ph_corey_1d",
         description=