diff --git a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.air b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.air
index 40de4981f5..fe69aba1ea 100644
--- a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.air
+++ b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.air
@@ -19,10 +19,16 @@ simulationType  RAS;
 
 RAS
 {
-    RASModel continuousGasKEpsilon;
+    RASModel mixtureKEpsilon; //continuousGasKEpsilon;
 
     turbulence      on;
     printCoeffs     on;
+
+    // mixtureKEpsilonCoeffs
+    // {
+    //     Cp  1;
+    //     C3  1;
+    // }
 }
 
 // ************************************************************************* //
diff --git a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.water b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.water
index 2c0b129e11..b82e81c0cb 100644
--- a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.water
+++ b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/constant/turbulenceProperties.water
@@ -19,7 +19,7 @@ simulationType  RAS;
 
 RAS
 {
-    RASModel LaheyKEpsilon;
+    RASModel mixtureKEpsilon; //LaheyKEpsilon;
 
     turbulence      on;
     printCoeffs     on;
diff --git a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/system/fvSchemes b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/system/fvSchemes
index f55c52ee56..d3f22a16e8 100644
--- a/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/system/fvSchemes
+++ b/tutorials/multiphase/compressibleTwoPhaseEulerFoam/RAS/bubbleColumn/system/fvSchemes
@@ -40,6 +40,7 @@ divSchemes
     "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
 
     "div\(alphaPhi.*,(k|epsilon).*\)"  Gauss limitedLinear 1;
+    "div\(phim,(k|epsilon)m\)"    Gauss limitedLinear 1;
 
     "div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
 }