[Reactor] Fix net surface mass flux calc in FlowReactor (Cantera#2069)

speth · speth · commit 9666a779571c · 2026-02-20T22:29:52.000-05:00
diff --git a/src/zeroD/FlowReactor.cpp b/src/zeroD/FlowReactor.cpp
@@ -251,7 +251,7 @@ void FlowReactor::evalDae(double time, double* y, double* ydot, double* residual
     const vector<double>& mw = m_thermo->molecularWeights();
     double sk_wk = 0;
     for (size_t i = 0; i < m_nsp; ++i) {
-        sk_wk = m_sdot[i] * mw[i];
+        sk_wk += m_sdot[i] * mw[i];
     }
     m_thermo->getPartialMolarEnthalpies(m_hk.data());
     // get net production

Original file line number	Diff line number	Diff line change
`@@ -251,7 +251,7 @@ void FlowReactor::evalDae(double time, double* y, double* ydot, double* residual`
`251`	`251`	`const vector<double>& mw = m_thermo->molecularWeights();`
`252`	`252`	`double sk_wk = 0;`
`253`	`253`	`for (size_t i = 0; i < m_nsp; ++i) {`
`254`		`- sk_wk = m_sdot[i] * mw[i];`
	`254`	`+ sk_wk += m_sdot[i] * mw[i];`
`255`	`255`	`}`
`256`	`256`	`m_thermo->getPartialMolarEnthalpies(m_hk.data());`
`257`	`257`	`// get net production`