From 5225ebf39fde9320697459e97f49762d0fadbea7 Mon Sep 17 00:00:00 2001
From: Sergio Ferraris <sergio@alex.opencfd.co.uk>
Date: Mon, 15 Apr 2013 15:46:13 +0100
Subject: [PATCH] Correcting k(grad U) expression
---
.../incompressible/LES/dynLagrangian/dynLagrangian.H | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
index 60f86c53cc5..4809ee902ad 100644
--- a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
+++ b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
@@ -154,8 +154,9 @@ public:
tmp<volScalarField> k(const tmp<volTensorField>& gradU) const
{
return
- 2.0
- * pow((2.0/sqrt(2.0))*sqrt(ce_)*flm_/fmm_, 2.0/3.0)/ce_
+ pow(2.0, 4.0/3.0)
+ * pow(flm_/fmm_, 2.0/3.0)
+ * pow(ce_, -2.0/3.0)
* sqr(delta())*magSqr(dev(symm(gradU)));
}
--
GitLab