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::NonlinearKEShih
29 Shih's quadratic non-linear k-epsilon turbulence model for
35 \*---------------------------------------------------------------------------*/
37 #ifndef NonlinearKEShih_H
38 #define NonlinearKEShih_H
42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
46 namespace incompressible
51 /*---------------------------------------------------------------------------*\
52 Class NonlinearKEShih Declaration
53 \*---------------------------------------------------------------------------*/
63 dimensionedScalar C1_;
64 dimensionedScalar C2_;
65 dimensionedScalar sigmak_;
66 dimensionedScalar sigmaEps_;
67 dimensionedScalar A1_;
68 dimensionedScalar A2_;
69 dimensionedScalar Ctau1_;
70 dimensionedScalar Ctau2_;
71 dimensionedScalar Ctau3_;
72 dimensionedScalar alphaKsi_;
74 dimensionedScalar kappa_;
81 volScalarField epsilon_;
83 volTensorField gradU_;
91 volSymmTensorField nonlinearStress_;
96 //- Runtime type information
97 TypeName("NonlinearKEShih");
101 //- Construct from components
104 const volVectorField& U,
105 const surfaceScalarField& phi,
106 transportModel& transport
111 virtual ~NonlinearKEShih()
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_/sigmak_ + 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 //- Return the turbulence kinetic energy dissipation rate
148 virtual tmp<volScalarField> epsilon() const
153 //- Return the Reynolds stress tensor
154 virtual tmp<volSymmTensorField> R() const;
156 //- Return the effective stress tensor including the laminar stress
157 virtual tmp<volSymmTensorField> devReff() const;
159 //- Return the source term for the momentum equation
160 virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
162 //- Solve the turbulence equations and correct the turbulence viscosity
163 virtual void correct();
165 //- Read RASProperties dictionary
170 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
172 } // End namespace RASModels
173 } // End namespace incompressible
174 } // End namespace Foam
176 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
180 // ************************************************************************* //