ENH: Update icoReactingMultiPhaseInterFoam tutorials

tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution

@@ -39,6 +39,14 @@ solvers
         cYi             1;
     }
 
+    "pcorr.*"
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-5;
+        relTol          0;
+    }
+
     "(U.*|epsilon.*|k.*)"
     {
         solver          smoothSolver;

tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/controlDict

@@ -29,7 +29,7 @@ deltaT          1e-5;
 
 writeControl    adjustable;
 
-writeInterval   0.02;
+writeInterval   0.01;
 
 purgeWrite      0;
 

tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSchemes

@@ -29,7 +29,7 @@ divSchemes
 {
     default         none;
 
-    div(rhoPhi,U)           Gauss linearUpwind gradU;
+    div(rhoPhi,U)           Gauss limitedLinearV 1;
 
     "div\(phi,alpha.*\)"        Gauss vanLeer;
     "div\(phir,alpha.*\)"       Gauss linear;

tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSolution

@@ -39,7 +39,7 @@ solvers
                         );
 
         nAlphaCorr      2;
-        nAlphaSubCycles 2;
+        nAlphaSubCycles 3;
 
         // Compressiion factor for species in each alpha phase
         // NOTE: It should be similar to cAlpha
@@ -74,19 +74,16 @@ solvers
 
     p_rgh
     {
-        solver          PCG;
-        preconditioner  DIC;
-        tolerance       1e-9;
-        relTol          0.01;
+        solver              GAMG;
+        smoother            DICGaussSeidel;
+        tolerance           1e-11;
+        relTol              0.01;
 
-        nPreSweeps      0;
-        nPostSweeps     2;
-
-        cacheAgglomeration true;
+        cacheAgglomeration  true;
 
         nCellsInCoarsestLevel 10;
-        agglomerator    faceAreaPair;
-        mergeLevels     1;
+        agglomerator        faceAreaPair;
+        mergeLevels         1;
     };
 
     p_rghFinal
@@ -94,10 +91,19 @@ solvers
         $p_rgh;
         relTol          0;
     };
+
+    "pcorr.*"
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-5;
+        relTol          0;
+    }
 }
 
 PIMPLE
 {
+    correctPhi          false;
     momentumPredictor   false;
     nOuterCorrectors    1;
     nCorrectors         4;

tutorials/multiphase/icoReactingMultiPhaseInterFoam/poolEvaporation/system/fvSolution

@@ -49,6 +49,11 @@ solvers
         relTol          0;
     }
 
+    "pcorr.*"
+    {
+        $p_rgh;
+    }
+
     "U.*"
     {
         solver          smoothSolver;

tutorials/multiphase/icoReactingMultiPhaseInterFoam/solidMelting2D/system/fvSolution

@@ -36,7 +36,7 @@ solvers
         relTol          0;
     }
 
-    pcorr
+    "pcorr.*"
     {
         solver          PCG;
         preconditioner  DIC;