1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
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
13 the Free Software Foundation, either version 3 of the License, or
14 (at your 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, see <http://www.gnu.org/licenses/>.
25 Foam::incompressible::RASModels::LienCubicKELowRe
28 Lien cubic non-linear low-Reynolds k-epsilon turbulence models for
33 Lien, F.S., Chen, W.L., Leschziner, M.A.,
34 "Low-Reynolds-number eddy-viscosity modeling based on non-linear
35 stress-strain/vorticity relations"
36 Engineering Turbulence Modelling and Experiments 3
37 (Edited by Rodi, W. and Bergeles, G.), 91-100. 1996.
38 Elsevier Science Publishers.
41 "Turbulent flow and heat transfer in a square-sectioned U bend"
42 Progress in compuational fluid dynamics 6, 89-100. 2006.
48 \*---------------------------------------------------------------------------*/
50 #ifndef LienCubicKELowRe_H
51 #define LienCubicKELowRe_H
56 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
60 namespace incompressible
65 /*---------------------------------------------------------------------------*\
66 Class LienCubicKELowRe Declaration
67 \*---------------------------------------------------------------------------*/
69 class LienCubicKELowRe
80 dimensionedScalar C1_;
81 dimensionedScalar C2_;
82 dimensionedScalar sigmak_;
83 dimensionedScalar sigmaEps_;
84 dimensionedScalar A1_;
85 dimensionedScalar A2_;
86 dimensionedScalar Ctau1_;
87 dimensionedScalar Ctau2_;
88 dimensionedScalar Ctau3_;
89 dimensionedScalar alphaKsi_;
91 dimensionedScalar CmuWall_;
92 dimensionedScalar kappa_;
94 dimensionedScalar Am_;
95 dimensionedScalar Aepsilon_;
96 dimensionedScalar Amu_;
102 volScalarField epsilon_;
106 volTensorField gradU_;
110 volScalarField fEta_;
111 volScalarField C5viscosity_;
113 volScalarField yStar_;
117 volSymmTensorField nonlinearStress_;
122 //- Runtime type information
123 TypeName("LienCubicKELowRe");
127 //- Construct from components
130 const volVectorField& U,
131 const surfaceScalarField& phi,
132 transportModel& transport,
133 const word& turbulenceModelName = turbulenceModel::typeName,
134 const word& modelName = typeName
139 virtual ~LienCubicKELowRe()
145 //- Return the turbulence viscosity
146 virtual tmp<volScalarField> nut() const
151 //- Return the effective diffusivity for k
152 tmp<volScalarField> DkEff() const
154 return tmp<volScalarField>
156 new volScalarField("DkEff", nut_/sigmak_ + nu())
160 //- Return the effective diffusivity for epsilon
161 tmp<volScalarField> DepsilonEff() const
163 return tmp<volScalarField>
165 new volScalarField("DepsilonEff", nut_/sigmaEps_ + nu())
169 //- Return the turbulence kinetic energy
170 virtual tmp<volScalarField> k() const
175 //- Return the turbulence kinetic energy dissipation rate
176 virtual tmp<volScalarField> epsilon() const
181 //- Return the Reynolds stress tensor
182 virtual tmp<volSymmTensorField> R() const;
184 //- Return the effective stress tensor including the laminar stress
185 virtual tmp<volSymmTensorField> devReff() const;
187 //- Return the source term for the momentum equation
188 virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
190 //- Solve the turbulence equations and correct the turbulence viscosity
191 virtual void correct();
193 //- Read RASProperties dictionary
198 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
200 } // End namespace RASModels
201 } // Edn namespace incompressible
202 } // End namespace Foam
204 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
208 // ************************************************************************* //