From df79844443154d50475044abbc7a569bf9d68131 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Tue, 15 Dec 2015 15:41:01 -0800
Subject: [PATCH] ENH: Adding kMin to k in dTilda function as dTilda can not be
 zero in case k = 0

---
 .../turbulenceModels/DES/kOmegaSSTDES/kOmegaSSTDES.C       | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/src/TurbulenceModels/turbulenceModels/DES/kOmegaSSTDES/kOmegaSSTDES.C b/src/TurbulenceModels/turbulenceModels/DES/kOmegaSSTDES/kOmegaSSTDES.C
index 449eb03e15..ef36104006 100644
--- a/src/TurbulenceModels/turbulenceModels/DES/kOmegaSSTDES/kOmegaSSTDES.C
+++ b/src/TurbulenceModels/turbulenceModels/DES/kOmegaSSTDES/kOmegaSSTDES.C
@@ -62,7 +62,12 @@ tmp<volScalarField> kOmegaSSTDES<BasicTurbulenceModel>::dTilda
     const volScalarField& k = this->k_;
     const volScalarField& omega = this->omega_;
 
-    return min(CDES*this->delta(), sqrt(k)/(this->betaStar_*omega));
+    return
+        min
+        (
+            CDES*this->delta(),
+            sqrt(k + this->kMin_)/(this->betaStar_*omega)
+        );
 }