Feature thermal single-phase mixed mimetic

.integrated_tests.yaml

@@ -1,6 +1,6 @@
 baselines:
   bucket: geosx
-  baseline: integratedTests/baseline_integratedTests-pr4169-17867-5d0c7de
+  baseline: integratedTests/baseline_integratedTests-pr4139-17884-0f6d03e
 allow_fail:
   all: ''
   streak: ''

BASELINE_NOTES.md

@@ -5,6 +5,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 #4139 (2026-10-31) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4139-17884-0f6d03e.tar.gz>
+=====================
+Add the mixed mimetic discretization: new `NumericalMethods/MixedMimetic` group in every restart file and new mixed MFD tests.
+
 PR #4169 (2026-10-30) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4169-17867-5d0c7de.tar.gz>
 =====================
 Fix surface generator splitting for contact mechanics case with wells

inputFiles/singlePhaseFlow/polyhedral/lapwoodConvectionFVM_Ra30.xml

@@ -0,0 +1,64 @@
+<?xml version="1.0" ?>
+
+<!-- Ra = 30, below the threshold 4.5 pi^2 of the seeded mode: the disturbance decays, n = -8.01e-10 1/s -->
+
+<Problem>
+  <Included>
+    <File name="./lapwoodConvectionFVM_base.xml"/>
+  </Included>
+
+  <!-- k = Ra / 7.848e12 -->
+  <FieldSpecifications>
+    <FieldSpecification
+      name="permX"
+      initialCondition="1"
+      component="0"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="3.82263e-12"/>
+
+    <FieldSpecification
+      name="permY"
+      initialCondition="1"
+      component="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="3.82263e-12"/>
+
+    <FieldSpecification
+      name="permZ"
+      initialCondition="1"
+      component="2"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="3.82263e-12"/>
+  </FieldSpecifications>
+
+  <!-- 2.0e10 s is the diffusion time h^2 / kappa -->
+  <Events
+    maxTime="2.0e10">
+    <PeriodicEvent
+      name="solverApplications"
+      forceDt="1.0e8"
+      target="/Solvers/flow"/>
+
+    <PeriodicEvent
+      name="probeCollection"
+      timeFrequency="1.0e8"
+      target="/Tasks/probeTemperature"/>
+
+    <PeriodicEvent
+      name="probeOutput"
+      timeFrequency="1.0e9"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput"/>
+
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="1.0e9"
+      target="/Outputs/vtkOutput"/>
+  </Events>
+</Problem>

inputFiles/singlePhaseFlow/polyhedral/lapwoodConvectionFVM_Ra50.xml

@@ -0,0 +1,64 @@
+<?xml version="1.0" ?>
+
+<!-- Ra = 50, above the threshold 4.5 pi^2 of the seeded mode: the disturbance grows, n = 3.10e-10 1/s -->
+
+<Problem>
+  <Included>
+    <File name="./lapwoodConvectionFVM_base.xml"/>
+  </Included>
+
+  <!-- k = Ra / 7.848e12 -->
+  <FieldSpecifications>
+    <FieldSpecification
+      name="permX"
+      initialCondition="1"
+      component="0"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="6.37105e-12"/>
+
+    <FieldSpecification
+      name="permY"
+      initialCondition="1"
+      component="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="6.37105e-12"/>
+
+    <FieldSpecification
+      name="permZ"
+      initialCondition="1"
+      component="2"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="rockPerm_permeability"
+      scale="6.37105e-12"/>
+  </FieldSpecifications>
+
+  <!-- 2.0e10 s is the diffusion time h^2 / kappa -->
+  <Events
+    maxTime="2.0e10">
+    <PeriodicEvent
+      name="solverApplications"
+      forceDt="1.0e8"
+      target="/Solvers/flow"/>
+
+    <PeriodicEvent
+      name="probeCollection"
+      timeFrequency="1.0e8"
+      target="/Tasks/probeTemperature"/>
+
+    <PeriodicEvent
+      name="probeOutput"
+      timeFrequency="1.0e9"
+      targetExactTimestep="0"
+      target="/Outputs/timeHistoryOutput"/>
+
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="1.0e9"
+      target="/Outputs/vtkOutput"/>
+  </Events>
+</Problem>

