diff --git a/AMR_tidy/dynamicRefineFvMeshNew/multiCritRefinement.C b/AMR_tidy/dynamicRefineFvMeshNew/multiCritRefinement.C
index 556763c1..9ac82af2 100644
--- a/AMR_tidy/dynamicRefineFvMeshNew/multiCritRefinement.C
+++ b/AMR_tidy/dynamicRefineFvMeshNew/multiCritRefinement.C
@@ -317,7 +317,7 @@ void Foam::multiCritRefinement::updateRefinementField()
                 IOobject::NO_WRITE,
                 false
             )
-        ).subDict("dynamicRefineFvMeshCoeffs")
+        ).subDict("dynamicRefineFvMeshNewCoeffs")
     );
     scalar globalMaxRefLevel = readScalar(refineDict.lookup("maxRefinement"));
 
diff --git a/AMR_tidy/hexRef/hexRef2/hexRef2.C b/AMR_tidy/hexRef/hexRef2/hexRef2.C
index fbcc9ea2..c1aaef5a 100644
--- a/AMR_tidy/hexRef/hexRef2/hexRef2.C
+++ b/AMR_tidy/hexRef/hexRef2/hexRef2.C
@@ -336,23 +336,6 @@ Foam::labelListList Foam::hexRef2::setRefinement
         isDivisibleEdge[edgeI] = false;
     }
 
