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

[Kinetics] Fixed a NaN result in falloff reaction rates

Browse source 
A NaN result could be caused by an underflow in the high-pressure rate
expression.
This commit is contained in:
Ray Speth 2013-05-06 20:05:57 +00:00
parent 259ceaa9a6
commit 2e26abd4f9
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/kinetics/GasKinetics.cpp
View file

@ -340,7 +340,7 @@ void GasKinetics::processFalloffReactions()
vector_fp& pr = m_ropr;
for (size_t i = 0; i < m_nfall; i++) {
pr[i] = concm_falloff_values[i] * m_rfn_low[i] / m_rfn_high[i];
pr[i] = concm_falloff_values[i] * m_rfn_low[i] / (m_rfn_high[i] + SmallNumber);
AssertFinite(pr[i], "GasKinetics::processFalloffReactions",
"pr[" + int2str(i) + "] is not finite.");
}