Skip to content

Commit e959180

Browse files
committed
Delete comments and fix typo in Srad Disproportionation and R_Recombination reactions
Found bug after getting Ea = -50, causing system to be too stiff and crash. We found that the alpha and n kinetics values in Disproportionation were switched due to a typo. The form of the kinetics should not be Evans-Polanyi. We also found that the MRH comments in the kinetics descriptions are no longer relevant as it corresponds to a previous version of the rate rule. The rate rule in the database was constructed by AG Vandeputte.
1 parent 7a4207d commit e959180

2 files changed

Lines changed: 5 additions & 39 deletions

File tree

input/kinetics/families/Disproportionation/rules.py

Lines changed: 5 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -855,8 +855,8 @@
855855
label = "S_rad/NonDeS;C/H2/Nd_Csrad/H/Cd",
856856
kinetics = ArrheniusEP(
857857
A = (6.44e+08, 'cm^3/(mol*s)'),
858-
n = 0,
859-
alpha = 1.19,
858+
n = 1.19,
859+
alpha = 0.0,
860860
E0 = (0.51, 'kcal/mol'),
861861
Tmin = (300, 'K'),
862862
Tmax = (1500, 'K'),
@@ -865,15 +865,7 @@
865865
shortDesc = u"""Very rough based on R_Recomb #491""",
866866
longDesc =
867867
u"""
868-
[91] Tsang, W.; Journal of Physical and Chemical Reference Data (1988), 17(2), 887-951.
869-
Literature review. C2H3 + n-C3H7 --> C3H6 + C2H4
870-
871-
pg. 922: Discussion on evaluated data
872-
873-
Entry 41,19 (a): No data available at the time. Author estimates the rate coefficient
874-
875-
based on the rxn C2H5+n-C3H7=C3H6=C2H6.
876-
MRH 30-Aug-2009
868+
877869
""",
878870
)
879871

@@ -910,8 +902,8 @@
910902
label = "S_rad/NonDeS;S_Csrad",
911903
kinetics = ArrheniusEP(
912904
A = (6.44e+08, 'cm^3/(mol*s)'),
913-
n = 0,
914-
alpha = 1.19,
905+
n = 1.19,
906+
alpha = 0.0,
915907
E0 = (0.51, 'kcal/mol'),
916908
Tmin = (300, 'K'),
917909
Tmax = (1500, 'K'),
@@ -920,16 +912,6 @@
920912
shortDesc = u"""Very rough based on R_Recomb #491""",
921913
longDesc =
922914
u"""
923-
[91] Tsang, W.; Journal of Physical and Chemical Reference Data (1988), 17(2), 887-951.
924-
Literature review. C2H + n-C3H7 --> C3H6 + C2H2
925-
926-
pg. 923: Discussion on evaluated data
927-
928-
Entry 41,21 (a): No data available at the time. Author notes that the rxn is more exothermic
929-
930-
than the rxn CH3+n-C3H7=C3H6+CH4 and suggests a rate coefficient 3x larger,
931-
namely 1.0x10^-11 cm3/molecule/s.
932-
MRH 30-Aug-2009
933915
""",
934916
)
935917

input/kinetics/families/R_Recombination/rules.py

Lines changed: 0 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -2156,22 +2156,6 @@
21562156
shortDesc = u"""A.G. Vandeputte""",
21572157
longDesc =
21582158
u"""
2159-
MRH estimate
2160-
2161-
A reasonable estimate for the total k_inf(T) for the recombination of H radical with a heavy atom
2162-
is a temperature-independent 1e+14 cm3 mol-1 s-1. HOWEVER, the value I choose to store in the database
2163-
is 1e+13 cm3 mol-1 s-1, because this is the single-event value. Not knowing what species RMG will find
2164-
in trying this estimate, I want to fail on the low side.
2165-
2166-
Using a hydrogen on a primary carbon as an example: Ethane has a multiplicity of six, iso-butane has a
2167-
multiplicity of nine, and neo-pentane has a multiplicity of 12. Using the 1e+13 cm3 mol-1 s-1 with any
2168-
of these multiplicities will not result in ridiculously fast kinetics. The purpose of adding entries
2169-
491 and 492 are to reduce the chance RMG sends ridiculously fast high-P-limit kinetics to fame, thereby
2170-
giving us ridiculously fast k(T,P) in our chem.inp files and causing stiffness issues in flame solvers.
2171-
2172-
NOTE TO RMG USERS: If your model proves to be sensitive to the kinetics of the H+R(+M)=H-R(+M), I would
2173-
encourage you to run a fame job separately, with the 1e+14 cm3 mol-1 s-1 as the total k(T) (if no better
2174-
estimate is known).
21752159
""",
21762160
)
21772161

0 commit comments

Comments
 (0)