Wall-functions: Corrected conversion of muw->nuw
This commit is contained in:
parent
a5f0ea6867
commit
187f928ff6
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -92,7 +92,7 @@ void epsilonLowReWallFunctionFvPatchScalarField::calculate
|
||||
{
|
||||
label cellI = patch.faceCells()[faceI];
|
||||
|
||||
scalar yPlus = Cmu25*sqrt(k[cellI])*y[faceI]/muw[faceI]/rhow[faceI];
|
||||
scalar yPlus = Cmu25*sqrt(k[cellI])*y[faceI]/(muw[faceI]/rhow[faceI]);
|
||||
|
||||
scalar w = cornerWeights[faceI];
|
||||
|
||||
|
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -198,7 +198,7 @@ void fWallFunctionFvPatchScalarField::updateCoeffs()
|
||||
|
||||
scalar uTau = Cmu25*sqrt(k[faceCellI]);
|
||||
|
||||
scalar yPlus = uTau*y[faceI]/muw[faceI]/rhow[faceI];
|
||||
scalar yPlus = uTau*y[faceI]/(muw[faceI]/rhow[faceI]);
|
||||
|
||||
if (yPlus > yPlusLam_)
|
||||
{
|
||||
|
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -195,7 +195,7 @@ void kLowReWallFunctionFvPatchScalarField::updateCoeffs()
|
||||
|
||||
scalar uTau = Cmu25*sqrt(k[faceCellI]);
|
||||
|
||||
scalar yPlus = uTau*y[faceI]/muw[faceI]/rhow[faceI];
|
||||
scalar yPlus = uTau*y[faceI]/(muw[faceI]/rhow[faceI]);
|
||||
|
||||
if (yPlus > yPlusLam_)
|
||||
{
|
||||
|
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -190,7 +190,7 @@ void v2WallFunctionFvPatchScalarField::updateCoeffs()
|
||||
|
||||
scalar uTau = Cmu25*sqrt(k[faceCellI]);
|
||||
|
||||
scalar yPlus = uTau*y[faceI]/muw[faceI]/rhow[faceI];
|
||||
scalar yPlus = uTau*y[faceI]/(muw[faceI]/rhow[faceI]);
|
||||
|
||||
if (yPlus > yPlusLam_)
|
||||
{
|
||||
|
Loading…
Reference in New Issue
Block a user