FreeFOAM The Cross-Platform CFD Toolkit
LienCubicKELowRe.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
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.
15 
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
19  for more details.
20 
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/>.
23 
24 Class
25  Foam::incompressible::RASModels::LienCubicKELowRe
26 
27 Description
28  Lien cubic non-linear low-Reynolds k-epsilon turbulence models for
29  incompressible flows.
30 
31  References:
32  @verbatim
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.
39 
40  Etemad, S., et al.,
41  "Turbulent flow and heat transfer in a square-sectioned U bend"
42  Progress in compuational fluid dynamics 6, 89-100. 2006.
43  @endverbatim
44 
45 SourceFiles
46  LienCubicKELowRe.C
47 
48 \*---------------------------------------------------------------------------*/
49 
50 #ifndef LienCubicKELowRe_H
51 #define LienCubicKELowRe_H
52 
54 #include <finiteVolume/wallDist.H>
55 
56 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
57 
58 namespace Foam
59 {
60 namespace incompressible
61 {
62 namespace RASModels
63 {
64 
65 /*---------------------------------------------------------------------------*\
66  Class LienCubicKELowRe Declaration
67 \*---------------------------------------------------------------------------*/
68 
70 :
71  public RASModel
72 {
73  // Private data
74 
75  // Model coefficients
76 
79  dimensionedScalar sigmak_;
80  dimensionedScalar sigmaEps_;
83  dimensionedScalar Ctau1_;
84  dimensionedScalar Ctau2_;
85  dimensionedScalar Ctau3_;
86  dimensionedScalar alphaKsi_;
87 
88  dimensionedScalar CmuWall_;
89  dimensionedScalar kappa_;
90 
92  dimensionedScalar Aepsilon_;
93  dimensionedScalar Amu_;
94 
95 
96  // Fields
97 
98  volScalarField k_;
99  volScalarField epsilon_;
100 
101  wallDist y_;
102 
103  volTensorField gradU_;
104  volScalarField eta_;
105  volScalarField ksi_;
106  volScalarField Cmu_;
107  volScalarField fEta_;
108  volScalarField C5viscosity_;
109 
110  volScalarField yStar_;
111 
112  volScalarField nut_;
113 
114  volSymmTensorField nonlinearStress_;
115 
116 
117 public:
118 
119  //- Runtime type information
120  TypeName("LienCubicKELowRe");
121 
122  // Constructors
123 
124  //- Construct from components
126  (
127  const volVectorField& U,
128  const surfaceScalarField& phi,
130  );
131 
132 
133  //- Destructor
135  {}
136 
137 
138  // Member Functions
139 
140  //- Return the turbulence viscosity
141  virtual tmp<volScalarField> nut() const
142  {
143  return nut_;
144  }
145 
146  //- Return the effective diffusivity for k
148  {
149  return tmp<volScalarField>
150  (
151  new volScalarField("DkEff", nut_/sigmak_ + nu())
152  );
153  }
154 
155  //- Return the effective diffusivity for epsilon
157  {
158  return tmp<volScalarField>
159  (
160  new volScalarField("DepsilonEff", nut_/sigmaEps_ + nu())
161  );
162  }
163 
164  //- Return the turbulence kinetic energy
165  virtual tmp<volScalarField> k() const
166  {
167  return k_;
168  }
169 
170  //- Return the turbulence kinetic energy dissipation rate
171  virtual tmp<volScalarField> epsilon() const
172  {
173  return epsilon_;
174  }
175 
176  //- Return the Reynolds stress tensor
177  virtual tmp<volSymmTensorField> R() const;
178 
179  //- Return the effective stress tensor including the laminar stress
180  virtual tmp<volSymmTensorField> devReff() const;
181 
182  //- Return the source term for the momentum equation
184 
185  //- Solve the turbulence equations and correct the turbulence viscosity
186  virtual void correct();
187 
188  //- Read RASProperties dictionary
189  virtual bool read();
190 };
191 
192 
193 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
194 
195 } // End namespace RASModels
196 } // Edn namespace incompressible
197 } // End namespace Foam
198 
199 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
200 
201 #endif
202 
203 // ************************ vim: set sw=4 sts=4 et: ************************ //