diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/createFluidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/createFluidFields.H
index 885b1f080ede729a943812d13bc89b0f86e65c9b..6a0d52b1c25e3bb6e7bc6318e85c0e3184d6c7c8 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/createFluidFields.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/createFluidFields.H
@@ -175,8 +175,35 @@
 
         simpleDict.readIfPresent("frozenFlow", frozenFlowFluid[i]);
 
-        rhoMax.set(i, new dimensionedScalar(simpleDict.lookup("rhoMax")));
-        rhoMin.set(i, new dimensionedScalar(simpleDict.lookup("rhoMin")));
+        rhoMax.set
+        (
+            i,
+            new dimensionedScalar
+            (
+                dimensionedScalar::lookupOrDefault
+                (
+                    "rhoMax",
+                    simpleDict,
+                    GREAT,
+                    dimDensity
+                )
+            )
+        );
+
+        rhoMin.set
+        (
+            i,
+            new dimensionedScalar
+            (
+                dimensionedScalar::lookupOrDefault
+                (
+                    "rhoMin",
+                    simpleDict,
+                    0,
+                    dimDensity
+                )
+            )
+        );
 
         Info<< "    Adding fvOptions\n" << endl;
         fluidFvOptions.set
@@ -185,5 +212,3 @@
             new fv::IOoptionList(fluidRegions[i])
         );
     }
-
-