Capillary interface conditions

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/immiscibleTwoPhase_GravitySegregation_1d.ats

@@ -0,0 +1,18 @@
+from geos.ats.test_builder import TestDeck, RestartcheckParameters, generate_geos_tests
+
+restartcheck_params = {}
+restartcheck_params['atol'] = 1.0E-8
+restartcheck_params['rtol'] = 1.0E-8
+
+decks = [
+    TestDeck(
+        name="immiscibleTwoPhase_GravitySegregation_1d",
+        description=
+        'Test that 2 fluids can be seperated based on gravity.',
+        partitions=((1, 1, 1), ),
+        restart_step=0,
+        check_step=45,
+        restartcheck_params=RestartcheckParameters(**restartcheck_params))
+]
+
+generate_geos_tests(decks)

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/immiscibleTwoPhase_GravitySegregation_1d.xml

@@ -0,0 +1,203 @@
+<?xml version="1.0"?>
+
+<Problem>
+  <Solvers
+    gravityVector="{ 0.0, 0.0, -10 }">
+    <ImmiscibleMultiphaseFlow
+      name="FlowSolver"
+      discretization="TPFA"
+      targetRegions="{ Domain }"
+      writeLinearSystem="0"
+      temperature="300"
+      initialDt="0.001">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-6"
+        lineSearchAction="None"
+        newtonMaxIter="10" />
+      <LinearSolverParameters
+        directParallel="0" />
+    </ImmiscibleMultiphaseFlow>
+  </Solvers>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0, 1 }"
+      yCoords="{ 0, 1 }"
+      zCoords="{ 0, 10 }"
+      nx="{ 1 }"
+      ny="{ 1 }"
+      nz="{ 10 }"
+      cellBlockNames="{ block1 }" />
+  </Mesh>
+
+  <Geometry>
+    <Box
+      name="sink"
+      xMin="{ -0.01, -0.01, -0.01 }"
+      xMax="{ 1.01, 1.01, 1.01 }" />
+  </Geometry>
+
+  <Events
+    maxTime="0.4">
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="0.01"
+      target="/Outputs/vtkOutput" />
+
+    <PeriodicEvent
+      name="solverApplications"
+      target="/Solvers/FlowSolver" />
+
+    <PeriodicEvent
+      name="restarts"
+      timeFrequency="0.1"
+      targetExactTimestep="0"
+      target="/Outputs/restartOutput"/>
+  </Events>
+
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="TPFA" />
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="Domain"
+      cellBlocks="{ block1 }"
+      materialList="{ fluid, rock, relperm }" />
+  </ElementRegions>
+
+  <Functions>
+    <TableFunction
+      name="densityTablePhase1"
+      coordinates="{ 0.0}"
+      values="{ 10 }" />
+
+    <TableFunction
+      name="densityTablePhase2"
+      coordinates="{ 0.0 }"
+      values="{ 1 }" />
+
+    <TableFunction
+      name="viscosityTablePhase1"
+      coordinates="{ 0.0 }"
+      values="{ 1.0 }" />
+
+    <TableFunction
+      name="viscosityTablePhase2"
+      coordinates="{ 0.0 }"
+      values="{ 1.0 }" />
+
+    <TableFunction
+      name="saturationTable1"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ x.txt, y.txt, z.txt }"
+      voxelFile="initialSaturation1.txt"
+      interpolation="nearest" />
+
+    <TableFunction
+      name="saturationTable2"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ x.txt, y.txt, z.txt }"
+      voxelFile="initialSaturation2.txt"
+      interpolation="nearest" />
+
+    <TableFunction
+      name="pressureTable"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ x.txt, y.txt, z.txt }"
+      voxelFile="initialPressure.txt"
+      interpolation="nearest" />
+  </Functions>
+
+  <Constitutive>
+    <TwoPhaseImmiscibleFluid
+      name="fluid"
+      phaseNames="{ water, gas }"
+      densityTableNames="{ densityTablePhase1, densityTablePhase2 }"
+      viscosityTableNames="{ viscosityTablePhase1, viscosityTablePhase2 }" />
+
+    <CompressibleSolidConstantPermeability
+      name="rock"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm" />
+
+    <NullModel
+      name="nullSolid" />
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="1.0"
+      referencePressure="0.0"
+      compressibility="1e-9" />
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1, 1, 1 }" />
+
+    <BrooksCoreyRelativePermeability
+      name="relperm"
+      phaseNames="{ water, oil }"
+      phaseMinVolumeFraction="{ 0.0, 0.0 }"
+      phaseRelPermExponent="{ 1.0, 1.0 }"
+      phaseRelPermMaxValue="{ 1.0, 1.0 }" />
+  </Constitutive>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="Porosity"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Domain/block1"
+      fieldName="rockPorosity_referencePorosity"
+      scale="1.0" />
+
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Domain/block1"
+      fieldName="pressure"
+      scale="1.0"
+      functionName="pressureTable" />
+
+    <FieldSpecification
+      name="initialSat1"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Domain/block1"
+      fieldName="phaseVolumeFraction"
+      component="0"
+      scale="1.0"
+      functionName="saturationTable1" />
+
+    <FieldSpecification
+      name="initialSat2"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Domain/block1"
+      fieldName="phaseVolumeFraction"
+      component="1"
+      scale="1.0"
+      functionName="saturationTable2" />
+
+    <FieldSpecification
+      name="sinkTermP"
+      objectPath="ElementRegions/Domain/block1"
+      fieldName="pressure"
+      scale="500.0"
+      setNames="{ sink }" />
+  </FieldSpecifications>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput" />
+    <Restart
+      name="restartOutput"/>      
+  </Outputs>
+</Problem>

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/initialPressure.txt

