Loading...
Searching...
No Matches
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]/specieData[i].W();
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}
thermo
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
p
volScalarField & p
Definition createFieldRefs.H:8
Y
PtrList< volScalarField > & Y
Definition createFieldRefs.H:7
specieData
PtrList< gasHThermoPhysics > specieData(Y.size())
rho0
scalar rho0
Definition readInitialConditions.H:89
constProp
word constProp(initialConditions.get< word >("constantProperty"))
forAll
#define forAll(list, i)
Loop across all elements in list.
Definition stdFoam.H:299