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::kOmega
29 Standard high Reynolds-number k-omega turbulence model for
34 "Turbulence Modeling for CFD"
36 DCW Industries, Inc., La Canada,
40 http://www.cfd-online.com/Wiki/Wilcox's_k-omega_model
43 The default model coefficients correspond to the following:
47 Cmu 0.09; // Equivalent to betaStar
58 \*---------------------------------------------------------------------------*/
63 #include <incompressibleRASModels/RASModel.H>
65 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
69 namespace incompressible
74 /*---------------------------------------------------------------------------*\
75 Class kOmega Declaration
76 \*---------------------------------------------------------------------------*/
86 dimensionedScalar Cmu_;
87 dimensionedScalar beta_;
88 dimensionedScalar alpha_;
89 dimensionedScalar alphaK_;
90 dimensionedScalar alphaOmega_;
96 volScalarField omega_;
102 //- Runtime type information
107 //- Construct from components
110 const volVectorField& U,
111 const surfaceScalarField& phi,
112 transportModel& transport
123 //- Return the turbulence viscosity
124 virtual tmp<volScalarField> nut() const
129 //- Return the effective diffusivity for k
130 tmp<volScalarField> DkEff() const
132 return tmp<volScalarField>
134 new volScalarField("DkEff", alphaK_*nut_ + nu())
138 //- Return the effective diffusivity for omega
139 tmp<volScalarField> DomegaEff() const
141 return tmp<volScalarField>
143 new volScalarField("DomegaEff", alphaOmega_*nut_ + nu())
147 //- Return the turbulence kinetic energy
148 virtual tmp<volScalarField> k() const
153 //- Return the turbulence specific dissipation rate
154 virtual tmp<volScalarField> omega() const
159 //- Return the turbulence kinetic energy dissipation rate
160 virtual tmp<volScalarField> epsilon() const
162 return tmp<volScalarField>
169 mesh_.time().timeName(),
173 omega_.boundaryField().types()
178 //- Return the Reynolds stress tensor
179 virtual tmp<volSymmTensorField> R() const;
181 //- Return the effective stress tensor including the laminar stress
182 virtual tmp<volSymmTensorField> devReff() const;
184 //- Return the source term for the momentum equation
185 virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
187 //- Solve the turbulence equations and correct the turbulence viscosity
188 virtual void correct();
190 //- Read RASProperties dictionary
195 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
197 } // End namespace RASModels
198 } // End namespace incompressible
199 } // End namespace Foam
201 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
205 // ************************ vim: set sw=4 sts=4 et: ************************ //