diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/0/omega.gas b/tutorials/modules/multiphaseEuler/wallBoiling/0/omega.gas
index 840b72f0c8..52493f420f 100644
--- a/tutorials/modules/multiphaseEuler/wallBoiling/0/omega.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/0/omega.gas
@@ -10,7 +10,7 @@ FoamFile
     format      ascii;
     class       volScalarField;
     location    "0";
-    object      omega.liquid;
+    object      omega.gas;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/Allrun b/tutorials/modules/multiphaseEuler/wallBoiling/Allrun
index 2b3eeca408..aa7c72cb87 100755
--- a/tutorials/modules/multiphaseEuler/wallBoiling/Allrun
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/Allrun
@@ -31,4 +31,16 @@ runApplication foamPostProcess -latestTime -func "
 
 ./validation/createGraphs
 
+runApplication -append foamPostProcess -latestTime -func "
+    patchSurface
+    (
+        funcName=patchWallBoilingProperties,
+        patch=wall,
+        surfaceFormat=raw,
+        interpolate=false,
+        fields=(dDeparture.liquid fDeparture.liquid nucleationSiteDensity.liquid fLiquid.liquid quenchingHeatFlux.liquid evaporativeHeatFlux.liquid)
+    )"
+
+./validation/createWallBoilingPropertiesGraphs
+
 #------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/constant/fvModels b/tutorials/modules/multiphaseEuler/wallBoiling/constant/fvModels
new file mode 100644
index 0000000000..966ca382a8
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/constant/fvModels
@@ -0,0 +1,38 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      fvModels;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phaseTurbulenceStabilisationGas
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas;
+
+    alphaInversion  0.1;
+}
+
+phaseTurbulenceStabilisationLiquid
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   liquid;
+
+    alphaInversion  0.1;
+}
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/constant/momentumTransport.gas b/tutorials/modules/multiphaseEuler/wallBoiling/constant/momentumTransport.gas
index bceba9af18..d0b2ee96f5 100644
--- a/tutorials/modules/multiphaseEuler/wallBoiling/constant/momentumTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/constant/momentumTransport.gas
@@ -14,6 +14,15 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  laminar;
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+    printCoeffs     on;
+}
+
 
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/constant/thermophysicalTransport.gas b/tutorials/modules/multiphaseEuler/wallBoiling/constant/thermophysicalTransport.gas
index 6126ac6692..058612ec9b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoiling/constant/thermophysicalTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/constant/thermophysicalTransport.gas
@@ -19,4 +19,11 @@ laminar
     model Fourier;
 }
 
