BUG: corrected pressure non-ortho loops when using correction(ddt(p))

applications/solvers/combustion/rhoReactingFoam/pEqn.H

@@ -24,12 +24,17 @@
             fvc::interpolate(thermo.psi())*phiv
         );
 
+        fvScalarMatrix pDDtEqn
+        (
+            fvc::ddt(rho) + fvc::div(phi)
+          + correction(fvm::ddt(psi, p) + fvm::div(phid, p))
+        );
+
         for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
         {
             fvScalarMatrix pEqn
             (
-                fvc::ddt(rho) + fvc::div(phi)
-              + correction(fvm::ddt(psi, p) + fvm::div(phid, p))
+                pDDtEqn
               - fvm::laplacian(rho*rAU, p)
             );
 
@@ -63,12 +68,17 @@
               + fvc::ddtPhiCorr(rAU, rho, U, phi)
             );
 
+        fvScalarMatrix pDDtEqn
+        (
+            fvc::ddt(rho) + psi*correction(fvm::ddt(p))
+          + fvc::div(phi)
+        );
+
         for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
         {
             fvScalarMatrix pEqn
             (
-                fvc::ddt(rho) + psi*correction(fvm::ddt(p))
-              + fvc::div(phi)
+                pDDtEqn
               - fvm::laplacian(rho*rAU, p)
             );
 

applications/solvers/heatTransfer/buoyantPimpleFoam/pEqn.H

@@ -19,12 +19,17 @@
     surfaceScalarField buoyancyPhi(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf());
     phi += buoyancyPhi;
 
+    fvScalarMatrix p_rghDDtEqn
+    (
+        fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh))
+      + fvc::div(phi)
+    );
+
     for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     {
         fvScalarMatrix p_rghEqn
         (
-            fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh))
-          + fvc::div(phi)
+            p_rghDDtEqn
           - fvm::laplacian(rhorAUf, p_rgh)
         );
 

applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H

@@ -23,7 +23,7 @@
             );
     }
 
-    fvScalarMatrix DrhoDtEqn
+    fvScalarMatrix pDDtEqn
     (
         fvc::ddt(rho) + psi*correction(fvm::ddt(p))
       + fvc::div(phi)
@@ -36,7 +36,7 @@
     {
         fvScalarMatrix pEqn
         (
-            DrhoDtEqn
+            pDDtEqn
           - fvm::laplacian(rho*rAU, p)
         );
 

applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/pEqn.H

@@ -23,16 +23,21 @@
             );
     }
 
+    fvScalarMatrix pDDtEqn
+    (
+        fvc::ddt(rho) + psi*correction(fvm::ddt(p))
+      + fvc::div(phi)
+     ==
+        parcels.Srho()
+      + massSource.SuTot()
+    );
+
     for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     {
         fvScalarMatrix pEqn
         (
-            fvc::ddt(rho) + psi*correction(fvm::ddt(p))
-          + fvc::div(phi)
+            pDDtEqn
           - fvm::laplacian(rho*rAU, p)
-        ==
-            parcels.Srho()
-          + massSource.SuTot()
         );
 
         if (corr == nCorr-1 && nonOrth == nNonOrthCorr)