diff --git a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C index e62de890f917fc7ec391c5bb0da8d65e786009a5..4e39627742b28b4f44a2142029650b8e4fa7ecaa 100644 --- a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C +++ b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C @@ -69,7 +69,7 @@ dimensionedScalar homogeneousDynSmagorinsky::cD tmp<volSymmTensorField> LL = dev(filter_(sqr(U())) - (sqr(filter_(U())))); - return average(LL && MM)/MMMM; + return 0.5*average(LL && MM)/MMMM; } else { diff --git a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H index 2a25f48fe05f4bff04e94426e491f12b6510bce2..d25f0e19f1e6e69993505e68a7767bd961d2fe92 100644 --- a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H +++ b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H @@ -48,7 +48,7 @@ Description and - cD=<L.M>/<M.M>, + cD=1/2*<L.M>/<M.M>, where