diff --git a/src/turbulenceModels/compressible/LES/oneEqEddy/oneEqEddy.H b/src/turbulenceModels/compressible/LES/oneEqEddy/oneEqEddy.H
index aa00969499aacb503e55d1c7cdedf6cd06f8bb63..0b463fb5c997439ebc7bc3ea51088ee34c050863 100644
--- a/src/turbulenceModels/compressible/LES/oneEqEddy/oneEqEddy.H
+++ b/src/turbulenceModels/compressible/LES/oneEqEddy/oneEqEddy.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -25,14 +25,14 @@ Class
Foam::compressible::LESModels::oneEqEddy
Description
- One Equation Eddy Viscosity Model for incompressible flows
+ One Equation Eddy Viscosity Model for compressible flows
Eddy viscosity SGS model using a modeled balance equation to simulate the
behaviour of k, hence,
\verbatim
d/dt(rho*k) + div(rho*U*k) - div(muEff*grad(k))
=
- -rho*D:B - ce*rho*k^3/2/delta
+ -rho*D:B - ce*rho*k^(3/2)/delta
and
@@ -41,10 +41,11 @@ Description
where
D = symm(grad(U));
- muSgs = ck*rho*sqrt(k)*delta
+ nuSgs = ck*sqrt(k)*delta
+ muSgs = rho*nuSgs
+ muEff = muSgs + mu
\endverbatim
-
SourceFiles
oneEqEddy.C
diff --git a/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H b/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H
index 0f040d829a3cdc7164a7ec8600fa016cc79892ee..143329c8c62a37131450b3c9f6f76ed9f1a35540 100644
--- a/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H
+++ b/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -32,7 +32,7 @@ Description
\verbatim
d/dt(k) + div(U*k) - div(nuEff*grad(k))
=
- -B*L - ce*k^3/2/delta
+ -D:B - ce*k^(3/2)/delta
and