feat: Geothermal Gradient for Single Phase flow

.integrated_tests.yaml

@@ -1,6 +1,6 @@
 baselines:
   bucket: geosx
-  baseline: integratedTests/baseline_integratedTests-pr3359-14024-3c9ebb4
+  baseline: integratedTests/baseline_integratedTests-pr3193-14118-8ee1c34
 
 allow_fail:
   all: ''

BASELINE_NOTES.md

@@ -6,6 +6,10 @@ This file is designed to track changes to the integrated test baselines.
 Any developer who updates the baseline ID in the .integrated_tests.yaml file is expected to create an entry in this file with the pull request number, date, and their justification for rebaselining.
 These notes should be in reverse-chronological order, and use the following time format: (YYYY-MM-DD).
 
+PR #3193 (2025-10-10) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr3193-14118-8ee1c34.tar.gz>
+=====================
+Enable geothermal gradient in HydrostaticEquilibrium for single-phase flow.
+
 PR #3359 (2025-10-07) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr3359-14024-3c9ebb4.tar.gz>
 =====================
 Add functions to connect well perforation to surface elements.

inputFiles/poromechanics/PoroElastic_hybridHexPrism_co2_base.xml

@@ -44,7 +44,7 @@
 
     <PeriodicEvent
       name="outputs"
-      timeFrequency="5e3"
+      timeFrequency="1e4"
       target="/Outputs/vtkOutput"/>
     <PeriodicEvent
       name="linearElasticityStatistics"

inputFiles/poromechanics/poroElasticCoupling.ats

@@ -265,7 +265,7 @@ def _build_PoroElasticPEBICO2FIM_cases():
         "Francois Hamon", True)
     restartcheck_params = {"atol": 1.0e-5, "rtol": 2.0e-7}
     return _build_poro_elastic_coupling_case(
-        "PoroElastic_hybridHexPrism_co2_fim_3d.xml", (10, 20),
+        "PoroElastic_hybridHexPrism_co2_fim_3d.xml", (5, 10),
         ((1, 1, 1), (2, 2, 1)), description, restartcheck_params)
 
 
@@ -275,7 +275,7 @@ def _build_PoroElasticPEBICO2Sequential_cases():
         "auto", "Francois Hamon", True)
     restartcheck_params = {"atol": 1.0e-5, "rtol": 2.0e-7}
     return _build_poro_elastic_coupling_case(
-        "PoroElastic_hybridHexPrism_co2_sequential_3d.xml", (10, 20),
+        "PoroElastic_hybridHexPrism_co2_sequential_3d.xml", (5, 10),
         ((1, 1, 1), (2, 2, 1)), description, restartcheck_params)
 
 

inputFiles/singlePhaseFlow/FieldCaseTutorial3_Isothermal_base.xml

