Update rho = thermo.rho(), in the pEq for reactingFoam.

Updating fvSolution's for closed domains for chtMultiRegionFoam cases

applications/solvers/combustion/reactingFoam/pEqn.H

@@ -92,9 +92,10 @@ K = 0.5*magSqr(U);
 if (pressureControl.limit(p))
 {
     p.correctBoundaryConditions();
-    rho = thermo.rho();
 }
 
+rho = thermo.rho();
+
 if (thermo.dpdt())
 {
     dpdt = fvc::ddt(p);

tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/fvSolution

@@ -62,6 +62,8 @@ PIMPLE
     momentumPredictor   yes;
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors

tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/fvSolution

@@ -62,6 +62,8 @@ PIMPLE
     momentumPredictor   yes;
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors

tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/cabin/fvSolution

@@ -66,6 +66,8 @@ PIMPLE
     momentumPredictor no;
     nCorrectors     1;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors

tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/ice/fvSolution

@@ -82,6 +82,8 @@ PIMPLE
     momentumPredictor no;
     nCorrectors     1;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors