diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowSchillerNaumann/GidaspowSchillerNaumann.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowSchillerNaumann/GidaspowSchillerNaumann.C
index b918d9ac5d..0571242d55 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowSchillerNaumann/GidaspowSchillerNaumann.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowSchillerNaumann/GidaspowSchillerNaumann.C
@@ -68,7 +68,7 @@ Foam::dragModels::GidaspowSchillerNaumann::CdRe() const
     volScalarField Re(alpha2*pair_.Re());
     volScalarField CdsRe
     (
-        neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687))
+        neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687))/alpha2
       + pos(Re - 1000)*0.44*max(Re, residualRe_)
     );
 
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties b/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties
index 111a730d93..07252b178c 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties
@@ -40,7 +40,7 @@ blending
     default
     {
         type            none;
-        residualAlpha   1e-6;
+        residualAlpha   1e-3;
         continuousPhase air;
     }
 }
@@ -73,7 +73,7 @@ virtualMass
     (particles in air)
     {
         type            constantCoefficient;
-        Cvm             0.5;
+        Cvm             0;
     }
 );
 
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties b/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties
index 111a730d93..a1aec908d9 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties
@@ -40,7 +40,7 @@ blending
     default
     {
         type            none;
-        residualAlpha   1e-6;
+        residualAlpha   1e-3;
         continuousPhase air;
     }
 }