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
LES
spectEddyVisc
spectEddyVisc.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::LESModels::spectEddyVisc
26
27
Description
28
The Isochoric spectral Eddy Viscosity Model for incompressible flows.
29
30
Algebraic eddy viscosity SGS model founded on the assumption that
31
local equilibrium prevail and that viscous effects may be of importance.
32
33
Thus,
34
@verbatim
35
B = 2/3*k*I - 2*nuSgs*dev(D)
36
Beff = 2/3*k*I - 2*nuEff*dev(D)
37
38
where
39
40
k = cK1*delta^(2/3)*eps^(2/3)*exp(-cK2*delta^(4/3)*nu*eps^(-1/3)) -
41
ck3*sqrt(eps*nu)*erfc(ck4*delta^(-2/3)*sqrt(nu)*eps^(-1/6)))
42
43
nuEff = nu/(1 - exp(-cB*pow((nu/nuEff), 1.0/3.0)*pow(Re, -2.0/3.0)))
44
nuSgs = nuEff - nu
45
46
Re = delta^2*mag(D)/nu
47
@endverbatim
48
49
SourceFiles
50
spectEddyVisc.C
51
52
\*---------------------------------------------------------------------------*/
53
54
#ifndef spectEddyVisc_H
55
#define spectEddyVisc_H
56
57
#include <
incompressibleLESModels/GenEddyVisc.H
>
58
59
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
60
61
namespace
Foam
62
{
63
namespace
incompressible
64
{
65
namespace
LESModels
66
{
67
68
/*---------------------------------------------------------------------------*\
69
Class spectEddyVisc Declaration
70
\*---------------------------------------------------------------------------*/
71
72
class
spectEddyVisc
73
:
74
public
GenEddyVisc
75
{
76
// Private data
77
78
dimensionedScalar
cB_;
79
dimensionedScalar
cK1_;
80
dimensionedScalar
cK2_;
81
dimensionedScalar
cK3_;
82
dimensionedScalar
cK4_;
83
84
85
// Private Member Functions
86
87
//- Update sub-grid scale fields
88
void
updateSubGridScaleFields(
const
volTensorField
& gradU);
89
90
// Disallow default bitwise copy construct and assignment
91
spectEddyVisc
(
const
spectEddyVisc
&);
92
spectEddyVisc
& operator=(
const
spectEddyVisc
&);
93
94
95
public
:
96
97
//- Runtime type information
98
TypeName
(
"spectEddyVisc"
);
99
100
// Constructors
101
102
// Construct from components
103
spectEddyVisc
104
(
105
const
volVectorField
&
U
,
106
const
surfaceScalarField
&
phi
,
107
transportModel
&
transport
108
);
109
110
111
//- Destructor
112
virtual
~spectEddyVisc
()
113
{}
114
115
116
// Member Functions
117
118
//- Return SGS kinetic energy
119
virtual
tmp<volScalarField>
k
()
const
;
120
121
//- Correct Eddy-Viscosity and related properties
122
virtual
void
correct
(
const
tmp<volTensorField>
&);
123
124
//- Read LESProperties dictionary
125
virtual
bool
read
();
126
};
127
128
129
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
130
131
}
// End namespace LESModels
132
}
// End namespace incompressible
133
}
// End namespace Foam
134
135
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
136
137
#endif
138
139
// ************************ vim: set sw=4 sts=4 et: ************************ //