research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

kLowReWallFunction: Limit k to avoid failure of the turbulence model due to division by k

Browse source
This commit is contained in:
Henry 2015-02-04 22:11:15 +00:00
parent 287e970b18
commit c481b625ef
1 changed files with 3 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C
View file

@ -215,6 +215,9 @@ void kLowReWallFunctionFvPatchScalarField::updateCoeffs()
kw[faceI] *= sqr(uTau);
}
// Limit kw to avoid failure of the turbulence model due to division by kw
kw = max(kw, SMALL);
fixedValueFvPatchField<scalar>::updateCoeffs();
// TODO: perform averaging for cells sharing more than one boundary face