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

[Transport] replace sumxw with mmw in GasTransport::getMixDiffCoeffs

Browse source 
The sumxw is redundant, so I replace it with mmw.
This commit is contained in:
BangShiuh 2016-10-21 23:17:11 -04:00 • committed by Ray Speth
parent 3f3b96587b
commit 04f9596cc5
1 changed files with 1 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/transport/GasTransport.cpp
View file

@ -255,14 +255,10 @@ void GasTransport::getMixDiffCoeffs(doublereal* const d)
}
doublereal mmw = m_thermo->meanMolecularWeight();
doublereal sumxw = 0.0;
doublereal p = m_thermo->pressure();
if (m_nsp == 1) {
d[0] = m_bdiff(0,0) / p;
} else {
for (size_t k = 0; k < m_nsp; k++) {
sumxw += m_molefracs[k] * m_mw[k];
}
for (size_t k = 0; k < m_nsp; k++) {
double sum2 = 0.0;
for (size_t j = 0; j < m_nsp; j++) {
@ -273,7 +269,7 @@ void GasTransport::getMixDiffCoeffs(doublereal* const d)
if (sum2 <= 0.0) {
d[k] = m_bdiff(k,k) / p;
} else {
d[k] = (sumxw - m_molefracs[k] * m_mw[k])/(p * mmw * sum2);
d[k] = (mmw - m_molefracs[k] * m_mw[k])/(p * mmw * sum2);
}
}
}