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

[Transport] fix temperature update in IonGasTransport::getMobilities

Browse source
This commit is contained in:
bangshiuh 2018-10-18 20:20:40 -04:00 • committed by Ray Speth
parent 87b042a231
commit 8e7e27ad62
2 changed files with 9 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
interfaces/cython/cantera/test/test_onedim.py
View file

@ -961,7 +961,7 @@ class TestIonFreeFlame(utilities.CanteraTest):
self.sim.solve(loglevel=0, stage=2, enable_energy=True)
# Regression test
self.assertNear(max(self.sim.E), 136.1234, 1e-3)
self.assertNear(max(self.sim.E), 142.2677, 1e-3)
class TestIonBurnerFlame(utilities.CanteraTest):
@ -986,4 +986,4 @@ class TestIonBurnerFlame(utilities.CanteraTest):
self.sim.solve(loglevel=0, stage=2, enable_energy=True)
# Regression test
self.assertNear(max(self.sim.E), 539.59, 1e-2)
self.assertNear(max(self.sim.E), 552.33, 1e-2)

7
src/transport/IonGasTransport.cpp
View file

@ -355,6 +355,13 @@ void IonGasTransport::getMixDiffCoeffs(double* const d)
void IonGasTransport::getMobilities(double* const mobi)
{
update_T();
update_C();
// update the binary diffusion coefficients if necessary
if (!m_bindiff_ok) {
updateDiff_T();
}
double p = m_thermo->pressure();
for (size_t k = 0; k < m_nsp; k++) {
if (k == m_kElectron) {