From 473fbfb8a3a78056c9b612b604cf80b306df1012 Mon Sep 17 00:00:00 2001
From: henry <henry@dm.opencfd.co.uk>
Date: Thu, 1 May 2008 18:43:49 +0100
Subject: [PATCH] Updates
---
src/LESmodels/compressible/Smagorinsky/Smagorinsky.C | 2 +-
src/LESmodels/compressible/Smagorinsky/Smagorinsky.H | 4 ++--
src/LESmodels/compressible/oneEqEddy/oneEqEddy.C | 1 -
src/LESmodels/compressible/oneEqEddy/oneEqEddy.H | 6 +++---
4 files changed, 6 insertions(+), 7 deletions(-)
diff --git a/src/LESmodels/compressible/Smagorinsky/Smagorinsky.C b/src/LESmodels/compressible/Smagorinsky/Smagorinsky.C
index 3cf16914963..0bee1a50ceb 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 ef2bbf72c9e..ba951aa86e1 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 6a9484b6774..9788ad61158 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 06f24988147..a4f0d150f78 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())
);
--
GitLab