diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict
index d1fc63c252..3719b30a03 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict
@@ -51,8 +51,6 @@ maxCo           0.5;
 
 maxDeltaT       1;
 
-alphaTauSpecie  1e10;
-
 alphaTauTemp    0.005;
 
 // ************************************************************************* //
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes
index 45f26de8bb..3f08dd0d6e 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes
@@ -17,7 +17,7 @@ FoamFile
 
 ddtSchemes
 {
-    default         localEuler invTauFlow;
+    default         localEuler invTau;
 }
 
 gradSchemes
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution
index 31c2340be4..e5a074276a 100755
--- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution
@@ -24,6 +24,13 @@ solvers
         tolerance       0;
         relTol          0.1;
     }
+    rhoFinal
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-4;
+        relTol          0;
+    }
     p
     {
         solver           GAMG;
@@ -44,6 +51,14 @@ solvers
         maxIter          50;
     };
 
+    pFinal
+    {
+        $p;
+
+        tolerance        1e-4;
+        relTol           0;
+    };
+
     "(U|Yi|hs|k|epsilon)"
     {
         solver          smoothSolver;
@@ -51,9 +66,17 @@ solvers
         tolerance       0;
         relTol          0.1;
     }
+
+    "(U|Yi|hs|k|epsilon)Final"
+    {
+        solver          smoothSolver;
+        smoother        GaussSeidel;
+        tolerance       1e-4;
+        relTol          0;
+    }
 }
 
-PISO
+PIMPLE
 {
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
@@ -64,8 +87,8 @@ PISO
 
 additional
 {
-    eWork           true;
-    hWork           true;
+    pressureWork    true;
+    pressureWorkTimeDerivative true;
     solveSpecies    true;
 }
 
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties
index b415ea939b..2215e348f7 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties
@@ -122,7 +122,7 @@ subModels
         massFlowRate    0.8e-03;
         parcelBasisType mass;
         patchName       inletCentral;
-        parcelsPerSecond 400;
+        parcelsPerSecond 100;
         duration        1; // NOTE: set to 1 for steady state
         U0              (0 40 0);
         flowRateProfile constant 1;
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict
index 7bf404c0a6..8569d05e8d 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict
@@ -45,17 +45,13 @@ timePrecision   6;
 
 runTimeModifiable yes;
 
-maxCo               5;
+maxCo           5;
 
-alphaTauSmooth      1e10;
+alphaTauSmooth  1;
 
-alphaTauTemp        1e10;
+alphaTauTemp    1;
 
-alphaTauRho         1e10;
-
-alphaTauSpecie      1e10;
-
-maxDeltaT           1e10;
+maxDeltaT       1;
 
 functions
 {
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes
index da93c451d2..0a99ca6ddf 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes
@@ -17,7 +17,7 @@ FoamFile
 
 ddtSchemes
 {
-    default         localEuler invTauFlow;
+    default         localEuler invTau;
 }
 
 gradSchemes
diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution
index 1d4e95dae0..7184f4b370 100644
--- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution
+++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution
@@ -1,14 +1,14 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  dev                                   |
+|  \\    /   O peration     | Version:  1.6                                   |
 |   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
     version     2.0;
-    format      ascii;
+    format      binary;
     class       dictionary;
     location    "system";
     object      fvSolution;
@@ -17,11 +17,21 @@ FoamFile
 
 solvers
 {
+    "(rho|.*Final)"
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       0;
+        relTol          0.1;
+    }
+
+//    "(p|.*Final)"
     p
     {
         solver           GAMG;
+
         tolerance        0;
-        relTol           0.1;
+        relTol           0.05;
 
         smoother         DICGaussSeidel;
         nPreSweeps       0;
@@ -33,47 +43,39 @@ solvers
         agglomerator     faceAreaPair;
         mergeLevels      1;
 
-        maxIter          20;
+        maxIter          50;
     };
-    "(rho|G)"
+
+    pFinal
     {
-        solver           PCG;
-        preconditioner   DIC;
-        tolerance        0;
-        relTol           0.1;
-        maxIter          20;
-    };
-    "(Yi|hs)"
+        $p;
+        tolerance        1e-4;
+    }
+
+    "(U|Yi|hs|k|omega|.*Final)"
     {
-        solver           smoothSolver;
-        smoother         GaussSeidel;
-        tolerance        0;
-        relTol           0.1;
-        maxIter          20;
-    };
-    "(U|k|omega)"
-    {
-        solver           smoothSolver;
-        smoother         GaussSeidel;
-        tolerance        0;
-        relTol           0.1;
-        maxIter          20;
-    };
+        solver          smoothSolver;
+        smoother        GaussSeidel;
+        tolerance       0;
+        relTol          0.1;
+        maxIter         20;
+    }
 }
 
-PISO
+PIMPLE
 {
-    nCorrectors          2;
-    nNonOrthogonalCorrectors 0; // 10;
+    nOuterCorrectors    1;
+    nCorrectors         2;
+    nNonOrthogonalCorrectors 0;
     momentumPredictor yes;
-    rhoMin rhoMin [ 1 -3 0 0 0 ] 0; // 0.5;
-    rhoMax rhoMax [ 1 -3 0 0 0 ] 100; // 1.5;
+    rhoMin rhoMin [1 -3 0 0 0] 0.1;
+    rhoMax rhoMax [1 -3 0 0 0] 1.5;
 }
 
 additional
 {
-    eWork           true;
-    hWork           true;
+    pressureWork    true;
+    pressureWorkTimeDerivative false; // true;
     solveSpecies    true;
 }