Fixed calculation of some multicomponent transport properties

For some properties, the internal temperature-dependent properties
were being updated in the wrong order.
This commit is contained in:
Ray Speth 2012-07-10 16:23:32 +00:00
parent b83a3b7e96
commit 7194bf96d4
2 changed files with 2 additions and 1 deletions

View file

@ -212,6 +212,7 @@ void GasTransport::updateViscosity_T()
void GasTransport::updateSpeciesViscosities()
{
update_T();
if (m_mode == CK_Mode) {
for (size_t k = 0; k < m_nsp; k++) {
m_visc[k] = exp(dot4(m_polytempvec, m_visccoeffs[k]));
@ -229,6 +230,7 @@ void GasTransport::updateSpeciesViscosities()
void GasTransport::updateDiff_T()
{
update_T();
// evaluate binary diffusion coefficients at unit pressure
size_t ic = 0;
if (m_mode == CK_Mode) {

View file

@ -655,7 +655,6 @@ void MultiTransport::updateThermal_T()
}
// we need species viscosities and binary diffusion coefficients
updateSpeciesViscosities();
update_T();
updateDiff_T();
// evaluate polynomial fits for A*, B*, C*