diff --git a/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSchemes b/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSchemes
index f020e7c415..98ffbd757f 100644
--- a/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSchemes
+++ b/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSchemes
@@ -31,10 +31,12 @@ divSchemes
     default             none;
 
     div(phi,U)          bounded Gauss linearUpwind limited;
+    div(phi,h)          bounded Gauss linearUpwind limited;
     div(phi,e)          bounded Gauss linearUpwind limited;
     div(phi,epsilon)    bounded Gauss linearUpwind limited;
     div(phi,k)          bounded Gauss linearUpwind limited;
     div(phi,Ekp)        bounded Gauss linearUpwind limited;
+    div(phi,K)          bounded Gauss linearUpwind limited;
 
     div(((rho*nuEff)*dev2(T(grad(U)))))      Gauss linear;
 }
diff --git a/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSolution b/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSolution
index d518c654f5..911b47bac9 100644
--- a/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSolution
+++ b/tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSolution
@@ -25,7 +25,7 @@ solvers
         relTol          0.1;
     }
 
-    "(U|e|k|epsilon)"
+    "(U|h|e|k|epsilon)"
     {
         solver          smoothSolver;
         smoother        symGaussSeidel;
@@ -46,7 +46,7 @@ SIMPLE
     {
         p               1e-4;
         U               1e-4;
-        "(k|omega|epsilon|e|h)" 1e-3;
+        "(k|omega|epsilon|h|e)" 1e-3;
     }
 }
 
@@ -59,6 +59,7 @@ relaxationFactors
     equations
     {
         U               0.7;
+        h               0.7;
         e               0.7;
         k               0.7;
         epsilon         0.7;