tutorials/multiphase/reactingTwoPhaseEulerFoam: Initialize by first running without boiling

then restart with boiling.

Contributed by Juho Peltola, VTT

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/0/T.gas

@@ -17,21 +17,21 @@ FoamFile
 
 dimensions      [0 0 0 1 0 0 0];
 
-internalField   uniform 373.55;
+internalField   uniform 368;
 
 boundaryField
 {
     inlet
     {
         type            fixedValue;
-        value           uniform 373.55;
+        value           uniform 368;
     }
     outlet
     {
         type            inletOutlet;
         phi             phi.gas;
-        inletValue      uniform 373.55;
+        inletValue      uniform 368;
-        value           uniform 373.55;
+        value           uniform 368;
     }
     wall1
     {
@@ -41,7 +41,7 @@ boundaryField
     {
         type            copiedFixedValue;
         sourceFieldName T.liquid;
-        value           uniform 373.55;
+        value           uniform 368;
     }
     defaultFaces
     {

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/0/T.liquid

@@ -17,21 +17,21 @@ FoamFile
 
 dimensions      [0 0 0 1 0 0 0];
 
-internalField   uniform 370;
+internalField   uniform 368;
 
 boundaryField
 {
     inlet
     {
         type            fixedValue;
-        value           uniform 370;
+        value           uniform 368;
     }
     outlet
     {
         type            inletOutlet;
         phi             phi.liquid;
-        inletValue      uniform 370;
+        inletValue      uniform 368;
-        value           uniform 370;
+        value           uniform 368;
     }
     wall1
     {
@@ -41,9 +41,9 @@ boundaryField
     {
         type            fixedMultiPhaseHeatFlux;
         relax           0.5;
-        q               uniform 100000;
+        q               uniform 0;
         phase           "liquid";
-        value           uniform 370;
+        value           uniform 368;
     }
     defaultFaces
     {

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/0/alpha.gas

@@ -10,7 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       volScalarField;
-    location    "5";
+    location    "0";
     object      alpha.gas;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/Allclean

@@ -0,0 +1,10 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+cp ./system/controlDict.org ./system/controlDict
+
+#------------------------------------------------------------------------------

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/Allrun

@@ -0,0 +1,15 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application=`getApplication`
+
+runApplication blockMesh
+runApplication $application
+foamDictionary system/controlDict -entry endTime -set 5
+foamDictionary 2/T.liquid -entry boundaryField.wall2.q -set 1e5
+runApplication -a $application
+
+#------------------------------------------------------------------------------

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/system/controlDict

@@ -17,7 +17,7 @@ FoamFile
 
 application     reactingTwoPhaseEulerFoam;
 
-startFrom       startTime;
+startFrom       latestTime;
 
 startTime       0;
 
@@ -25,11 +25,11 @@ stopAt          endTime;
 
 endTime         2;
 
-deltaT          1e-5;
+deltaT          1e-4;
 
 writeControl    adjustableRunTime;
 
-writeInterval   0.1;
+writeInterval   0.5;
 
 purgeWrite      0;
 
@@ -37,7 +37,7 @@ writeFormat     ascii;
 
 writePrecision  9;
 
-writeCompression compressed;
+writeCompression uncompressed;
 
 timeFormat      general;
 

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/system/controlDict.org

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     reactingTwoPhaseEulerFoam;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         2;
+
+deltaT          1e-4;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.5;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  9;
+
+writeCompression uncompressed;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  yes;
+
+maxCo           0.05;
+
+maxDeltaT       0.001;
+
+// ************************************************************************* //

tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling/system/fvSolution

@@ -80,10 +80,10 @@ solvers
 
 PIMPLE
 {
-    nOuterCorrectors    3;
+    nOuterCorrectors    6;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;
-    nEnergyCorrectors   3;
+    nEnergyCorrectors   2;
     faceMomentum        yes;
 }
 
@@ -97,7 +97,7 @@ relaxationFactors
     equations
     {
         ".*"            1;
-        "h.*"           0.4;
+        "h.*"           0.3;
     }
 }