diff --git a/tutorials/multicomponentFluid/aachenBomb/0/N2 b/tutorials/multicomponentFluid/aachenBomb/0/N2
index 729adeb35b..76973f4a3f 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/N2
+++ b/tutorials/multicomponentFluid/aachenBomb/0/N2
@@ -19,6 +19,8 @@ internalField       uniform 0.766;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type                zeroGradient;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/O2 b/tutorials/multicomponentFluid/aachenBomb/0/O2
index f03d59f6be..9b9983e9d5 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/O2
+++ b/tutorials/multicomponentFluid/aachenBomb/0/O2
@@ -19,6 +19,8 @@ internalField       uniform 0.234;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type                zeroGradient;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/T b/tutorials/multicomponentFluid/aachenBomb/0/T
index 2d1e607b83..388de1c14a 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/T
+++ b/tutorials/multicomponentFluid/aachenBomb/0/T
@@ -19,6 +19,8 @@ internalField       uniform 800;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type                zeroGradient;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/U b/tutorials/multicomponentFluid/aachenBomb/0/U
index a0e3ef8f01..76358880f4 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/U
+++ b/tutorials/multicomponentFluid/aachenBomb/0/U
@@ -19,6 +19,8 @@ internalField       uniform (0 0 0);
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type                fixedValue;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/Ydefault b/tutorials/multicomponentFluid/aachenBomb/0/Ydefault
index 06d65c2d3c..2415ecf101 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/Ydefault
+++ b/tutorials/multicomponentFluid/aachenBomb/0/Ydefault
@@ -19,6 +19,8 @@ internalField       uniform 0;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type                zeroGradient;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/alphat b/tutorials/multicomponentFluid/aachenBomb/0/alphat
index c90e2650a9..84e7b611e4 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/alphat
+++ b/tutorials/multicomponentFluid/aachenBomb/0/alphat
@@ -20,6 +20,8 @@ internalField   uniform 0;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type            compressible::alphatWallFunction;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/epsilon b/tutorials/multicomponentFluid/aachenBomb/0/epsilon
index c5a1ee33ae..07bb3ac83a 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/epsilon
+++ b/tutorials/multicomponentFluid/aachenBomb/0/epsilon
@@ -20,6 +20,8 @@ internalField   uniform 90;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type            epsilonWallFunction;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/k b/tutorials/multicomponentFluid/aachenBomb/0/k
index 1d9f8c67e5..7618a42a34 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/k
+++ b/tutorials/multicomponentFluid/aachenBomb/0/k
@@ -20,6 +20,8 @@ internalField   uniform 1;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type            kqRWallFunction;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/nut b/tutorials/multicomponentFluid/aachenBomb/0/nut
index 43c4d9c817..ce9f6c1098 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/nut
+++ b/tutorials/multicomponentFluid/aachenBomb/0/nut
@@ -20,6 +20,8 @@ internalField   uniform 0;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type            nutkWallFunction;
diff --git a/tutorials/multicomponentFluid/aachenBomb/0/p b/tutorials/multicomponentFluid/aachenBomb/0/p
index 8b64f298e6..e6cccbcb28 100644
--- a/tutorials/multicomponentFluid/aachenBomb/0/p
+++ b/tutorials/multicomponentFluid/aachenBomb/0/p
@@ -19,6 +19,8 @@ internalField       uniform 5e+06;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     walls
     {
         type            zeroGradient;
diff --git a/tutorials/multicomponentFluid/aachenBomb/Allrun b/tutorials/multicomponentFluid/aachenBomb/Allrun
index ed4eb7ef27..ebb7bfe67c 100755
--- a/tutorials/multicomponentFluid/aachenBomb/Allrun
+++ b/tutorials/multicomponentFluid/aachenBomb/Allrun
@@ -4,14 +4,11 @@ cd ${0%/*} || exit 1    # Run from this directory
 # Source tutorial run functions
 . $WM_PROJECT_DIR/bin/tools/RunFunctions
 
-# Set application name
-application=$(getApplication)
-
 runApplication chemkinToFoam \
                chemkin/chem.inp chemkin/therm.dat chemkin/transportProperties \
                constant/reactions constant/speciesThermo
 
 runApplication blockMesh
-runApplication $application
+runApplication $(getApplication)
 
 #------------------------------------------------------------------------------
diff --git a/tutorials/multicomponentFluid/aachenBomb/Allrun-parallel b/tutorials/multicomponentFluid/aachenBomb/Allrun-parallel
new file mode 100755
index 0000000000..8c242e3f84
--- /dev/null
+++ b/tutorials/multicomponentFluid/aachenBomb/Allrun-parallel
@@ -0,0 +1,17 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+runApplication chemkinToFoam \
+               chemkin/chem.inp chemkin/therm.dat chemkin/transportProperties \
+               constant/reactions constant/speciesThermo
+
+runApplication blockMesh
+
+runApplication decomposePar
+runParallel $(getApplication)
+runApplication reconstructPar
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/multicomponentFluid/aachenBomb/constant/chemistryProperties b/tutorials/multicomponentFluid/aachenBomb/constant/chemistryProperties
index cf75029281..2bda592323 100644
--- a/tutorials/multicomponentFluid/aachenBomb/constant/chemistryProperties
+++ b/tutorials/multicomponentFluid/aachenBomb/constant/chemistryProperties
@@ -21,6 +21,8 @@ chemistryType
 
 chemistry       on;
 
+cpuLoad         true;
+
 initialChemicalTimeStep 1e-07;
 
 odeCoeffs
diff --git a/tutorials/multicomponentFluid/aachenBomb/constant/cloudProperties b/tutorials/multicomponentFluid/aachenBomb/constant/cloudProperties
index d28023a0bb..b182787035 100644
--- a/tutorials/multicomponentFluid/aachenBomb/constant/cloudProperties
+++ b/tutorials/multicomponentFluid/aachenBomb/constant/cloudProperties
@@ -16,6 +16,8 @@ FoamFile
 
 type        sprayCloud;
 
+cpuLoad     true;
+
 solution
 {
     coupled         true;
diff --git a/tutorials/multicomponentFluid/aachenBomb/constant/dynamicMeshDict b/tutorials/multicomponentFluid/aachenBomb/constant/dynamicMeshDict
new file mode 100644
index 0000000000..a07e40541b
--- /dev/null
+++ b/tutorials/multicomponentFluid/aachenBomb/constant/dynamicMeshDict
@@ -0,0 +1,27 @@
+/*--------------------------------*- 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      dynamicMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+distributor
+{
+    type            loadBalancer;
+
+    libs            ("libfvMeshDistributors.so");
+
+    multiConstraint true;
+    redistributionInterval  20;
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/multicomponentFluid/aachenBomb/system/blockMeshDict b/tutorials/multicomponentFluid/aachenBomb/system/blockMeshDict
index 4e9d971585..9083f9ad86 100644
--- a/tutorials/multicomponentFluid/aachenBomb/system/blockMeshDict
+++ b/tutorials/multicomponentFluid/aachenBomb/system/blockMeshDict
@@ -47,6 +47,12 @@ boundary
             (3 2 1 0)
         );
     }
+
+    internalFaces
+    {
+        type internal;
+        faces ();
+    }
 );
 
 
diff --git a/tutorials/multicomponentFluid/aachenBomb/system/decomposeParDict b/tutorials/multicomponentFluid/aachenBomb/system/decomposeParDict
new file mode 100644
index 0000000000..0b0a5b4976
--- /dev/null
+++ b/tutorials/multicomponentFluid/aachenBomb/system/decomposeParDict
@@ -0,0 +1,37 @@
+/*--------------------------------*- 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;
+    location    "system";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 12;
+
+decomposer      simple;
+
+distributor     zoltan;
+libs            ("libzoltanDecomp.so");
+
+// distributor     parMetis;
+// libs            ("libparMetisDecomp.so");
+
+simpleCoeffs
+{
+    n           (2 2 3);
+}
+
+zoltanCoeffs
+{
+    lb_method   rcb;
+}
+
+// ************************************************************************* //