From 3e08747b29f75def403560a76833b391338e9dab Mon Sep 17 00:00:00 2001
From: Sergio Ferraris <sergio@alex.opencfd.co.uk>
Date: Mon, 15 Apr 2013 15:18:41 +0100
Subject: [PATCH] Updating copyright years

---
 .../incompressible/LES/dynLagrangian/dynLagrangian.C       | 2 +-
 .../incompressible/LES/dynLagrangian/dynLagrangian.H       | 7 +++++--
 2 files changed, 6 insertions(+), 3 deletions(-)

diff --git a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C
index 32d23f4eab9..bfadebd7554 100644
--- a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C
+++ b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C
@@ -47,7 +47,7 @@ void dynLagrangian::updateSubGridScaleFields
     const tmp<volTensorField>& gradU
 )
 {
-    nuSgs_ = (flm_/fmm_)*delta()*sqrt(k(gradU));
+    nuSgs_ = 2.0*(flm_/fmm_)*sqr(delta())*mag(dev(symm(gradU)));
     nuSgs_.correctBoundaryConditions();
 }
 
diff --git a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
index 004efeaaed7..60f86c53cc5 100644
--- a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
+++ b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -153,7 +153,10 @@ public:
         //- Return SGS kinetic energy
         tmp<volScalarField> k(const tmp<volTensorField>& gradU) const
         {
-            return 2.0*sqr(delta())*magSqr(dev(symm(gradU)));
+            return
+                2.0
+              * pow((2.0/sqrt(2.0))*sqrt(ce_)*flm_/fmm_, 2.0/3.0)/ce_
+              * sqr(delta())*magSqr(dev(symm(gradU)));
         }
 
         //- Return SGS kinetic energy
-- 
GitLab