diff --git a/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H b/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H
index b75db346f5..632432ff8c 100644
--- a/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H
+++ b/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H
@@ -13,6 +13,7 @@
         {
             scalarField& mutw = mut.boundaryField()[patchi];
             const scalarField& mulw = mul.boundaryField()[patchi];
+            const scalarField& rhow = rho.boundaryField()[patchi];
 
             forAll(curPatch, facei)
             {
@@ -20,7 +21,7 @@
 
                 scalar yPlus =
                     Cmu25*y[patchi][facei]*::sqrt(k[faceCelli])
-                   /(mulw[facei]/rho[facei]);
+                   /(mulw[facei]/rhow[facei]);
 
                 if (yPlus > 11.6)
                 {