+RAS
+{
+    model   eddyDiffusivity;
+
+    Prt     1;
+}
+
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/system/controlDict.orig b/tutorials/modules/multiphaseEuler/wallBoiling/system/controlDict.orig
index cc905e581f..284437de28 100644
--- a/tutorials/modules/multiphaseEuler/wallBoiling/system/controlDict.orig
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/system/controlDict.orig
@@ -54,6 +54,13 @@ maxDeltaT       0.001;
 
 functions
 {
+    writeWallBoilingProperties
+    {
+        type            wallBoilingProperties;
+        functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" );
+        writeControl    writeTime;
+        phase           liquid;
+    }
     outflow
     {
         type            surfaceFieldValue;
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/system/fvSolution b/tutorials/modules/multiphaseEuler/wallBoiling/system/fvSolution
index 33cf3e5667..c4a1580136 100644
--- a/tutorials/modules/multiphaseEuler/wallBoiling/system/fvSolution
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/system/fvSolution
@@ -79,7 +79,7 @@ relaxationFactors
     equations
     {
         ".*"            1;
-        "h\..*"         1.0;
+        "h\..*"         1;
     }
 }
 
diff --git a/tutorials/modules/multiphaseEuler/wallBoiling/validation/createWallBoilingPropertiesGraphs b/tutorials/modules/multiphaseEuler/wallBoiling/validation/createWallBoilingPropertiesGraphs
new file mode 100755
index 0000000000..dcb8b8e97f
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoiling/validation/createWallBoilingPropertiesGraphs
@@ -0,0 +1,49 @@
+#!/bin/sh
+
+if ! which gnuplot > /dev/null 2>&1
+then
+    echo 'gnuplot not found - skipping graph creation' >&2
+    exit 1
+fi
+
+graphFile=$(foamListTimes -latestTime)/patch.xy
+
+gnuplot<<EOF
+
+    set terminal postscript eps size 8,9 color enhanced font "Helvetica,20"
+    set output "./validation/$(basename "$PWD")Properties.eps"
+    set multiplot layout 2,2
+    set decimalsign '.'
+    set grid
+
+    set key at graph 0.65,0.95
+    set ylabel 'Vertical coordinate (m)'
+    set yrange [0:4]
+
+    set xlabel 'Departure diameter (mm)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u (\$4*1000):1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Departure frequency (Hz)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Nucleation Site Density (1/m^2)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Heat flux partitioning function, liquid (-)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 7:1 w lp lt 1 t 'Simulation' \
+
+    set key at graph 0.99,0.95
+
+    unset multiplot
+
+EOF
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.gas b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.gas
index cae32b7186..284cc88317 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.gas
@@ -16,7 +16,7 @@ FoamFile
 
 dimensions      [0 1 -1 0 0 0 0];
 
-internalField   uniform (1.96244 0 0);
+internalField   uniform (1.75175 0 0);
 
 boundaryField
 {
@@ -24,13 +24,13 @@ boundaryField
     {
         type            mappedInternalValue;
         interpolationScheme cell;
-        value           uniform (1.96244 0 0);
+        value           uniform (1.75175 0 0);
     }
     outlet
     {
         type            pressureInletOutletVelocity;
         phi             phi.gas;
-        value           uniform (1.96244 0 0);
+        value           uniform (1.75175 0 0);
     }
     wall
     {
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.liquid b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.liquid
index 04de7ede06..7f580af892 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.liquid
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/U.liquid
@@ -16,22 +16,22 @@ FoamFile
 
 dimensions      [0 1 -1 0 0 0 0];
 
-internalField   uniform (1.743138395 0 0);
+internalField   uniform (1.75175 0 0);
 
 boundaryField
 {
     inlet
     {
         type            mappedInternalValue;
-        average         (1.743138395 0 0);
+        average         (1.75175 0 0);
         interpolationScheme cell;
-        value           uniform (1.743138395 0 0);
+        value           uniform (1.75175 0 0);
     }
     outlet
     {
         type            pressureInletOutletVelocity;
         phi             phi.liquid;
-        value           uniform (1.743138395 0 0);
+        value           uniform (1.75175 0 0);
     }
     wall
     {
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/omega.gas b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/omega.gas
new file mode 100644
index 0000000000..840b72f0c8
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/0/omega.gas
@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 0.01;
+
+boundaryField
+{
+    inlet
+    {
+        type            mappedInternalValue;
+        interpolationScheme cell;
+        value           uniform 0.01;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      uniform 0.0015;
+        value           uniform 0.01;
+    }
+    wall
+    {
+        type            omegaWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0.001;
+    }
+    front
+    {
+        type            wedge;
+    }
+    back
+    {
+        type            wedge;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/Allrun b/tutorials/modules/multiphaseEuler/wallBoilingIATE/Allrun
index e706c9ea8d..10cc7832ee 100755
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/Allrun
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/Allrun
@@ -31,4 +31,16 @@ runApplication foamPostProcess -latestTime -func "
 
 ./validation/createGraphs
 
+runApplication -append foamPostProcess -latestTime -func "
+    patchSurface
+    (
+        funcName=patchWallBoilingProperties,
+        patch=wall,
+        surfaceFormat=raw,
+        interpolate=false,
+        fields=(dDeparture.liquid fDeparture.liquid nucleationSiteDensity.liquid fLiquid.liquid quenchingHeatFlux.liquid evaporativeHeatFlux.liquid)
+    )"
+
+./validation/createWallBoilingPropertiesGraphs
+
 #------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/fvModels b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/fvModels
new file mode 100644
index 0000000000..966ca382a8
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/fvModels
@@ -0,0 +1,38 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      fvModels;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phaseTurbulenceStabilisationGas
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas;
+
+    alphaInversion  0.1;
+}
+
+phaseTurbulenceStabilisationLiquid
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   liquid;
+
+    alphaInversion  0.1;
+}
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/momentumTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/momentumTransport.gas
index bceba9af18..d0b2ee96f5 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/momentumTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/momentumTransport.gas
@@ -14,6 +14,15 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  laminar;
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+    printCoeffs     on;
+}
+
 
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/phaseProperties b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/phaseProperties
index 00c2373502..a512ec8673 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/phaseProperties
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/phaseProperties
@@ -98,15 +98,12 @@ heatTransfer
 {
     gas_dispersedIn_liquid_inThe_gas
     {
-        type constantNu;
-        Nu  1e1;
-        residualAlpha 1e-4;
+        type spherical;
     }
 
     gas_dispersedIn_liquid_inThe_liquid
     {
         type RanzMarshall;
-        residualAlpha 1e-4;
     }
 }
 
@@ -139,7 +136,7 @@ lift
 
 wallLubrication
 {
-   gas_dispersedIn_liquid
+    gas_dispersedIn_liquid
     {
         type            Antal;
         Cw1             -0.01;
@@ -203,7 +200,4 @@ surfaceTension
 phaseTransfer
 {}
 
-interfaceCompression
-{}
-
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/thermophysicalTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/thermophysicalTransport.gas
index 6126ac6692..058612ec9b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/thermophysicalTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/constant/thermophysicalTransport.gas
@@ -19,4 +19,11 @@ laminar
     model Fourier;
 }
 
+RAS
+{
+    model   eddyDiffusivity;
+
+    Prt     1;
+}
+
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/system/controlDict.orig b/tutorials/modules/multiphaseEuler/wallBoilingIATE/system/controlDict.orig
index 368b03ab22..e3da7f558f 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingIATE/system/controlDict.orig
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/system/controlDict.orig
@@ -54,6 +54,13 @@ maxDeltaT       0.001;
 
 functions
 {
+    writeWallBoilingProperties
+    {
+        type            wallBoilingProperties;
+        functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" );
+        writeControl    writeTime;
+        phase           liquid;
+    }
     outflow
     {
         type            surfaceFieldValue;
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingIATE/validation/createWallBoilingPropertiesGraphs b/tutorials/modules/multiphaseEuler/wallBoilingIATE/validation/createWallBoilingPropertiesGraphs
new file mode 100755
index 0000000000..dcb8b8e97f
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingIATE/validation/createWallBoilingPropertiesGraphs
@@ -0,0 +1,49 @@
+#!/bin/sh
+
+if ! which gnuplot > /dev/null 2>&1
+then
+    echo 'gnuplot not found - skipping graph creation' >&2
+    exit 1
+fi
+
+graphFile=$(foamListTimes -latestTime)/patch.xy
+
+gnuplot<<EOF
+
+    set terminal postscript eps size 8,9 color enhanced font "Helvetica,20"
+    set output "./validation/$(basename "$PWD")Properties.eps"
+    set multiplot layout 2,2
+    set decimalsign '.'
+    set grid
+
+    set key at graph 0.65,0.95
+    set ylabel 'Vertical coordinate (m)'
+    set yrange [0:4]
+
+    set xlabel 'Departure diameter (mm)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u (\$4*1000):1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Departure frequency (Hz)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Nucleation Site Density (1/m^2)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Heat flux partitioning function, liquid (-)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 7:1 w lp lt 1 t 'Simulation' \
+
+    set key at graph 0.99,0.95
+
+    unset multiplot
+
+EOF
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/Allrun b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/Allrun
index 887f4530a2..05a5c3504b 100755
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/Allrun
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/Allrun
@@ -32,4 +32,16 @@ runApplication foamPostProcess -latestTime -func "
 
 ./validation/createGraphs
 
+runApplication -append foamPostProcess -latestTime -func "
+    patchSurface
+    (
+        funcName=patchWallBoilingProperties,
+        patch=wall,
+        surfaceFormat=raw,
+        interpolate=false,
+        fields=(dDeparture.liquid fDeparture.liquid nucleationSiteDensity.liquid fLiquid.liquid quenchingHeatFlux.liquid evaporativeHeatFlux.liquid)
+    )"
+
+./validation/createWallBoilingPropertiesGraphs
+
 #------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/fvModels b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/fvModels
new file mode 100644
index 0000000000..966ca382a8
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/fvModels
@@ -0,0 +1,38 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      fvModels;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phaseTurbulenceStabilisationGas
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas;
+
+    alphaInversion  0.1;
+}
+
+phaseTurbulenceStabilisationLiquid
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   liquid;
+
+    alphaInversion  0.1;
+}
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/momentumTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/momentumTransport.gas
index bceba9af18..d0b2ee96f5 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/momentumTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/momentumTransport.gas
@@ -14,6 +14,15 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  laminar;
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+    printCoeffs     on;
+}
+
 
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/phaseProperties b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/phaseProperties
index 26a0bcc4e4..50d4c8cdbe 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/phaseProperties
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/phaseProperties
@@ -246,7 +246,4 @@ surfaceTension
 phaseTransfer
 {}
 
-interfaceCompression
-{}
-
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/thermophysicalTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/thermophysicalTransport.gas
index 6126ac6692..058612ec9b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/thermophysicalTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/constant/thermophysicalTransport.gas
@@ -19,4 +19,11 @@ laminar
     model Fourier;
 }
 
+RAS
+{
+    model   eddyDiffusivity;
+
+    Prt     1;
+}
+
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/system/controlDict.orig b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/system/controlDict.orig
index 49d0831464..aa211b5811 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/system/controlDict.orig
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/system/controlDict.orig
@@ -54,6 +54,13 @@ maxDeltaT       0.001;
 
 functions
 {
+    writeWallBoilingProperties
+    {
+        type            wallBoilingProperties;
+        functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" );
+        writeControl    writeTime;
+        phase           liquid;
+    }
     volumeDensity.diameter.bubbles
     {
         type            populationBalanceSizeDistribution;
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/validation/createWallBoilingPropertiesGraphs b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/validation/createWallBoilingPropertiesGraphs
new file mode 100755
index 0000000000..dcb8b8e97f
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperse/validation/createWallBoilingPropertiesGraphs
@@ -0,0 +1,49 @@
+#!/bin/sh
+
+if ! which gnuplot > /dev/null 2>&1
+then
+    echo 'gnuplot not found - skipping graph creation' >&2
+    exit 1
+fi
+
+graphFile=$(foamListTimes -latestTime)/patch.xy
+
+gnuplot<<EOF
+
+    set terminal postscript eps size 8,9 color enhanced font "Helvetica,20"
+    set output "./validation/$(basename "$PWD")Properties.eps"
+    set multiplot layout 2,2
+    set decimalsign '.'
+    set grid
+
+    set key at graph 0.65,0.95
+    set ylabel 'Vertical coordinate (m)'
+    set yrange [0:4]
+
+    set xlabel 'Departure diameter (mm)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u (\$4*1000):1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Departure frequency (Hz)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Nucleation Site Density (1/m^2)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Heat flux partitioning function, liquid (-)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 7:1 w lp lt 1 t 'Simulation' \
+
+    set key at graph 0.99,0.95
+
+    unset multiplot
+
+EOF
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/Allrun b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/Allrun
index e49b6660a9..421afa9512 100755
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/Allrun
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/Allrun
@@ -32,4 +32,16 @@ runApplication foamPostProcess -latestTime -func "
 
 ./validation/createGraphs
 
+runApplication -append foamPostProcess -latestTime -func "
+    patchSurface
+    (
+        funcName=patchWallBoilingProperties,
+        patch=wall,
+        surfaceFormat=raw,
+        interpolate=false,
+        fields=(dDeparture.liquid fDeparture.liquid nucleationSiteDensity.liquid fLiquid.liquid quenchingHeatFlux.liquid evaporativeHeatFlux.liquid)
+    )"
+
+./validation/createWallBoilingPropertiesGraphs
+
 #------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/fvModels b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/fvModels
new file mode 100644
index 0000000000..6c57595ddd
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/fvModels
@@ -0,0 +1,49 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      fvModels;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phaseTurbulenceStabilisationGas
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas;
+
+    alphaInversion  0.1;
+}
+
+phaseTurbulenceStabilisationGas2
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas2;
+
+    alphaInversion  0.1;
+}
+
+phaseTurbulenceStabilisationLiquid
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   liquid;
+
+    alphaInversion  0.1;
+}
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas
index bceba9af18..d0b2ee96f5 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas
@@ -14,6 +14,15 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  laminar;
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+    printCoeffs     on;
+}
+
 
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas2 b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas2
index bceba9af18..12414285ea 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas2
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/momentumTransport.gas2
@@ -14,6 +14,14 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  laminar;
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+    printCoeffs     on;
+}
 
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/phaseProperties b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/phaseProperties
index 67628ac2f1..1c595d384b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/phaseProperties
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/phaseProperties
@@ -192,13 +192,11 @@ heatTransfer
     gas_dispersedIn_liquid_inThe_liquid
     {
         type RanzMarshall;
-        residualAlpha 1e-4;
     }
 
     gas2_dispersedIn_liquid_inThe_liquid
     {
         type RanzMarshall;
-        residualAlpha 1e-4;
     }
 }
 
