1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software; you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
26 Foam::compressible::LESModels::oneEqEddy
29 One Equation Eddy Viscosity Model for incompressible flows
31 Eddy viscosity SGS model using a modeled balance equation to simulate the
32 behaviour of k, hence,
34 d/dt(rho*k) + div(rho*U*k) - div(muEff*grad(k))
36 -rho*D:B - ce*rho*k^3/2/delta
40 B = 2/3*k*I - 2*nuSgs*dev(D)
45 muSgs = ck*rho*sqrt(k)*delta
52 \*---------------------------------------------------------------------------*/
54 #ifndef compressibleOneEqEddy_H
55 #define compressibleOneEqEddy_H
57 #include "GenEddyVisc.H"
59 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
63 namespace compressible
68 /*---------------------------------------------------------------------------*\
69 Class oneEqEddy Declaration
70 \*---------------------------------------------------------------------------*/
78 dimensionedScalar ck_;
81 // Private Member Functions
83 //- Update sub-grid scale fields
84 void updateSubGridScaleFields();
86 // Disallow default bitwise copy construct and assignment
87 oneEqEddy(const oneEqEddy&);
88 oneEqEddy& operator=(const oneEqEddy&);
93 //- Runtime type information
94 TypeName("oneEqEddy");
99 //- Constructor from components
102 const volScalarField& rho,
103 const volVectorField& U,
104 const surfaceScalarField& phi,
105 const basicThermo& thermoPhysicalModel
116 //- Return the effective diffusivity for k
117 tmp<volScalarField> DkEff() const
119 return tmp<volScalarField>
121 new volScalarField("DkEff", muSgs_ + mu())
125 //- Correct Eddy-Viscosity and related properties
126 virtual void correct(const tmp<volTensorField>& gradU);
128 //- Read LESProperties dictionary
133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
135 } // End namespace LESModels
136 } // End namespace compressible
137 } // End namespace Foam
139 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
143 // ************************************************************************* //