OpenFOAM
4.1
The OpenFOAM Foundation
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
applications
solvers
heatTransfer
buoyantBoussinesqPimpleFoam
buoyantBoussinesqPimpleFoam.C
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) 2011-2016 OpenFOAM Foundation
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
Application
25
buoyantBoussinesqPimpleFoam
26
27
Description
28
Transient solver for buoyant, turbulent flow of incompressible fluids.
29
30
Uses the Boussinesq approximation:
31
\f[
32
rho_{k} = 1 - beta(T - T_{ref})
33
\f]
34
35
where:
36
\f$ rho_{k} \f$ = the effective (driving) kinematic density
37
beta = thermal expansion coefficient [1/K]
38
T = temperature [K]
39
\f$ T_{ref} \f$ = reference temperature [K]
40
41
Valid when:
42
\f[
43
\frac{beta(T - T_{ref})}{rho_{ref}} << 1
44
\f]
45
46
\*---------------------------------------------------------------------------*/
47
48
#include "
fvCFD.H
"
49
#include "
singlePhaseTransportModel.H
"
50
#include "
turbulentTransportModel.H
"
51
#include "
radiationModel.H
"
52
#include "
fvOptions.H
"
53
#include "
pimpleControl.H
"
54
55
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
56
57
int
main(
int
argc,
char
*argv[])
58
{
59
#include "
postProcess.H
"
60
61
#include "
setRootCase.H
"
62
#include "
createTime.H
"
63
#include "
createMesh.H
"
64
#include "createControl.H"
65
#include "createFields.H"
66
#include "
createFvOptions.H
"
67
#include "
createTimeControls.H
"
68
#include "CourantNo.H"
69
#include "setInitialDeltaT.H"
70
#include "initContinuityErrs.H"
71
72
turbulence
->validate();
73
74
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
75
76
Info
<<
"\nStarting time loop\n"
<<
endl
;
77
78
while
(runTime.run())
79
{
80
#include "
readTimeControls.H
"
81
#include "CourantNo.H"
82
#include "setDeltaT.H"
83
84
runTime++;
85
86
Info
<<
"Time = "
<< runTime.timeName() <<
nl
<<
endl
;
87
88
// --- Pressure-velocity PIMPLE corrector loop
89
while
(
pimple
.loop())
90
{
91
#include "UEqn.H"
92
#include "TEqn.H"
93
94
// --- Pressure corrector loop
95
while
(
pimple
.correct())
96
{
97
#include "pEqn.H"
98
}
99
100
if
(
pimple
.turbCorr())
101
{
102
laminarTransport
.correct();
103
turbulence
->correct();
104
}
105
}
106
107
runTime.write();
108
109
Info
<<
"ExecutionTime = "
<< runTime.elapsedCpuTime() <<
" s"
110
<<
" ClockTime = "
<< runTime.elapsedClockTime() <<
" s"
111
<<
nl
<<
endl
;
112
}
113
114
Info
<<
"End\n"
<<
endl
;
115
116
return
0;
117
}
118
119
120
// ************************************************************************* //
turbulence
autoPtr< compressible::turbulenceModel > turbulence
Definition:
createFields.H:23
fvCFD.H
createFvOptions.H
radiationModel.H
pimple
const dictionary & pimple
Definition:
readFluidMultiRegionPIMPLEControls.H:1
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition:
Ostream.H:253
readTimeControls.H
Read the control parameters used by setDeltaT.
createTime.H
turbulentTransportModel.H
laminarTransport
singlePhaseTransportModel laminarTransport(U, phi)
Foam::nl
static const char nl
Definition:
Ostream.H:262
createMesh.H
fvOptions.H
setRootCase.H
pimpleControl.H
Foam::Info
messageStream Info
postProcess.H
Execute application functionObjects to post-process existing results.
createTimeControls.H
Read the control parameters used by setDeltaT.
singlePhaseTransportModel.H
Generated by
1.8.11