diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/TEqns.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/TEqns.H
index 0af1e17946..861590b235 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/TEqns.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/TEqns.H
@@ -11,7 +11,7 @@
      ==
         heatTransferCoeff*T2/Cp1/rho1
       - fvm::Sp(heatTransferCoeff/Cp1/rho1, T1)
-      + alpha1*(dpdt - (fvc::ddt(rho1, K1) + fvc::div(phi1, K1)))/Cp1/rho1
+      + alpha1*(dpdt/rho1 - (fvc::ddt(K1) + fvc::div(phi1, K1)))/Cp1
     );
 
     fvScalarMatrix T2Eqn
@@ -23,7 +23,7 @@
      ==
         heatTransferCoeff*T1/Cp2/rho2
       - fvm::Sp(heatTransferCoeff/Cp2/rho2, T2)
-      + alpha2*(dpdt - (fvc::ddt(rho2, K2) + fvc::div(phi2, K2)))/Cp2/rho2
+      + alpha2*(dpdt/rho2 - (fvc::ddt(K2) + fvc::div(phi2, K2)))/Cp2
     );
 
     T1Eqn.relax();
diff --git a/tutorials/compressible/rhoPorousMRFSimpleFoam/angledDuctImplicit/constant/polyMesh/boundary b/tutorials/compressible/rhoPorousMRFSimpleFoam/angledDuctImplicit/constant/polyMesh/boundary
index a664893343..0abd1608ab 100644
--- a/tutorials/compressible/rhoPorousMRFSimpleFoam/angledDuctImplicit/constant/polyMesh/boundary
+++ b/tutorials/compressible/rhoPorousMRFSimpleFoam/angledDuctImplicit/constant/polyMesh/boundary
@@ -8,7 +8,7 @@
 FoamFile
 {
     version     2.0;
-    format      binary;
+    format      ascii;
     class       polyBoundaryMesh;
     location    "constant/polyMesh";
     object      boundary;
diff --git a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/bubbleColumn/system/fvSchemes b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/bubbleColumn/system/fvSchemes
index 93b115565d..34b14edfcc 100644
--- a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/bubbleColumn/system/fvSchemes
+++ b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/bubbleColumn/system/fvSchemes
@@ -44,8 +44,8 @@ divSchemes
     div(phi,Theta)      Gauss limitedLinear 1;
     div(phid1,p)        Gauss upwind;
     div(phid2,p)        Gauss upwind;
-    div(phi1,K1)        Gauss limitedLinearV 1;
-    div(phi2,K2)        Gauss limitedLinearV 1;
+    div(phi1,K1)        Gauss limitedLinear 1;
+    div(phi2,K2)        Gauss limitedLinear 1;
 }
 
 laplacianSchemes
diff --git a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/fluidisedBed/system/fvSchemes b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/fluidisedBed/system/fvSchemes
index 27ba67b86d..4412903f93 100644
--- a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/fluidisedBed/system/fvSchemes
+++ b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/fluidisedBed/system/fvSchemes
@@ -44,8 +44,8 @@ divSchemes
     div((alpha2*Rc2))        Gauss linear;
     div(phid1,p)        Gauss upwind;
     div(phid2,p)        Gauss upwind;
-    div(phi1,K1)        Gauss limitedLinearV 1;
-    div(phi2,K2)        Gauss limitedLinearV 1;
+    div(phi1,K1)        Gauss limitedLinear 1;
+    div(phi2,K2)        Gauss limitedLinear 1;
 }
 
 laplacianSchemes