diff --git a/src/LESmodels/compressible/Smagorinsky/Smagorinsky.C b/src/LESmodels/compressible/Smagorinsky/Smagorinsky.C
index 3cf1691496384cad8833baf0ef77c80495aec022..0bee1a50ceb0bcbd4d5307fdff00082e56c4a277 100644
--- a/src/LESmodels/compressible/Smagorinsky/Smagorinsky.C
+++ b/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();
diff --git a/src/LESmodels/compressible/Smagorinsky/Smagorinsky.H b/src/LESmodels/compressible/Smagorinsky/Smagorinsky.H
index ef2bbf72c9e03cbc360607fa82cbfb5f009c7e02..ba951aa86e15a830a98366d5449b768b45730469 100644
--- a/src/LESmodels/compressible/Smagorinsky/Smagorinsky.H
+++ b/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
diff --git a/src/LESmodels/compressible/oneEqEddy/oneEqEddy.C b/src/LESmodels/compressible/oneEqEddy/oneEqEddy.C
index 6a9484b6774aec63870d5085294876ed202785e9..9788ad6115802a812f8da418c6765df788d24f27 100644
--- a/src/LESmodels/compressible/oneEqEddy/oneEqEddy.C
+++ b/src/LESmodels/compressible/oneEqEddy/oneEqEddy.C
@@ -43,7 +43,6 @@ addToRunTimeSelectionTable(LESmodel, oneEqEddy, dictionary);
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-// from components
 oneEqEddy::oneEqEddy
 (
     const volScalarField& rho,
diff --git a/src/LESmodels/compressible/oneEqEddy/oneEqEddy.H b/src/LESmodels/compressible/oneEqEddy/oneEqEddy.H
index 06f249881479b1536479fe073378d1fdaf20532c..a4f0d150f789662bc41c6e6d12b4851e9b065818 100644
--- a/src/LESmodels/compressible/oneEqEddy/oneEqEddy.H
+++ b/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())
             );