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::incompressible::RASModels::LaunderSharmaKE
29 Launder and Sharma low-Reynolds k-epsilon turbulence model for
32 The default model coefficients correspond to the following:
47 \*---------------------------------------------------------------------------*/
49 #ifndef LaunderSharmaKE_H
50 #define LaunderSharmaKE_H
54 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
58 namespace incompressible
63 /*---------------------------------------------------------------------------*\
64 Class LaunderSharmaKE Declaration
65 \*---------------------------------------------------------------------------*/
75 dimensionedScalar Cmu_;
76 dimensionedScalar C1_;
77 dimensionedScalar C2_;
78 dimensionedScalar sigmaEps_;
84 volScalarField epsilonTilda_;
88 // Private member functions
90 tmp<volScalarField> fMu() const;
91 tmp<volScalarField> f2() const;
96 //- Runtime type information
97 TypeName("LaunderSharmaKE");
101 //- Construct from components
104 const volVectorField& U,
105 const surfaceScalarField& phi,
106 transportModel& transport
111 virtual ~LaunderSharmaKE()
117 //- Return the turbulence viscosity
118 virtual tmp<volScalarField> nut() const
123 //- Return the effective diffusivity for k
124 tmp<volScalarField> DkEff() const
126 return tmp<volScalarField>
128 new volScalarField("DkEff", nut_ + nu())
132 //- Return the effective diffusivity for epsilon
133 tmp<volScalarField> DepsilonEff() const
135 return tmp<volScalarField>
137 new volScalarField("DepsilonEff", nut_/sigmaEps_ + nu())
141 //- Return the turbulence kinetic energy
142 virtual tmp<volScalarField> k() const
147 //- Note that epsilonTilda is returned as epsilon.
148 // This is the appropriate variable for most purposes.
149 virtual tmp<volScalarField> epsilon() const
151 return epsilonTilda_;
154 //- Return the Reynolds stress tensor
155 virtual tmp<volSymmTensorField> R() const;
157 //- Return the effective stress tensor including the laminar stress
158 virtual tmp<volSymmTensorField> devReff() const;
160 //- Return the source term for the momentum equation
161 virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
163 //- Solve the turbulence equations and correct the turbulence viscosity
164 virtual void correct();
166 //- Read RASProperties dictionary
171 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
173 } // End namespace RASModels
174 } // End namespace incompressible
175 } // End namespace Foam
177 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
181 // ************************************************************************* //