@@ -0,0 +1,10 @@
+495.0
+485.0
+475.0
+465.0
+455.0
+400.0
+300.0
+200.0
+100.0
+0.0

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/initialSaturation1.txt

@@ -0,0 +1,10 @@
+0
+0
+0
+0
+0
+1
+1
+1
+1
+1

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/initialSaturation2.txt

@@ -0,0 +1,10 @@
+1
+1
+1
+1
+1
+0
+0
+0
+0
+0

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/x.txt

@@ -0,0 +1 @@
+0.5

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/y.txt

@@ -0,0 +1 @@
+0.5

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_GravitySegregation_1d/z.txt

@@ -0,0 +1,10 @@
+0.5
+1.5
+2.5
+3.5
+4.5
+5.5
+6.5
+7.5
+8.5
+9.5

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/initialPressure.txt

@@ -0,0 +1,10 @@
+1e7
+1e7
+1e7
+1e7
+1e7
+1e7
+1e7
+1e7
+1e7
+1e7

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/initialSaturation1.txt

@@ -0,0 +1,10 @@
+0.95
+0.95
+0.95
+0.95
+0.95
+0.95
+0.95
+0.05
+0.05
+0.05

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/initialSaturation2.txt

@@ -0,0 +1,10 @@
+0.05
+0.05
+0.05
+0.05
+0.05
+0.05
+0.05
+0.95
+0.95
+0.95

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/jFunction_linear.txt

@@ -0,0 +1,2 @@
+4.33172935918785
+1.66604975353379

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/permx.geos

@@ -0,0 +1,10 @@
+50
+50
+50
+50
+50
+200
+200
+200
+200
+200

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/permy.geos

@@ -0,0 +1,10 @@
+50
+50
+50
+50
+50
+200
+200
+200
+200
+200

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/permz.geos

@@ -0,0 +1,10 @@
+50
+50
+50
+50
+50
+200
+200
+200
+200
+200

inputFiles/immiscibleMultiphaseFlow/immiscibleTwoPhase_interface_conditions/1D_case/phaseVolumeFraction_water_linear.txt

@@ -0,0 +1,2 @@
+0
+1