FreeFOAM The Cross-Platform CFD Toolkit
LienLeschzinerLowRe.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::LienLeschzinerLowRe
26 
27 Description
28  Lien and Leschziner low-Reynolds k-epsilon turbulence model for
29  incompressible flows.
30 
31 SourceFiles
32  LienLeschzinerLowRe.C
33 
34 \*---------------------------------------------------------------------------*/
35 
36 #ifndef LienLeschzinerLowRe_H
37 #define LienLeschzinerLowRe_H
38 
40 #include <finiteVolume/wallDist.H>
41 
42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43 
44 namespace Foam
45 {
46 namespace incompressible
47 {
48 namespace RASModels
49 {
50 
51 /*---------------------------------------------------------------------------*\
52  Class LienLeschzinerLowRe Declaration
53 \*---------------------------------------------------------------------------*/
54 
56 :
57  public RASModel
58 {
59  // Private data
60 
61  // Model coefficients
62 
65  dimensionedScalar sigmak_;
66  dimensionedScalar sigmaEps_;
67  dimensionedScalar Cmu_;
68  dimensionedScalar kappa_;
69 
71  dimensionedScalar Aepsilon_;
72  dimensionedScalar Amu_;
73 
74 
75  // Fields
76 
77  volScalarField k_;
78  volScalarField epsilon_;
79 
80  wallDist y_;
81 
82  volScalarField yStar_;
83 
84  volScalarField nut_;
85 
86 
87 public:
88 
89  TypeName("LienLeschzinerLowRe");
90 
91  // Constructors
92 
93  //- Construct from components
95  (
96  const volVectorField& U,
97  const surfaceScalarField& phi,
99  );
100 
101 
102  //- Destructor
104  {}
105 
106 
107  // Member Functions
108 
109  //- Return the turbulence viscosity
110  virtual tmp<volScalarField> nut() const
111  {
112  return nut_;
113  }
114 
115  //- Return the effective diffusivity for k
117  {
118  return tmp<volScalarField>
119  (
120  new volScalarField("DkEff", nut_/sigmak_ + nu())
121  );
122  }
123 
124  //- Return the effective diffusivity for epsilon
126  {
127  return tmp<volScalarField>
128  (
129  new volScalarField("DepsilonEff", nut_/sigmaEps_ + nu())
130  );
131  }
132 
133  //- Return the turbulence kinetic energy
134  virtual tmp<volScalarField> k() const
135  {
136  return k_;
137  }
138 
139  //- Return the turbulence kinetic energy dissipation rate
140  virtual tmp<volScalarField> epsilon() const
141  {
142  return epsilon_;
143  }
144 
145  //- Return the Reynolds stress tensor
146  virtual tmp<volSymmTensorField> R() const;
147 
148  //- Return the effective stress tensor including the laminar stress
149  virtual tmp<volSymmTensorField> devReff() const;
150 
151  //- Return the source term for the momentum equation
153 
154  //- Solve the turbulence equations and correct the turbulence viscosity
155  virtual void correct();
156 
157  //- Read RASProperties dictionary
158  virtual bool read();
159 };
160 
161 
162 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
163 
164 } // End namespace RASModels
165 } // End namespace incompressible
166 } // End namespace Foam
167 
168 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
169 
170 #endif
171 
172 // ************************ vim: set sw=4 sts=4 et: ************************ //