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 }
rho
rho
Definition: pEqn.H:1
p
p
Definition: pEqn.H:50
thermo
rhoReactionThermo & thermo
Definition: createFields.H:28
forAll
forAll(phases, phasei)
Definition: pEqn.H:3
rho0
volScalarField rho0(rho - psi *p)
constProp
const word constProp(initialConditions.lookup("constantProperty"))
R
#define R(A, B, C, D, E, F, K, M)
W
const scalarList W(::W(thermo))
Y
PtrList< volScalarField > & Y
Definition: createFieldRefs.H:11
Rspecific
Rspecific[0]
Definition: readInitialConditions.H:95