diff --git a/tutorials/incompressible/simpleFoam/bump2D/0/epsilon b/tutorials/incompressible/simpleFoam/bump2D/0/epsilon
new file mode 100644
index 0000000000..8726020aa2
--- /dev/null
+++ b/tutorials/incompressible/simpleFoam/bump2D/0/epsilon
@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      epsilon;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 0 2 -3 0 0 0 0 ];
+
+internalField   uniform 469.872; // computed from omega and k
+
+boundaryField
+{
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+
+    top
+    {
+        type            symmetryPlane;
+    }
+
+    "(symUp|symDown)"
+    {
+        type            symmetryPlane;
+    }
+
+    bump
+    {
+        type            epsilonWallFunction;
+        value           $internalField;
+        lowReCorrection on;
+    }
+
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/incompressible/simpleFoam/bump2D/0/f b/tutorials/incompressible/simpleFoam/bump2D/0/f
new file mode 100644
index 0000000000..60bc4e10d1
--- /dev/null
+++ b/tutorials/incompressible/simpleFoam/bump2D/0/f
@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      f;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 0 0 -1 0 0 0 0 ];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    inlet
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    top
+    {
+        type            symmetryPlane;
+    }
+
+    "(symUp|symDown)"
+    {
+        type            symmetryPlane;
+    }
+
+    bump
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/incompressible/simpleFoam/bump2D/0/nuTilda b/tutorials/incompressible/simpleFoam/bump2D/0/nuTilda
index 0e45430030..35a8d6c028 100644
--- a/tutorials/incompressible/simpleFoam/bump2D/0/nuTilda
+++ b/tutorials/incompressible/simpleFoam/bump2D/0/nuTilda
@@ -45,7 +45,7 @@ boundaryField
 
     bump
     {
-        type            fixedValue; 
+        type            fixedValue;
         value           uniform 0;
     }
 
diff --git a/tutorials/incompressible/simpleFoam/bump2D/0/p b/tutorials/incompressible/simpleFoam/bump2D/0/p
index 14d7280e5b..e50954ce2e 100644
--- a/tutorials/incompressible/simpleFoam/bump2D/0/p
+++ b/tutorials/incompressible/simpleFoam/bump2D/0/p
@@ -43,7 +43,7 @@ boundaryField
 
     bump
     {
-        type            zeroGradient; 
+        type            zeroGradient;
     }
 
     frontAndBack
diff --git a/tutorials/incompressible/simpleFoam/bump2D/0/phit b/tutorials/incompressible/simpleFoam/bump2D/0/phit
new file mode 100644
index 0000000000..b152668137
--- /dev/null
+++ b/tutorials/incompressible/simpleFoam/bump2D/0/phit
@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      phit;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 0 0 0 0 0 0 0 ];
+
+internalField   uniform 0.66;
+
+boundaryField
+{
+    inlet
+    {
+        type            zeroGradient;
+    }
+
+    outlet
+    {
+        type            zeroGradient;
+    }
+
+    top
+    {
+        type            symmetryPlane;
+    }
+
+    "(symUp|symDown)"
+    {
+        type            symmetryPlane;
+    }
+
+    bump
+    {
+        type            fixedValue;
+        value           uniform 1e-10;
+    }
+
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/incompressible/simpleFoam/bump2D/Allrun b/tutorials/incompressible/simpleFoam/bump2D/Allrun
index 32ad6c2947..962a15de40 100755
--- a/tutorials/incompressible/simpleFoam/bump2D/Allrun
+++ b/tutorials/incompressible/simpleFoam/bump2D/Allrun
@@ -10,6 +10,7 @@ runApplication blockMesh
 models="
 kOmegaSST
 SpalartAllmaras
+kEpsilonPhitF
 "
 
 # Extract a value (Eg, from boundaryField/bump/value)
diff --git a/tutorials/incompressible/simpleFoam/bump2D/plot b/tutorials/incompressible/simpleFoam/bump2D/plot
index a400d328e3..52f8966b27 100755
--- a/tutorials/incompressible/simpleFoam/bump2D/plot
+++ b/tutorials/incompressible/simpleFoam/bump2D/plot
@@ -4,7 +4,7 @@ cd "${0%/*}" || exit                                # Run from this directory
 #------------------------------------------------------------------------------
 
 # Turbulence closure models
-models=('kOmegaSST' 'SpalartAllmaras')
+models=('kOmegaSST' 'SpalartAllmaras' 'kEpsilonPhitF')
 
 # Note: CFL3D data available from:
 # https://turbmodels.larc.nasa.gov/bump_sa.html
diff --git a/tutorials/incompressible/simpleFoam/bump2D/system/fvSchemes b/tutorials/incompressible/simpleFoam/bump2D/system/fvSchemes
index dfd865261c..a9cd9a7c3b 100644
--- a/tutorials/incompressible/simpleFoam/bump2D/system/fvSchemes
+++ b/tutorials/incompressible/simpleFoam/bump2D/system/fvSchemes
@@ -27,12 +27,14 @@ gradSchemes
 
 divSchemes
 {
-    default         none;
-    div(phi,U)      bounded Gauss linearUpwind grad(U);
+    default          none;
+    div(phi,U)       bounded Gauss linearUpwind grad(U);
     div(phi,epsilon) bounded Gauss upwind;
-    div(phi,omega) bounded Gauss upwind;
-    div(phi,k)     bounded Gauss upwind;
+    div(phi,omega)   bounded Gauss upwind;
+    div(phi,k)       bounded Gauss upwind;
     div(phi,nuTilda) bounded Gauss upwind;
+    div(phi,phit)    bounded Gauss upwind;
+    div(phi,f)    bounded Gauss upwind;
     div((nuEff*dev2(T(grad(U))))) Gauss linear;
 }
 
diff --git a/tutorials/incompressible/simpleFoam/bump2D/system/fvSolution b/tutorials/incompressible/simpleFoam/bump2D/system/fvSolution
index 643f69a990..eecd3ada6a 100644
--- a/tutorials/incompressible/simpleFoam/bump2D/system/fvSolution
+++ b/tutorials/incompressible/simpleFoam/bump2D/system/fvSolution
@@ -25,13 +25,21 @@ solvers
         relTol          0.1;
     }
 
-    "(U|k|epsilon|omega|nuTilda)"
+    "(U|k|epsilon|omega|nuTilda|phit)"
     {
         solver          PBiCGStab;
         preconditioner  DILU;
         tolerance       1e-8;
         relTol          0;
     }
+
+    f
+    {
+        solver          PBiCGStab;
+        preconditioner  DIC;
+        tolerance       1e-8;
+        relTol          0;
+    }
 }
 
 SIMPLE
@@ -42,7 +50,7 @@ SIMPLE
     {
         p               1e-5;
         U               1e-5;
-        "(k|epsilon|omega|nuTilda)" 1e-5;
+        "(k|epsilon|omega|nuTilda|phit|f)" 1e-5;
     }
 }
 
@@ -51,7 +59,7 @@ relaxationFactors
     equations
     {
         U               0.9;
-        "(k|epsilon|omega|nuTilda)" 0.7;
+        "(k|epsilon|omega|nuTilda|phit|f)" 0.7;
     }
 }
 
@@ -59,6 +67,8 @@ relaxationFactors
 cache
 {
     grad(U);
+    grad(k);
+    grad(phit);
 }
 
 // ************************************************************************* //