diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H
index b543c07dde..779b857294 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H
@@ -4,7 +4,7 @@ if (pimple.frozenFlow())
 }
 else
 {
-    if (!mesh.steady() && pimples.nCorrPimple() <= 1)
+    if (!mesh.steady() && pimples.firstPimpleIter())
     {
         #include "rhoEqn.H"
     }
diff --git a/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C b/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
index 9b81decdc5..0cebe103dd 100644
--- a/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
+++ b/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
@@ -94,10 +94,7 @@ int main(int argc, char *argv[])
 
         if (solvePrimaryRegion)
         {
-            if (pimple.nCorrPimple() <= 1)
-            {
-                #include "rhoEqn.H"
-            }
+            #include "rhoEqn.H"
 
             // --- PIMPLE loop
             while (pimple.loop())