Corrected rho used for yPlus

57b96804 · Henry · Andrew Heather · fcbca30a · 57b96804
Commit 57b96804 authored 11 years ago by Henry Committed by Andrew Heather 11 years ago
--- 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)
                {