From 661758d64dd786d0a1e880d643aef647487ee8be Mon Sep 17 00:00:00 2001 From: Steven DeCaluwe Date: Tue, 16 Dec 2014 13:32:08 +0000 Subject: [PATCH] Correcting Cp_ref calculation in MultiTransport, which previously assumed an ideal gas phase. --- src/transport/MultiTransport.cpp | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/src/transport/MultiTransport.cpp b/src/transport/MultiTransport.cpp index 88b41d243..b09cb9e00 100644 --- a/src/transport/MultiTransport.cpp +++ b/src/transport/MultiTransport.cpp @@ -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; }