Modifying TransportFactory::fitProperties to eliminate an assumption

that the phase of interest is an ideal gas.
This commit is contained in:
Steven DeCaluwe 2014-07-28 18:18:41 +00:00
parent 6c3867ca79
commit 8954090881

View file

@ -1145,7 +1145,10 @@ void TransportFactory::fitProperties(GasTransportParams& tr,
t = tr.tmin + dt*n;
tr.thermo->setTemperature(t);
cp_R = ((IdealGasPhase*)tr.thermo)->cp_R_ref()[k];
vector_fp cp_R_all(tr.thermo->nSpecies());
tr.thermo->getCp_R_ref(&cp_R_all[0]);
cp_R = cp_R_all[k];
tstar = Boltzmann * t/ tr.eps[k];
sqrt_T = sqrt(t);