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
applications
solvers
compressible
rhoPimpleFoam
createFields.H
Go to the documentation of this file.
1
Info
<<
"Reading thermophysical properties\n"
<<
endl
;
2
3
autoPtr<basicPsiThermo>
pThermo
4
(
5
basicPsiThermo::New(
mesh
)
6
);
7
basicPsiThermo&
thermo
=
pThermo
();
8
9
volScalarField
&
p
= thermo.p();
10
volScalarField
&
h
= thermo.h();
11
const
volScalarField
&
psi
= thermo.psi();
12
13
volScalarField
rho
14
(
15
IOobject
16
(
17
"rho"
,
18
runTime.timeName(),
19
mesh
,
20
IOobject::READ_IF_PRESENT,
21
IOobject::AUTO_WRITE
22
),
23
thermo.rho()
24
);
25
26
Info
<<
"Reading field U\n"
<<
endl
;
27
volVectorField
U
28
(
29
IOobject
30
(
31
"U"
,
32
runTime.timeName(),
33
mesh
,
34
IOobject::MUST_READ,
35
IOobject::AUTO_WRITE
36
),
37
mesh
38
);
39
40
# include <
finiteVolume/compressibleCreatePhi.H
>
41
42
dimensionedScalar
pMin
43
(
44
mesh
.solutionDict().subDict(
"PIMPLE"
).lookup(
"pMin"
)
45
);
46
47
Info
<<
"Creating turbulence model\n"
<<
endl
;
48
autoPtr<compressible::turbulenceModel>
turbulence
49
(
50
compressible::turbulenceModel::New
51
(
52
rho
,
53
U
,
54
phi
,
55
thermo
56
)
57
);
58
59
//dimensionedScalar initialMass = fvc::domainIntegrate(rho);
60
61
62
Info
<<
"Creating field DpDt\n"
<<
endl
;
63
volScalarField
DpDt
=
64
fvc::DDt
(
surfaceScalarField
(
"phiU"
,
phi
/
fvc::interpolate
(
rho
)), p);
65
66
// ************************ vim: set sw=4 sts=4 et: ************************ //