diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/system/fvSchemes b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/system/fvSchemes
index 63b95e2727..c06c406b95 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/system/fvSchemes
+++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/system/fvSchemes
@@ -44,7 +44,7 @@ divSchemes
 
     div((((alpha.gas*thermo:rho.gas)*nuEff.gas)*dev2(T(grad(U.gas))))) Gauss linear;
 
-     div((((thermo:rho.solids*nut.solids)*dev2(T(grad(U.solids))))+(((thermo:rho.solids*lambda.solids)*div(phi.solids))*I)))  Gauss linear;
+    divDevTau(U.solids)             Gauss linear;
 }
 
 laplacianSchemes
diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/system/fvSchemes b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/system/fvSchemes
index 05e5407e18..84d77f9e61 100644
--- a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/system/fvSchemes
+++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/system/fvSchemes
@@ -44,7 +44,7 @@ divSchemes
 
     div((((alpha.air*thermo:rho.air)*nuEff.air)*dev2(T(grad(U.air))))) Gauss linear;
 
-    div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I)))  Gauss linear;
+    divDevTau(U.particles)          Gauss linear;
 }
 
 laplacianSchemes