research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Correcting Cp_ref calculation in MultiTransport, which previously

Browse source 
assumed an ideal gas phase.
This commit is contained in:
Steven DeCaluwe 2014-12-16 13:32:08 +00:00
parent f3b66ec5f7
commit 661758d64d
1 changed files with 5 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
src/transport/MultiTransport.cpp
View file

@ -542,6 +542,11 @@ void MultiTransport::updateThermal_T()
* The original Dixon-Lewis paper subtracted 1.5 here.
*/
const vector_fp& cp = ((IdealGasPhase*)m_thermo)->cp_R_ref();
vector_fp cp_R_all(m_thermo->nSpecies());
m_thermo->getCp_R_ref(&cp_R_all[0]);
const vector_fp& cp = cp_R_all;
for (size_t k = 0; k < m_nsp; k++) {
m_cinternal[k] = cp[k] - 2.5;
}