Updates

src/LESmodels/compressible/Smagorinsky/Smagorinsky.C

@@ -71,7 +71,7 @@ void Smagorinsky::correct(const tmp<volTensorField>& gradU)
     volScalarField b = (2.0/3.0)*tr(D);
     volScalarField c = 2*ck_*delta()*(dev(D) && D);
 
-    k_ = sqr((2*b + sqrt(sqr(b) + 4*a*c))/(2*a));
+    k_ = sqr((-b + sqrt(sqr(b) + 4*a*c))/(2*a));
 
     muSgs_ = ck_*rho()*delta()*sqrt(k_);
     muSgs_.correctBoundaryConditions();

src/LESmodels/compressible/Smagorinsky/Smagorinsky.H

@@ -38,8 +38,8 @@ Description
     where
 
         D = symm(grad(U));
-        k = (2*ck/ce)*delta^2*grad(U):dev(D)
-        nuSgs = ck*sqrt(k)*delta
+        k from rho*D:B + ce*rho*k^3/2/delta = 0
+        muSgs = ck*rho*sqrt(k)*delta
     @endverbatim
 
 SourceFiles

src/LESmodels/compressible/oneEqEddy/oneEqEddy.C

@@ -43,7 +43,6 @@ addToRunTimeSelectionTable(LESmodel, oneEqEddy, dictionary);
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-// from components
 oneEqEddy::oneEqEddy
 (
     const volScalarField& rho,

src/LESmodels/compressible/oneEqEddy/oneEqEddy.H

@@ -33,7 +33,7 @@ Description
     @verbatim
         d/dt(rho*k) + div(rho*U*k) - div(muEff*grad(k))
         =
-        -rho*B*L - ce*rho*k^3/2/delta
+        -rho*D:B - ce*rho*k^3/2/delta
 
     and
 
@@ -42,7 +42,7 @@ Description
     where
 
         D = symm(grad(U));
-        nuSgs = ck*sqrt(k)*delta
+        muSgs = ck*rho*sqrt(k)*delta
     @endverbatim
 
 
@@ -114,7 +114,7 @@ public:
         //- Return the effective diffusivity for k
         tmp<volScalarField> DkEff() const
         {
-            return tmp<volScalarField> 
+            return tmp<volScalarField>
             (
                 new volScalarField("DkEff", muSgs_ + mu())
             );