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] = constant::physicoChemical::RR.value()*invW;
12 
13  p[0] = rho0*Rspecific[0]*thermo.T()[0];
14  rho[0] = rho0;
15  }
16 }
forAll
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
Foam::dimensioned::value
const Type & value() const
Return const reference to value.
Definition: dimensionedType.C:287
rho
rho
Definition: pEqn.H:1
Foam::constant::physicoChemical::RR
const dimensionedScalar RR
Universal gas constant: default SI units: [J/kmol/K].
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:4
Y
PtrList< volScalarField > & Y
Definition: createFieldRefs.H:3
thermo
fluidMulticomponentThermo & thermo
Definition: createFields.H:31