@@ -230,26 +228,7 @@ lift
 
     gas2_dispersedIn_liquid
     {
-        type wallDamped;
-
-        lift
-        {
-            type Tomiyama;
-            Cl 0.288;
-
-            aspectRatio
-            {
-                type constant;
-                E0 1;
-            }
-        }
-
-        wallDamping
-        {
-            type    cosine;
-            Cd      1.0;
-            zeroWallDist 0.0002;
-        }
+        $gas_dispersedIn_liquid
     }
 }
 
@@ -264,9 +243,7 @@ wallLubrication
 
     gas2_dispersedIn_liquid
     {
-        type            Antal;
-        Cw1             -0.01;
-        Cw2             0.05;
+        $gas_dispersedIn_liquid
     }
 }
 
@@ -322,24 +299,7 @@ saturation
 
     gas2_liquid
     {
-        type function1;
-
-        function scale;
-
-        xScale 1e-6;
-        scale   1;
-        value
-        {
-            type                tableFile;
-            format              csv;
-            nHeaderLine         1;
-            refColumn           1;
-            componentColumns    (0);
-            mergeSeparators     no;
-            file "$FOAM_TUTORIALS/resources/thermoData/wallBoiling-saturation.csv";
-            outOfBounds         clamp;
-            interpolationScheme linear;
-        }
+        $gas_liquid
     }
 }
 
