OpenFOAM
dev
The OpenFOAM Foundation
pEqn.H
Go to the documentation of this file.
1
{
2
rho
=
thermo
.rho();
3
if
(
constProp
==
"volume"
)
4
{
5
scalar invW = 0.0;
6
forAll
(
Y
, i)
7
{
8
invW +=
Y
[i][0]/
W
[i];
9
}
10
11
Rspecific
[0] = 1000.0*
constant::physicoChemical::R
.value()*invW;
12
13
p
[0] =
rho0
*
Rspecific
[0]*
thermo
.T()[0];
14
rho
[0] =
rho0
;
15
}
16
}
R
#define R(A, B, C, D, E, F, K, M)
forAll
#define forAll(list, i)
Loop across all elements in list.
Definition:
UList.H:434
rho
rho
Definition:
pEqn.H:1
rho0
scalar rho0
Definition:
readInitialConditions.H:92
Rspecific
Rspecific[0]
Definition:
readInitialConditions.H:95
W
const scalarList W(::W(thermo))
constProp
const word constProp(initialConditions.lookup("constantProperty"))
p
volScalarField & p
Definition:
createFieldRefs.H:5
Y
PtrList< volScalarField > & Y
Definition:
createFieldRefs.H:4
thermo
fluidMulticomponentThermo & thermo
Definition:
createFields.H:31
applications
solvers
chemFoam
pEqn.H
Generated by
1.9.1