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

dynOneEqEddy: Bound KK during construction

Browse source 
Resolves bug-report http://openfoam.org/mantisbt/view.php?id=1850
This commit is contained in:
Henry Weller 2015-09-09 19:45:00 +01:00
parent 51457432ee
commit e040f431ba
1 changed files with 3 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -133,7 +133,8 @@ dynOneEqEddy::dynOneEqEddy
{
bound(k_, kMin_);
const volScalarField KK(0.5*(filter_(magSqr(U)) - magSqr(filter_(U))));
volScalarField KK(0.5*(filter_(magSqr(U)) - magSqr(filter_(U))));
KK.max(dimensionedScalar("small", KK.dimensions(), SMALL));
updateSubGridScaleFields(symm(fvc::grad(U)), KK);
printCoeffs();