diff --git a/applications/solvers/combustion/reactingFoam/pEqn.H b/applications/solvers/combustion/reactingFoam/pEqn.H
index 81ab6424b4c1a7bfac43df5538ff4d089da5dca0..36f3e777530b11439c4fda6bc1c2730a32d1f19c 100644
--- a/applications/solvers/combustion/reactingFoam/pEqn.H
+++ b/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);
diff --git a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
index e56035f24a703420b00d2e539ee89f3a7af56090..3ad7cb746412436db03a7d75885c8fd559157828 100644
--- a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
@@ -113,3 +113,26 @@ volScalarField Qdot
 );
 
 #include "createMRF.H"
+
+dimensionedScalar rhoMax
+(
+    dimensionedScalar::lookupOrDefault
+    (
+        "rhoMax",
+        pimple.dict(),
+        dimDensity,
+        GREAT
+    )
+);
+
+dimensionedScalar rhoMin
+(
+    dimensionedScalar::lookupOrDefault
+    (
+        "rhoMin",
+        pimple.dict(),
+        dimDensity,
+        0
+    )
+);
+
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/fvSolution
index 1185c9d7c422484dd5a9fd14259d2a89820bab53..51011f702cbd0ad8e0737fa9e0d7f83230a23ae6 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/fvSolution
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/fvSolution
@@ -62,6 +62,8 @@ PIMPLE
     momentumPredictor   yes;
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/fvSolution
index 1185c9d7c422484dd5a9fd14259d2a89820bab53..51011f702cbd0ad8e0737fa9e0d7f83230a23ae6 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/fvSolution
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/fvSolution
@@ -62,6 +62,8 @@ PIMPLE
     momentumPredictor   yes;
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/cabin/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/cabin/fvSolution
index 0eabd4031a504817596bd53e25826e26cd44c36b..6dc2762c5a57d2e676482d8fbbef193dbcbc2814 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/cabin/fvSolution
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/cabin/fvSolution
@@ -66,6 +66,8 @@ PIMPLE
     momentumPredictor no;
     nCorrectors     1;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/ice/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/ice/fvSolution
index 0c2e7c9b8b95069b40121ccfc76b1d209bbf8ce9..afd8468ca6538fb2f1cbdef77c5707934fb9f84b 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/ice/fvSolution
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/windshieldDefrost/system/ice/fvSolution
@@ -82,6 +82,8 @@ PIMPLE
     momentumPredictor no;
     nCorrectors     1;
     nNonOrthogonalCorrectors 0;
+    pRefCell    0;
+    pRefValue   1e5;
 }
 
 relaxationFactors