Home
Downloads
Documentation
Installation
User Guide
man-pages
API Documentation
README
Release Notes
Changes
License
Support
SourceForge Project
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
src
turbulenceModels
incompressible
RAS
LienCubicKELowRe
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
53
#include <
incompressibleRASModels/RASModel.H
>
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
69
class
LienCubicKELowRe
70
:
71
public
RASModel
72
{
73
// Private data
74
75
// Model coefficients
76
77
dimensionedScalar
C1_;
78
dimensionedScalar
C2_;
79
dimensionedScalar
sigmak_;
80
dimensionedScalar
sigmaEps_;
81
dimensionedScalar
A1_;
82
dimensionedScalar
A2_;
83
dimensionedScalar
Ctau1_;
84
dimensionedScalar
Ctau2_;
85
dimensionedScalar
Ctau3_;
86
dimensionedScalar
alphaKsi_;
87
88
dimensionedScalar
CmuWall_;
89
dimensionedScalar
kappa_;
90
91
dimensionedScalar
Am_;
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
125
LienCubicKELowRe
126
(
127
const
volVectorField
&
U
,
128
const
surfaceScalarField
&
phi
,
129
transportModel
&
transport
130
);
131
132
133
//- Destructor
134
virtual
~LienCubicKELowRe
()
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
147
tmp<volScalarField>
DkEff
()
const
148
{
149
return
tmp<volScalarField>
150
(
151
new
volScalarField
(
"DkEff"
, nut_/sigmak_ +
nu
())
152
);
153
}
154
155
//- Return the effective diffusivity for epsilon
156
tmp<volScalarField>
DepsilonEff
()
const
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
183
virtual
tmp<fvVectorMatrix>
divDevReff
(
volVectorField
& U)
const
;
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: ************************ //