diff --git a/applications/solvers/multiphase/compressibleInterFoam/alphaEqnsSubCycle.H b/applications/solvers/multiphase/compressibleInterFoam/alphaEqnsSubCycle.H
index 5a74d1c456a70089f94697a4ad88f7a088881313..89e56191503730cf9634ba0d7d72169b0aac8634 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/alphaEqnsSubCycle.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/alphaEqnsSubCycle.H
@@ -18,7 +18,7 @@
                 mesh
             ),
             mesh,
-            dimensionedScalar("0", dimDensity, 0)
+            dimensionedScalar("0", rhoPhi.dimensions(), 0)
         );
 
         for
diff --git a/applications/solvers/multiphase/interFoam/LTSInterFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/interFoam/LTSInterFoam/alphaEqnSubCycle.H
index ac40a9bc9e1d163440129b15fc4eb53d49e85cf6..81c65c1caf737480956f72ff2d3952f7311436c5 100644
--- a/applications/solvers/multiphase/interFoam/LTSInterFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interFoam/LTSInterFoam/alphaEqnSubCycle.H
@@ -12,7 +12,7 @@ if (nAlphaSubCycles > 1)
             mesh
         ),
         mesh,
-        dimensionedScalar("0", dimDensity, 0)
+        dimensionedScalar("0", rhoPhi.dimensions(), 0)
     );
 
     for
diff --git a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
index ac40a9bc9e1d163440129b15fc4eb53d49e85cf6..81c65c1caf737480956f72ff2d3952f7311436c5 100644
--- a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
@@ -12,7 +12,7 @@ if (nAlphaSubCycles > 1)
             mesh
         ),
         mesh,
-        dimensionedScalar("0", dimDensity, 0)
+        dimensionedScalar("0", rhoPhi.dimensions(), 0)
     );
 
     for
diff --git a/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqnsSubCycle.H b/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqnsSubCycle.H
index 9682b30212c9af69eea682bb4f4c5e93400995ce..259ebb9cef6c986ee2bd2861ac1ff997a7dfb6ac 100644
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqnsSubCycle.H
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqnsSubCycle.H
@@ -12,7 +12,7 @@ if (nAlphaSubCycles > 1)
             mesh
         ),
         mesh,
-        dimensionedScalar("0", dimDensity, 0)
+        dimensionedScalar("0", rhoPhi.dimensions(), 0)
     );
 
     for
diff --git a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C
index 31e326793bfb8e1b3649e5293ef772605239ba2d..2a5d9a03b0a14c2576ef9485675fab3310eb731d 100644
--- a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C
+++ b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C
@@ -322,7 +322,7 @@ void Foam::multiphaseMixture::solve()
                 mesh_
             ),
             mesh_,
-            dimensionedScalar("0", dimDensity, 0)
+            dimensionedScalar("0", rhoPhi_.dimensions(), 0)
         );
 
         dimensionedScalar totalDeltaT = runTime.deltaT();
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqnSubCycle.H
index ac40a9bc9e1d163440129b15fc4eb53d49e85cf6..81c65c1caf737480956f72ff2d3952f7311436c5 100644
--- a/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqnSubCycle.H
@@ -12,7 +12,7 @@ if (nAlphaSubCycles > 1)
             mesh
         ),
         mesh,
-        dimensionedScalar("0", dimDensity, 0)
+        dimensionedScalar("0", rhoPhi.dimensions(), 0)
     );
 
     for