From 04f9596cc57c85505b09f6f33cb400c8f0f43c4f Mon Sep 17 00:00:00 2001 From: BangShiuh Date: Fri, 21 Oct 2016 23:17:11 -0400 Subject: [PATCH] [Transport] replace sumxw with mmw in GasTransport::getMixDiffCoeffs The sumxw is redundant, so I replace it with mmw. --- src/transport/GasTransport.cpp | 6 +----- 1 file changed, 1 insertion(+), 5 deletions(-) diff --git a/src/transport/GasTransport.cpp b/src/transport/GasTransport.cpp index 6122e5070..5375f9ab2 100644 --- a/src/transport/GasTransport.cpp +++ b/src/transport/GasTransport.cpp @@ -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); } } }