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

kLowReWallFunction: Limit k to avoid /0

Browse source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1724
This commit is contained in:
Henry 2015-05-31 21:18:34 +01:00
parent 25d7706560
commit 9b0d7415bb
3 changed files with 9 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
src/TurbulenceModels/turbulenceModels/RAS/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

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

@ -213,6 +213,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

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

@ -211,6 +211,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