diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/MRFZones b/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/MRFZones
new file mode 100644
index 0000000000..0034a58c22
--- /dev/null
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/MRFZones
@@ -0,0 +1,20 @@
+/*--------------------------------*- 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    "constant";
+    object      MRFZones;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+0()
+
+// ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/transportProperties b/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/transportProperties
index ef9a0fccf7..6d05c03c61 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/transportProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed/constant/transportProperties
@@ -17,25 +17,22 @@ FoamFile
 
 phase1
 {
-    rho             rho [ 1 -3 0 0 0 ] 2640;
-    nu              nu [ 0 2 -1 0 0 ] 1e-06;
-    d               d [ 0 1 0 0 0 0 0 ] 0.00048;
+    rho             2640;
+    nu              1e-06;
+    d               0.00048;
 }
 
 phase2
 {
-    rho             rho [ 1 -3 0 0 0 ] 1.28;
-    nu              nu [ 0 2 -1 0 0 ] 1.328e-05;
-    d               d [ 0 1 0 0 0 0 0 ] 1;
+    rho             1.28;
+    nu              1.328e-05;
+    d               1;
 }
 
-Cvm             Cvm [ 0 0 0 0 0 ] 0;
+Cvm             0;
 
-Cl              Cl [ 0 0 0 0 0 ] 0;
-
-Ct              Ct [ 0 0 0 0 0 ] 0;
-
-alpha1Alpha      alpha1Alpha [ 0 0 0 0 0 ] 0;
+Cl              0;
 
+Ct              0;
 
 // ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed/system/fvSolution b/tutorials/multiphase/twoPhaseEulerFoam/bed/system/fvSolution
index bdb02666e3..5800e2a4b4 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bed/system/fvSolution
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed/system/fvSolution
@@ -49,25 +49,21 @@ solvers
 
     "(k|epsilon)Final"
     {
-        solver          PBiCG;
-        preconditioner  DILU;
-        tolerance       1e-05;
+        $k;
         relTol          0;
     }
 
     alpha
     {
-        solver          PBiCG;
-        preconditioner  DILU;
+        solver          PCG;
+        preconditioner  DIC;
         tolerance       1e-10;
         relTol          0.1;
     }
 
     alpha1Final
     {
-        solver          PBiCG;
-        preconditioner  DILU;
-        tolerance       1e-10;
+        $alpha;
         relTol          0;
     }
 }
@@ -77,6 +73,7 @@ PIMPLE
     nCorrectors     2;
     nNonOrthogonalCorrectors 0;
     nAlphaCorr      2;
+    nAlphaSubCycles 3;
     correctAlpha    yes;
     pRefCell        0;
     pRefValue       0;
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/MRFZones b/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/MRFZones
new file mode 100644
index 0000000000..0034a58c22
--- /dev/null
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/MRFZones
@@ -0,0 +1,20 @@
+/*--------------------------------*- 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    "constant";
+    object      MRFZones;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+0()
+
+// ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/transportProperties b/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/transportProperties
index 859bd94263..1afd6e92c5 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/transportProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed2/constant/transportProperties
@@ -17,25 +17,20 @@ FoamFile
 
 phase1
 {
-    rho             rho [ 1 -3 0 0 0 ] 2500;
-    nu              nu [ 0 2 -1 0 0 ] 1e-06;
-    d               d [ 0 1 0 0 0 0 0 ] 0.0003;
+    rho             2500;
+    nu              1e-06;
+    d               0.0003;
 }
 
 phase2
 {
-    rho             rho [ 1 -3 0 0 0 ] 1.2;
-    nu              nu [ 0 2 -1 0 0 ] 1.5e-05;
-    d               d [ 0 1 0 0 0 0 0 ] 1;
+    rho             1.2;
+    nu              1.5e-05;
+    d               1;
 }
 
-Cvm             Cvm [ 0 0 0 0 0 ] 0;
-
-Cl              Cl [ 0 0 0 0 0 ] 0;
-
-Ct              Ct [ 0 0 0 0 0 ] 0;
-
-alpha1Alpha      alpha1Alpha [ 0 0 0 0 0 ] 0;
-
+Cvm             0;
+Cl              0;
+Ct              0;
 
 // ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bed2/system/fvSolution b/tutorials/multiphase/twoPhaseEulerFoam/bed2/system/fvSolution
index bdb02666e3..2e3cdfa4aa 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bed2/system/fvSolution
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bed2/system/fvSolution
@@ -76,7 +76,8 @@ PIMPLE
 {
     nCorrectors     2;
     nNonOrthogonalCorrectors 0;
-    nAlphaCorr      2;
+    nAlphaCorr      1;
+    nAlphaSubCycles 3;
     correctAlpha    yes;
     pRefCell        0;
     pRefValue       0;
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/MRFZones b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/MRFZones
new file mode 100644
index 0000000000..0034a58c22
--- /dev/null
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/MRFZones
@@ -0,0 +1,20 @@
+/*--------------------------------*- 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    "constant";
+    object      MRFZones;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+0()
+
+// ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/transportProperties b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/transportProperties
index 276c521ad3..e013730cb5 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/transportProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/constant/transportProperties
@@ -17,23 +17,21 @@ FoamFile
 
 phase1
 {
-    nu              nu [ 0 2 -1 0 0 0 0 ] 1.6e-05;
-    rho             rho [ 1 -3 0 0 0 0 0 ] 1;
-    d               d [ 0 1 0 0 0 0 0 ] 0.003;
+    nu              1.6e-05;
+    rho             1;
+    d               0.003;
 }
 
 phase2
 {
-    nu              nu [ 0 2 -1 0 0 0 0 ] 1e-06;
-    rho             rho [ 1 -3 0 0 0 0 0 ] 1000;
-    d               d [ 0 1 0 0 0 0 0 ] 0.0001;
+    nu              1e-06;
+    rho             1000;
+    d               0.0001;
 }
 
-Cvm             Cvm [ 0 0 0 0 0 0 0 ] 0.5;
-
-Cl              Cl [ 0 0 0 0 0 0 0 ] 0;
-
-Ct              Ct [ 0 0 0 0 0 0 0 ] 1;
+Cvm             0.5;
+Cl              0;
+Ct              1;
 
 
 // ************************************************************************* //
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/system/fvSolution b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/system/fvSolution
index d0491f6d65..76a95c8327 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/system/fvSolution
+++ b/tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/system/fvSolution
@@ -76,7 +76,8 @@ PIMPLE
 {
     nCorrectors     2;
     nNonOrthogonalCorrectors 0;
-    nAlphaCorr      2;
+    nAlphaCorr      1;
+    nAlphaSubCycles 2;
     correctAlpha    yes;
     pRefCell        0;
     pRefValue       0;
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/controlDict b/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/controlDict
index dd1a00f171..7b2bf5a816 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/controlDict
+++ b/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/controlDict
@@ -47,7 +47,7 @@ runTimeModifiable yes;
 
 adjustTimeStep  yes;
 
-maxCo           0.2;
+maxCo           0.5;
 
 maxDeltaT       1;
 
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/fvSolution b/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/fvSolution
index 891746fcf3..69208cbe95 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/fvSolution
+++ b/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/fvSolution
@@ -76,6 +76,7 @@ PIMPLE
     nCorrectors     3;
     nNonOrthogonalCorrectors 0;
     nAlphaCorr      1;
+    nAlphaSubCycles 3;
     correctAlpha    yes;
     pRefCell        0;
     pRefValue       0;