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

Fixed update logic for temperature-dependent properties in MultiTransport

Browse source
This commit is contained in:
Ray Speth 2012-05-30 16:02:35 +00:00
parent cddd2ec724
commit f24f8f6d3b
2 changed files with 10 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
src/transport/GasTransport.cpp
View file

@ -136,7 +136,12 @@ bool GasTransport::initGas(GasTransportParams& tr)
}
void GasTransport::update_T(void) {
m_temp = m_thermo->temperature();
double T = m_thermo->temperature();
if (T == m_temp) {
return;
}
m_temp = T;
m_kbt = Boltzmann * m_temp;
m_sqrt_kbt = sqrt(Boltzmann*m_temp);
m_logt = log(m_temp);

4
src/transport/MultiTransport.cpp
View file

@ -307,6 +307,7 @@ void MultiTransport::getSpeciesFluxes(size_t ndim, const doublereal* const grad_
size_t ldf, doublereal* const fluxes)
{
// update the binary diffusion coefficients if necessary
update_T();
updateDiff_T();
// If any component of grad_T is non-zero, then get the
@ -456,6 +457,7 @@ void MultiTransport::getMassFluxes(const doublereal* state1, const doublereal* s
m_thermo->getMoleFractions(DATA_PTR(m_molefracs));
// update the binary diffusion coefficients if necessary
update_T();
updateDiff_T();
// If there is a temperature gradient, then get the
@ -576,6 +578,7 @@ void MultiTransport::getMultiDiffCoeffs(const size_t ld, doublereal* const d)
update_C();
// update the binary diffusion coefficients
update_T();
updateDiff_T();
// evaluate L0000 if the temperature or concentrations have
@ -652,6 +655,7 @@ void MultiTransport::updateThermal_T()
}
// we need species viscosities and binary diffusion coefficients
updateSpeciesViscosities();
update_T();
updateDiff_T();
// evaluate polynomial fits for A*, B*, C*