BUG: Corrected kappat Jayatilleke wall function

This commit is contained in:
andy
2010-11-23 15:37:00 +00:00
parent 49331bf917
commit d21d66ef8c

View File

@ -235,7 +235,7 @@ void kappatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs()
if (yPlus > yPlusTherm)
{
scalar nu = nuw[faceI];
scalar kt = nu*(yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P) - 1/Pr);
scalar kt = nu*(yPlus/(Prt_*(log(E_*yPlus)/kappa_ + P)) - 1/Pr)
kappatw[faceI] = max(0.0, kt);
}
else