diff --git a/tutorials/lagrangian/sprayFoam/aachenBomb/0/N2 b/tutorials/lagrangian/sprayFoam/aachenBomb/0/N2
new file mode 100644
index 0000000000000000000000000000000000000000..291efd9c59f9fe4625e0de4859a8c6c195dff97e
--- /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 0000000000000000000000000000000000000000..bacf89cb28e2548f3a353903146f2f4dff99ce19
--- /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 0000000000000000000000000000000000000000..b6ffeed49c6cbe61c3465714fbd096c662a4dc6a
--- /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 0000000000000000000000000000000000000000..98260b8fb4b6cfb5e9283c77d5647c8ffb7cfe6e
--- /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 0000000000000000000000000000000000000000..1039e373fe58da0263bd3e7dc4439da0ffb1ac8c
--- /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 0000000000000000000000000000000000000000..e4dedc2b7c2fe0a64c3e3cdb9dc260bd187f0577
--- /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 0000000000000000000000000000000000000000..6f3bc5bfc6a0ff483120434adad63b6dfa808e02
--- /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 0000000000000000000000000000000000000000..29b4aafe910a5f5559176cc729f1784d8f7e410b
--- /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 0000000000000000000000000000000000000000..c57432b38a161254d4fa0bdb71d66fddb58d4ad4
--- /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 0000000000000000000000000000000000000000..ec8fdc7ad188a793c76821d36107a1bd5e87a313
--- /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 0000000000000000000000000000000000000000..0d914724259f8f13476daa8bf94cc2304487e4c7
--- /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 0000000000000000000000000000000000000000..b94225f096c9a656b5d69defbc0e0cac621e51b7
--- /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 0000000000000000000000000000000000000000..3a7d1a555bbec1ec28710b42867483c893df9741
--- /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 0000000000000000000000000000000000000000..24985ce80296be655fc60df1f60b9dadbf9ca99c
--- /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 0000000000000000000000000000000000000000..822baede79395fb5541c76bdcb2ff22244ca9b06
--- /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 0000000000000000000000000000000000000000..278bfcb94fe598d9e37c952bf646646fd8ab6fe3
--- /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 0000000000000000000000000000000000000000..3e3cafa0a5f90bc02a31ee6676df8c1eb996d376
--- /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 0000000000000000000000000000000000000000..4a033313dc5b1ad83f29af80f418dff055a8db47
--- /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 0000000000000000000000000000000000000000..8c8ddf3c6dfc00ac41ee6fa01f8a88be1aa4a9ce
--- /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 0000000000000000000000000000000000000000..8750aee18e18ed1ab7b36f50670044cf6778cbc7
--- /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 0000000000000000000000000000000000000000..38dd94a7dcb3f6ae340dae9b13919bcbdf5698b5
--- /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 0000000000000000000000000000000000000000..6c5f4f9422d34484d90295bc77bca7b61bcfcd23
--- /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 0000000000000000000000000000000000000000..5391299bfc561010d4bf041f97da2b1da59f1463
--- /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;
+}
+
+// ************************************************************************* //