diff --git a/tutorials/cfdemSolverMultiphase/damBreak/CFD/system/fvSchemes b/tutorials/cfdemSolverMultiphase/damBreak/CFD/system/fvSchemes
index e15a3682..4d2ff463 100644
--- a/tutorials/cfdemSolverMultiphase/damBreak/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverMultiphase/damBreak/CFD/system/fvSchemes
@@ -38,6 +38,7 @@ divSchemes
     div(((dev(grad(U).T())*rho)*dev(grad(U).T()))) Gauss linear;
     div((((((alpha.water*rho)*nu.water)|(alpha.water*rho))*rho)*dev(grad(U).T()))) Gauss linear;
     div(((nu*rho)*dev(grad(U).T()))) Gauss linear;
+	default	Gauss linear;
 }
 
 laplacianSchemes
diff --git a/tutorials/cfdemSolverMultiphase/particlesWaterEntry/CFD/system/fvSchemes b/tutorials/cfdemSolverMultiphase/particlesWaterEntry/CFD/system/fvSchemes
index 2f12a845..764a3064 100644
--- a/tutorials/cfdemSolverMultiphase/particlesWaterEntry/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverMultiphase/particlesWaterEntry/CFD/system/fvSchemes
@@ -39,6 +39,7 @@ divSchemes
     div(((dev(grad(U).T())*rho)*dev(grad(U).T()))) Gauss linear;
     div((((((alpha.water*rho)*nu.water)|(alpha.water*rho))*rho)*dev(grad(U).T()))) Gauss linear;
     div(((nu*rho)*dev(grad(U).T()))) Gauss linear;
+	default	Gauss linear;
 }
 
 laplacianSchemes
diff --git a/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/constant/couplingProperties b/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/constant/couplingProperties
index fd787c1e..b1764bf6 100644
--- a/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/constant/couplingProperties
+++ b/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/constant/couplingProperties
@@ -304,8 +304,8 @@ gradPForceSmoothProps
     {
         lowerLimit 0.1;
         upperLimit 1e10;
-		referenceField "p_rgh";
-		alpha 0.01;
+		refField "p_rgh";
+		gamma 0.01;
 	}
 
 }
diff --git a/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/system/fvSchemes b/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/system/fvSchemes
index e15a3682..4d2ff463 100644
--- a/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverMultiphase/small_scale_hearth_drainage/run/CFD/system/fvSchemes
@@ -38,6 +38,7 @@ divSchemes
     div(((dev(grad(U).T())*rho)*dev(grad(U).T()))) Gauss linear;
     div((((((alpha.water*rho)*nu.water)|(alpha.water*rho))*rho)*dev(grad(U).T()))) Gauss linear;
     div(((nu*rho)*dev(grad(U).T()))) Gauss linear;
+	default	Gauss linear;
 }
 
 laplacianSchemes
diff --git a/tutorials/cfdemSolverMultiphase/tankDrainage/CFD/system/fvSchemes b/tutorials/cfdemSolverMultiphase/tankDrainage/CFD/system/fvSchemes
index e15a3682..d0a271ce 100644
--- a/tutorials/cfdemSolverMultiphase/tankDrainage/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverMultiphase/tankDrainage/CFD/system/fvSchemes
@@ -38,6 +38,7 @@ divSchemes
     div(((dev(grad(U).T())*rho)*dev(grad(U).T()))) Gauss linear;
     div((((((alpha.water*rho)*nu.water)|(alpha.water*rho))*rho)*dev(grad(U).T()))) Gauss linear;
     div(((nu*rho)*dev(grad(U).T()))) Gauss linear;
+	default Gauss linear;
 }
 
 laplacianSchemes