From 332b72d5619675a0ce24abc8605affb89854ebd9 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Wed, 31 Jul 2019 15:57:18 +0100
Subject: [PATCH] tutorials: reactingTwoPhaseEulerFoam: Updates to
 steamInjection

Added limiters for the phase temperatures to prevent divergence, and
monitors to report the minimum and maximum values. Removed the
setTimeStep functionObject as the temperature limiters make this
unnecessary. Dereased the number of energy correctors and set a higher
Courant number limit to reduce the execution-time of the case.

Patch contributed by Juho Peltola, VTT.
---
 .../laminar/steamInjection/constant/fvOptions | 20 +++++++++++++++++++
 .../laminar/steamInjection/system/controlDict | 17 +++++-----------
 .../steamInjection/system/deltaTvalues        |  9 ---------
 .../laminar/steamInjection/system/fvSolution  |  2 +-
 4 files changed, 26 insertions(+), 22 deletions(-)
 delete mode 100644 tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/deltaTvalues

diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/fvOptions b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/fvOptions
index 8628ed295f..e32427e29f 100644
--- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/fvOptions
+++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/fvOptions
@@ -66,6 +66,26 @@ options
             h.steam      (3700 0); // kg*m^2/s^3
         }
     }
+
+    limitTsteam
+    {
+        type            limitTemperature;
+        active          yes;
+        selectionMode   all;
+        min             270;
+        max             2000;
+        phase           steam;
+    }
+
+    limitTwater
+    {
+        type            limitTemperature;
+        active          yes;
+        selectionMode   all;
+        min             270;
+        max             2000;
+        phase           water;
+    }
 }
 
 
diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/controlDict b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/controlDict
index cf3faac2f5..000d6ed513 100644
--- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/controlDict
+++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/controlDict
@@ -47,28 +47,21 @@ runTimeModifiable yes;
 
 adjustTimeStep  yes;
 
-maxCo           0.1;
+maxCo           0.25;
 
 maxDeltaT       1e-2;
 
 functions
 {
-    timeStepping
-    {
-        type            setTimeStep;
-        functionObjectLibs ("libutilityFunctionObjects.so");
-        enabled         yes;
-        deltaT          tableFile;
-        file            "system/deltaTvalues";
-    }
-
-    minMaxp
+    minMax
     {
         type        fieldMinMax;
         functionObjectLibs ("libfieldFunctionObjects.so");
         fields
         (
-             p
+            T.steam
+            T.water
+            p
         );
         location        no;
         writeControl    timeStep;
diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/deltaTvalues b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/deltaTvalues
deleted file mode 100644
index 03c2defea0..0000000000
--- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/deltaTvalues
+++ /dev/null
@@ -1,9 +0,0 @@
-(
-    (0 1e-3)
-    (0.99 1e-3)
-    (0.999 1e-4)
-    (0.9999 1e-5)
-    (1.001 1e-5)
-    (1.01 1e-4)
-    (1.1 1e-3)
-);
diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/fvSolution b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/fvSolution
index 5e353fee62..4e76f88020 100644
--- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/fvSolution
+++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/system/fvSolution
@@ -81,7 +81,7 @@ PIMPLE
     nOuterCorrectors 3;
     nCorrectors      1;
     nNonOrthogonalCorrectors 0;
-    nEnergyCorrectors 2;
+    nEnergyCorrectors 1;
     faceMomentum    yes;
 }