From 9c1bb8d707510c00b3b197989c20cb7494935039 Mon Sep 17 00:00:00 2001
From: Andrew Heather <a.heather@opencfd.co.uk>
Date: Fri, 15 Dec 2017 17:25:22 +0000
Subject: [PATCH] ENH: maxDeltaxyz - simplified the delta calc and reset the
default coeff to 2
---
.../LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C b/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
index 351d0fb11f3..284af6ecf2c 100644
--- a/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
+++ b/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
@@ -64,14 +64,14 @@ void Foam::LESModels::maxDeltaxyz::calcDelta()
const point& fc = faceC[facei];
const vector& n = faceN[facei];
- scalar tmp = magSqr(n*(n & (fc - cc)));
+ scalar tmp = mag(n & (fc - cc));
if (tmp > deltaMaxTmp)
{
deltaMaxTmp = tmp;
}
}
- hmax[celli] = deltaCoeff_*Foam::sqrt(deltaMaxTmp);
+ hmax[celli] = deltaCoeff_*deltaMaxTmp;
}
if (nD == 3)
@@ -113,7 +113,7 @@ Foam::LESModels::maxDeltaxyz::maxDeltaxyz
dict.optionalSubDict(type() + "Coeffs").lookupOrDefault<scalar>
(
"deltaCoeff",
- 1
+ 2
)
)
{
--
GitLab