@@ -353,8 +313,7 @@ surfaceTension
 
     gas2_liquid
     {
-        type            constant;
-        sigma           0.00176574;
+        $gas_liquid
     }
 }
 
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas
index 6126ac6692..058612ec9b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas
@@ -19,4 +19,11 @@ laminar
     model Fourier;
 }
 
+RAS
+{
+    model   eddyDiffusivity;
+
+    Prt     1;
+}
+
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas2 b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas2
index 6126ac6692..058612ec9b 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas2
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/constant/thermophysicalTransport.gas2
@@ -19,4 +19,11 @@ laminar
     model Fourier;
 }
 
+RAS
+{
+    model   eddyDiffusivity;
+
+    Prt     1;
+}
+
 // ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/system/controlDict.orig b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/system/controlDict.orig
index 93f0e68dde..d593c66f8d 100644
--- a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/system/controlDict.orig
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/system/controlDict.orig
@@ -54,6 +54,13 @@ maxDeltaT       0.001;
 
 functions
 {
+    writeWallBoilingProperties
+    {
+        type            wallBoilingProperties;
+        functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" );
+        writeControl    writeTime;
+        phase           liquid;
+    }
     volumeDensity.diameter.bubbles
     {
         type            populationBalanceSizeDistribution;
diff --git a/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/validation/createWallBoilingPropertiesGraphs b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/validation/createWallBoilingPropertiesGraphs
new file mode 100755
index 0000000000..dcb8b8e97f
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/wallBoilingPolydisperseTwoGroups/validation/createWallBoilingPropertiesGraphs
@@ -0,0 +1,49 @@
+#!/bin/sh
+
+if ! which gnuplot > /dev/null 2>&1
+then
+    echo 'gnuplot not found - skipping graph creation' >&2
+    exit 1
+fi
+
+graphFile=$(foamListTimes -latestTime)/patch.xy
+
+gnuplot<<EOF
+
+    set terminal postscript eps size 8,9 color enhanced font "Helvetica,20"
+    set output "./validation/$(basename "$PWD")Properties.eps"
+    set multiplot layout 2,2
+    set decimalsign '.'
+    set grid
+
+    set key at graph 0.65,0.95
+    set ylabel 'Vertical coordinate (m)'
+    set yrange [0:4]
+
+    set xlabel 'Departure diameter (mm)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u (\$4*1000):1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Departure frequency (Hz)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Nucleation Site Density (1/m^2)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 5:1 w lp lt 1 t 'Simulation' \
+
+    set xlabel 'Heat flux partitioning function, liquid (-)'
+    plot \
+        "postProcessing/patchWallBoilingProperties/$graphFile" \
+        u 7:1 w lp lt 1 t 'Simulation' \
+
+    set key at graph 0.99,0.95
+
+    unset multiplot
+
+EOF
+
+#------------------------------------------------------------------------------