inputFiles/singlePhaseFlow/polyhedral/lapwoodConvectionFVM_base.xml

@@ -0,0 +1,208 @@
+<?xml version="1.0" ?>
+
+<!-- The problem of lapwoodConvection_base.xml with the cell-centred finite volume solver and the two-point flux
+     approximation: same domain, medium, fluid, initial state and events.
+
+     Difference in the thermal boundary condition. In SinglePhaseFVM a face Dirichlet condition prescribes the pressure
+     and the temperature together, so an impervious face of given temperature is not available. The conductive
+     temperature 320 - 0.2 z is imposed instead in the cells of the bottom and top layers, at their centres
+     z = dz_b / 2 and h - dz_b / 2: a condition of first order in the thickness dz_b of those layers.
+
+     The mesh is therefore graded towards z = 0 and z = h: dz = 10/3 m in [5, 95] as in the mixed deck, 5/6 m in the
+     next 10/3 m and dz_b = 5/24 m in the last 5/3 m. With dz_b = h / 480, n = 3.10e-10 1/s at Ra = 50 and
+     -7.86e-10 1/s at Ra = 30 (theory 3.10e-10 and -8.01e-10). With uniform layers dz_b = h / 30 the rates were
+     2.52e-10 and -8.39e-10. -->
+
+<Problem>
+
+  <Solvers
+    gravityVector="{ 0.0, 0.0, -9.81 }">
+    <SinglePhaseFVM
+      name="flow"
+      logLevel="1"
+      isThermal="1"
+      discretization="tpfa"
+      temperature="310.0"
+      targetRegions="{ domain }">
+      <NonlinearSolverParameters
+        newtonTol="1.0e-5"
+        newtonMaxIter="10"/>
+      <LinearSolverParameters
+        solverType="gmres"
+        preconditionerType="amg"
+        krylovTol="1.0e-7"
+        logLevel="1"/>
+    </SinglePhaseFVM>
+  </Solvers>
+
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="tpfa"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh"
+      elementTypes="{ C3D4 }"
+      xCoords="{ 0, 200 }"
+      yCoords="{ 0, 200 }"
+      zCoords="{ 0, 1.6666666667, 5, 95, 98.3333333333, 100 }"
+      nx="{ 60 }"
+      ny="{ 30 }"
+      nz="{ 8, 4, 27, 4, 8 }"
+      cellBlockNames="{ block }"/>
+  </Mesh>
+
+  <Geometry>
+    <!-- cells of the bottom and top layers, dz = 5/24 m -->
+    <Box
+      name="bottomCells"
+      xMin="{ -0.01, -0.01, -0.01 }"
+      xMax="{ 200.01, 200.01, 0.21 }"/>
+
+    <Box
+      name="topCells"
+      xMin="{ -0.01, -0.01, 99.79 }"
+      xMax="{ 200.01, 200.01, 100.01 }"/>
+
+    <!-- cell at mid-height ( z = 50 ), at the maximum of the initial perturbation -->
+    <Box
+      name="probe"
+      xMin="{ -0.01, -0.01, 48.32 }"
+      xMax="{ 3.34, 6.68, 51.68 }"/>
+  </Geometry>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="domain"
+      cellBlocks="{ * }"
+      materialList="{ rock, fluid, thermalCond }"/>
+  </ElementRegions>
+
+  <Constitutive>
+    <CompressibleSolidConstantPermeability
+      name="rock"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"
+      solidInternalEnergyModelName="rockInternalEnergy"/>
+
+    <NullModel
+      name="nullSolid"/>
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.2"
+      referencePressure="1.0e7"
+      compressibility="0.0"/>
+
+    <SolidInternalEnergy
+      name="rockInternalEnergy"
+      referenceVolumetricHeatCapacity="2.0e6"
+      referenceTemperature="310.0"
+      referenceInternalEnergy="0.0"/>
+
+    <!-- k at Ra = 4 pi^2; the variants overwrite it -->
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 5.03038e-12, 5.03038e-12, 5.03038e-12 }"/>
+
+    <ThermalCompressibleSinglePhaseFluid
+      name="fluid"
+      defaultDensity="1000"
+      defaultViscosity="0.001"
+      referencePressure="1.0e7"
+      referenceTemperature="310.0"
+      compressibility="4.5e-10"
+      thermalExpansionCoeff="2.0e-4"
+      viscosibility="0.0"
+      temperatureViscosityCoefficient="0.0"
+      specificHeatCapacity="4000"
+      referenceInternalEnergy="1.0e6"/>
+
+    <SinglePhaseThermalConductivity
+      name="thermalCond"
+      defaultThermalConductivityComponents="{ 2.0, 2.0, 2.0 }"
+      thermalConductivityGradientComponents="{ 0, 0, 0 }"
+      referenceTemperature="310.0"/>
+  </Constitutive>
+
+  <!-- tabulated over the whole domain, independent of the mesh -->
+  <Functions>
+    <!-- 1.0e7 + rho g ( h - z ) -->
+    <TableFunction
+      name="hydrostaticPressure"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ lapwoodConvection_tables/xMid.geos, lapwoodConvection_tables/y.geos, lapwoodConvection_tables/z.geos }"
+      voxelFile="lapwoodConvection_tables/pressure.geos"/>
+
+    <!-- 320 - 0.2 z + 0.01 sin( pi z / h ) cos( pi x / h ) cos( pi y / h ), tabulated every 5 m -->
+    <TableFunction
+      name="disturbedTemperature"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ lapwoodConvection_tables/x3D.geos, lapwoodConvection_tables/y3D.geos, lapwoodConvection_tables/z3D.geos }"
+      voxelFile="lapwoodConvection_tables/temperature3D.geos"/>
+
+    <!-- 320 - 0.2 z -->
+    <TableFunction
+      name="conductiveTemperature"
+      inputVarNames="{ elementCenter }"
+      coordinateFiles="{ lapwoodConvection_tables/xMid.geos, lapwoodConvection_tables/y.geos, lapwoodConvection_tables/z.geos }"
+      voxelFile="lapwoodConvection_tables/conduction.geos"/>
+  </Functions>
+
+  <FieldSpecifications>
+    <FieldSpecification
+      name="initialPressure"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="pressure"
+      functionName="hydrostaticPressure"
+      scale="1.0"/>
+
+    <FieldSpecification
+      name="initialTemperature"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="temperature"
+      functionName="disturbedTemperature"
+      scale="1.0"/>
+
+    <!-- conductive temperature in the bottom and top cell layers; no flow condition: every boundary is impervious -->
+    <FieldSpecification
+      name="boundaryTemperature"
+      setNames="{ bottomCells, topCells }"
+      objectPath="ElementRegions"
+      fieldName="temperature"
+      functionName="conductiveTemperature"
+      scale="1.0"/>
+  </FieldSpecifications>
+
+  <!-- the disturbance is temperature - ( 320 - 0.2 z ) in the probe cell, z stored with the history -->
+  <Tasks>
+    <PackCollection
+      name="probeTemperature"
+      objectPath="ElementRegions/domain/block"
+      fieldName="temperature"
+      setNames="{ probe }"/>
+  </Tasks>
+
+  <Outputs>
+    <VTK
+      name="vtkOutput"
+      fieldNames="{ pressure, temperature }"/>
+
+    <TimeHistory
+      name="timeHistoryOutput"
+      sources="{ /Tasks/probeTemperature }"
+      filename="lapwoodProbe"/>
+
+    <Restart
+      name="restartOutput"/>
+  </Outputs>
+
+</Problem>