ENH: (re-)added under-relaxation to pEqn in LTSReactingParcelFoam solver

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"

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;
 }