1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2008 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::LamBremhorstKE
29 Lam and Bremhorst low-Reynolds number k-epsilon turbulence model
30 for incompressible flows
35 \*---------------------------------------------------------------------------*/
37 #ifndef LamBremhorstKE_H
38 #define LamBremhorstKE_H
43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
47 namespace incompressible
52 /*---------------------------------------------------------------------------*\
53 Class LamBremhorstKE Declaration
54 \*---------------------------------------------------------------------------*/
62 dimensionedScalar Cmu_;
63 dimensionedScalar C1_;
64 dimensionedScalar C2_;
65 dimensionedScalar alphaEps_;
68 volScalarField epsilon_;
79 //- Runtime type information
80 TypeName("LamBremhorstKE");
88 const volVectorField& U,
89 const surfaceScalarField& phi,
90 transportModel& transport
102 //- Return the turbulence viscosity
103 tmp<volScalarField> nut() const
108 //- Return the effective diffusivity for k
109 tmp<volScalarField> DkEff() const
111 return tmp<volScalarField>
113 new volScalarField("DkEff", nut_ + nu())
117 //- Return the effective diffusivity for epsilon
118 tmp<volScalarField> DepsilonEff() const
120 return tmp<volScalarField>
122 new volScalarField("DepsilonEff", alphaEps_*nut_ + nu())
126 //- Return the turbulence kinetic energy
127 tmp<volScalarField> k() const
132 //- Return the turbulence kinetic energy dissipation rate
133 tmp<volScalarField> epsilon() const
138 //- Return the Reynolds stress tensor
139 tmp<volSymmTensorField> R() const;
141 //- Return the effective stress tensor including the laminar stress
142 tmp<volSymmTensorField> devReff() const;
144 //- Return the source term for the momentum equation
145 tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
147 //- Solve the turbulence equations and correct the turbulence viscosity
150 //- Read turbulenceProperties dictionary
155 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
157 } // End namespace RASModels
158 } // End namespace incompressible
159 } // End namespace Foam
161 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
165 // ************************************************************************* //