Correcting k(grad U) expression

5225ebf3 · Sergio Ferraris · 3e08747b · 5225ebf3
Commit 5225ebf3 authored 12 years ago by Sergio Ferraris
--- 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)));
        }