@@ -0,0 +1,183 @@
+<?xml version="1.0" ?>
+
+<Problem>
+<!-- SPHINX_FIELD_CASE_REGION -->
+  <ElementRegions>
+    <CellElementRegion
+      name="Reservoir"
+      cellBlocks="{ 2 }"
+      materialList="{ water, rock }"/>
+
+    <CellElementRegion
+      name="Burden"
+      cellBlocks="{ 1, 3 }"
+      materialList="{ water, rock }"/>
+  </ElementRegions>
+<!-- SPHINX_FIELD_CASE_REGION_END -->
+
+<!-- SPHINX_FIELD_CASE_CONSTITUTIVE -->
+  <Constitutive>
+    <CompressibleSinglePhaseFluid
+      name="water"
+      defaultDensity="1000"
+      defaultViscosity="0.001"
+      referencePressure="0.0"
+      compressibility="1e-9"
+      viscosibility="0.0"/>
+
+    <CompressibleSolidConstantPermeability
+      name="rock"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"/>
+
+    <NullModel
+      name="nullSolid"/>
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.05"
+      referencePressure="10e7"
+      compressibility="1.0e-9"/>
+
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-13, 1.0e-13, 1.0e-16 }"/>
+  </Constitutive>
+<!-- SPHINX_FIELD_CASE_CONSTITUTIVE_END -->
+
+<!-- SPHINX_FIELD_CASE_NUMERICAL -->
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="singlePhaseTPFA"
+      />
+    </FiniteVolume>
+  </NumericalMethods>
+<!-- SPHINX_FIELD_CASE_NUMERICAL_END -->
+
+<!-- SPHINX_FIELD_CASE_GEOMETRY -->
+  <Geometry>
+    <Box
+      name="all"
+      xMin="{ -1e9, -1e9, -1e9 }"
+      xMax="{ 1e9, 1e9, 1e9 }"/>
+
+    <Box
+      name="source"
+      xMin="{ 15500, 7000, -5000 }"
+      xMax="{ 16000, 7500, 0 }"/>
+
+    <Box
+      name="sink"
+      xMin="{ 6500, 1500, -5000 }"
+      xMax="{ 7000, 2000, 0 }"/>
+  </Geometry>
+<!-- SPHINX_FIELD_CASE_GEOMETRY_END -->
+
+<!-- SPHINX_FIELD_CASE_FIELD -->
+  <FieldSpecifications>
+
+    <FieldSpecification
+       name="permx"
+       initialCondition="1"
+       component="0"
+       setNames="{ all }"
+       objectPath="ElementRegions/Reservoir"
+       fieldName="rockPerm_permeability"
+       scale="1e-15"
+       functionName="permxFunc"/>
+
+     <FieldSpecification
+       name="permy"
+       initialCondition="1"
+       component="1"
+       setNames="{ all }"
+       objectPath="ElementRegions/Reservoir"
+       fieldName="rockPerm_permeability"
+       scale="1e-15"
+       functionName="permyFunc"/>
+
+     <FieldSpecification
+       name="permz"
+       initialCondition="1"
+       component="2"
+       setNames="{ all }"
+       objectPath="ElementRegions/Reservoir"
+       fieldName="rockPerm_permeability"
+       scale="3e-15"
+       functionName="permzFunc"/>
+
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/Reservoir/2_tetrahedra"
+      fieldName="pressure"
+      scale="1e7"
+      />
+
+    <FieldSpecification
+      name="sourceTerm"
+      objectPath="ElementRegions/Reservoir/2_tetrahedra"
+      fieldName="pressure"
+      scale="15e7"
+      setNames="{ source }"
+      />
+
+    <FieldSpecification
+      name="sinkTerm"
+      objectPath="ElementRegions/Reservoir/2_tetrahedra"
+      fieldName="pressure"
+      scale="5e7"
+      setNames="{ sink }"/>
+  </FieldSpecifications>
+<!-- SPHINX_FIELD_CASE_FIELD_END -->
+
+<!-- SPHINX_FIELD_CASE_OUTPUT -->
+  <Outputs>
+    <VTK
+      name="reservoir_with_properties"/>
+    <Restart
+      name="restartOutput"/>
+  </Outputs>
+<!-- SPHINX_FIELD_CASE_OUTPUT_END -->
+
+<!-- SPHINX_FIELD_CASE_TFUNC -->
+  <Functions>
+    <TableFunction
+      name="timeInj"
+      inputVarNames="{ time }"
+      coordinates="{ 1e6, 10e6, 50e6 }"
+      values="{ 1, 0.01, 0.00001 }"/>
+
+    <TableFunction
+      name="initialPressureFunc"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ tables_FieldCaseTuto/xlin.geos, tables_FieldCaseTuto/ylin.geos, tables_FieldCaseTuto/zlin.geos }"
+      voxelFile="tables_FieldCaseTuto/pressure.geos"/>
+
+    <TableFunction
+      name="permxFunc"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ tables_FieldCaseTuto/xlin.geos, tables_FieldCaseTuto/ylin.geos, tables_FieldCaseTuto/zlin.geos }"
+      voxelFile="tables_FieldCaseTuto/permx.geos"
+      interpolation="nearest"/>
+
+    <TableFunction
+      name="permyFunc"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ tables_FieldCaseTuto/xlin.geos, tables_FieldCaseTuto/ylin.geos, tables_FieldCaseTuto/zlin.geos }"
+      voxelFile="tables_FieldCaseTuto/permy.geos"
+      interpolation="nearest"/>
+
+    <TableFunction
+      name="permzFunc"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ tables_FieldCaseTuto/xlin.geos, tables_FieldCaseTuto/ylin.geos, tables_FieldCaseTuto/zlin.geos }"
+      voxelFile="tables_FieldCaseTuto/permz.geos"
+      interpolation="nearest"/>
+  </Functions>
+<!-- SPHINX_FIELD_CASE_TFUNC_END -->
+
+</Problem>

inputFiles/singlePhaseFlow/FieldCaseTutorial3_Isothermal_smoke.xml

@@ -0,0 +1,47 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Included>
+    <File name="./FieldCaseTutorial3_Isothermal_base.xml" />
+  </Included>
+
+<!-- SPHINX_FIELD_CASE_SOLVER -->
+  <Solvers>
+    <SinglePhaseFVM
+      name="SinglePhaseFlow"
+      discretization="singlePhaseTPFA"
+      targetRegions="{ Reservoir }"
+      logLevel="1">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-6"
+        newtonMaxIter="8"
+        lineSearchAction="None"
+        logLevel="1"/>
+      <LinearSolverParameters
+        directParallel="0"/>
+    </SinglePhaseFVM>
+  </Solvers>
+<!-- SPHINX_FIELD_CASE_SOLVER_END -->
+
+  <!-- SPHINX_FIELD_CASE_MESH -->
+  <Mesh>
+    <VTKMesh name="SyntheticMesh"
+                file="synthetic.vtu" />
+  </Mesh>
+  <!-- SPHINX_FIELD_CASE_MESH_END -->
+
+  <!-- SPHINX_FIELD_CASE_EVENTS -->
+  <Events maxTime="100.0e6">
+    <PeriodicEvent name="solverApplications"
+                   forceDt="50.0e6"
+                   target="/Solvers/SinglePhaseFlow" />
+
+
+    <PeriodicEvent name="restarts"
+                   timeFrequency="50.0e6"
+                   target="/Outputs/restartOutput"/>
+
+  </Events>
+  <!-- SPHINX_FIELD_CASE_EVENTS_END -->
+
+</Problem>