From 4af28562c5d7177a9c523eab8184cfacf36367dd Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 22 Jun 2021 22:04:05 +0100
Subject: [PATCH] 
 tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D: Changed to
 enthalpy

Solving for enthalpy provides better convergence and stability than internal
energy.  Also correctPhi is now off pending the addition of compressibility
effects to the pcorr equation.
---
 .../constant/thermophysicalProperties.air             |  2 +-
 .../constant/thermophysicalProperties.mercury         |  2 +-
 .../constant/thermophysicalProperties.oil             |  2 +-
 .../constant/thermophysicalProperties.water           |  2 +-
 .../laminar/mixerVesselAMI2D/system/fvSolution        | 11 ++++-------
 5 files changed, 8 insertions(+), 11 deletions(-)

diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.air b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.air
index aca85e96f4..d04ab6ddce 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.air
+++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.air
@@ -23,7 +23,7 @@ thermoType
     thermo          hConst;
     equationOfState perfectGas;
     specie          specie;
-    energy          sensibleInternalEnergy;
+    energy          sensibleEnthalpy;
 }
 
 mixture
diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.mercury b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.mercury
index 18091e495f..98932db82b 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.mercury
+++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.mercury
@@ -23,7 +23,7 @@ thermoType
     thermo          hConst;
     equationOfState rhoConst;
     specie          specie;
-    energy          sensibleInternalEnergy;
+    energy          sensibleEnthalpy;
 }
 
 mixture
diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.oil b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.oil
index 2a84f2b931..97cb162f8d 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.oil
+++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.oil
@@ -23,7 +23,7 @@ thermoType
     thermo          hConst;
     equationOfState rhoConst;
     specie          specie;
-    energy          sensibleInternalEnergy;
+    energy          sensibleEnthalpy;
 }
 
 mixture
diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.water b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.water
index f745ebe9c5..48318e4da5 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.water
+++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/constant/thermophysicalProperties.water
@@ -23,7 +23,7 @@ thermoType
     thermo          eConst;
     equationOfState rPolynomial;
     specie          specie;
-    energy          sensibleInternalEnergy;
+    energy          sensibleEnthalpy;
 }
 
 mixture
diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/system/fvSolution b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/system/fvSolution
index 3141a6212d..6030d58bef 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/system/fvSolution
+++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVesselAMI2D/system/fvSolution
@@ -63,7 +63,7 @@ solvers
         minIter         1;
     }
 
-    "e.*"
+    "(h|e).*"
     {
         solver          smoothSolver;
         smoother        symGaussSeidel;
@@ -76,14 +76,11 @@ solvers
 
 PIMPLE
 {
-    nOuterCorrectors     1;
-    nCorrectors     3;
+    nOuterCorrectors    1;
+    nCorrectors         3;
     nNonOrthogonalCorrectors 0;
 
-    correctPhi      yes;
-
-    pRefCell        0;
-    pRefValue       0;
+    correctPhi          no;
 }
 
 relaxationFactors