From 57b968044ed23ec904121e74e4537a2ba13d0fe3 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 26 Mar 2014 15:19:39 +0000
Subject: [PATCH] Corrected rho used for yPlus
---
applications/solvers/multiphase/driftFluxFoam/wallViscosity.H | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H b/applications/solvers/multiphase/driftFluxFoam/wallViscosity.H
index b75db346f50..632432ff8c1 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)
{
--
GitLab