From e040f431bac28ce3aac3eff31c918a8c98135231 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Wed, 9 Sep 2015 19:45:00 +0100
Subject: [PATCH] dynOneEqEddy: Bound KK during construction Resolves
 bug-report http://openfoam.org/mantisbt/view.php?id=1850

---
 .../incompressible/LES/dynOneEqEddy/dynOneEqEddy.C           | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
index 5e125d53..aa0c1f21 100644
--- a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
+++ b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
@@ -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();