From 9fa8b85056229c7e6e9f2f53d231a665172c2f7f Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Tue, 15 Nov 2022 15:59:44 +0000
Subject: [PATCH] tutorials: multiphaseEuler: Added hydrofoil tutorial

This serves as an example of cavitation modelling with the
multiphaseEuler module. This case also contains validation of the
pressure profile along the hydrofoil against experimental data.

Based on a case contributed by Petteri Peltonen, VTT.
---
 .../modules/multiphaseEuler/hydrofoil/0/T.gas |  49 ++++
 .../multiphaseEuler/hydrofoil/0/T.liquid      |  49 ++++
 .../multiphaseEuler/hydrofoil/0/U.gas.orig    |  46 ++++
 .../multiphaseEuler/hydrofoil/0/U.liquid.orig |  46 ++++
 .../multiphaseEuler/hydrofoil/0/alpha.gas     |  49 ++++
 .../multiphaseEuler/hydrofoil/0/alpha.liquid  |  49 ++++
 .../multiphaseEuler/hydrofoil/0/alphat.gas    |  57 ++++
 .../multiphaseEuler/hydrofoil/0/alphat.liquid |  57 ++++
 .../modules/multiphaseEuler/hydrofoil/0/k.gas |  51 ++++
 .../multiphaseEuler/hydrofoil/0/k.liquid      |  52 ++++
 .../multiphaseEuler/hydrofoil/0/nut.gas       |  49 ++++
 .../multiphaseEuler/hydrofoil/0/nut.liquid    |  49 ++++
 .../multiphaseEuler/hydrofoil/0/omega.gas     |  60 ++++
 .../multiphaseEuler/hydrofoil/0/omega.liquid  |  59 ++++
 .../modules/multiphaseEuler/hydrofoil/0/p     |  49 ++++
 .../modules/multiphaseEuler/hydrofoil/0/p_rgh |  47 ++++
 .../multiphaseEuler/hydrofoil/Allclean        |   7 +
 .../modules/multiphaseEuler/hydrofoil/Allrun  |  20 ++
 .../hydrofoil/constant/fvModels               |  28 ++
 .../multiphaseEuler/hydrofoil/constant/g      |  20 ++
 .../geometry/NACAMOD66_6deg_lower.stl.gz      | Bin 0 -> 3026 bytes
 .../geometry/NACAMOD66_6deg_tip.stl.gz        | Bin 0 -> 6831 bytes
 .../geometry/NACAMOD66_6deg_upper.stl.gz      | Bin 0 -> 4108 bytes
 .../hydrofoil/constant/momentumTransport.gas  |  28 ++
 .../constant/momentumTransport.liquid         |  28 ++
 .../hydrofoil/constant/phaseProperties        | 118 ++++++++
 .../hydrofoil/constant/physicalProperties.gas |  46 ++++
 .../constant/physicalProperties.liquid        |  50 ++++
 .../hydrofoil/system/blockMeshDict            | 257 ++++++++++++++++++
 .../hydrofoil/system/controlDict              |  63 +++++
 .../hydrofoil/system/decomposeParDict         |  21 ++
 .../hydrofoil/system/fvConstraints            |  24 ++
 .../hydrofoil/system/fvSchemes                |  68 +++++
 .../hydrofoil/system/fvSolution               | 100 +++++++
 .../hydrofoil/system/hydrofoilPressure        | 149 ++++++++++
 .../hydrofoil/validation/createGraph          |  48 ++++
 36 files changed, 1893 insertions(+)
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/T.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/T.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/U.gas.orig
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/U.liquid.orig
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/k.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/k.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/nut.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/nut.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/omega.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/omega.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/p
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/0/p_rgh
 create mode 100755 tutorials/modules/multiphaseEuler/hydrofoil/Allclean
 create mode 100755 tutorials/modules/multiphaseEuler/hydrofoil/Allrun
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/fvModels
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/g
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_lower.stl.gz
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_tip.stl.gz
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_upper.stl.gz
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/phaseProperties
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.gas
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.liquid
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/blockMeshDict
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/controlDict
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/decomposeParDict
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/fvConstraints
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/fvSchemes
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/fvSolution
 create mode 100644 tutorials/modules/multiphaseEuler/hydrofoil/system/hydrofoilPressure
 create mode 100755 tutorials/modules/multiphaseEuler/hydrofoil/validation/createGraph

diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/T.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/T.gas
new file mode 100644
index 0000000000..ec5afe5875
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/T.gas
@@ -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       volScalarField;
+    location    "0";
+    object      T.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 293.15;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.gas;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            zeroGradient;
+    }
+    walls
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/T.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/T.liquid
new file mode 100644
index 0000000000..0407bf043a
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/T.liquid
@@ -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       volScalarField;
+    location    "0";
+    object      T.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 293.15;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            zeroGradient;
+    }
+    walls
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/U.gas.orig b/tutorials/modules/multiphaseEuler/hydrofoil/0/U.gas.orig
new file mode 100644
index 0000000000..cdb30a35ad
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/U.gas.orig
@@ -0,0 +1,46 @@
+/*--------------------------------*- 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       volVectorField;
+    location    "0";
+    object      U.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (5.33 0 0);
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           uniform (5.33 0 0);
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+    hydrofoil
+    {
+        type            noSlip;
+    }
+    walls
+    {
+        type            noSlip;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/U.liquid.orig b/tutorials/modules/multiphaseEuler/hydrofoil/0/U.liquid.orig
new file mode 100644
index 0000000000..ee78804820
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/U.liquid.orig
@@ -0,0 +1,46 @@
+/*--------------------------------*- 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       volVectorField;
+    location    "0";
+    object      U.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (5.33 0 0);
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           uniform (5.33 0 0);
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+    hydrofoil
+    {
+        type            noSlip;
+    }
+    walls
+    {
+        type            noSlip;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.gas
new file mode 100644
index 0000000000..c850ac1c2c
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.gas
@@ -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       volScalarField;
+    location    "0";
+    object      alpha.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.gas;
+        inletValue      uniform 0;
+        value           uniform 0;
+    }
+    hydrofoil
+    {
+        type            zeroGradient;
+    }
+    walls
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.liquid
new file mode 100644
index 0000000000..bc54b10e67
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/alpha.liquid
@@ -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       volScalarField;
+    location    "0";
+    object      alpha.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1.0;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      uniform 1;
+        value           uniform 1;
+    }
+    hydrofoil
+    {
+        type            zeroGradient;
+    }
+    walls
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.gas
new file mode 100644
index 0000000000..75775e3f26
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.gas
@@ -0,0 +1,57 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      alphat.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -1 0 0 0 0];
+
+internalField   uniform 1e-8;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            calculated;
+        value           uniform 1e-8;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 1e-8;
+    }
+    hydrofoil
+    {
+        type            compressible::alphatPhaseJayatillekeWallFunction;
+        Prt             0.85;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    walls
+    {
+        type            compressible::alphatPhaseJayatillekeWallFunction;
+        Prt             0.85;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.liquid
new file mode 100644
index 0000000000..1ca935c585
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/alphat.liquid
@@ -0,0 +1,57 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      alphat.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 1 -1 -1 0 0 0 0 ];
+
+internalField   uniform 1e-08;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 1e-08;
+    }
+    hydrofoil
+    {
+        type            compressible::alphatPhaseJayatillekeWallFunction;
+        Prt             0.85;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    walls
+    {
+        type            compressible::alphatPhaseJayatillekeWallFunction;
+        Prt             0.85;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/k.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/k.gas
new file mode 100644
index 0000000000..7d977934c2
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/k.gas
@@ -0,0 +1,51 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      k.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0.0042613;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.gas;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            kqRWallFunction;
+        value           $internalField;
+    }
+    walls
+    {
+        type            kqRWallFunction;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/k.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/k.liquid
new file mode 100644
index 0000000000..9ec37cc5b2
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/k.liquid
@@ -0,0 +1,52 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      k.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0.0042613;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            kqRWallFunction;
+        value           $internalField;
+    }
+    walls
+    {
+        type            kqRWallFunction;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.gas
new file mode 100644
index 0000000000..9f3b3c397e
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.gas
@@ -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       volScalarField;
+    location    "0";
+    object      nut.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 9.993764000E-07;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+          type          nutUSpaldingWallFunction;
+          value         $internalField;
+    }
+    walls
+    {
+          type          nutUSpaldingWallFunction;
+          value         $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.liquid
new file mode 100644
index 0000000000..05a29fa91f
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/nut.liquid
@@ -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       volScalarField;
+    location    "0";
+    object      nut.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 9.993764000E-07;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+          type          nutUSpaldingWallFunction;
+          value         $internalField;
+    }
+    walls
+    {
+          type          nutUSpaldingWallFunction;
+          value         $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.gas b/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.gas
new file mode 100644
index 0000000000..ad361c0bad
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.gas
@@ -0,0 +1,60 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      omega.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 355.3333333;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            omegaWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        blended         on;
+        value           $internalField;
+    }
+    walls
+    {
+        type            omegaWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        blended         on;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.liquid
new file mode 100644
index 0000000000..0529839074
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/omega.liquid
@@ -0,0 +1,59 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      omega.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 355.3333333;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        phi             phi.liquid;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            omegaWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        blended         on;
+        value           $internalField;
+    }
+    walls
+    {
+        type            omegaWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        blended         on;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/p b/tutorials/modules/multiphaseEuler/hydrofoil/0/p
new file mode 100644
index 0000000000..360ab5ce19
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/p
@@ -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       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 1e5;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    outlet
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            calculated;
+        value           $internalField;
+    }
+    walls
+    {
+        type            calculated;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/0/p_rgh b/tutorials/modules/multiphaseEuler/hydrofoil/0/p_rgh
new file mode 100644
index 0000000000..f2ec3cab33
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/0/p_rgh
@@ -0,0 +1,47 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      p_rgh;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 1e5;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedFluxPressure;
+    }
+    outlet
+    {
+        type            prghPressure;
+        p               $internalField;
+        value           $internalField;
+    }
+    hydrofoil
+    {
+        type            fixedFluxPressure;
+    }
+    walls
+    {
+        type            fixedFluxPressure;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/Allclean b/tutorials/modules/multiphaseEuler/hydrofoil/Allclean
new file mode 100755
index 0000000000..77bdf0f3ff
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/Allclean
@@ -0,0 +1,7 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase && rm -rf 0/phi 0/U 0/uniform 0/yPlus.liquid validation/*.eps
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/Allrun b/tutorials/modules/multiphaseEuler/hydrofoil/Allrun
new file mode 100755
index 0000000000..4adb9f2fd4
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/Allrun
@@ -0,0 +1,20 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+runApplication blockMesh
+
+cp 0/U.liquid.orig 0/U
+runApplication potentialFoam -pName p_rgh
+cp 0/U 0/U.gas
+cp 0/U 0/U.liquid
+
+runApplication decomposePar
+runParallel $(getApplication)
+runApplication reconstructPar
+
+validation/createGraph
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/fvModels b/tutorials/modules/multiphaseEuler/hydrofoil/constant/fvModels
new file mode 100644
index 0000000000..b0c49f3e93
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/fvModels
@@ -0,0 +1,28 @@
+/*--------------------------------*- 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       dictionary;
+    object      fvOptions;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phaseTurbulenceStabilisation
+{
+    type    phaseTurbulenceStabilisation;
+
+    libs    ("libmultiphaseEulerFoamFvModels.so");
+
+    phase   gas;
+
+    alphaInversion  0.1;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/g b/tutorials/modules/multiphaseEuler/hydrofoil/constant/g
new file mode 100644
index 0000000000..a11e931c9a
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/g
@@ -0,0 +1,20 @@
+/*--------------------------------*- 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       uniformDimensionedVectorField;
+    object      g;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -2 0 0 0 0];
+value           (0 -9.81 0);
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_lower.stl.gz b/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_lower.stl.gz
new file mode 100644
index 0000000000000000000000000000000000000000..96ab117bba7d84b01dfc9aecf19220ef59ac4d58
GIT binary patch
literal 3026
zcmYLLc_0)1|CVFq7)cCs6>TfaQDlyoE9S@?<%1+w%#kypKCVoTNsi^LvD|!?QljRF
z&3#5<j<UIvT$Llf+ow;z-`}tI@jS2Rc|WhWOfnBo#j{{HHV#Xq0n*}<p}P7tb$5>&
z*L-h=c?6va4)JB{E%@9g1p#6Yew+&8sAh|+Q7RqxEJlraj@-HuxxI`CD7sYwqRCnB
zmO{FADXYOQ%AV1yb^TBN-d}tpD6_v!*S)a6|M^jW)BbY%C7H%cy0t*@Z##xL%Uei8
znG1LJ`&VV6etlcqn2&1OhaUX;K_@ju?M)|a8ldXQ`w7M9Q==D+lTVE%MZtSaY#eJg
zb@!5GqPD#w(k>}212~Viz!^V$l8O;$zdqJyM5{p>Qk2i?O}jC=s^=eV84AnbK~aTH
zhhgX_zK1@)^3<^PX1~xs$&U_|4^_g-L>RtzpgmHQRL3T$na|cTvM;D=(l67W-%9`h
zXncHxkE-ekCjPjPRo+PLBGEOHYt-w|5wKi-Tl(Py>t2|f__>~#@CDUlz;je7I+QY$
zSFv&2wj`Zf|Fo?yyoOg?p^VKY5y3?*(`qxtzO2Iy2|Sc+5*8S_H!?2VO$bc2TqOZ`
zILcRnL3p<nf8wUAo|zl8)Bz@GL)G-kCeEwG3wz){PftsSr4hkNZ*W7&=(0K_a#RR|
zNe%7MFh2dF^yHa`IiGQKMe6v3@nT7ujO6=w0O}qYC|>%>x<pv(yPn(LVY(t^oMyZj
z-vQrJfn*9-JILTX29&v5w8~3z`F0)Yz@SWnKqfO4;VcBa&8g~AyQWuu+>oNH$9N9=
zN+36)yXICMx<TdX6&&Y%i~YlHk&1Cx#Ouw`rVI%iuiAz<@T}E@U#4YB7yL&evY@hB
z=L3=sE*Q5}KYq}&bm$Mc`@7$l>f|T8H3qM?r?_TvT9sNSz19n!(6Ds%no~!eH~J@B
z<b_mmLQby1VTpMc5J&h4|B-81x_ZB5aZUja0|1>MjNbz$`h=dnX;e>4)&zG48c3Nm
zUZm(Twli(WRBM~aDkeHqr*c;DgVg;{*^3B~Fv(&ugLM#YR=P9h)tADy6gQ#Xx2KLF
z32OeJ8!-G-NMT-~t;EKWJ~9n2G^#=r3V<I7TyK@=pogVoh$PkHhRo2jr8|VzcyV7`
z4NzjP%*qqBk~4OV>w;ct5cI1$N(t-Uq=}DkpYglHRmLe0Y+|PqdCD4j;8^+xAGPCK
z1trQQJfeDqdBVKxW7x6EI~FnDU&rLn+P|Se$)Mu7!D@|{R2N?sKki%E3gRW>9P|+>
z1~g%YWe)eYwd1izA+^;GFk~xOici49fY5w9lpE1i_?iZt1T9J!D5d1%k}sorlmVv3
z1zqW8eCOJqAW~hS%=Ri#)OOFr6{*o_K?+{kagx1>mMc13OQ~A_^{<W1-;F46qr%->
zdN~`-Z^CnygmU<<3;wqQ=dUh!q-ow#bTW~XmuedpWv;ldc=SMyvm$1;m78kHOJTXM
z1MxuG?&&`-$$56MSOo9rnAsr*K5g;;Gw)#G0}tslWDnY#znBn+cx_0&QKLlE)a~fs
z26E;%%$?+QhrUz*vg$7pt*|t9pUZ-&M37H#!vNplZN*LHjs<UUC0>rmuS15MuNUK$
zPnD#Y@=%!m#->@GlfIdjKKv?zC4rUFS>jwiM)u|e@;{A?{fG-k`7th))GPuzH8<dk
zxC#8wjRWa)?haF6FUv#GeHx<Sn7-cMCy0QTE9|~-iVDhrFRFZ)3n`R5zNMOeE#5z$
z@sFFa9Hw0bW81KEivafrL0}@mlEsgTnoUSfHK0me1y7IxH6XvsiK81E#>D@xhaMi;
zF^O7#??TflV&YTAZ;y1?ONVQG{bSJ+fuB{heAnZs-|e|^LABQ5U!=--KC*}J0fW*F
zf>|#RRkF5Tg{)Ca5I&nS7|^HT6KCbKC(ewX@w_u3sv@vFXBzDua!HPsiIF<X+v44r
zSb0qNCwIjkiHZr%?T$YFWXb*qr;$Xk{mw_YtViG+;}TskEd!j(Oc>n5u8;6;+*o#d
z7U4SZ>GrEaPtWxSwt{~GXiyHfrQPWb8!M?t4L~pUtM#+%T@X9KjBqgI6XsYpbuq{X
zUAc?F2r@O~xISTiJZsIzu7C3XGi;ED`kZ~n-e1k7%FS41K~rW!xNG1EJ8?vSE+0Q*
zA=V%+pVupGM16l_ISmbn9Wz2=bDGCqX#FR1e3|kp$94C<R~Fh#1PPhG^eS}hxV2y{
zRq77}HhoEu&WUkRU={(Hnvd$Ia|w0)e8O(jF2#^dSYm+`0c3;NUS%*Gm>m9|R>I4^
ztUtLMDholJV`haUF*8v1aioVT?h32kxw^6!|1|nrlXO&cK4XBHbj~RgCppz=$iLw3
zne2qka$g;_1G7>|<~!(R<*dm}95=wj=1>_W#nm<^+tMR3lYTxaL&(RcqR&EBt>6Cs
zSKo>9XF%~2&DvqR#XTZpli#B?+6mz_Xn|19p^E~|@*94RW>BU=<o-ZMPyOLvt?%RL
z@>o|(u?RVyqbU)K<Lwx;$>nSnp2ZKr+qT$%>{474kPeQ=Fv(emptP-}N8BrSYne9J
z*FL2z?x*}LXQ?Gn*g&O}Z&GOSm+_~DknRtl<!t1AV?XIBJgLmurF%Y(E{k3AhinBA
zdKOM{7x|sel2A+0*IBbSe`^^n7Zc^U$uS2__7AAsdS2VLVJ8s!%CB;$73(egBV6tE
z=0vXJr@0yKUNmU(ljV|8xehM5F<XP?45&c2LIuJZaUmmIuz;SqmLL(Mt-6c`Q{m6I
z*ZN0m6m2TKm@90(8;R+iR*6nAYc7bSq8*wkS>9?o>loxY{CYW>2tBs6OR_0H9O1GP
zdIOPIw9Be<%@?=!<X9v_Fayem?;=+2NKN^~P^z@L(842AmP=fn8xB%)`d$)^!vu8t
z<<ECH*0&}NziKoR^W_IuR-XE&<<Cc|15hP{W;h9>mXK;+es8W=3GYnFZc~xTY0D*r
zhYae*LsE6NdR`e5;;}AmPPoo(>AOEscp|i-sbJXjxNBSZ-?5VJHsm6#zze5Q(f)?Q
zEHL!)W5tW;b5YE(b8n&HZg8}H7`rz*PL}NS=)U>v)@iG~;$KH+WfQ;NgoI`=3|9+5
zk*Yc{%tz&9PSOj!&nQ*;09y~~L5q&CI=L9XbH}1rs=}Lumr(E0&$+!H`mVX=9i-D#
zh*D6oq3=Rg=QHPd<+sipt8y+kH>RTRjjm4dGSa!J2;PxN%cWVbpU$anC!ZCVTD2Pg
zr4iN>{8JCM^HFyt;$6P`pupFs&Y37&Bf(1DE|rOaPq^t~+GCX=ML{OX+$QFyB?z9;
zSJjbY4s%YH;OLxTpMK|n`lLyLbRvW%r;Hb8qxFN`bZLK^RZE;rYN_t`NcAB`@R`Tm
z1VdELu~Je0WHh<IQ=3?BZ%YIqQMojdVJ8>6Sc0x8Ft+SFbMkW^8QMkMEYaE<E<+2M
zZ@i-z9_QH`FJdw5WnM4`9T;+wobeCuEd$_!l-JBjs&IJBwDuYsqO-7`H9G5O1K}Uy
zk?vh`4#ToM&2Lp~Gl%S8Rxfl-e&3*!FkFA7MicqkQFLMC8Cn>;HDETzF2=LR#cJZ?
zGn+0>x>*E{t9+%{SX8oTZPZooJtO<UZTCA`Y`vq#i(G_Zo1!q?=U^l2H~<2J4L`HF
zl)GWLGR<xHvJS0@cy6w@zFU9NX8paZu?pK1>23K_{I5&Y^Ruk9uuNgs7HcD+)g9-s
z4n0ug%B@K>Nx-)6-dgp)adXr_71UHW9I)cvV#VAJ?WHTIF3CpL4c_(nJ(Zesn^NjY
q!JCo%3Z5bcUGy*q3zAKAL~uvK^o6LM?b6kKL73R1=%V{<Z2tppI>PP%

literal 0
HcmV?d00001

diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_tip.stl.gz b/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_tip.stl.gz
new file mode 100644
index 0000000000000000000000000000000000000000..237c31f20a2602f181e87044fe810b15b2bfa189
GIT binary patch
literal 6831
zcmX|GcRbbq_fJAuMMieq%#08+x>mLi+1;$_8X@c2Tq`T&a_M%-Dp`?Dc1Fg%MwdSJ
z4!KsYy~po;`FwwWpT~K<p68tB^SsXcyt!hpUM(o~fL<ciQ&v^}NB@zuw6(O2t-ZCU
zgPW*_=c`NY=}rjSYGn9X0AnMG)V}8T_t3XeAtZ592~uH@?+0>u8io`VJeJ|lnuV0z
zUE!#9nWTBe;L{~S@KHDV?06;kcz*>Qe0<y;e6o`#rdU(*rG|c-d9R9`0NtNG-Tc#y
zJKI``Z-O5z;WN+b&j6p34YYOeaZ|C|&%A^9_fKQz^4{yeKK-kwLc@J(#oOE1`gqFq
zaDyw@Joh9Ddq9{y{fUnv@Q^0QpFGV(G^D{a2gxoSvuU06$IdIBnpVvhbswCaEN61f
zO|u3nwq06KZ5pC@WA4Ozc+kjwcq(?%d9(5+!7;V8cwnR)p1CAnN9gwNcXRKqYiu&h
zRy1u1Dlyb~Xa?5S1PKh{-+|011e2#@45Q88nZF5^mdKJVxHZqJ^Y!MajDZvPz4Cie
z^VQnrg%)XIN}#mQF)v(l5KTT_ytXtf`?nIHsYup>8>+YGjz@bE7V=*tm}^Ny2XgdI
z(y<&oN$`LRc@QkaskZQvpbwv8lz!^HU{7~;sJ((>xdlo7ylTkTSU5imCuRN3ZbDr_
zFe%ZqvkvS}4n;<c=A5kPbcH%)x?Wsb-1vlg;e+9?Lk1h0XS6x1lp<yEYFMRUx%+JZ
z+=bg%T&qFjmeq50b*6=RPSkZ2+uNxS^#>Y^?GvWGLOH%wwW@&ScDZrcj)5r4A})$(
z`@tU2E#*I(=(Lk1P6{HeaIJDSaeu{ja?p1JpZ_Acy+kkwTZ?EqMzXL#lDj1rV_JVg
zFYr9q1LukxEl>IiYC*}>nd}_uV`DzW;4cS~5JO?`AV9`!-k&PQ5yKo-X@6Wh^(!&-
z6#g>Y!&tlr;+Sy~GylGdTTr6>UVSMOuQ}l(eD?o+dtJmk*((8~#J&1o-rcgM!Pu1F
zcvV#(_s>NeJ#`n7d7V@ZZET+F;TEMenI!s{Khpg)XuoGRJjO2Ae&a(I^;_W{aYD5F
zRoMi|A5?u;@%dTFI$USe3kf?xWI!Z)3Ev7WXl?9$&%siABuNjq#kz_*Qq^PL*srM>
zw*8F`FPWEE!W(h<{*-uK_ZIs_cllO>N?f|`h??X0&KO?B#m<tX^HB`qG9i6an7{Y|
zicM~`*$Ck>>^4o|`U)eK0Oe8Z4B6wf`4-w}wLH^@CFgX4&)q{m_>wSv``l6N-zI(P
zI|AJ@6#Q(y>%Uo2b{sT7ItMf`-R#9%>)@6TnTvrzQT>=tsJIEOaao=s<~Eg%V@f%7
zLr(_3RHLsJynI!aO8hr{ivrh6SB^Amxm8aOW^OGCLrbC8JT`8l-+d@w2q7Y#Fj_rK
z{-?$!;_ghcX219h{cO*bu0TwVEe*%lSMy7~f-AQ7D9ZO5UL%kU9BNjkotYIqjmwg<
zmE}mC7}wibbGgbs818;Zhg;qhDn3JjKhTXTPIsW6V1pA~nUfZsU@qmxLDUV*6wO#+
z%l}$ei9gcgTLqt_Dh>|uN7w|wNZB?3dVA_g7)8L&qU64*0E9tjRQ?9vWRPH6CYL4$
zLh_+W;*H(%j|(^D+m!I88Xx+jp1We9taw{cF0x@;&pt7V&j6waelm(F4DPZ@Z2kLv
znnXiMgHesc^rsRcF>#2N2SnweRUXyuH6N9rx=uB63y5{1%4f4;e4`Q6|3T%IEeJ=*
zuR(~wV7@Ri4v?T6qTwf4{;1&3I#=4O8No^!+D3bjhM){a9GL}eHh<}ET&V!NQYm!}
z@I=o{t@UW}m>Cs+gx*opNHD|W%O5%By!7)XmT%jZ&@~`>r~81BSUrLD<b|{%{h+YA
z-N>R>^DR2#9OeWX?4vpP#amka2J?->*7H`|c63iiF&ZBrlXdK*GK_U$#N?~41)wGT
zH^A=2%OO+%p#rwn2lmH6tuz+J7lfjass7QFOP^+?`Uw@SP<=Z9`*{enP;il_n2F9`
ze6_JAYzr^$viZfGI@|sbsIusDJx7cjzK{PYShY|o7Xx5YJGQDR9+ZO3!y?$}(`~7J
zJhwgvZdx?J_Q}{%@;%0Kitl?fFOo~O$`%ML&~FKICsvI80Rai_I`}e^qk&kJ5_<-{
zt=8!li}P5P-_<A|+In`Et=9PLEvL3i9Q6UC=UX5Tn+*ux9Po=6M;vL<+xe5T45jQI
z=^39s?!+u8Cdqt-Ch(A{F!4IPv>e%~aR<tk7_-ml((FNR&@<_8`j7P^`HcuTb+htu
zMfR}$41#7mPZi3krtuv8V6Qk4?Z-6tmIep~tO%j5+Skr#Zn#L{l%jkIxz<m77H!c-
zR2)1z>+gZ}^DRe0f@)U9+ux!tG<LSal1}E8eb~sqK4mLOO0c)Eh{K}?{B$Ntju}oa
zzY<?koFK1XPy6sLZ?0>um80A^Sb95Dqu-$Yvp$pEM*dXnw`&OC;z0J}fxl@weCV$i
zoFHb{!tgNJxG-s?OLJ|qg5<r?3usUTI~_5>C23nOn%GbGonh2OO-?-w1&i_PfBUbC
zlz3b{GBI>kX6pt@ckp2Da4`$YT=d_$==<`{SZJITD9ThFmT8^N6@j8UWb}Lk^oh;)
zLnu3*#DNGIfTIMXj%xFjOdfOeT_Vxw=AG|E_s-G>q#LaW5Vd?)Uow?7#UdbzSaJrR
zCwG3~>6T}&c$Po!l$osZai4eiBR`X%*>5oOYa!@c_`892BJ<-1o9ZSMWc2=*UX90+
zdzKA%{z|Ndp^ysm_Mu-o;}lLku5s`P%Ec&zrx=KwfU=mFU)qAM^EQMg%wrOR&czl9
z&xTWnxVvWh$C7@NNQ_#St`L_i*J=u*1&D8A@ifhxr>s4H(ijG?;b^mXEvH#}$&A=`
z)KQGqcE!5vvCLyZ+yGdLy;tB|_mb-njBAK0kD5j(-s`i{s!?557zm)?xY^;bJMWF^
z$N&m<5)KoFojPOw-OG;$t+@7+6IA_!f|pdUmdn;E(fBcNgmIjVHN-$1P0)fLQPCDp
z1JQfgouR^$GJV%IV-P$g?xJ?;T9AK~noZR=_g&aEX>F|j@)Gkd0E6gfwGEt|K6P&~
zmEIi7aI$mcDl2pDKIi={*tIQz<>6au1167O+1}A7B>(fmW8&I^&f*QC{L2?8#P1JR
zx^S7ZU{8#8Tc5=KQ3)brzXX8wa&{lmYogr~8km+TzlSb2sy#|WHn<*;-v@;FVBwzG
ztt>G}<7*0*lw@d-&e&CQwq;yrq*4sxx8A3o`&jU}$TIdaaHWsykZtBieprDY#d%pQ
z>`s>dKGRuVrv#eHmYm6`32HK_CX~iakB3qgSN}|ItIJ<eGI+Z4GW)bDcbTs~t+uF=
zDp70kIO^38i625kh=iZoW2#=oKm~mI#EA;_(A?q-?Ni={kXA}0;y1}vSlPDaFJh#(
z-&MbrCpYku^C=5<&PqJaUBR<0;+XQmCAk8D<yHJ0X9d`C=<Cfu(9Mb?(}j*8YL->1
zksAQ_CK>~9N#g+}-im!)>Y@P$qKCTRD(@(^^|VI}9(T7qm2PY}=H6K^yOg&-w0iJ_
zoEaEpA{Al~EuvO-?UZFNalbgz-p>4X)>hSdRS^huH$gc{g(9Tk)$fb{+y-#g`%pX_
zm`Xv4H@zkjPffIt^cuHZqc;!G#)VxCjRj_F!ko>BPr&Eq!k(LcEjzWLa={@01Ge*F
zcS_pouF!@cgd*$YWHd~_4+<^cW4e`+JXhZ|!Y{TOvQ932#n983dx@iu9kVNDD5m!V
zkgS+%N>abD+n*y!xkrHYIZ{p%4PRn3a;65}_J-=OxS0oS&BYPc%AEk~K-PF!e>7p>
zs}oR1Lea{UB{jEAIicZPCK}&G*Vg9M2E5=eZSUX+Pp#HNCo3)#q-M2Hbjv!hKzzBK
z48nrV3^$c2cI&oD4LWB4LZRn~8SZ!}at~-toKwH>qZ&<OND#3jx&%$u&mE;b`5T>p
zj2pw}O&o|T|DFvDW4jDIkn%mx|Mi6${Bzb0#Icu)-~L3AO6hriWC843t;)urPX4ka
z^+=cc+RD}|=#F1yr$s%Cl=!$7)U+wsGfhQ(7!;xey8rZ=z@jY*PB^6}CV~CHU#2a9
zF3<a<uS+vyTJ`;8hUJgO^BC1KGISj`xC@eAF4{iQ|5LPugnwLMl2eiEez2a-GXx*(
ztifO-MW&Yn60{etT9QN&r!bt$!_+78)l+AxVt4BTRsy&V0WtD=w`DIr*gh$?_~)iJ
zhyp0#Vy4Lr#i2(d`O`{KM1_mp6&s1Gw;mW@Ek)CvjH70X_hRN(&*<fCfQiK4YsFJ{
zckl~A>$w=HE?B<0j6p^Hg!TeS`NFk1@JN0b>0QuXz|b#7=8>{)4p=yF-)7xrV9vm9
zsPd&viL7{Nh&h__VNmT4+7t$kh^bb1(a-kN=Q}RStF~&vTaD>gYa4v2IBuugvin--
ztG0s^2jjJQ?%u2p$W&SLP}2b4@rz<VpHwPGzmxG}3u!>x<ylQW74n%v<N4MvG<0JS
zs?)<{+dmr@$LspK-&l=L@rr`p%uCGidSkK|HMbwEI1n5gf*vvYL;5ql@~1Sf;JkCu
zBXmhwNB8{s2O%x!v!k!FkC`W?LQh1(vq0T#842|giI{;5iaN;eCdTsQ#lTvp?Kj%}
z$aVQY_M+;Gjn^e4%bU_iF<~ENL_}Lu1=KWl58VSj-UPh}tQ<c}Pl=nX=Xb7%N3jJ=
z7U+NZS(OI<#fkGdc8e5)jOBKT@cUkS<>xxGT*un?ZHLR20e*FsU%9JFnuD`GWoEc5
zip_ODp}J`L@)@{a>qonM=MM7ST+OuA4{}ol4^MsCM@uwU_858?UtaMvHE2P&q(}Yp
zCEM@E4hx$<r1|7!LbZ^;otqa8JCC4*=mKO2%C5Eifj48cEcM#d4#|#7S0{?xA{7?P
zp&}?DwEhe4!Up~tTk>;#r1$dl%3(#IQO8P`#f(@PR6wnR4E${{Ixl%n;lLwYu_k3^
zrJBT(3+q0<I6T7rj%V(&b(eO%_8i@vsyfaI-SwRJ)i)4@iZ#vFFQOjDyrqQ4I2YCB
z&SY7=EKP}iC{xWPEbv{XZi3Va-N&88^SL#n^KW+2q*Sz!ytHIl-X7)2tSWbAQByd;
z6c)3wabMC`rAHgAjg;||uhw1l>^D^GzE43lfRpbMMvLl%0+UgLf>SNMSEli3S7Afm
zOv9#j+D{BD%=_Rl9jfo!OI3iU?z%WUiIIEyO8vQwasO9TRYhG_eyzV(y9?j`jXxsc
zz3$cKme32tUUxZkVl)Rdpd&qPYyT#NnY`Aqcs(QHiOI&)lty_uT6Z-d)uEkse7oR+
zAlsHCT#A)Sly1yhsv}=d$JRY9>;B(msbt`|wxWmVYvRnsDYb2gb4Z4^Kk8HZgKZR(
zM@J^*Yp9|FPJXk>D#;b*C&wqzY=L-eCpEu(1W8i){ry~5XOhxjc)%oUV`yVBu*Ldg
z7wnJ?a#{p1SsmU#L+-@8sh!Zv0ao3S)D@IROXs`oh=H)wYhbjGvci4on+Htt;ehK5
z?^oLxzkXDmnzvg;OA~qrwFs(a@{8wC2AHuv+CHcm`2aCfxZc-$C3DK2Pn#Omkh^s5
zYCJQ1(@$|#W^>z{9^fQzQ$U6RrW)n)36xtA6I-|ckg#ZN?CJT^j>m`a?qO~h$Co-X
z+Jv7%qVu+jtnB6onoU=Fh3CZ=>`7ltIAeG6(HLUEv%tkDi&ZIf7VbeQ0G7dzoVg4q
zTU#D<&PBQPKcT}ueAiv1eZ^M#^6&l%&{(_lTgv-N4yLAguUz76yUc8Soy2$1?)9Dn
zp+?J9x~5^&eE!msK$Z&UQ50V3nT;a+qmhag(<*9R(YD`y?anx_OL#U0lTiJw|LC?Z
z*VYS}n+-|$FAF@07<&hoGC;T;WZ!G5$P}mgo&0&*U6_1XeD19>Dr@z6a4hqBQ|DkX
z&HQL4xP#gxCAemE<|r8M7<?@yRnYRBhLrG3haUa5kzX1kpHiY5ziqx1=|Bkhgm0bK
z#$u$iD(U8TnMAF>X|yiV`wUhPC|7Z92jU}<qF%91T#x-qz9A{Z*y@~ctXi8E)&9{(
zM@kxW%xk7*uZ_1qE!LX!b%$7PW{BqAo;Y3Nb-d)sV9=shy!Hw0ge`O$=<9yXvB?$+
zeu#hXE9MTmqt+3j&{@lOB2{@Ksm|ZiO|k`lDC*7cFmEJXeIeJ(U6J-rKBxX`GcCR3
znN#d%;8Y(QzozyJNFRofk&knDy(ZUY?85d4U)~ScFWW?4SQ+<M$Dg{;ZV9K-UI#R_
zG#+VJuo}?w-??0S1|whx%s+&T(o~w$afOEuK}(967c>I}I=93-MH+j*w1e*9P9F2_
zo#}M|?0o?cbi^FO{JA&pU>i&BE!lz_SQMekn3$+JS1QYX4Q%CyZ=_nL_cIwbco#0a
zOolBbRpP(YA8Tr-1|Qxl{ySWF*Sn%cV^28LJ~|GbPhgUNO{9JRhvlFYpBeblb)b0=
zR)Tr?4<TLVt5>y=?xAAbMW<MbJg!HwKuc?CBzXekn{KL&KqLbMUatMx7BLe*$Sz4!
zk-X{Wzqiq8dLJ?Zj0${e9cqn{(u(vQ8<E#bRRK)YIQ5_)J?*mt;;{yyB|(30?lWz(
zZ-k*~CKdI)QB%K>SS>V$E|F*M_4iEOEvY{RV(!l*Zw9=?nRM7%D7BKE(8{XF71kuj
znnUkdYxleAHrw%S&RUzCMZ<O1=w*pH8c6y|^K6<Qcev3SmP!N08l7r;h6=A#kI)et
zaHfwM2isMyXHj<2yROk6=hk#L`|&!9Z?ure=lbsGl;qdo-_3xc&N6-rMeDNEMADuI
zx#eEVH8uRl*wY43@KW558qP?nLta=LAeQ`LuZe<+sCXd4n+>i3xh}WbEqEw2HRu>O
z!2^~uc}F(_O~s7E0hMj83I9tlnfP=H?lX8qKCS^7!LC#bsEASR9$7qPS&EtetnX6;
z1P@Z>9CA=nw0bjK-3ssnye9Z!2ZFg<?TQTMS`wI55p}}McRs|+NxPLRot>61E2guO
z2ukox1lQOc>7Mql<SF_y?r{O*2(|5Gpjlwi7+dN{ccc|;N@uH9OJ^Px<})Imlm+dx
zq}{UV86JSs|LxUSwN>TX8^8>;YfkaGMRQ03Gq5GxlKGW%v|2zJYd}h0p~oMQ5O7kw
z*~D|jdxCP88}}maCJEF(8Q5zn1wNV^xDVbb;xW~|YZtqNs0P@YT~NgC0N-x;rzn&D
z&dKK48j?9)*dtG4@sX;p105x`@uY+J0%Q`I+)$+n4l=BWW(Af$43$2GA=?`VnrRA>
zq`>-M(?M)-{?uchFH@#eJRu8ObCm4sY`UG%biZh&#Y%cxc<0Ns6lKa1>;iP74>#4r
zg@d3u77qW>zHZ+|=n9;-<Vk%S0+^=4GCZuPq9L!1oCJZd+v;_j@9~A)*hd}zZCP4M
z;IYDu%!WtsMOzhDmqwB^FRKC^n94uuq#5{f-hOwyyb9A6e5s%MT#U{CdYhG4H?NF~
z*%~mP%5@94;fI;=t6_&R!j%&48d~}&F~LM;gK2K!h|Y~J&y?;7bVQoK7Hw70{FKP_
z9qY>!O8j)B98OPC=xg8HFWRC|C!b{>KN_b%WBL!s9I1hi#<fa4t8GZ@1IeOis=$Kg
z%TCdp-sRDc7Uy$6_jSn7CzM6(Q{#_GG!57-b`4I4_<77Rxhx+ri~7E;|KX~aDBaU)
z#pFrI*sf2D<+rA<Q9(o_gviI^n@(FfV^!^6E!u?AI*0|nE#r_CZBUN$MF>%jmHeC4
z9^nW1l|o#_wksBl@QnIfpKpeLOeiMFcxq*N%JyL|r*~P&_Hrj`^p!jYWqgha`B{YW
zsvY&cue~z)DT%BOlVNR-#_B9zSYIQ&Og6xxn^usvC_Sddg=(Q47U^LTt550*8jPft
zIvG(9JC(W|j(Gd9m;PE#s);aiTNl<V)A#{^<W>vMx_z6TjKU)T*90HZ4|Q+Zf-U}Y
zg^l;r1|vI&8(bb-i5Oy{u7M=Z4Ivcc2X`!54NoESjwA=%8n*txD3LuM8dZovlFN0T
zP^uwW2iM4CAs<&h2*gK8dhk;1&^WdM>*l^j^&+f=mO+d(UPX_5u{8c}t#B1GuhsRj
zDiG-Fwe&95%v`xGNPHg#dzquvMJeZs1jR9>)GV-9e~Hy&@+^5b#wqPtuB1b4j>>jp
zaWc!N!x^Y4*VoZTyqhX?I-vr#*aYL<;+Xa_<o{OO?pz>2Uf#T37bG&Gc0uD0TP?x&
zV$1I#P@YGUv*3;&Vqty2@)7?qSm#fNa-h&J=6r7Au}6V>s!O2UEM*Ws)wWR%l#4Cl
z;VjL|5D$8B3^-e=QJ^<^c$xd2#0|opGOXv7&<2f-5@2g9UgMwDiG`#OUIvl6d&O(h
z`-!#C*a7=?gPxWHle4F&f}JT0_L4}6fX$}vF%7VZr>Dbw(ekR;Di$5eLeV_4xsAXv
z@W}1fm9q5!EAto0@_E)65m{%^SfCg<uf*DTa!jg?%f;^qgaZxIF*C1#?GcbP%kIx(
z^`XI{_nT#-Z_U&Hc$)s@`<~EPl`#%suX)tvEtYr|vK7$ZL^2V)kH;yVoVW%bFykFQ
KkJ_q}T>3vKR9<HQ

literal 0
HcmV?d00001

diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_upper.stl.gz b/tutorials/modules/multiphaseEuler/hydrofoil/constant/geometry/NACAMOD66_6deg_upper.stl.gz
new file mode 100644
index 0000000000000000000000000000000000000000..7d13f2274b8b0dbfff86a61cdedb0d50a25e3428
GIT binary patch
literal 4108
zcmX|^c_38l8^>ofm@FA(nI_c?%5ce&J?ogVWhab%A4Zm}C0nM!m55P8c4LH5vZP#P
z3oXVLBh#hDO=L?LBKe(hyT5<WIq&m+pXd8|p7%Y+Goa9#`UoG;K1)L*Lklacy1IwD
zuU~*iR9KkbrBe}+!5~J_{LN!6f~4Iqa%>*NEH=HZurxdw6yho!IQ%-x(&OOZDJ^MT
zE^m1KzINwlgJJ&CHmN_>Uw<h{__?&W=%N2}^-XA~e&UbM@Aa=HY#VPL!{Ogfb!||*
zkyrWMR@v;&?d|bfpBGo}DDLL(P9*LuzflApB#vK)dHpgpAopcK?bT{>lq!FgsffCM
zLi|B)|J84$zx&v2t=(eTO4-%>JG@OxxBh_C=wqcIx$1jL^`mWF(QFa!7?Y?7Bm{D7
zQ?ALs!dm0{p$gj}Kclp{p3?8Hj;PBgivAKr$3|#i;#9mVG?lhaAB@@{!``Id!a%p@
zN`IYU<%VZtlwzb*W%ik0yKaNN#bfm6TrcQ`i_=z=n#{C3H(}`<-trvb38Up=x0#uM
z9!*@TH_#jE;#iOyPQIqWSy|EU^Th$vFI>Y1?Awc1-Bn|6H&48`4M1szHVhP}X9vVU
zAHK(#LDi|bsJxt(`^2EZZbJyOfPT6PY-X|ggl4T=4d2YRBuO}yJeuFCohl^dLg{MG
z0hUJ9t97iyS3<4FR$!PbhkJXDG`vl`+ZaYp3P1U(U+YEYS4jgu{yv_yABxt~%yk!5
zVRFPTXCGkKWP`w;ps&nG$W^Qi*FOz@lnh<&FCa~y^8UhPGAN8|lpZ|tHcc-H5%w}C
zpg(L#Kj1t${}7C@-;c}M-eKQnlI3m_EA4)Y|8W_KwZX%=Q8NNWvx|`j&@<(PzRfgC
zBI#S}Ng|u{0STTz&z#7Gu@)ZN(I3t&N$F8@!O9C6Ij75X`FVw%2n_=Tr1Qf5i@CO;
z#gT3V;B|7CM)uPd?~q$r2g3%i`C7b`@aTE;=NqltX_>QJFTe!1)*Pm4$9qUp%*pwF
zEqqjSK3R@TGlXHZ4}U6Y#W8d%Dt_0e5+2Zf8=X~!>+DF4RH2S_FKa{bA)+@W>zTeq
zx`N7mm;d-|AGvIgqUM%usCa0b%$ypr`LFbh#*|wq$e5Ecbj~-49FxU;FaM8w#(R<r
zNh2eO)h|q?+C5|V0%M$mz;QBj(p&@!X8xOOS;kqx(Nuji9!|DKliUokYfmp~rgz|S
zqu8oz;fX~8Pre?*>y;l$4r2K+cy6YC+^bW5ZuYrtb>Sea!%<u4oDW*Lt9Y)<Q}Lf7
z?1#FS+7Srv#fzi99qNd>I37STs(IX+Z_$p=c)w-<;9E-4#b~MhSK3{SJI<|vMPl;~
zq02A3FG_<r21F={X;>0n4c6uPQAGl3<2#pTCq*Nix6vHI&Tocz3suD{VIa9>X^!Fu
z7mLs6MM@w4?`(M>&nFtdw3-?a6csN3W0{ZV&G>(m<_hH`JA|mq@0xmktEN$b7^(74
zeU+PclTh)$-s$psDdf)Uwp@sW1Y?hb$PJ1naucpBz5N~905TGzI^rL#6$I_G8eqfe
zLu}Q*XFo`<@@?fp=%)d?xjJxUOtZIrlp--(kA(PI(-M8Oz?dPk-)S%1sEtU_b6RWW
zW|H3SVf@$+A8=E9#=wt&OrTOcD#M7cO;eWDW{MaS86$qL`BLe)yH^h;H=mqe2uy2K
zm=}1OQC1);g4u%OOyV3a&~GhRbQwaXI8BGj(~RRA<(p-xaf&>yAL;8M5G}%4+8>zK
z{CuLmxgaHoW(`n*Pjn?eI`b{@$iK*-kk{9-8Y2n|oUlry)_3NRhL~FiB)|O&+xXZs
zVXSBhY~&?Kn0>a)+oK|WFDgASZ8WZbHN7_>;GR33qc`j6DKkW4CsS48tmmHbuvhYY
zS#p#m%TH2UFnln^rB#olq9Mq;KQ^FYppeK^q93V-OcLvW*amNz*6Nia{X5iuo50E_
z^!`lK<@U<R%F8HfhSz;4E+Q))<ayliI2K0dP!%M1`8K#--QnmtzB!*+DYVmxfvMe=
z@n)+%;aY;$Y(42<DoarIdZfySV?!s?d=uf8zgx4?))^w=u<c1qJCy_}TX@$%3ohQ$
zhjyId*z0Q+eqWxdB;N8~io=c1_2n9N)$N5Ug~C0a;}G7xJ+KI-J%|eUd+_^Gz2{B#
zxxG#x=SO6J*QO<*{l>!8F#tJ_+$>DUi`%i5T=fIKMOwVV_zOkpmk@P#xsrwEe3CAV
z=Oj^5d|W=7{<hGwFl$B=jwbZ}vNbv0U9$d5!n~n{k5@Prm^a0w?C#b)zGVN~SRi4%
z0$p}^rFDRPYjw~6MjFc6jKX{hQK-EF$DBV7k?g$U@I`>FBjNYqDNruSdLfr+#|L6V
zA|NW|RGR7?e+9SXp!*>uzOXjVYn(!3?l=l(Et*<M(M!|&#JA0c^j*49A<(7)FpcdQ
z4O6yv@UjJVenH^vqBK{tcOLyldS59IA#hAH&T%zJ6w?Yddsj8v^LI`B;SN(S0L;dK
zG9iC<oM`P=O(9vDQ27jwU6k&bC{={@tV>sk5Hnj|tYcipmIb~7VNArqqy}+@pPA=<
zpF-)q;URWO)@eKwd*WV!I*KNA4KX$uh=Q#C>iSPuyPA8mX+i4i_eM}<H&{H(AzI_v
z+V0CKFWnc&SinszR-U2`sR{D2J%1c~vnj^?_OCx5wv5Oej^&(UY1XSt%j;ioZ^p8{
zw~$EM+DO?>aL#|R6Wy4OWtVg&5PYb6S>=_Rt-8YrH>aF=n0L{qKYJk!AWr6GeLlKN
zH{y8oX#B~@%C;zm=JGe;t`tZcaL!p69XF#j+XP2$o8)n$X#C}C{?YMAQ8!OV_FsOi
zqLYSn%8g>`P)Fa_HNg@0ig`HIY08<d-OjoApi&54yhj`1*mrbv-6iEr(X;{`sB%0^
zVbRA$<?#4IyumkVP{4WH{L@MC_D^=IYF>km4EwZuTpQOBNZ3g}acB3-!fLN?k=njW
zrh{wQ!)G$JD#+)NhHD!U6L*y&Cf#0ydd0#F>rc+pFVDtj+-7czcv(WHFebnr+#Y{{
zP8+D}3p$cy_Sh>(E$9W~Gf@2}7#lEvF07I3GM%!$-J?GF4cwgz<)q6nI`UrZZ@ZJF
zL2zpJDH2Q3_nm0*+j1>IaU$=22&dt%PkFf4mR8?y><C3s0m>v^6B<+3BL{_G#nPTd
z=_-VSu<wKeE%zCxHlSAl77Y61QVm@mY=%n4!Z?Q}-L$k#Nan=7vt#lmKw}CXS|o7w
znwD%pthnlXvH8l}%FhvE2Ax!M@bTe&Ll~BaCE|=ue^#gH2d&*}$^gcl6r&PIzJY^5
z4lwC#lN**_C1kLe=bstJ=Iq6HhM!r#(C$tDXzynT9Aru1u%Dwn7b<?(JVxr6a)tv_
zo{GIZmSUG)Gk!~C<YoQFP7f<MqIh{nKWNv-<=Ts@E6Q~W-Kv95Hf~o>7I$mY9#<ty
zQ-YK|emr34&lI(sn5T>DkGH*6*fGV9^@w9y;bDaq&!VOqob193W3vsyfFzyKeu@D`
za`tM@*?`I`iaSv+LK9_N_|G{KW%dSV0k_G;^%`M2hvpt&O6(im9ElQQYWIkK1o|zX
z-soS(!xz`6UOn@Mnx+R4xX56Wyt-neKb`@WLjp7~5ePJ0s9GSR2Jr~~Ce_d<Swh|m
zth5xHV6`Gm{AckTo^h7om2jl;ovERNfvAtR@(2a#(zUWEIX^m(_sIs@`7U#)fPbOx
zo%V>7n%j!Wjk)JZQ&R(F*u<E%`M?mFQzBUy=u`-i>AD7%QTP+?5bvl=%`MMs5VIQ;
zEbLvcSw<f7gV86UnfWHGj|*C&AUK;S&mM}UK+v`w*ejT!9P;pr)e4&U0njqQ-H8xt
z4wo;L6SqAt-n+^E?&>ujCCvP5#szuIHybzz2$J5fo8Vt@cfxp4q%~!7o1X6BEwz{U
zx<xnPV(AY>z!ze+WfNn|gFzR@p7a1C8^<)uRNDH%`03vIe_Q`)vU+F<i_E$P8%|BC
zw+!8Ht!uhj@WB)`MWHv&ltNflY|=`nyTw5wfOdidb3zd8UDtyPc`B)6g46`&j6AaK
z=GcQ8_geT!-%#ksA6aE<q!FFTu(SW61)|Ky_Ku&x4b!_}MqBV#ESmbW2e&I$T&^tt
zUnjfe3L(&8LmCQ>hbV2FWmtu1A0s~asUKL<eBu@TPW$zzX=OCC*3_h|SuT8_1(;a4
zzyFN*eUQ-B7tZy84^cOqC|HUk%K(M&bbBu<P%TtWKf(N?GZ#cxw(+;<L<?GhI_NAd
z$+MAQSRF$bwQh6aq1c57zatR%&=0IuF-K+Kcoifz5=DLN7_bhO_m7*KC0yeLLx3wd
zCFSov+4*3n7Yel`Jb<53ms?l5sH}7Df0_0><C~&`jQVRn6#=ENX5Ug?eN!-_)$&JI
zE~E{(A2BelRf%?GF#BaG*>a(>N0SZ`9v(am(V3^9RtKR205_^s`_~Ioc{T5hPxM84
zBNqozHx8J;?p-#0))!00sKivMR;A!YQ_S5DO>ZE*tz}F@1G@{>5G5DxZf~@wPumE*
z{2L(jW{7WsOHiSBZOu$Nsz~ICu}9aa$B`N1X#&Yum~;OqAlZ6aApC<Z9C|p7CimH4
zT1ZRlE1Wa)AOsBzK@v%8ekMImEYQMh(S$09j+al9e|=IfR{T-?$rSn$XoJEGp4Cc2
zt%|StFe7bQ9P5k@Ce*U?s$@6!CmWtl7VyFT=MEAto<Q8>s1(Z^T6guY!^_t-XE%Pg
z2+cK8`9^ie%IZT}j;i?w3>N3XCqP6N{-t!l#f8A?);#E-N&5M8#B7+CQSs2Aqdc&w
zCoLrH(Va}5AhLw<(PApBZ4zAvll>$ii@afH;2?6<M8C&fot4+DX(_Fj(RueTn*v#u
z_OBz`cQMWU=6rKi(7&i>pN3l8<ef(Igf5*6tb9?)u=A?#aE+#yc^b9K2Bqx0!avl0
zwh|#?q!d%pWJ`umw^nkAG9HP@QcEPe-7VZ@HT49ph^1)6RH)j)7<DC-?NX~4^A*&j
zh`~h<3a&H3V9ADW>aKi~3W%(A_*w<7*2hz-RXFJFB3_|$<>-0^tgYoprKq<=qw1NN
z>h}oFFJy^i-p<U@1HM6{a1l=ACflYVCR(;<f*go`o9+0~7QC3{e48ly`AntB(C+v7
TJA?W=o0H2p@%oOMHW284b!(p1

literal 0
HcmV?d00001

diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.gas b/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.gas
new file mode 100644
index 0000000000..a200e5c3b1
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.gas
@@ -0,0 +1,28 @@
+/*--------------------------------*- 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       dictionary;
+    object      momentumTransport.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+
+    printCoeffs     on;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.liquid
new file mode 100644
index 0000000000..62775f7b79
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/momentumTransport.liquid
@@ -0,0 +1,28 @@
+/*--------------------------------*- 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       dictionary;
+    object      momentumTransport.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RAS;
+
+RAS
+{
+    model    kOmegaSST;
+
+    turbulence      on;
+
+    printCoeffs     on;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/phaseProperties b/tutorials/modules/multiphaseEuler/hydrofoil/constant/phaseProperties
new file mode 100644
index 0000000000..fd474a285d
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/phaseProperties
@@ -0,0 +1,118 @@
+/*--------------------------------*- 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      phaseProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+type    basicMultiphaseSystem;
+
+phases (gas liquid);
+
+gas
+{
+    type            purePhaseModel;
+
+    diameterModel   constant;
+    constantCoeffs
+    {
+        d               0.0002;
+    }
+
+    residualAlpha   1e-6;
+}
+
+liquid
+{
+    type            purePhaseModel;
+
+    diameterModel   none;
+
+    residualAlpha   1e-6;
+}
+
+blending
+{
+    default
+    {
+        type        continuous;
+        phase       liquid;
+    }
+}
+
+drag
+{
+    gas_dispersedIn_liquid
+    {
+        type        IshiiZuber;
+    }
+}
+
+virtualMass
+{}
+
+heatTransfer
+{
+    gas_dispersedIn_liquid
+    {
+        type        RanzMarshall;
+    }
+}
+
+lift
+{}
+
+wallLubrication
+{}
+
+turbulentDispersion
+{}
+
+aspectRatio
+{}
+
+surfaceTension
+{
+    gas_liquid
+    {
+        type        constant;
+        sigma       0.071;
+    }
+}
+
+saturation
+{}
+
+phaseTransfer
+{
+    gas_dispersedIn_liquid
+    {
+        type        cavitation;
+
+        model       Kunz;
+
+        liquid      liquid;
+
+        pSat        79995.75943;
+
+        UInf        5.33;
+        tInf        0.028142589;
+        Cc          100;
+        Cv          100;
+    }
+}
+
+interfaceCompression
+{}
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.gas b/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.gas
new file mode 100644
index 0000000000..ab6f58b8e4
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.gas
@@ -0,0 +1,46 @@
+/*--------------------------------*- 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       dictionary;
+    object      thermophysicalProperties.gas;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+thermoType
+{
+    type            heRhoThermo;
+    mixture         pureMixture;
+    transport       const;
+    thermo          hConst;
+    equationOfState perfectGas;
+    specie          specie;
+    energy          sensibleEnthalpy;
+}
+
+mixture
+{
+    specie
+    {
+        molWeight   18.015;
+    }
+    thermodynamics
+    {
+        Cp          1006.3;
+        Hf          6e3;
+    }
+    transport
+    {
+        mu          9.86e-06;
+        Pr          0.71;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.liquid b/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.liquid
new file mode 100644
index 0000000000..f01a6ad7d4
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/constant/physicalProperties.liquid
@@ -0,0 +1,50 @@
+/*--------------------------------*- 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       dictionary;
+    object      thermophysicalProperties.liquid;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+thermoType
+{
+    type            heRhoThermo;
+    mixture         pureMixture;
+    transport       const;
+    thermo          eConst;
+    equationOfState rhoConst;
+    specie          specie;
+    energy          sensibleInternalEnergy;
+}
+
+mixture
+{
+    specie
+    {
+        molWeight   18.015;
+    }
+    equationOfState
+    {
+        rho         998.4;
+    }
+    thermodynamics
+    {
+        Cv          4180;
+        Hf          40e3;
+    }
+    transport
+    {
+        mu          9.9768e-04;
+        Pr          2.289;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/blockMeshDict b/tutorials/modules/multiphaseEuler/hydrofoil/system/blockMeshDict
new file mode 100644
index 0000000000..75d26cfdaa
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/blockMeshDict
@@ -0,0 +1,257 @@
+/*--------------------------------*- 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       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 1.0;
+
+x0 -0.381;
+x1 -0.081;
+x2 -0.02;
+xTip 0.0015;
+xTail 0.149571;
+x3 0.21;
+x4 0.96;
+
+y0 -0.096;
+y1 -0.03;
+y1Tip -0.022;
+y1Tail -0.038;
+yTipL 0.00625802;
+yTipU 0.0095993;
+yTail -0.00801537;
+y2 0;
+y3 0.03;
+y3Tip 0.038;
+y3Tail 0.022;
+y4 0.096;
+
+z 0.01;
+
+nX0 24;
+nX1 20;
+nX2 80;
+nX3 20;
+nX4 60;
+
+nY0 11;
+nY1 11;
+nY2 11;
+
+nO 16;
+
+geometry
+{
+    hydrofoilTip
+    {
+        type triSurfaceMesh;
+        file "NACAMOD66_6deg_tip.stl";
+    }
+    hydrofoilLower
+    {
+        type triSurfaceMesh;
+        file "NACAMOD66_6deg_lower.stl";
+    }
+    hydrofoilUpper
+    {
+        type triSurfaceMesh;
+        file "NACAMOD66_6deg_upper.stl";
+    }
+}
+
+vertices
+(
+    ($x0    $y0     #neg $z)
+    ($x0    $y1     #neg $z)
+    ($x0    $y3     #neg $z)
+    ($x0    $y4     #neg $z)
+
+    ($x1    $y0     #neg $z)
+    ($x1    $y1     #neg $z)
+    ($x1    $y3     #neg $z)
+    ($x1    $y4     #neg $z)
+
+    ($x2    $y0     #neg $z)
+    ($x2    $y1Tip  #neg $z)
+    ($x2    $y3Tip  #neg $z)
+    ($x2    $y4     #neg $z)
+
+    ($xTip  $yTipL  #neg $z)
+    ($xTip  $yTipU  #neg $z)
+
+    ($xTail $y0     #neg $z)
+    ($xTail $y1Tail #neg $z)
+    ($xTail $yTail  #neg $z)
+    ($xTail $y3Tail #neg $z)
+    ($xTail $y4     #neg $z)
+
+    ($x3    $y0     #neg $z)
+    ($x3    $y1     #neg $z)
+    ($x3    $y2     #neg $z)
+    ($x3    $y3     #neg $z)
+    ($x3    $y4     #neg $z)
+
+    ($x4    $y0     #neg $z)
+    ($x4    $y1     #neg $z)
+    ($x4    $y2     #neg $z)
+    ($x4    $y3     #neg $z)
+    ($x4    $y4     #neg $z)
+
+
+    ($x0    $y0     $z)
+    ($x0    $y1     $z)
+    ($x0    $y3     $z)
+    ($x0    $y4     $z)
+
+    ($x1    $y0     $z)
+    ($x1    $y1     $z)
+    ($x1    $y3     $z)
+    ($x1    $y4     $z)
+
+    ($x2    $y0     $z)
+    ($x2    $y1Tip  $z)
+    ($x2    $y3Tip  $z)
+    ($x2    $y4     $z)
+
+    ($xTip  $yTipL  $z)
+    ($xTip  $yTipU  $z)
+
+    ($xTail $y0     $z)
+    ($xTail $y1Tail $z)
+    ($xTail $yTail  $z)
+    ($xTail $y3Tail $z)
+    ($xTail $y4     $z)
+
+    ($x3    $y0     $z)
+    ($x3    $y1     $z)
+    ($x3    $y2     $z)
+    ($x3    $y3     $z)
+    ($x3    $y4     $z)
+
+    ($x4    $y0     $z)
+    ($x4    $y1     $z)
+    ($x4    $y2     $z)
+    ($x4    $y3     $z)
+    ($x4    $y4     $z)
+);
+
+blocks
+(
+    hex ( 0  4  5  1 29 33 34 30) ($nX0 $nY0 1) simpleGrading (0.1 1 1)
+    hex ( 1  5  6  2 30 34 35 31) ($nX0 $nY1 1) simpleGrading (0.1 1 1)
+    hex ( 2  6  7  3 31 35 36 32) ($nX0 $nY2 1) simpleGrading (0.1 1 1)
+
+    hex ( 4  8  9  5 33 37 38 34) ($nX1 $nY0 1) simpleGrading (1 1 1)
+    hex ( 5  9 10  6 34 38 39 35) ($nX1 $nY1 1) simpleGrading (1 1 1)
+    hex ( 6 10 11  7 35 39 40 36) ($nX1 $nY2 1) simpleGrading (1 1 1)
+
+    hex ( 8 14 15  9 37 43 44 38) ($nX2 $nY0 1) simpleGrading (1 1 1)
+    hex ( 9 15 16 12 38 44 45 41) ($nX2 $nO 1) edgeGrading (1 2 2 1 0.2 0.2 0.2 0.2 1 1 1 1)
+    hex ( 9 12 13 10 38 41 42 39) ($nO $nY1 1) simpleGrading (0.2 1 1)
+    hex (13 16 17 10 42 45 46 39) ($nX2 $nO 1) edgeGrading (2 1 1 2 5 5 5 5 1 1 1 1)
+    hex (10 17 18 11 39 46 47 40) ($nX2 $nY2 1) simpleGrading (1 1 1)
+
+    hex (14 19 20 15 43 48 49 44) ($nX3 $nY0 1) simpleGrading (1 1 1)
+    hex (15 20 21 16 44 49 50 45) ($nX3 $nO 1) edgeGrading (1 1 1 1 0.2 1 1 0.2 1 1 1 1)
+    hex (16 21 22 17 45 50 51 46) ($nX3 $nO 1) edgeGrading (1 1 1 1 5 1 1 5 1 1 1 1)
+    hex (17 22 23 18 46 51 52 47) ($nX3 $nY2 1) simpleGrading (1 1 1)
+
+    hex (19 24 25 20 48 53 54 49) ($nX4 $nY0 1) simpleGrading (10 1 1)
+    hex (20 25 26 21 49 54 55 50) ($nX4 $nO 1) simpleGrading (10 1 1)
+    hex (21 26 27 22 50 55 56 51) ($nX4 $nO 1) simpleGrading (10 1 1)
+    hex (22 27 28 23 51 56 57 52) ($nX4 $nY2 1) simpleGrading (10 1 1)
+);
+
+edges
+(
+    project 12 13 (hydrofoilTip)
+    project 12 16 (hydrofoilLower)
+    project 13 16 (hydrofoilUpper)
+
+    project 41 42 (hydrofoilTip)
+    project 41 45 (hydrofoilLower)
+    project 42 45 (hydrofoilUpper)
+
+    arc 9 10 30 (0 0 -1)
+    arc 9 15 20 (0 0 1)
+    arc 10 17 30 (0 0 -1)
+
+    arc 38 39 30 (0 0 -1)
+    arc 38 44 20 (0 0 1)
+    arc 39 46 30 (0 0 -1)
+);
+
+defaultPatch
+{
+    name frontAndBack;
+    type empty;
+}
+
+boundary
+(
+    inlet
+    {
+        type patch;
+        faces
+        (
+            (0 1 30 29)
+            (1 2 31 30)
+            (2 3 32 31)
+        );
+    }
+    outlet
+    {
+        type patch;
+        faces
+        (
+            (24 25 54 53)
+            (25 26 55 54)
+            (26 27 56 55)
+            (27 28 57 56)
+        );
+    }
+    hydrofoil
+    {
+        type wall;
+        faces
+        (
+            (12 13 42 41)
+            (12 16 45 41)
+            (13 16 45 42)
+        );
+    }
+    walls
+    {
+        type wall;
+        faces
+        (
+            ( 0  4 33 29)
+            ( 4  8 37 33)
+            ( 8 14 43 37)
+            (14 19 48 43)
+            (19 24 53 48)
+
+            ( 3  7 36 32)
+            ( 7 11 40 36)
+            (11 18 47 40)
+            (18 23 52 47)
+            (23 28 57 52)
+        );
+    }
+);
+
+mergePatchPairs
+();
+
+// ************************************************************************* //
+
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/controlDict b/tutorials/modules/multiphaseEuler/hydrofoil/system/controlDict
new file mode 100644
index 0000000000..e012a01ca6
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/controlDict
@@ -0,0 +1,63 @@
+/*--------------------------------*- 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       dictionary;
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     foamRun;
+
+solver          multiphaseEuler;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         0.2;
+
+deltaT          2e-6;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.01;
+
+purgeWrite      0;
+
+writeFormat     binary;
+
+writePrecision  9;
+
+writeCompression uncompressed;
+
+timePrecision   14;
+
+timeFormat      general;
+
+runTimeModifiable yes;
+
+adjustTimeStep  yes;
+
+maxCo           0.8;
+
+maxDeltaT       1e-4;
+
+functions
+{
+    #includeFunc yPlus(phase=liquid)
+    #includeFunc cellMin(funcName=min, alpha.liquid, alpha.gas, p, p_rgh, U.liquid, U.gas, T.gas, T.liquid)
+    #includeFunc cellMax(funcName=max, alpha.liquid, alpha.gas, p, p_rgh, U.liquid, U.gas, T.gas, T.liquid)
+    #includeFunc hydrofoilPressure
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/decomposeParDict b/tutorials/modules/multiphaseEuler/hydrofoil/system/decomposeParDict
new file mode 100644
index 0000000000..7d3802d855
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/decomposeParDict
@@ -0,0 +1,21 @@
+/*--------------------------------*- 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       dictionary;
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 4;
+
+method          scotch;
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/fvConstraints b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvConstraints
new file mode 100644
index 0000000000..cd06a71b91
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvConstraints
@@ -0,0 +1,24 @@
+/*--------------------------------*- 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       dictionary;
+    object      fvConstraints;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+limitp
+{
+    type       limitPressure;
+    min        1;
+    max        1e12;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSchemes b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSchemes
new file mode 100644
index 0000000000..867f0ab5a7
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSchemes
@@ -0,0 +1,68 @@
+/*--------------------------------*- 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       dictionary;
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+    limited         cellLimited Gauss linear 1;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+divSchemes
+{
+    default         none;
+
+    div(phi,alpha)  Gauss vanLeer;
+
+    div(alphaRhoPhi,U) Gauss linearUpwind limited;
+    div(alphaRhoPhi,h) Gauss linearUpwind limited;
+    div(alphaRhoPhi,e) Gauss linearUpwind limited;
+    div(alphaRhoPhi,f) Gauss upwind;
+    div(alphaRhoPhi,K) Gauss linearUpwind limited;
+    div(alphaRhoPhi,k) Gauss upwind;
+    div(alphaRhoPhi,omega) Gauss upwind;
+
+    div(alphaRhoPhi,(p|rho)) Gauss linearUpwind limited;
+
+    div((((alpha*rho)*nuEff)*dev2(T(grad(U))))) Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear corrected;
+}
+
+snGradSchemes
+{
+    default        corrected;
+}
+
+wallDist
+{
+    method          meshWave;
+    correctWalls    true;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSolution b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSolution
new file mode 100644
index 0000000000..46fe1f368b
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/fvSolution
@@ -0,0 +1,100 @@
+/*--------------------------------*- 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;
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "alpha.*"
+    {
+        nAlphaCorr      1;
+        nAlphaSubCycles 1;
+    };
+
+    p_rgh
+    {
+        solver          GAMG;
+        smoother        DICGaussSeidel;
+        tolerance       1e-8;
+        relTol          0.1;
+    };
+
+    p_rghFinal
+    {
+        $p_rgh;
+        relTol          0;
+    }
+
+    "U.*"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-6;
+        relTol          0;
+    }
+
+    "(e|h).*"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-8;
+        relTol          0.001;
+    }
+
+    "(k|omega).*"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-8;
+        relTol          0;
+    }
+
+    Phi
+    {
+        $p_rgh;
+        relTol          0;
+    };
+}
+
+potentialFlow
+{
+    nNonOrthogonalCorrectors 10;
+}
+
+PIMPLE
+{
+    nOuterCorrectors    3;
+    nCorrectors         3;
+    nEnergyCorrectors   1;
+
+    transportCorrectionFinal yes;
+
+    faceMomentum        no;
+}
+
+cache
+{
+    grad(U);
+}
+
+relaxationFactors
+{
+    equations
+    {
+        ".*"            1;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/system/hydrofoilPressure b/tutorials/modules/multiphaseEuler/hydrofoil/system/hydrofoilPressure
new file mode 100644
index 0000000000..2821f57fdf
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/system/hydrofoilPressure
@@ -0,0 +1,149 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            sets;
+libs            ("libsampling.so");
+
+writeControl    writeTime;
+
+interpolationScheme cellPoint;
+
+setFormat       raw;
+
+fields          (p);
+
+axis            xyz;
+
+sets
+(
+    lower
+    {
+        type            boundaryPoints;
+        axis            xyz;
+        points
+        (
+            (0      0      0)
+            (0.003 -0.0003 0)
+            (0.006 -0.0006 0)
+            (0.009 -0.0009 0)
+            (0.012 -0.0012 0)
+            (0.015 -0.0015 0)
+            (0.018 -0.0018 0)
+            (0.021 -0.0021 0)
+            (0.024 -0.0024 0)
+            (0.027 -0.0027 0)
+            (0.03  -0.003  0)
+            (0.033 -0.0033 0)
+            (0.036 -0.0036 0)
+            (0.039 -0.0039 0)
+            (0.042 -0.0042 0)
+            (0.045 -0.0045 0)
+            (0.048 -0.0048 0)
+            (0.051 -0.0051 0)
+            (0.054 -0.0054 0)
+            (0.057 -0.0057 0)
+            (0.06  -0.006  0)
+            (0.063 -0.0063 0)
+            (0.066 -0.0066 0)
+            (0.069 -0.0069 0)
+            (0.072 -0.0072 0)
+            (0.075 -0.0075 0)
+            (0.078 -0.0078 0)
+            (0.081 -0.0081 0)
+            (0.084 -0.0084 0)
+            (0.087 -0.0087 0)
+            (0.09  -0.009  0)
+            (0.093 -0.0093 0)
+            (0.096 -0.0096 0)
+            (0.099 -0.0099 0)
+            (0.102 -0.0102 0)
+            (0.105 -0.0105 0)
+            (0.108 -0.0108 0)
+            (0.111 -0.0111 0)
+            (0.114 -0.0114 0)
+            (0.117 -0.0117 0)
+            (0.12  -0.012  0)
+            (0.123 -0.0123 0)
+            (0.126 -0.0126 0)
+            (0.129 -0.0129 0)
+            (0.132 -0.0132 0)
+            (0.135 -0.0135 0)
+            (0.138 -0.0138 0)
+            (0.141 -0.0141 0)
+            (0.144 -0.0144 0)
+            (0.147 -0.0147 0)
+            (0.15  -0.015  0)
+        );
+        patches         (hydrofoil);
+        maxDistance     1;
+    }
+    upper
+    {
+        type            boundaryPoints;
+        axis            xyz;
+        points
+        (
+            (0     0.015  0)
+            (0.003 0.0147 0)
+            (0.006 0.0144 0)
+            (0.009 0.0141 0)
+            (0.012 0.0138 0)
+            (0.015 0.0135 0)
+            (0.018 0.0132 0)
+            (0.021 0.0129 0)
+            (0.024 0.0126 0)
+            (0.027 0.0123 0)
+            (0.03  0.012  0)
+            (0.033 0.0117 0)
+            (0.036 0.0114 0)
+            (0.039 0.0111 0)
+            (0.042 0.0108 0)
+            (0.045 0.0105 0)
+            (0.048 0.0102 0)
+            (0.051 0.0099 0)
+            (0.054 0.0096 0)
+            (0.057 0.0093 0)
+            (0.06  0.009  0)
+            (0.063 0.0087 0)
+            (0.066 0.0084 0)
+            (0.069 0.0081 0)
+            (0.072 0.0078 0)
+            (0.075 0.0075 0)
+            (0.078 0.0072 0)
+            (0.081 0.0069 0)
+            (0.084 0.0066 0)
+            (0.087 0.0063 0)
+            (0.09  0.006  0)
+            (0.093 0.0057 0)
+            (0.096 0.0054 0)
+            (0.099 0.0051 0)
+            (0.102 0.0048 0)
+            (0.105 0.0045 0)
+            (0.108 0.0042 0)
+            (0.111 0.0039 0)
+            (0.114 0.0036 0)
+            (0.117 0.0033 0)
+            (0.12  0.003  0)
+            (0.123 0.0027 0)
+            (0.126 0.0024 0)
+            (0.129 0.0021 0)
+            (0.132 0.0018 0)
+            (0.135 0.0015 0)
+            (0.138 0.0012 0)
+            (0.141 0.0009 0)
+            (0.144 0.0006 0)
+            (0.147 0.0003 0)
+            (0.15  0      0)
+        );
+        patches         (hydrofoil);
+        maxDistance     1;
+    }
+);
+
+
+// ************************************************************************* //
diff --git a/tutorials/modules/multiphaseEuler/hydrofoil/validation/createGraph b/tutorials/modules/multiphaseEuler/hydrofoil/validation/createGraph
new file mode 100755
index 0000000000..172a705585
--- /dev/null
+++ b/tutorials/modules/multiphaseEuler/hydrofoil/validation/createGraph
@@ -0,0 +1,48 @@
+#!/bin/sh
+
+if ! which gnuplot > /dev/null 2>&1
+then
+    echo 'gnuplot not found - skipping graph creation' >&2
+    exit 1
+fi
+
+lowerFile=postProcessing/hydrofoilPressure/$(foamListTimes -latestTime)/lower.xy
+upperFile=postProcessing/hydrofoilPressure/$(foamListTimes -latestTime)/upper.xy
+
+gnuplot<<EOF
+
+\$experiment << end
+0.05 -1.416
+0.10 -1.400
+0.20 -1.375
+0.30 -1.148
+0.40 -0.744
+0.50 -0.695
+0.60 -0.582
+0.70 -0.520
+0.80 -0.341
+0.90 -0.108
+0.50  0.040
+end
+
+set terminal postscript eps size 8,9 color enhanced font "Helvetica,20"
+set output "./validation/$(basename "$PWD").eps"
+
+set grid
+
+c = 0.15
+UInf = 5.33
+pInf = 1e5
+rhoInf = 998.4
+
+set xlabel 'x/L'
+set ylabel 'Cp'
+plot [0:1] \
+    '$lowerFile' u (\$1/0.15):(2*(\$4-pInf)/(rhoInf*UInf*UInf)) w l lt 1 t 'Simulation Lower', \
+    '$upperFile' u (\$1/0.15):(2*(\$4-pInf)/(rhoInf*UInf*UInf)) w l lt 2 t 'Simulation Upper', \
+    '\$experiment' u 1:2 w p lt rgb "black" pt 4 t 'Experiment'
+
+EOF
+
+#------------------------------------------------------------------------------
+