diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H.gz
new file mode 100644
index 000000000..1dca86ce9
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2.gz
new file mode 100644
index 000000000..16d49e23b
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O.gz
new file mode 100644
index 000000000..fd61cd404
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O2.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O2.gz
new file mode 100644
index 000000000..d5e52d995
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/H2O2.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/HO2.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/HO2.gz
new file mode 100644
index 000000000..92fd7d76c
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/HO2.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/N2.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/N2.gz
new file mode 100644
index 000000000..ed9d83f2c
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/N2.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O.gz
new file mode 100644
index 000000000..e92aa33fd
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O2.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O2.gz
new file mode 100644
index 000000000..3aed4f3b7
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/O2.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/OH.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/OH.gz
new file mode 100644
index 000000000..f27e9dd3b
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/OH.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/T.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/T.gz
new file mode 100644
index 000000000..49e66201d
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/T.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/U.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/U.gz
new file mode 100644
index 000000000..0f020faa1
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/U.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/p.gz b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/p.gz
new file mode 100644
index 000000000..bef352a51
Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/0/p.gz differ
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allclean b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allclean
new file mode 100644
index 000000000..2cd9e2246
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allclean
@@ -0,0 +1,9 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+echo "Cleaning log.*"
+rm log.*
+echo "Cleaning processor*"
+rm -r processor*
+echo "Cleaning polyMesh/"
+rm -r constant/polyMesh
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allrun b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allrun
new file mode 100644
index 000000000..e4baa099e
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Allrun
@@ -0,0 +1,11 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application=dfLowMachFoam
+
+runApplication blockMesh
+runApplication decomposePar
+runApplication mpirun -np 8 --allow-run-as-root $application -parallel
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Burke2012_s9r23.yaml b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Burke2012_s9r23.yaml
new file mode 100644
index 000000000..d77a77bbf
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/Burke2012_s9r23.yaml
@@ -0,0 +1,242 @@
+generator: cti2yaml
+cantera-version: 2.5.1
+date: Tue, 24 Oct 2023 17:09:01 +0800
+input-files: [Burke2012_s9s.cti]
+
+units: {length: cm, quantity: mol, activation-energy: cal/mol}
+
+phases:
+- name: gas
+  thermo: ideal-gas
+  elements: [H, O, N]
+  species: [H, H2, O, OH, H2O, O2, HO2, H2O2, N2]
+  kinetics: gas
+  reactions: all
+  transport: mixture-averaged
+  state:
+    T: 300.0
+    P: 1.01325e+05
+
+species:
+- name: H
+  composition: {H: 1}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [2.5, 0.0, 0.0, 0.0, 0.0, 2.547163e+04, -0.4601176]
+    - [2.5, 0.0, 0.0, 0.0, 0.0, 2.547163e+04, -0.4601176]
+  transport:
+    model: gas
+    geometry: atom
+    diameter: 2.05
+    well-depth: 145.0
+  note: '120186'
+- name: H2
+  composition: {H: 2}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [3.298124, 8.249442e-04, -8.143015e-07, -9.475434e-11, 4.134872e-13,
+      -1012.521, -3.294094]
+    - [2.991423, 7.000644e-04, -5.633829e-08, -9.231578e-12, 1.582752e-15,
+      -835.034, -1.35511]
+  transport:
+    model: gas
+    geometry: linear
+    diameter: 2.92
+    well-depth: 38.0
+    polarizability: 0.79
+    rotational-relaxation: 280.0
+  note: '121286'
+- name: O
+  composition: {O: 1}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [2.946429, -1.638166e-03, 2.421032e-06, -1.602843e-09, 3.890696e-13,
+      2.914764e+04, 2.963995]
+    - [2.54206, -2.755062e-05, -3.102803e-09, 4.551067e-12, -4.368052e-16,
+      2.92308e+04, 4.920308]
+  transport:
+    model: gas
+    geometry: atom
+    diameter: 2.75
+    well-depth: 80.0
+  note: '120186'
+- name: OH
+  composition: {H: 1, O: 1}
+  thermo:
+    model: NASA7
+    temperature-ranges: [200.0, 1000.0, 6000.0]
+    data:
+    - [4.12530561, -3.22544939e-03, 6.52764691e-06, -5.79853643e-09, 2.06237379e-12,
+      3346.30913, -0.69043296]
+    - [2.86472886, 1.05650448e-03, -2.59082758e-07, 3.05218674e-11, -1.33195876e-15,
+      3683.62875, 5.70164073]
+  transport:
+    model: gas
+    geometry: linear
+    diameter: 2.75
+    well-depth: 80.0
+  note: S9/01
+- name: H2O
+  composition: {H: 2, O: 1}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [3.386842, 3.474982e-03, -6.354696e-06, 6.968581e-09, -2.506588e-12,
+      -3.020811e+04, 2.590233]
+    - [2.672146, 3.056293e-03, -8.73026e-07, 1.200996e-10, -6.391618e-15,
+      -2.989921e+04, 6.862817]
+  transport:
+    model: gas
+    geometry: nonlinear
+    diameter: 2.605
+    well-depth: 572.4
+    dipole: 1.844
+    rotational-relaxation: 4.0
+  note: '20387'
+- name: O2
+  composition: {O: 2}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [3.212936, 1.127486e-03, -5.75615e-07, 1.313877e-09, -8.768554e-13,
+      -1005.249, 6.034738]
+    - [3.697578, 6.135197e-04, -1.258842e-07, 1.775281e-11, -1.136435e-15,
+      -1233.93, 3.189166]
+  transport:
+    model: gas
+    geometry: linear
+    diameter: 3.458
+    well-depth: 107.4
+    polarizability: 1.6
+    rotational-relaxation: 3.8
+  note: '121386'
+- name: HO2
+  composition: {H: 1, O: 2}
+  thermo:
+    model: NASA7
+    temperature-ranges: [200.0, 1000.0, 3500.0]
+    data:
+    - [4.30179801, -4.74912051e-03, 2.11582891e-05, -2.42763894e-08, 9.29225124e-12,
+      294.80804, 3.71666245]
+    - [4.0172109, 2.23982013e-03, -6.3365815e-07, 1.1424637e-10, -1.07908535e-14,
+      111.856713, 3.78510215]
+  transport:
+    model: gas
+    geometry: nonlinear
+    diameter: 3.458
+    well-depth: 107.4
+    rotational-relaxation: 1.0
+  note: L5/89
+- name: H2O2
+  composition: {H: 2, O: 2}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [3.388754, 6.569226e-03, -1.485013e-07, -4.625806e-09, 2.471515e-12,
+      -1.766315e+04, 6.785363]
+    - [4.573167, 4.336136e-03, -1.474689e-06, 2.348904e-10, -1.431654e-14,
+      -1.800696e+04, 0.501137]
+  transport:
+    model: gas
+    geometry: nonlinear
+    diameter: 3.458
+    well-depth: 107.4
+    rotational-relaxation: 3.8
+  note: '120186'
+- name: N2
+  composition: {N: 2}
+  thermo:
+    model: NASA7
+    temperature-ranges: [300.0, 1000.0, 5000.0]
+    data:
+    - [3.298677, 1.40824e-03, -3.963222e-06, 5.641515e-09, -2.444855e-12,
+      -1020.9, 3.950372]
+    - [2.92664, 1.487977e-03, -5.684761e-07, 1.009704e-10, -6.753351e-15,
+      -922.7977, 5.980528]
+  transport:
+    model: gas
+    geometry: linear
+    diameter: 3.621
+    well-depth: 97.53
+    polarizability: 1.76
+    rotational-relaxation: 4.0
+  note: '121286'
+
+reactions:
+- equation: H + O2 <=> O + OH  # Reaction 1
+  rate-constant: {A: 1.04e+14, b: 0.0, Ea: 1.5286e+04}
+- equation: O + H2 <=> H + OH  # Reaction 2
+  rate-constant: {A: 3.818e+12, b: 0.0, Ea: 7948.0}
+  duplicate: true
+- equation: O + H2 <=> H + OH  # Reaction 3
+  rate-constant: {A: 8.792e+14, b: 0.0, Ea: 1.917e+04}
+  duplicate: true
+- equation: H2 + OH <=> H2O + H  # Reaction 4
+  rate-constant: {A: 2.16e+08, b: 1.51, Ea: 3430.0}
+- equation: OH + OH <=> O + H2O  # Reaction 5
+  rate-constant: {A: 3.34e+04, b: 2.42, Ea: -1930.0}
+- equation: H2 + M <=> H + H + M  # Reaction 6
+  type: three-body
+  rate-constant: {A: 4.577e+19, b: -1.4, Ea: 1.0438e+05}
+  efficiencies: {H2: 2.5, H2O: 12.0}
+- equation: O + O + M <=> O2 + M  # Reaction 7
+  type: three-body
+  rate-constant: {A: 6.165e+15, b: -0.5, Ea: 0.0}
+  efficiencies: {H2: 2.5, H2O: 12.0}
+- equation: O + H + M <=> OH + M  # Reaction 8
+  type: three-body
+  rate-constant: {A: 4.714e+18, b: -1.0, Ea: 0.0}
+  efficiencies: {H2: 2.5, H2O: 12.0}
+- equation: H2O + M <=> H + OH + M  # Reaction 9
+  type: three-body
+  rate-constant: {A: 6.064e+27, b: -3.322, Ea: 1.2079e+05}
+  efficiencies: {H2: 3.0, H2O: 0.0, N2: 2.0, O2: 1.5}
+- equation: H2O + H2O <=> H + OH + H2O  # Reaction 10
+  rate-constant: {A: 1.006e+26, b: -2.44, Ea: 1.2018e+05}
+- equation: H + O2 (+ M) <=> HO2 (+ M)  # Reaction 11
+  type: falloff
+  low-P-rate-constant: {A: 6.366e+20, b: -1.72, Ea: 524.8}
+  high-P-rate-constant: {A: 4.65084e+12, b: 0.44, Ea: 0.0}
+  Troe: {A: 0.5, T3: 1.0e-30, T1: 1.0e+30}
+  efficiencies: {H2: 2.0, H2O: 14.0, O2: 0.78}
+- equation: HO2 + H <=> H2 + O2  # Reaction 12
+  rate-constant: {A: 2.75e+06, b: 2.09, Ea: -1451.0}
+- equation: HO2 + H <=> OH + OH  # Reaction 13
+  rate-constant: {A: 7.079e+13, b: 0.0, Ea: 295.0}
+- equation: HO2 + O <=> O2 + OH  # Reaction 14
+  rate-constant: {A: 2.85e+10, b: 1.0, Ea: -723.93}
+- equation: HO2 + OH <=> H2O + O2  # Reaction 15
+  rate-constant: {A: 2.89e+13, b: 0.0, Ea: -497.0}
+- equation: HO2 + HO2 <=> H2O2 + O2  # Reaction 16
+  rate-constant: {A: 4.2e+14, b: 0.0, Ea: 1.1982e+04}
+  duplicate: true
+- equation: HO2 + HO2 <=> H2O2 + O2  # Reaction 17
+  rate-constant: {A: 1.3e+11, b: 0.0, Ea: -1629.3}
+  duplicate: true
+- equation: H2O2 (+ M) <=> OH + OH (+ M)  # Reaction 18
+  type: falloff
+  low-P-rate-constant: {A: 2.49e+24, b: -2.3, Ea: 4.8749e+04}
+  high-P-rate-constant: {A: 2.0e+12, b: 0.9, Ea: 4.8749e+04}
+  Troe: {A: 0.43, T3: 1.0e-30, T1: 1.0e+30}
+  efficiencies: {H2: 3.7, H2O: 7.5, H2O2: 7.7, N2: 1.5, O2: 1.2}
+- equation: H2O2 + H <=> H2O + OH  # Reaction 19
+  rate-constant: {A: 2.41e+13, b: 0.0, Ea: 3970.0}
+- equation: H2O2 + H <=> HO2 + H2  # Reaction 20
+  rate-constant: {A: 4.82e+13, b: 0.0, Ea: 7950.0}
+- equation: H2O2 + O <=> OH + HO2  # Reaction 21
+  rate-constant: {A: 9.55e+06, b: 2.0, Ea: 3970.0}
+- equation: H2O2 + OH <=> HO2 + H2O  # Reaction 22
+  rate-constant: {A: 1.74e+12, b: 0.0, Ea: 318.0}
+  duplicate: true
+- equation: H2O2 + OH <=> HO2 + H2O  # Reaction 23
+  rate-constant: {A: 7.59e+13, b: 0.0, Ea: 7270.0}
+  duplicate: true
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/CanteraTorchProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/CanteraTorchProperties
new file mode 100644
index 000000000..9bb575dc2
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/CanteraTorchProperties
@@ -0,0 +1,47 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version          2.0;
+    format           ascii;
+    class            dictionary;
+    location         "constant";
+    object           CanteraTorchProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+chemistry            on;
+
+CanteraMechanismFile "Burke2012_s9r23.yaml";
+
+transportModel       "Mix";
+
+odeCoeffs
+{
+    "relTol"         1e-6;
+    "absTol"         1e-10;
+}
+
+inertSpecie          "N2";
+
+splittingStrategy    off;
+
+TorchSettings
+{
+    torch            off;
+}
+
+loadbalancing
+{
+    active           true;
+    log              false;
+    algorithm        allAverage;//headTail;
+}
+
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/combustionProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/combustionProperties
new file mode 100644
index 000000000..7418500d0
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/combustionProperties
@@ -0,0 +1,20 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      combustionProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+combustionModel  laminar;
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/g b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/g
new file mode 100644
index 000000000..28cfc2369
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/g
@@ -0,0 +1,22 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     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/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/sprayCloudProperties
new file mode 100644
index 000000000..2010a8858
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/sprayCloudProperties
@@ -0,0 +1,26 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      binary;
+    class       dictionary;
+    location    "constant";
+    object      SprayCloudProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solution
+{
+    active          false;
+    coupled         false;
+    transient       no;
+    cellValueSourceCorrection off;
+}
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/thermophysicalProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/thermophysicalProperties
new file mode 100644
index 000000000..f82ee9785
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/thermophysicalProperties
@@ -0,0 +1,18 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.4.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      thermophysicalProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/turbulenceProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/turbulenceProperties
new file mode 100644
index 000000000..ab70e94d3
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/constant/turbulenceProperties
@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  laminar;
+
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/blockMeshDict b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/blockMeshDict
new file mode 100644
index 000000000..38710de65
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/blockMeshDict
@@ -0,0 +1,126 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 0.001;
+
+vertices #codeStream
+{
+    codeInclude
+    #{
+        #include "pointField.H"
+    #};
+
+    code
+    #{
+        pointField points(8);
+        points[0]  = point(0, 0, -0.5);
+        points[1]  = point(16, 0, -0.5);
+        points[2]  = point(16, 7.5, -0.5);
+        points[3]  = point(0, 7.5, -0.5);
+        points[4]  = point(0, 8.5, -0.5);
+        points[5]  = point(16, 8.5, -0.5);
+        points[6]  = point(16, 16, -0.5);
+        points[7]  = point(0, 16, -0.5);
+
+        // Duplicate z points
+        label sz = points.size();
+        points.setSize(2*sz);
+        for (label i = 0; i < sz; i++)
+        {
+            const point& pt = points[i];
+            points[i+sz] = point(pt.x(), pt.y(), -pt.z());
+        }
+
+        os  << points;
+    #};
+};
+
+blocks
+(
+    hex (0 1 2 3 8 9 10 11) (800 250 1) 
+	simpleGrading
+    (
+        1
+        (
+            (0.6 0.4 0.25)
+            (0.4 0.6 1)
+        )
+        1
+    )
+    hex (3 2 5 4 11 10 13 12) (800 50 1) simpleGrading (1 1 1)
+    hex (4 5 6 7 12 13 14 15) (800 250 1) 
+	simpleGrading
+    (
+        1
+        (
+            (0.4 0.6 1)
+            (0.6 0.4 4)
+        )
+        1
+    )
+);
+
+edges
+(
+);
+
+boundary
+(
+    left
+    {
+        type cyclic;
+        neighbourPatch right;
+        faces
+        (
+            (0 3 11 8)
+	    (3 4 12 11)
+	    (4 7 15 12)
+        );
+    }
+    right
+    {
+        type cyclic;
+        neighbourPatch left;
+        faces
+        (
+            (1 2 10 9)
+	    (2 5 13 10)
+	    (5 6 14 13)
+        );
+    }
+    top
+    {
+        type patch;
+        faces
+        (
+            (7 6 14 15)
+        );
+    }
+    down
+    {
+        type patch;
+        faces
+        (
+            (0 1 9 8)
+        );
+    }
+);
+
+// mergePatchPairs
+// (
+// );
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/controlDict b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/controlDict
new file mode 100644
index 000000000..07e45fce7
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/controlDict
@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     dfLowMachFoam;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         1.0e-4;
+
+deltaT          1e-8;
+
+writeControl    adjustableRunTime;
+
+writeInterval   1e-5;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+// writePrecision    6;
+
+writeCompression on;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable true;
+
+adjustTimeStep no;
+
+maxCo  0.3;
+
+maxDeltaT 1e-4;
+
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/decomposeParDict b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/decomposeParDict
new file mode 100644
index 000000000..7e67de0fd
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/decomposeParDict
@@ -0,0 +1,45 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 8;
+
+method          scotch;
+
+simpleCoeffs
+{
+    n               (8 4 1);
+    delta           0.001;
+}
+
+hierarchicalCoeffs
+{
+    n               (1 1 1);
+    delta           0.001;
+    order           xyz;
+}
+
+manualCoeffs
+{
+    dataFile        "";
+}
+
+distributed     no;
+
+roots           ( );
+
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSchemes b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSchemes
new file mode 100644
index 000000000..2f2cf3dd9
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSchemes
@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+}
+
+divSchemes
+{
+    default             none;
+
+    div(phi,U)          Gauss linear;
+    div(phi,Yi)         Gauss limitedLinear01 1;
+    div(phi,h)          Gauss limitedLinear 1;
+    div(phi,ha)          Gauss limitedLinear 1;
+    div(phi,K)          Gauss limitedLinear 1;
+    div(phid,p)         Gauss limitedLinear 1;
+    div(phi,epsilon)    Gauss limitedLinear 1;
+    div(phi,Yi_h)       Gauss limitedLinear01 1;
+    div(phi,k)          Gauss limitedLinear 1;
+    div(hDiffCorrFlux)             Gauss cubic;
+    div(((rho*nuEff)*dev2(T(grad(U)))))     Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear orthogonal;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         orthogonal;
+}
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSolution b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSolution
new file mode 100644
index 000000000..15823f3de
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/fvSolution
@@ -0,0 +1,89 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "rho.*"
+    {
+        solver          diagonal;
+    }
+
+    p
+    {
+        solver          GAMG;
+        tolerance       1e-7;
+        relTol          0.01;
+        smoother        GaussSeidel;
+        cacheAgglomeration true;
+        nCellsInCoarsestLevel 20;
+        agglomerator    faceAreaPair;
+        mergeLevels     1;
+    }
+
+    pFinal
+    {
+        solver          GAMG;
+        tolerance       1e-7;
+        relTol          0.01;
+        smoother        GaussSeidel;
+        cacheAgglomeration true;
+        nCellsInCoarsestLevel 20;
+        agglomerator    faceAreaPair;
+        mergeLevels     1;
+    }
+
+    "(U|ha|h|k|epsilon)"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-7;
+        relTol          0.01;
+    }
+
+    "(U|ha|h|k|epsilon)Final"
+    {
+        $U;
+        tolerance       1e-7;
+        relTol          0.01;
+    }
+
+    "Yi"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-7;
+        relTol          0.01;
+    }
+    "YiFinal"
+    {
+        solver          PBiCGStab;
+        preconditioner  DILU;
+        tolerance       1e-7;
+        relTol          0.01;
+    }
+}
+
+PIMPLE
+{
+    momentumPredictor yes;
+    nOuterCorrectors  1;
+    nCorrectors     2;
+    nNonOrthogonalCorrectors 0;
+}
+
+
+// ************************************************************************* //
diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/setFieldsDict b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/setFieldsDict
new file mode 100644
index 000000000..688953686
--- /dev/null
+++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/H2/system/setFieldsDict
@@ -0,0 +1,52 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  7
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      setFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+defaultFieldValues
+(
+    volScalarFieldValue T      1428.5
+    volScalarFieldValue H      0.0
+    volScalarFieldValue H2     0.0
+    volScalarFieldValue O      0.0
+    volScalarFieldValue OH     3.3e-05
+    volScalarFieldValue H2O    0.103708
+    volScalarFieldValue O2     0.138161
+    volScalarFieldValue HO2    0.0
+    volScalarFieldValue H2O2   0.0
+    volScalarFieldValue N2     0.758098
+);
+
+regions
+(
+    boxToCell
+    {
+        box (0 0.0075 -0.0005) (0.016 0.0085 0.0005);
+        fieldValues
+        (
+            volScalarFieldValue T      300
+            volScalarFieldValue H      0.0
+            volScalarFieldValue H2     0.011608
+            volScalarFieldValue O      0.0
+            volScalarFieldValue OH     0.0
+            volScalarFieldValue H2O    0.0
+            volScalarFieldValue O2     0.230295
+            volScalarFieldValue HO2    0.0
+            volScalarFieldValue H2O2   0.0
+            volScalarFieldValue N2     0.758097
+        );
+    }
+);
+
+// ************************************************************************* //