From f947f0b46e55938bbdfdad0ce8b9f8bc9fdfa99b Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 15 Apr 2011 16:09:39 +0100
Subject: [PATCH] bubbleFoam: Updated the viscosity used to the effective
viscosity
in line with the rest of the wall-function implementations
---
applications/solvers/multiphase/bubbleFoam/wallFunctions.H | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/applications/solvers/multiphase/bubbleFoam/wallFunctions.H b/applications/solvers/multiphase/bubbleFoam/wallFunctions.H
index bc9e93c63df..0d047468408 100644
--- a/applications/solvers/multiphase/bubbleFoam/wallFunctions.H
+++ b/applications/solvers/multiphase/bubbleFoam/wallFunctions.H
@@ -4,6 +4,7 @@
scalar Cmu25 = ::pow(Cmu.value(), 0.25);
scalar Cmu75 = ::pow(Cmu.value(), 0.75);
scalar kappa_ = kappa.value();
+ scalar nub_ = nub.value();
const fvPatchList& patches = mesh.boundary();
@@ -32,7 +33,7 @@
if (isA<wallFvPatch>(currPatch))
{
- const scalarField& nuw = nutb.boundaryField()[patchi];
+ const scalarField& nutbw = nutb.boundaryField()[patchi];
scalarField magFaceGradU(mag(Ub.boundaryField()[patchi].snGrad()));
@@ -43,7 +44,7 @@
scalar yPlus =
Cmu25*y[patchi][facei]
*::sqrt(k[faceCelli])
- /nub.value();
+ /nub_;
// For corner cells (with two boundary or more faces),
@@ -59,7 +60,7 @@
if (yPlus > 11.6)
{
G[faceCelli] +=
- nuw[facei]*magFaceGradU[facei]
+ (nutbw[facei] + nub_)*magFaceGradU[facei]
*Cmu25*::sqrt(k[faceCelli])
/(kappa_*y[patchi][facei]);
}
--
GitLab