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

[Kinetics] Ensure that equilibrium constants are not zero

Browse source 
This prevents reverse rate constants from becoming NaN when the forward rate
constant underflows.
This commit is contained in:
Ray Speth 2013-02-15 17:31:15 +00:00
parent 3aa84e15e9
commit 8847b5abac
1 changed files with 2 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
src/kinetics/GasKinetics.cpp
View file

@ -232,7 +232,8 @@ void GasKinetics::updateKc()
doublereal rrt = 1.0/(GasConstant * thermo().temperature());
for (size_t i = 0; i < m_nrev; i++) {
size_t irxn = m_revindex[i];
m_rkcn[irxn] = exp(m_rkcn[irxn]*rrt - m_dn[irxn]*m_logStandConc);
m_rkcn[irxn] = std::min(exp(m_rkcn[irxn]*rrt - m_dn[irxn]*m_logStandConc),
BigNumber);
}
for (size_t i = 0; i != m_nirrev; ++i) {