zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

bubbleFoam: Updated the viscosity used to the effective viscosity

Browse Source 
in line with the rest of the wall-function implementations
This commit is contained in:
Henry
2011-04-15 16:09:39 +01:00
parent 022df6cffd
commit f947f0b46e
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
applications/solvers/multiphase/bubbleFoam/wallFunctions.H
View File

@ -4,6 +4,7 @@
scalar Cmu25 = ::pow(Cmu.value(), 0.25); scalar Cmu25 = ::pow(Cmu.value(), 0.25);
scalar Cmu75 = ::pow(Cmu.value(), 0.75); scalar Cmu75 = ::pow(Cmu.value(), 0.75);
scalar kappa_ = kappa.value(); scalar kappa_ = kappa.value();
scalar nub_ = nub.value();
const fvPatchList& patches = mesh.boundary(); const fvPatchList& patches = mesh.boundary();
@ -32,7 +33,7 @@
if (isA<wallFvPatch>(currPatch)) if (isA<wallFvPatch>(currPatch))
{ {
const scalarField& nuw = nutb.boundaryField()[patchi]; const scalarField& nutbw = nutb.boundaryField()[patchi];
scalarField magFaceGradU(mag(Ub.boundaryField()[patchi].snGrad())); scalarField magFaceGradU(mag(Ub.boundaryField()[patchi].snGrad()));
@ -43,7 +44,7 @@
scalar yPlus = scalar yPlus =
Cmu25*y[patchi][facei] Cmu25*y[patchi][facei]
*::sqrt(k[faceCelli]) *::sqrt(k[faceCelli])
/nub.value(); /nub_;
// For corner cells (with two boundary or more faces), // For corner cells (with two boundary or more faces),
@ -59,7 +60,7 @@
if (yPlus > 11.6) if (yPlus > 11.6)
{ {
G[faceCelli] += G[faceCelli] +=
nuw[facei]*magFaceGradU[facei] (nutbw[facei] + nub_)*magFaceGradU[facei]
*Cmu25*::sqrt(k[faceCelli]) *Cmu25*::sqrt(k[faceCelli])
/(kappa_*y[patchi][facei]); /(kappa_*y[patchi][facei]);
} }