From f2b5b2527cee44ee60191751145d76c01dfb2ae5 Mon Sep 17 00:00:00 2001
From: Andrew Heather <>
Date: Fri, 1 Mar 2019 11:30:20 +0000
Subject: [PATCH] ENH: DEShybrid - added robustness for case that nuEff goes
negative (typically from negative denisty in compressible cals). See #1224
---
.../schemes/DEShybrid/DEShybrid.H | 15 ++++++++++++++-
1 file changed, 14 insertions(+), 1 deletion(-)
diff --git a/src/TurbulenceModels/schemes/DEShybrid/DEShybrid.H b/src/TurbulenceModels/schemes/DEShybrid/DEShybrid.H
index c4218e2295..8a17edf4ef 100644
--- a/src/TurbulenceModels/schemes/DEShybrid/DEShybrid.H
+++ b/src/TurbulenceModels/schemes/DEShybrid/DEShybrid.H
@@ -197,13 +197,26 @@ class DEShybrid
CH3_*Omega*max(S, Omega)
/max(0.5*(sqr(S) + sqr(Omega)), sqr(OmegaLim_/tau0_))
);
+
const volScalarField K
(
max(Foam::sqrt(0.5*(sqr(S) + sqr(Omega))), 0.1/tau0_)
);
- const volScalarField lTurb(Foam::sqrt(nuEff/(pow(0.09, 1.5)*K)));
+ const volScalarField lTurb
+ (
+ Foam::sqrt
+ (
+ max
+ (
+ nuEff/(pow(0.09, 1.5)*K),
+ dimensionedScalar("l0", sqr(dimLength), 0)
+ )
+ )
+ );
+
const volScalarField g(tanh(pow4(B)));
+
const volScalarField A
(
CH2_*max
--
GitLab