From 592a0c2ad1048510ad77863c2ed831b80b581e64 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 27 May 2011 12:28:59 +0100
Subject: [PATCH] ENH: Added aachenBomb tutorial case for sprayFoam

---
 .../lagrangian/sprayFoam/aachenBomb/0/N2      |  29 +++
 .../lagrangian/sprayFoam/aachenBomb/0/O2      |  29 +++
 tutorials/lagrangian/sprayFoam/aachenBomb/0/T |  29 +++
 tutorials/lagrangian/sprayFoam/aachenBomb/0/U |  30 +++
 .../sprayFoam/aachenBomb/0/Ydefault           |  29 +++
 .../lagrangian/sprayFoam/aachenBomb/0/alphat  |  32 +++
 .../lagrangian/sprayFoam/aachenBomb/0/epsilon |  32 +++
 tutorials/lagrangian/sprayFoam/aachenBomb/0/k |  32 +++
 .../lagrangian/sprayFoam/aachenBomb/0/mut     |  32 +++
 tutorials/lagrangian/sprayFoam/aachenBomb/0/p |  29 +++
 .../sprayFoam/aachenBomb/chemkin/chem.inp     |  11 +
 .../aachenBomb/constant/RASProperties         |  25 ++
 .../aachenBomb/constant/chemistryProperties   |  37 +++
 .../sprayFoam/aachenBomb/constant/g           |  22 ++
 .../constant/polyMesh/blockMeshDict           |  57 +++++
 .../aachenBomb/constant/polyMesh/boundary     |  28 +++
 .../aachenBomb/constant/radiationProperties   |  67 ++++++
 .../aachenBomb/constant/sprayCloudProperties  | 215 ++++++++++++++++++
 .../constant/thermophysicalProperties         |  41 ++++
 .../aachenBomb/constant/turbulenceProperties  |  21 ++
 .../sprayFoam/aachenBomb/system/controlDict   |  53 +++++
 .../sprayFoam/aachenBomb/system/fvSchemes     |  72 ++++++
 .../sprayFoam/aachenBomb/system/fvSolution    |  94 ++++++++
 23 files changed, 1046 insertions(+)
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/N2
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/O2
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/T
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/U
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/Ydefault
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/alphat
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/epsilon
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/k
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/mut
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/0/p
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/chemkin/chem.inp
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/RASProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/chemistryProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/g
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/blockMeshDict
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/boundary
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/radiationProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/sprayCloudProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/thermophysicalProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/constant/turbulenceProperties
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/system/controlDict
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSchemes
 create mode 100644 tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSolution

diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/N2 b/tutorials/lagrangian/sprayFoam/aachenBomb/0/N2
new file mode 100644
index 00000000000..291efd9c59f
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/N2
@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      N2;
+}
+// ************************************************************************* //
+
+dimensions          [ 0 0 0 0 0 0 0 ];
+
+internalField       uniform 0.766;
+
+boundaryField
+{
+    walls
+    {
+        type                zeroGradient;
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/O2 b/tutorials/lagrangian/sprayFoam/aachenBomb/0/O2
new file mode 100644
index 00000000000..bacf89cb28e
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/O2
@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      O2;
+}
+// ************************************************************************* //
+
+dimensions          [ 0 0 0 0 0 0 0 ];
+
+internalField       uniform 0.234;
+
+boundaryField
+{
+    walls
+    {
+        type                zeroGradient;
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/T b/tutorials/lagrangian/sprayFoam/aachenBomb/0/T
new file mode 100644
index 00000000000..b6ffeed49c6
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/T
@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      T;
+}
+// ************************************************************************* //
+
+dimensions          [ 0 0 0 1 0 0 0 ];
+
+internalField       uniform 800;
+
+boundaryField
+{
+    walls
+    {
+        type                zeroGradient;
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/U b/tutorials/lagrangian/sprayFoam/aachenBomb/0/U
new file mode 100644
index 00000000000..98260b8fb4b
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/U
@@ -0,0 +1,30 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// ************************************************************************* //
+
+dimensions          [ 0 1 -1 0 0 0 0 ];
+
+internalField       uniform ( 0 0 0 );
+
+boundaryField
+{
+    walls
+    {
+        type                fixedValue;
+        value               uniform ( 0 0 0 );
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/Ydefault b/tutorials/lagrangian/sprayFoam/aachenBomb/0/Ydefault
new file mode 100644
index 00000000000..1039e373fe5
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/Ydefault
@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      Ydefault;
+}
+// ************************************************************************* //
+
+dimensions          [ 0 0 0 0 0 0 0 ];
+
+internalField       uniform 0;
+
+boundaryField
+{
+    walls
+    {
+        type                zeroGradient;
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/alphat b/tutorials/lagrangian/sprayFoam/aachenBomb/0/alphat
new file mode 100644
index 00000000000..e4dedc2b7c2
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/alphat
@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       volScalarField;
+    location    "0";
+    object      alphat;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    walls
+    {
+        type            alphatWallFunction;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/epsilon b/tutorials/lagrangian/sprayFoam/aachenBomb/0/epsilon
new file mode 100644
index 00000000000..6f3bc5bfc6a
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/epsilon
@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       volScalarField;
+    location    "0";
+    object      epsilon;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -3 0 0 0 0];
+
+internalField   uniform 90;
+
+boundaryField
+{
+    walls
+    {
+        type            compressible::epsilonWallFunction;
+        value           uniform 90;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/k b/tutorials/lagrangian/sprayFoam/aachenBomb/0/k
new file mode 100644
index 00000000000..29b4aafe910
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/k
@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    walls
+    {
+        type            compressible::kqRWallFunction;
+        value           uniform 1;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/mut b/tutorials/lagrangian/sprayFoam/aachenBomb/0/mut
new file mode 100644
index 00000000000..c57432b38a1
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/mut
@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       volScalarField;
+    location    "0";
+    object      mut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    walls
+    {
+        type            mutkWallFunction;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/p b/tutorials/lagrangian/sprayFoam/aachenBomb/0/p
new file mode 100644
index 00000000000..ec8fdc7ad18
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/0/p
@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// ************************************************************************* //
+
+dimensions          [ 1 -1 -2 0 0 0 0 ];
+
+internalField       uniform 5e+06;
+
+boundaryField
+{
+    walls
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/chemkin/chem.inp b/tutorials/lagrangian/sprayFoam/aachenBomb/chemkin/chem.inp
new file mode 100644
index 00000000000..0d914724259
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/chemkin/chem.inp
@@ -0,0 +1,11 @@
+ELEMENTS
+ H   O    C   N   AR
+END
+SPECIE
+C7H16 O2 N2 CO2 H2O
+END
+REACTIONS
+ C7H16 + 11O2            => 7CO2 + 8H2O        5.00E+8  0.0   15780.0! 1
+	FORD	/ C7H16	0.25 /
+	FORD	/ O2 1.5 /
+END
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/RASProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/RASProperties
new file mode 100644
index 00000000000..b94225f096c
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/RASProperties
@@ -0,0 +1,25 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       dictionary;
+    location    "constant";
+    object      RASProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+RASModel        kEpsilon;
+
+turbulence      on;
+
+printCoeffs     on;
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/chemistryProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/chemistryProperties
new file mode 100644
index 00000000000..3a7d1a555bb
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/chemistryProperties
@@ -0,0 +1,37 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      chemistryProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+psiChemistryModel ODEChemistryModel<gasThermoPhysics>;
+
+chemistry       off;
+
+chemistrySolver ode;
+
+initialChemicalTimeStep 1e-07;
+
+turbulentReaction yes;
+
+Cmix Cmix [0 0 0 0 0] 1;
+
+odeCoeffs
+{
+    ODESolver       SIBS;
+    eps             0.05;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/g b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/g
new file mode 100644
index 00000000000..24985ce8029
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/g
@@ -0,0 +1,22 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       uniformDimensionedVectorField;
+    location    "constant";
+    object      g;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -2 0 0 0 0];
+value           ( 0 -9.81 0 );
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/blockMeshDict b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/blockMeshDict
new file mode 100644
index 00000000000..822baede793
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/blockMeshDict
@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 0.001;
+
+vertices
+(
+    (-10 0 -10)
+    (-10 0 10)
+    (10 0 10)
+    (10 0 -10)
+    (-10 100 -10)
+    (-10 100 10)
+    (10 100 10)
+    (10 100 -10)
+);
+
+blocks
+(
+    hex (0 1 2 3 4 5 6 7) (41 41 100) simpleGrading (1 1 1)
+);
+
+edges
+(
+);
+
+patches
+(
+    wall walls
+    (
+        (2 6 5 1)
+        (0 4 7 3)
+        (0 1 5 4)
+        (4 5 6 7)
+        (7 6 2 3)
+        (3 2 1 0)
+    )
+);
+
+mergePatchPairs
+(
+);
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/boundary b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/boundary
new file mode 100644
index 00000000000..278bfcb94fe
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/polyMesh/boundary
@@ -0,0 +1,28 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       polyBoundaryMesh;
+    location    "constant/polyMesh";
+    object      boundary;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+1
+(
+    walls
+    {
+        type            wall;
+        nFaces          19762;
+        startFace       494419;
+    }
+)
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/radiationProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/radiationProperties
new file mode 100644
index 00000000000..3e3cafa0a5f
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/radiationProperties
@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      radiationProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+radiation       off;
+
+radiationModel  none;
+
+solverFreq      10;
+
+P1Coeffs
+{
+    C               C [ 0 0 0 0 0 0 0 ] 0;
+}
+
+absorptionEmissionModel binaryAbsorptionEmission;
+
+binaryAbsorptionEmissionCoeffs
+{
+    model1
+    {
+        absorptionEmissionModel constantAbsorptionEmission;
+        constantAbsorptionEmissionCoeffs
+        {
+            a               a [ 0 -1 0 0 0 0 0 ] 0.5;
+            e               e [ 0 -1 0 0 0 0 0 ] 0.5;
+            E               E [ 1 -1 -3 0 0 0 0 ] 0;
+        }
+    }
+    model2
+    {
+        absorptionEmissionModel cloudAbsorptionEmission;
+        cloudAbsorptionEmissionCoeffs
+        {
+            cloudNames
+            (
+                sprayCloud
+            );
+        }
+    }
+}
+
+scatterModel    cloudScatter;
+
+cloudScatterCoeffs
+{
+    cloudNames
+    (
+        sprayCloud
+    );
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/sprayCloudProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/sprayCloudProperties
new file mode 100644
index 00000000000..4a033313dc5
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/sprayCloudProperties
@@ -0,0 +1,215 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       dictionary;
+    location    "constant";
+    object      SprayCloudProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solution
+{
+    active          true;
+    coupled         true;
+    transient       yes;
+    cellValueSourceCorrection on;
+
+    sourceTerms
+    {
+        schemes
+        {
+            rho             explicit 1;
+            U               explicit 1;
+            Yi              explicit 1;
+            hs              explicit 1;
+        }
+    }
+
+    interpolationSchemes
+    {
+        rho             cell;
+        U               cellPoint;
+        mu              cell;
+        T               cell;
+        Cp              cell;
+        p               cell;
+    }
+
+    integrationSchemes
+    {
+        U               Euler;
+        T               analytical;
+    }
+}
+
+
+constantProperties
+{
+    parcelTypeId    1;
+
+    rhoMin          1e-15;
+    TMin            200;
+    pMin            1000;
+    minParticleMass 1e-15;
+
+    rho0            1000;
+    T0              300;
+    Cp0             4187;
+
+    youngsModulus   1e9;
+    poissonsRatio   0.35;
+
+    epsilon0        1;
+    f0              0.5;
+    Pr              0.7;
+    Tvap            273;
+    Tbp             373;
+
+    constantVolume  false;
+}
+
+
+subModels
+{
+    particleForces
+    {
+        sphereDrag;
+    }
+
+    injectionModel  coneNozzleInjection;
+
+    dispersionModel none;
+
+    patchInteractionModel standardWallInteraction;
+
+    heatTransferModel RanzMarshall;
+
+    compositionModel singlePhaseMixture;
+
+    phaseChangeModel liquidEvaporation;
+
+    surfaceFilmModel none;
+
+    atomizationModel none;
+
+    breakupModel    PilchErdman;
+
+    stochasticCollisionModel none;
+
+    radiation       off;
+
+    coneNozzleInjectionCoeffs
+    {
+        SOI             0;
+        massTotal       6.0e-6;
+        parcelBasisType mass;
+        injectionMethod disc;
+        outerNozzleDiameter 1.9e-4;
+        innerNozzleDiameter 0;
+        duration        1.25e-3;
+        position        ( 0 0.0995 0 );
+        direction       ( 0 -1 0 );
+        parcelsPerSecond 100000000;
+        volumeFlowRate  table
+        (
+            (0 0.1272)
+            (4.16667e-05 6.1634)
+            (8.33333e-05 9.4778)
+            (0.000125 9.5806)
+            (0.000166667 9.4184)
+            (0.000208333 9.0926)
+            (0.00025 8.7011)
+            (0.000291667 8.2239)
+            (0.000333333 8.0401)
+            (0.000375 8.845)
+            (0.000416667 8.9174)
+            (0.000458333 8.8688)
+            (0.0005 8.8882)
+            (0.000541667 8.6923)
+            (0.000583333 8.0014)
+            (0.000625 7.2582)
+            (0.000666667 7.2757)
+            (0.000708333 6.968)
+            (0.00075 6.7608)
+            (0.000791667 6.6502)
+            (0.000833333 6.7695)
+            (0.000875 5.5774)
+            (0.000916667 4.8649)
+            (0.000958333 5.0805)
+            (0.001 4.9547)
+            (0.00104167 4.5613)
+            (0.00108333 4.4536)
+            (0.001125 5.2651)
+            (0.00116667 5.256)
+            (0.00120833 5.1737)
+            (0.00125 3.9213)
+        );
+
+        Cd              constant 0.85;
+
+        thetaInner      constant 0.0;
+        thetaOuter      constant 10.0;
+
+        sizeDistribution
+        {
+            type            uniform;
+
+            uniformDistribution
+            {
+                minValue        1.0e-4;
+                maxValue        1.1e-4;
+            }
+
+        }
+    }
+
+    standardWallInteractionCoeffs
+    {
+        type            rebound;
+    }
+
+    RanzMarshallCoeffs
+    {
+        BirdCorrection  true;
+    }
+
+    singlePhaseMixtureCoeffs
+    {
+        phases
+        (
+            liquid
+            {
+                C7H16               1;
+            }
+        );
+    }
+
+    liquidEvaporationCoeffs
+    {
+        enthalpyTransfer enthalpyDifference;
+
+        activeLiquids    ( C7H16 );
+    }
+
+    PilchErdmanCoeffs
+    {
+        solveOscillationEq false;
+        B1              0.375;
+        B2              0.236;
+    }
+}
+
+
+cloudFunctions
+{}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/thermophysicalProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/thermophysicalProperties
new file mode 100644
index 00000000000..8c8ddf3c6df
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/thermophysicalProperties
@@ -0,0 +1,41 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      thermophysicalProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+thermoType      hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
+
+CHEMKINFile     "$FOAM_CASE/chemkin/chem.inp";
+
+CHEMKINThermoFile "~OpenFOAM/thermoData/therm.dat";
+
+inertSpecie     N2;
+
+
+liquids
+{
+    liquidComponents ( C7H16 );
+    C7H16
+    {
+        defaultCoeffs   yes;
+    }
+}
+
+solids
+{
+    solidComponents ( );
+}
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/constant/turbulenceProperties b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/turbulenceProperties
new file mode 100644
index 00000000000..8750aee18e1
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/constant/turbulenceProperties
@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RASModel;
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/system/controlDict b/tutorials/lagrangian/sprayFoam/aachenBomb/system/controlDict
new file mode 100644
index 00000000000..38dd94a7dcb
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/system/controlDict
@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     sprayFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         0.01;
+
+deltaT          2.5e-06;
+
+writeControl    adjustableRunTime;
+
+writeInterval   5e-05;
+
+purgeWrite      0;
+
+writeFormat     binary;
+
+writePrecision  6;
+
+writeCompression uncompressed;
+
+timeFormat      general;
+
+timePrecision   6;
+
+adjustTimeStep  yes;
+
+maxCo           0.1;
+
+runTimeModifiable yes;
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSchemes b/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSchemes
new file mode 100644
index 00000000000..6c5f4f9422d
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSchemes
@@ -0,0 +1,72 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+    grad(p)         Gauss linear;
+}
+
+divSchemes
+{
+    default         none;
+    div(phi,U)      Gauss upwind;
+    div(phid,p)     Gauss upwind;
+    div(phiU,p)     Gauss linear;
+    div(phi,hs)     Gauss upwind;
+    div(phi,k)      Gauss upwind;
+    div(phi,epsilon) Gauss upwind;
+    div(U)          Gauss linear;
+    div((muEff*dev2(T(grad(U))))) Gauss linear;
+    div(phi,Yi_h)   Gauss upwind;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear corrected;
+    laplacian(muEff,U) Gauss linear corrected;
+    laplacian(mut,U) Gauss linear corrected;
+    laplacian(DkEff,k) Gauss linear corrected;
+    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
+    laplacian(DREff,R) Gauss linear corrected;
+    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
+    laplacian(alphaEff,hs) Gauss linear corrected;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         corrected;
+}
+
+fluxRequired
+{
+    default         no;
+    p               ;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSolution b/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSolution
new file mode 100644
index 00000000000..5391299bfc5
--- /dev/null
+++ b/tutorials/lagrangian/sprayFoam/aachenBomb/system/fvSolution
@@ -0,0 +1,94 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    rho
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-05;
+        relTol          0.1;
+    }
+
+    rhoFinal
+    {
+        $rho;
+        tolerance       1e-05;
+        relTol          0;
+    }
+
+    "(U|hs)"
+    {
+        solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-05;
+        relTol          0.1;
+    }
+
+    p
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-06;
+        relTol          0.1;
+    }
+
+    pFinal
+    {
+        $p;
+        tolerance       1e-06;
+        relTol          0;
+    }
+
+    "(U|hs)Final"
+    {
+        $U;
+        tolerance       1e-05;
+        relTol          0;
+    }
+
+    "(k|epsilon)"
+    {
+        $U;
+    }
+
+    "(k|epsilon)Final"
+    {
+        $UFinal;
+    }
+
+    "(Yi|O2|N2|H2O)"
+    {
+        $UFinal;
+    }
+}
+
+PIMPLE
+{
+    transonic       no;
+    nCorrectors     2;
+    nNonOrthogonalCorrectors 0;
+    momentumPredictor yes;
+}
+
+relaxationFactors
+{
+    ".*Final"       1;
+}
+
+// ************************************************************************* //
-- 
GitLab