diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
index dc1c09e7440fd597a996da4fe09e0149877ba030..9cb7908a1479660ed85b7b2412006360a637258a 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
@@ -93,6 +93,11 @@ int main(int argc, char *argv[])
             bool finalIter = oCorr == nOuterCorr-1;
             #include "addFinalIter.H"
 
+            if (nOuterCorr != 1)
+            {
+                p.storePrevIter();
+            }
+
             turbulence->correct();
 
             #include "UEqn.H"
diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
index c25c63df0b3c314414a6a3a85121a7ba0b22c8f6..97ba0c5c781a868bed51173f57420b1bfbb493c2 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
@@ -59,7 +59,7 @@
         }
     }
 
-    Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
+    p.relax();
 
     // Second part of thermodynamic density update
     thermo.rho() += psi*p;
@@ -75,4 +75,6 @@
     rho = min(rho, rhoMax);
 
     #include "setPressureWork.H"
+
+    Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
 }