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 }
p
Definition: pEqn.H:50
forAll(phases, phasei)
Definition: pEqn.H:5
psiReactionThermo & thermo
Definition: createFields.H:31
word constProp(initialConditions.lookup("constantProperty"))
volScalarField rho0(rho - psi *p)
PtrList< gasHThermoPhysics > specieData(Y.size())
#define R(A, B, C, D, E, F, K, M)
PtrList< volScalarField > & Y
Rspecific[0]
rho
Definition: pEqn.H:1