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

BUG: corrected laminar clause in kappt wall function

Browse Source
This commit is contained in:
andy
2010-10-06 17:41:40 +01:00
parent f2c3ba65a6
commit aadacc5812
1 changed files with 6 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C
View File

@ -231,19 +231,17 @@ void kappatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs()
scalar P = Psmooth(Prat); scalar P = Psmooth(Prat);
scalar yPlusTherm = this->yPlusTherm(P, Prat); scalar yPlusTherm = this->yPlusTherm(P, Prat);
// Evaluate new effective thermal diffusivity // Update turbulent thermal conductivity
scalar kappaEff = 0.0; if (yPlus > yPlusTherm)
if (yPlus < yPlusTherm)
{ {
kappaEff = Pr*yPlus; scalar nu = nuw[faceI];
scalar kt = nu*(yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P) - 1/Pr);
kappatw[faceI] = max(0.0, kt);
} }
else else
{ {
kappaEff = nuw[faceI]*yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P); kappatw[faceI] = 0.0;
} }
// Update turbulent thermal diffusivity
kappatw[faceI] = max(0.0, kappaEff - nuw[faceI]/Pr);
} }
fixedValueFvPatchField<scalar>::updateCoeffs(); fixedValueFvPatchField<scalar>::updateCoeffs();