-    // for (label facei = mesh_.nInternalFaces(); facei < mesh_.nFaces(); facei++)
-    // {
-    //     const label & patchID = mesh_.boundaryMesh().whichPatch(facei);
-
-    //     if (mesh_.boundaryMesh()[patchID].type() != "empty")
-    //     {
-    //         isDivisibleFace[facei] = false;
-    //         const labelList& fEdges = mesh_.faceEdges(facei);
-
-    //         forAll(fEdges, i)
-    //         {
-    //             label edgeI = fEdges[i];
-    //             isDivisibleEdge[edgeI] = false;
-    //         }
-    //     }
-    // }
-
     for (label facei = 0; facei < mesh_.nInternalFaces(); facei++)
     {
         isDivisibleFace[facei] = false;
@@ -898,11 +881,13 @@ Foam::labelListList Foam::hexRef2::setRefinement
                 face newFace(2);
 
                 // maybe need to adjust the point order
-                DynamicList<label> faceVerts(4);
-                
+                // DynamicList<label> faceVerts(4);              
 
                 forAll(f, fp)
                 {
+                    //maybe need to put inside (but nothing wrong to be put outside)
+                    DynamicList<label> faceVerts(4);
+
                     label pointi = f[fp];
 
                     label nextpointi = f[f.fcIndex(fp)];
diff --git a/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C b/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C
index a561db72..3088e778 100644
--- a/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C
+++ b/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C
@@ -60,6 +60,7 @@ int main(int argc, char *argv[])
     double time_monitor_flow;
     double time_monitor_chem;
     double time_monitor_Y;
+    double time_monitor_AMR;
     clock_t start, end;
 
     turbulence->validate();
@@ -69,6 +70,7 @@ int main(int argc, char *argv[])
     #include "readFluxScheme.H"
 
     dimensionedScalar v_zero("v_zero", dimVolume/dimTime, 0.0);
+    dimensionedScalar refCri("refCri", dimensionSet(1, -4, 0, 0, 0), 0.0);
 
     // Courant numbers used to adjust the time-step
     scalar CoNum = 0.0;
@@ -80,14 +82,16 @@ int main(int argc, char *argv[])
     {
         #include "readTimeControls.H"
 
-        tmp<volScalarField> tmagGradP = mag(fvc::grad(p));
-        volScalarField normalisedGradP
+        //used for AMR
+        refCri = max(mag(fvc::grad(rho)));
+        tmp<volScalarField> tmagGradrho = mag(fvc::grad(rho));
+        volScalarField normalisedGradrho
         (
-            "normalisedGradP",
-            tmagGradP()/max(tmagGradP())
+            "normalisedGradrho",
+            tmagGradrho()/refCri
         );
-        normalisedGradP.writeOpt() = IOobject::AUTO_WRITE;
-        tmagGradP.clear();
+        normalisedGradrho.writeOpt() = IOobject::AUTO_WRITE;
+        tmagGradrho.clear();
 
         if (!LTS)
         {
@@ -95,7 +99,11 @@ int main(int argc, char *argv[])
             runTime++;
 
             // Do any mesh changes
+            start = std::clock();
             mesh.update();
+            end = std::clock();
+            time_monitor_AMR += double(end - start) / double(CLOCKS_PER_SEC);
+            
         }
 
         // --- Directed interpolation of primitive fields onto faces
@@ -234,6 +242,7 @@ int main(int argc, char *argv[])
         Info<< "MonitorTime_chem = " << time_monitor_chem << " s" << nl << endl;
         Info<< "MonitorTime_Y = " << time_monitor_Y << " s" << nl << endl;
         Info<< "MonitorTime_flow = " << time_monitor_flow << " s" << nl << endl;
+        Info<< "MonitorTime_AMR = " << time_monitor_AMR << " s" << nl << endl;
 
         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/H2.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/H2.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/H2.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/H2.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/N2.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/N2.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/N2.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/N2.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/O2.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/O2.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/O2.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/O2.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/T.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/T.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/T.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/T.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/U.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/U.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/U.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/U.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/Ydefault.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/Ydefault.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/Ydefault.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/Ydefault.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/p.gz b/examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/p.gz
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/p.gz
rename to examples/dfHighSpeedFoam/oneD_detonationH2/0_orig/p.gz
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/Allclean b/examples/dfHighSpeedFoam/oneD_detonationH2/Allclean
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/Allclean
rename to examples/dfHighSpeedFoam/oneD_detonationH2/Allclean
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/Allrun b/examples/dfHighSpeedFoam/oneD_detonationH2/Allrun
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/Allrun
rename to examples/dfHighSpeedFoam/oneD_detonationH2/Allrun
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/H2_Li.xml b/examples/dfHighSpeedFoam/oneD_detonationH2/H2_Li.xml
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/H2_Li.xml
rename to examples/dfHighSpeedFoam/oneD_detonationH2/H2_Li.xml
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/CanteraTorchProperties
similarity index 91%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties
rename to examples/dfHighSpeedFoam/oneD_detonationH2/constant/CanteraTorchProperties
index 3436a18a..640b9e14 100755
--- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/CanteraTorchProperties
@@ -21,8 +21,8 @@ CanteraMechanismFile "H2_Li.xml";
 transportModel "UnityLewis";
 loadbalancing
 {
-        active  false;
-        //log   true;
+        active  true;
+        log   true;
 }
 odeCoeffs
 {}
@@ -32,9 +32,4 @@ zeroDReactor
 {
     //cvorcp "UV";
 }
-loadbalancing
-{
-        active  true;
-        //log   true;
-}
 // ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_detonationH2/constant/dynamicMeshDict b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/dynamicMeshDict
new file mode 100644
index 00000000..cf838491
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/dynamicMeshDict
@@ -0,0 +1,68 @@
+/*--------------------------------*- 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      dynamicMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dynamicFvMesh   dynamicRefineFvMeshNew;
+
+multiCritRefinementControls
+{
+    enableMultiCritRefinementControl  false;
+}
+
+dynamicRefineFvMeshNewCoeffs
+{
+    // How often to refine
+    refineInterval  1; 
+
+    // Field to be refinement on
+    field           normalisedGradrho;
+
+    // Refine field in between lower..upper
+    lowerRefineLevel 0.2;
+    upperRefineLevel 1;
+
+    // If value < unrefineLevel unrefine
+    unrefineLevel   0.2;
+
+    // Have slower than 2:1 refinement
+    nBufferLayers   3;
+
+    // Refine cells only up to maxRefinement levels
+    maxRefinement   3;
+
+    // Stop refinement if maxCells reached
+    maxCells        1000000;
+
+    // Flux field and corresponding velocity field. Fluxes on changed
+    // faces get recalculated by interpolating the velocity. Use 'none'
+    // on surfaceScalarFields that do not need to be reinterpolated.
+    correctFluxes
+    (
+        (phi none)
+        //(nHatf none) 
+        //(rhoPhi none)
+        //(alphaPhi0.water none)
+        //(ghf none)
+        (neg none)
+        (pos none)
+    );
+
+    // Write the refinement level as a volScalarField
+    dumpLevel       true;
+}
+
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/thermophysicalProperties b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/thermophysicalProperties
similarity index 97%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/thermophysicalProperties
rename to examples/dfHighSpeedFoam/oneD_detonationH2/constant/thermophysicalProperties
index e983a193..873b8e16 100755
--- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/thermophysicalProperties
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/thermophysicalProperties
@@ -14,5 +14,5 @@ FoamFile
     object      thermophysicalProperties;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-inviscid      true;
+inviscid      false;
 // ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/turbulenceProperties b/examples/dfHighSpeedFoam/oneD_detonationH2/constant/turbulenceProperties
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/turbulenceProperties
rename to examples/dfHighSpeedFoam/oneD_detonationH2/constant/turbulenceProperties
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict b/examples/dfHighSpeedFoam/oneD_detonationH2/system/blockMeshDict
similarity index 95%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/blockMeshDict
index 1a73e764..921bd986 100755
--- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/system/blockMeshDict
@@ -30,7 +30,7 @@ vertices
 
 blocks
 (
-    hex (0 1 2 3 4 5 6 7) (5000 1 1) simpleGrading (1 1 1)
+    hex (0 1 2 3 4 5 6 7) (625 1 1) simpleGrading (1 1 1)
 );
 
 edges
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/controlDict b/examples/dfHighSpeedFoam/oneD_detonationH2/system/controlDict
similarity index 87%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/controlDict
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/controlDict
index 57ed5645..d897bd2d 100755
--- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/controlDict
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/system/controlDict
@@ -14,7 +14,10 @@ FoamFile
     object      controlDict;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-libs ("libfieldFunctionObjects.so");
+libs (
+      "libfieldFunctionObjects.so"
+      "libnewdynamicFvMesh.so"
+      );
 
 application     dfHighSpeedFoam;
 
@@ -24,7 +27,7 @@ startTime       0;
 
 stopAt          endTime;
 
-endTime         2.4e-4;
+endTime         2.2e-4;
 
 deltaT          1e-09;
 
@@ -48,7 +51,7 @@ runTimeModifiable true;
 
 adjustTimeStep  yes;
 
-maxCo           0.02;
+maxCo           0.1;
 
 maxDeltaT       1;
 
@@ -58,7 +61,7 @@ functions
     {
       type          fieldMinMax;
       libs          ("libfieldFunctionObjects.so");
-      writeControl  timeStep; //writeTime;
+      writeControl  writeTime; //timeStep;
       fields        (p);
     }
 }
diff --git a/examples/dfHighSpeedFoam/oneD_detonationH2/system/decomposeParDict b/examples/dfHighSpeedFoam/oneD_detonationH2/system/decomposeParDict
new file mode 100755
index 00000000..a9715975
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_detonationH2/system/decomposeParDict
@@ -0,0 +1,32 @@
+/*--------------------------------*- 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  4;
+
+method          scotch;
+
+constraints
+{
+    refinementHistory
+    {
+        //- Decompose cells such that all cell originating from single cell
+        //  end up on same processor
+        type    refinementHistory;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes b/examples/dfHighSpeedFoam/oneD_detonationH2/system/fvSchemes
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/fvSchemes
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSolution b/examples/dfHighSpeedFoam/oneD_detonationH2/system/fvSolution
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSolution
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/fvSolution
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/sample b/examples/dfHighSpeedFoam/oneD_detonationH2/system/sample
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/sample
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/sample
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/setFieldsDict b/examples/dfHighSpeedFoam/oneD_detonationH2/system/setFieldsDict
similarity index 100%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/setFieldsDict
rename to examples/dfHighSpeedFoam/oneD_detonationH2/system/setFieldsDict
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/AR b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/AR
new file mode 100755
index 00000000..4a7914aa
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/AR
@@ -0,0 +1,42 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      AR;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.8861; 
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/H2 b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/H2
new file mode 100755
index 00000000..6c9c8e4c
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/H2
@@ -0,0 +1,41 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      H2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.01266;
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/O2 b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/O2
new file mode 100755
index 00000000..4d6f7000
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/O2
@@ -0,0 +1,41 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      O2;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.101266;
+
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/T b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/T
new file mode 100755
index 00000000..eecff994
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/T
@@ -0,0 +1,42 @@
+/*--------------------------------*- 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       volScalarField;
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 0.913683;
+//internalField   uniform 300;
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            inletOutlet;
+        inletValue      uniform 748.472;
+        value           uniform 748.472;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/U b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/U
new file mode 100755
index 00000000..4d76526e
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/U
@@ -0,0 +1,42 @@
+/*--------------------------------*- 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       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    Left
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    
+    Right
+    {
+        type            inletOutlet;
+        inletValue      uniform (-487.34 0 0);
+        value           uniform (-487.34 0 0);
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/Ydefault b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/Ydefault
new file mode 100755
index 00000000..7e583b1e
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/Ydefault
@@ -0,0 +1,42 @@
+/*--------------------------------*- 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       volScalarField;
+    location    "0";
+    object      Ydefault;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0.0;
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            zeroGradient;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/p b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/p
new file mode 100755
index 00000000..2b26b604
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/0_orig/p
@@ -0,0 +1,40 @@
+/*--------------------------------*- 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       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 35594;
+
+boundaryField
+{
+    Left
+    {
+        type            zeroGradient;
+    }
+    
+    Right
+    {
+        type            fixedValue;
+        value           uniform 35594;
+    }
+
+    empty
+    {
+        type            empty;
+    }
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allclean b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allclean
new file mode 100755
index 00000000..e40b9408
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allclean
@@ -0,0 +1,13 @@
+#!/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
+echo "Cleaning postProcessing/"
+rm -r postProcessing
+echo "Cleaning 0/"
+rm -r 0
\ No newline at end of file
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allrun b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allrun
new file mode 100755
index 00000000..f9647173
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/Allrun
@@ -0,0 +1,13 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application=dfHighSpeedFoam
+
+cp -r 0_orig/ 0/
+runApplication blockMesh
+runApplication setFields
+runApplication decomposePar
+runApplication mpirun -np 5 $application -parallel
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/H2_AR.xml b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/H2_AR.xml
new file mode 100755
index 00000000..4945c984
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/H2_AR.xml
@@ -0,0 +1,531 @@
+<?xml version="1.0"?>
+<ctml>
+  <validate species="yes" reactions="yes"/>
+
+  <!-- phase gas     -->
+  <phase id="gas" dim="3">
+    <elementArray datasrc="elements.xml">H O Ar</elementArray>
+    <speciesArray datasrc="#species_data">H2    O2    H     O     OH    HO2   H2O2  H2O   AR</speciesArray>
+    <reactionArray datasrc="#reaction_data"/>
+    <state>
+      <temperature units="K">300.0</temperature>
+      <pressure units="Pa">101325.0</pressure>
+    </state>
+    <thermo model="IdealGas"/>
+    <kinetics model="GasKinetics"/>
+    <transport model="Mix"/>
+  </phase>
+
+  <!-- species definitions     -->
+  <speciesData id="species_data">
+
+    <!-- species H2    -->
+    <species name="H2">
+      <atomArray>H:2 </atomArray>
+      <note>TPIS78</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.344331120E+00,   7.980520750E-03,  -1.947815100E-05,   2.015720940E-08, 
+             -7.376117610E-12,  -9.179351730E+02,   6.830102380E-01</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             3.337279200E+00,  -4.940247310E-05,   4.994567780E-07,  -1.795663940E-10, 
+             2.002553760E-14,  -9.501589220E+02,  -3.205023310E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">linear</string>
+        <LJ_welldepth units="K">38.000</LJ_welldepth>
+        <LJ_diameter units="A">2.920</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.790</polarizability>
+        <rotRelax>280.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species O2    -->
+    <species name="O2">
+      <atomArray>O:2 </atomArray>
+      <note>TPIS89</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             3.782456360E+00,  -2.996734160E-03,   9.847302010E-06,  -9.681295090E-09, 
+             3.243728370E-12,  -1.063943560E+03,   3.657675730E+00</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             3.282537840E+00,   1.483087540E-03,  -7.579666690E-07,   2.094705550E-10, 
+             -2.167177940E-14,  -1.088457720E+03,   5.453231290E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">linear</string>
+        <LJ_welldepth units="K">107.400</LJ_welldepth>
+        <LJ_diameter units="A">3.458</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">1.600</polarizability>
+        <rotRelax>3.800</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species H    -->
+    <species name="H">
+      <atomArray>H:1 </atomArray>
+      <note>L7/88</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.500000000E+00,   7.053328190E-13,  -1.995919640E-15,   2.300816320E-18, 
+             -9.277323320E-22,   2.547365990E+04,  -4.466828530E-01</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.500000010E+00,  -2.308429730E-11,   1.615619480E-14,  -4.735152350E-18, 
+             4.981973570E-22,   2.547365990E+04,  -4.466829140E-01</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">atom</string>
+        <LJ_welldepth units="K">145.000</LJ_welldepth>
+        <LJ_diameter units="A">2.050</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>0.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species O    -->
+    <species name="O">
+      <atomArray>O:1 </atomArray>
+      <note>L1/90</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             3.168267100E+00,  -3.279318840E-03,   6.643063960E-06,  -6.128066240E-09, 
+             2.112659710E-12,   2.912225920E+04,   2.051933460E+00</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.569420780E+00,  -8.597411370E-05,   4.194845890E-08,  -1.001777990E-11, 
+             1.228336910E-15,   2.921757910E+04,   4.784338640E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">atom</string>
+        <LJ_welldepth units="K">80.000</LJ_welldepth>
+        <LJ_diameter units="A">2.750</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>0.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species OH    -->
+    <species name="OH">
+      <atomArray>H:1 O:1 </atomArray>
+      <note>S9/01</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.125305610E+00,  -3.225449390E-03,   6.527646910E-06,  -5.798536430E-09, 
+             2.062373790E-12,   3.381538120E+03,  -6.904329600E-01</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="6000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.864728860E+00,   1.056504480E-03,  -2.590827580E-07,   3.052186740E-11, 
+             -1.331958760E-15,   3.718857740E+03,   5.701640730E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">linear</string>
+        <LJ_welldepth units="K">80.000</LJ_welldepth>
+        <LJ_diameter units="A">2.750</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>0.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species HO2    -->
+    <species name="HO2">
+      <atomArray>H:1 O:2 </atomArray>
+      <note>L5/89</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.301798010E+00,  -4.749120510E-03,   2.115828910E-05,  -2.427638940E-08, 
+             9.292251240E-12,   2.948080400E+02,   3.716662450E+00</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.017210900E+00,   2.239820130E-03,  -6.336581500E-07,   1.142463700E-10, 
+             -1.079085350E-14,   1.118567130E+02,   3.785102150E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">nonlinear</string>
+        <LJ_welldepth units="K">107.400</LJ_welldepth>
+        <LJ_diameter units="A">3.458</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>1.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species H2O2    -->
+    <species name="H2O2">
+      <atomArray>H:2 O:2 </atomArray>
+      <note>L7/88</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.276112690E+00,  -5.428224170E-04,   1.673357010E-05,  -2.157708130E-08, 
+             8.624543630E-12,  -1.770258210E+04,   3.435050740E+00</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.165002850E+00,   4.908316940E-03,  -1.901392250E-06,   3.711859860E-10, 
+             -2.879083050E-14,  -1.786178770E+04,   2.916156620E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">nonlinear</string>
+        <LJ_welldepth units="K">107.400</LJ_welldepth>
+        <LJ_diameter units="A">3.458</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>3.800</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species H2O    -->
+    <species name="H2O">
+      <atomArray>H:2 O:1 </atomArray>
+      <note>L8/89</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             4.198640560E+00,  -2.036434100E-03,   6.520402110E-06,  -5.487970620E-09, 
+             1.771978170E-12,  -3.029372670E+04,  -8.490322080E-01</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="3500.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             3.033992490E+00,   2.176918040E-03,  -1.640725180E-07,  -9.704198700E-11, 
+             1.682009920E-14,  -3.000429710E+04,   4.966770100E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">nonlinear</string>
+        <LJ_welldepth units="K">572.400</LJ_welldepth>
+        <LJ_diameter units="A">2.605</LJ_diameter>
+        <dipoleMoment units="Debye">1.844</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>4.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+
+    <!-- species AR    -->
+    <species name="AR">
+      <atomArray>Ar:1 </atomArray>
+      <note>120186</note>
+      <thermo>
+        <NASA Tmin="200.0" Tmax="1000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.500000000E+00,   0.000000000E+00,   0.000000000E+00,   0.000000000E+00, 
+             0.000000000E+00,  -7.453750000E+02,   4.366000000E+00</floatArray>
+        </NASA>
+        <NASA Tmin="1000.0" Tmax="5000.0" P0="100000.0">
+           <floatArray size="7" name="coeffs">
+             2.500000000E+00,   0.000000000E+00,   0.000000000E+00,   0.000000000E+00, 
+             0.000000000E+00,  -7.453750000E+02,   4.366000000E+00</floatArray>
+        </NASA>
+      </thermo>
+      <transport model="gas_transport">
+        <string title="geometry">atom</string>
+        <LJ_welldepth units="K">136.500</LJ_welldepth>
+        <LJ_diameter units="A">3.330</LJ_diameter>
+        <dipoleMoment units="Debye">0.000</dipoleMoment>
+        <polarizability units="A3">0.000</polarizability>
+        <rotRelax>0.000</rotRelax>
+        <dispersion_coefficient units="A5">0.000</dispersion_coefficient>
+        <quadrupole_polarizability units="A5">0.000</quadrupole_polarizability>
+      </transport>
+    </species>
+  </speciesData>
+  <reactionData id="reaction_data">
+
+    <!-- reaction 0001    -->
+    <reaction id="0001" reversible="yes">
+      <equation>H2 + O2 [=] 2 OH</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.700000E+10</A>
+           <b>0.0</b>
+           <E units="cal/mol">47780.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H2:1.0 O2:1</reactants>
+      <products>OH:2.0</products>
+    </reaction>
+
+    <!-- reaction 0002    -->
+    <reaction id="0002" reversible="yes">
+      <equation>OH + H2 [=] H2O + H</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.170000E+06</A>
+           <b>1.3</b>
+           <E units="cal/mol">3626.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>OH:1.0 H2:1</reactants>
+      <products>H2O:1.0 H:1</products>
+    </reaction>
+
+    <!-- reaction 0003    -->
+    <reaction id="0003" reversible="yes">
+      <equation>H + O2 [=] OH + O</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>5.130000E+13</A>
+           <b>-0.816</b>
+           <E units="cal/mol">16507.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H:1.0 O2:1</reactants>
+      <products>OH:1.0 O:1</products>
+    </reaction>
+
+    <!-- reaction 0004    -->
+    <reaction id="0004" reversible="yes">
+      <equation>O + H2 [=] OH + H</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.800000E+07</A>
+           <b>1.0</b>
+           <E units="cal/mol">8826.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>O:1.0 H2:1</reactants>
+      <products>OH:1.0 H:1</products>
+    </reaction>
+
+    <!-- reaction 0005    -->
+    <reaction id="0005" reversible="yes" type="threeBody">
+      <equation>H + O2 + M [=] HO2 + M</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>2.100000E+12</A>
+           <b>-1.0</b>
+           <E units="cal/mol">0.000000</E>
+        </Arrhenius>
+        <efficiencies default="1.0">H2:3.3 H2O:21.0 O2:0.0</efficiencies>
+      </rateCoeff>
+      <reactants>H:1.0 O2:1</reactants>
+      <products>HO2:1.0</products>
+    </reaction>
+
+    <!-- reaction 0006    -->
+    <reaction id="0006" reversible="yes">
+      <equation>H + O2 + O2 [=] HO2 + O2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>6.700000E+13</A>
+           <b>-1.42</b>
+           <E units="cal/mol">0.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H:1.0 O2:2</reactants>
+      <products>HO2:1.0 O2:1</products>
+    </reaction>
+
+    <!-- reaction 0007    -->
+    <reaction id="0007" reversible="yes">
+      <equation>OH + HO2 [=] H2O + O2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>5.000000E+10</A>
+           <b>0.0</b>
+           <E units="cal/mol">1000.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>OH:1.0 HO2:1</reactants>
+      <products>H2O:1.0 O2:1</products>
+    </reaction>
+
+    <!-- reaction 0008    -->
+    <reaction id="0008" reversible="yes">
+      <equation>H + HO2 [=] 2 OH</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>2.500000E+11</A>
+           <b>0.0</b>
+           <E units="cal/mol">1900.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H:1.0 HO2:1</reactants>
+      <products>OH:2.0</products>
+    </reaction>
+
+    <!-- reaction 0009    -->
+    <reaction id="0009" reversible="yes">
+      <equation>O + HO2 [=] O2 + OH</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>4.800000E+10</A>
+           <b>0.0</b>
+           <E units="cal/mol">1000.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>O:1.0 HO2:1</reactants>
+      <products>O2:1.0 OH:1</products>
+    </reaction>
+
+    <!-- reaction 0010    -->
+    <reaction id="0010" reversible="yes">
+      <equation>2 OH [=] O + H2O</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>6.000000E+05</A>
+           <b>1.3</b>
+           <E units="cal/mol">0.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>OH:2.0</reactants>
+      <products>O:1.0 H2O:1</products>
+    </reaction>
+
+    <!-- reaction 0011    -->
+    <reaction id="0011" reversible="yes" type="threeBody">
+      <equation>H2 + M [=] H + H + M</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>2.230000E+09</A>
+           <b>0.5</b>
+           <E units="cal/mol">92600.000000</E>
+        </Arrhenius>
+        <efficiencies default="1.0">H:2.0 H2:3.0 H2O:6.0</efficiencies>
+      </rateCoeff>
+      <reactants>H2:1.0</reactants>
+      <products>H:2.0</products>
+    </reaction>
+
+    <!-- reaction 0012    -->
+    <reaction id="0012" reversible="yes" type="threeBody">
+      <equation>O2 + M [=] O + O + M</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.850000E+08</A>
+           <b>0.5</b>
+           <E units="cal/mol">95560.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>O2:1.0</reactants>
+      <products>O:2.0</products>
+    </reaction>
+
+    <!-- reaction 0013    -->
+    <reaction id="0013" reversible="yes" type="threeBody">
+      <equation>H + OH + M [=] H2O + M</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>7.500000E+17</A>
+           <b>-2.6</b>
+           <E units="cal/mol">0.000000</E>
+        </Arrhenius>
+        <efficiencies default="1.0">H2O:20.0</efficiencies>
+      </rateCoeff>
+      <reactants>H:1.0 OH:1</reactants>
+      <products>H2O:1.0</products>
+    </reaction>
+
+    <!-- reaction 0014    -->
+    <reaction id="0014" reversible="yes">
+      <equation>H + HO2 [=] H2 + O2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>2.500000E+10</A>
+           <b>0.0</b>
+           <E units="cal/mol">700.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H:1.0 HO2:1</reactants>
+      <products>H2:1.0 O2:1</products>
+    </reaction>
+
+    <!-- reaction 0015    -->
+    <reaction id="0015" reversible="yes">
+      <equation>HO2 + HO2 [=] H2O2 + O2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>2.000000E+09</A>
+           <b>0.0</b>
+           <E units="cal/mol">0.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>HO2:2.0</reactants>
+      <products>H2O2:1.0 O2:1</products>
+    </reaction>
+
+    <!-- reaction 0016    -->
+    <reaction id="0016" reversible="yes" type="threeBody">
+      <equation>H2O2 + M [=] OH + OH + M</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.300000E+14</A>
+           <b>0.0</b>
+           <E units="cal/mol">45500.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H2O2:1.0</reactants>
+      <products>OH:2.0</products>
+    </reaction>
+
+    <!-- reaction 0017    -->
+    <reaction id="0017" reversible="yes">
+      <equation>H2O2 + H [=] HO2 + H2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.600000E+09</A>
+           <b>0.0</b>
+           <E units="cal/mol">3800.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H2O2:1.0 H:1</reactants>
+      <products>HO2:1.0 H2:1</products>
+    </reaction>
+
+    <!-- reaction 0018    -->
+    <reaction id="0018" reversible="yes">
+      <equation>H2O2 + OH [=] H2O + HO2</equation>
+      <rateCoeff>
+        <Arrhenius>
+           <A>1.000000E+10</A>
+           <b>0.0</b>
+           <E units="cal/mol">1800.000000</E>
+        </Arrhenius>
+      </rateCoeff>
+      <reactants>H2O2:1.0 OH:1</reactants>
+      <products>H2O:1.0 HO2:1</products>
+    </reaction>
+  </reactionData>
+</ctml>
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/CanteraTorchProperties b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/CanteraTorchProperties
new file mode 100755
index 00000000..4640968d
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/CanteraTorchProperties
@@ -0,0 +1,35 @@
+/*--------------------------------*- 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      chemistryProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+chemistry           on;
+CanteraMechanismFile "H2_AR.xml";
+transportModel "UnityLewis";
+loadbalancing
+{
+        active  true;
+        log   true;
+}
+odeCoeffs
+{}
+inertSpecie        "AR";
+
+zeroDReactor
+{
+    //cvorcp "UV";
+}
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/thermophysicalProperties b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/thermophysicalProperties
new file mode 100755
index 00000000..873b8e16
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/thermophysicalProperties
@@ -0,0 +1,18 @@
+/*--------------------------------*- 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      thermophysicalProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+inviscid      false;
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/turbulenceProperties b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/constant/turbulenceProperties
new file mode 100755
index 00000000..ab70e94d
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/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/dfHighSpeedFoam/oneD_reactiveShockTube/system/blockMeshDict b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/blockMeshDict
new file mode 100755
index 00000000..c4112e99
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/blockMeshDict
@@ -0,0 +1,75 @@
+/*--------------------------------*- 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.1;
+
+vertices
+(
+    (0 0 0)
+    (1.2 0 0)
+    (1.2 1 0)
+    (0 1 0)
+    (0 0 1)
+    (1.2 0 1)
+    (1.2 1 1)
+    (0 1 1)
+);
+
+blocks
+(
+    hex (0 1 2 3 4 5 6 7) (2400 1 1) simpleGrading (1 1 1)
+);
+
+edges
+(
+);
+
+boundary
+(
+    Left
+    {
+        type patch;
+        faces
+        (
+            (0 4 7 3)
+        );
+    }
+    Right
+    {
+        type patch;
+        faces
+        (
+            (1 2 6 5)
+        );
+    }
+    empty
+    {
+        type empty;
+        faces
+        (
+            (0 1 5 4)
+            (5 6 7 4)
+            (3 7 6 2)
+            (0 3 2 1)
+        );
+    }
+);
+
+mergePatchPairs
+(
+);
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/controlDict b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/controlDict
new file mode 100755
index 00000000..7567b959
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/controlDict
@@ -0,0 +1,59 @@
+/*--------------------------------*- 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     dfHighSpeedFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         23e-5;
+
+deltaT          1e-08;
+
+writeControl    adjustableRunTime;
+
+writeInterval   5e-6;
+
+cycleWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression on;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable true;
+
+adjustTimeStep  yes;
+
+maxCo           0.1;
+
+maxDeltaT       1;
+
+functions
+{
+    #includeFunc mag(U)
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/decomposeParDict b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/decomposeParDict
similarity index 96%
rename from examples/dfHighSpeedFoam/det_1D_H2_01mm/system/decomposeParDict
rename to examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/decomposeParDict
index 470777a2..dc84190c 100755
--- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/decomposeParDict
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/decomposeParDict
@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-numberOfSubdomains  4;
+numberOfSubdomains  5;
 
 method          scotch;
 
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSchemes b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSchemes
new file mode 100755
index 00000000..b8df2960
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSchemes
@@ -0,0 +1,58 @@
+/*--------------------------------*- 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;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+fluxScheme		Kurganov;
+
+ddtSchemes
+{
+    default               Euler;
+}
+
+gradSchemes
+{
+    default               Gauss linear;
+}
+
+divSchemes
+{
+    default               none;
+    div(tauMC)            Gauss linear;
+    div(hDiffCorrFlux)    Gauss cubic;
+    div(phi,Yi_h)         Gauss vanLeer;
+}
+
+laplacianSchemes
+{
+    default               Gauss linear uncorrected;	
+}
+
+interpolationSchemes
+{
+    default               linear;
+    reconstruct(rho)      Minmod;
+    reconstruct(U)        MinmodV;
+    reconstruct(T)        Minmod;
+    reconstruct(Yi)       Minmod;
+}
+
+snGradSchemes
+{
+    default               uncorrected;
+}
+
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSolution b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSolution
new file mode 100755
index 00000000..17c77b5d
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/fvSolution
@@ -0,0 +1,67 @@
+/*--------------------------------*- 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
+{
+	U
+	{
+		//solver		smoothSolver;
+		//smoother	GaussSeidel;
+		//nSweeps		2;
+		//tolerance	1e-17;
+		//relTol		0;
+
+		solver			PBiCGStab;
+		preconditioner	DIC;
+		tolerance		1e-11;
+		relTol			0;
+	}
+	
+	h
+	{
+		$U;
+		tolerance	1e-11;
+		relTol		0;
+	}
+
+	e
+	{
+		$U;
+		tolerance	1e-11;
+		relTol		0;
+	}
+
+	"rho.*"
+	{
+		solver		diagonal;
+	}
+
+	"(O2|N2|Yi)"
+	{
+		solver          PBiCGStab;
+        preconditioner	DILU;
+        tolerance       1e-11;
+        relTol          0;
+	}
+	
+}
+
+CENTRAL
+{
+}
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/sample b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/sample
new file mode 100755
index 00000000..f5ac29bb
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/sample
@@ -0,0 +1,39 @@
+/*--------------------------------*- 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      sample;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+type sets;
+libs            ("libsampling.so");
+
+interpolationScheme cellPoint;
+
+setFormat       raw;
+
+sets
+(
+    data
+    {
+        type    lineFace;
+        axis    x;
+        start   (-4.995 0 0);
+        end     (4.995 0 0);
+        nPoints 1000;
+    }
+);
+
+fields          (T mag(U) p);
+
+// ************************************************************************* //
diff --git a/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/setFieldsDict b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/setFieldsDict
new file mode 100755
index 00000000..4f686528
--- /dev/null
+++ b/examples/dfHighSpeedFoam/oneD_reactiveShockTube/system/setFieldsDict
@@ -0,0 +1,40 @@
+/*--------------------------------*- 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
+(
+    volVectorFieldValue U (0 0 0)
+    volScalarFieldValue T 378.656
+    volScalarFieldValue p 7173
+);
+
+regions
+(
+    boxToCell
+    {
+        box (0.06 -0.1 -0.1) (0.12 0.1 0.1);
+        fieldValues
+        (
+            volVectorFieldValue U (-487.34 0 0)
+            volScalarFieldValue T 748.472
+            volScalarFieldValue p 35594
+        );
+    }
+);
+
+
+// ************************************************************************* //