Rate coefficient parameter changes are noted on the individual reaction Web pages. The CH kinetics important to prompt NO formation were altered in light of new measurements. New expressions were also used for the H + O2 reactions, CH3 + O2, CH2O + H and CH2O decomposition. The methanol decomposition/chemical activation system was recomputed. The oxidation steps CH3 + O and CH2 + O2 feature new branching paths.
Two main additions were made to the kinetics mechanism, with the addition of 4 species. Acetaldehyde and vinoxy chemistry are included to better describe ethylene oxidation, and this path is now included among the Ox + C2Hy reaction products. Because natural gas contains propane (and some higher hydrocarbons that may be approximately represented by propane), a minimal set of propane kinetics is included to model this species, as a minor constituent only.
Other new target additions include a series of shock tube observations sensitive to the oxidation of the formaldehyde intermediate; a set of shock tube, low pressure flame, and flow reactor experiments concerning prompt NO formation and reburn; and a few targets concerning the shortening of methane shock tube ignition delays by small amounts of propane or ethane. The sets of shock tube ignition delays and laminar flame speeds were revised and expanded. Many of these changes reflect the acquisition of improved data.
The new mechanism contains 325 reactions (3 are duplicates because the sum of two rate parameter expressions is required) and 53 species (including argon). The final optimization to 77 targets altered 31 rate parameters.
The major results of the new mechanism optimization are:
1. Deviations from the target values are generally less than previously. (The chi square for the sum of common retained targets dropped from 0.96 to 0.81.)
2. Similar final values for the key rates CH3 + H, CH3 + OH, and CH3 + O2 were found.
3. Including formaldehyde target optimization did not require changes to the new expressions for the sensitive CH2O + M and CH2O + H rate constants. Other changes consistent with the rest of the targets sufficed.
4. Significant thermodynamics sensitivity was found only for HCN in some targets, leading to an alteration of the JANAF value.
5. A new, improved low pressure flame prompt NO target results in higher values for CH + N2, and predictions of increased prompt NO.
6. New lower experimental flame speed values remain overpredicted by the optimized mechanism. We have not examined transport uncertainties yet in this regard. Note also in making comparisons that some computational inaccuracies exist in the flame speed computations used for the GRI-Mech 1.2 target optimization.
Some details on the new optimization are provided.
grimech30.dat | A reaction mechanism and rate coefficient file, in Chemkin format |
thermo30.dat | A thermochemical data file to be used with grimech30.dat, as NASA polynomial coefficients |
transport.dat | A file containing the parameters needed for calculating transport coefficients to be used in the Sandia flame code |
bugfix.dat | A file containing selected user questions, comments, and suggestions related to the implementation of GRI-Mech |
SHOCK-TUBE IGNITION DELAY MEASUREMENTS | ||
---|---|---|
IG.1a | Seery, D.J., and Bowman, C.T. (1970) Combust.
Flame 14, 37.
CH4-O2-Ar
|
|
IG.1b | ||
IG.2 | Seery, D.J., and Bowman, C.T. (1970) Combust. Flame
14, 37.
CH4-O2-Ar
|
|
IG.6a | Hidaka, Y., Gardiner, W.C., and Eubank, C.S.
(1982) J. of Molec.Sci. (China) 2, p.141-53.
C2H6-O2-Ar
|
|
IG.6b | ||
IG.T1 | Frenklach, M., and Bornside, D.E., (1984) Combust. Flame
56, 1
CH4-C3H8-O2-Ar (9.5%-1.9%-19%-69.6%)
|
|
IG.T2 | Spadaccini, L.J., and Colket, M.B., (1994) Prog. Energy
Combust. Sci. 20, 431.
CH4-C3H8-O2-Ar (3.4%-0.1%-7%-89.5%)
|
|
IG.St1a | Spadaccini, L.J., and Colket, M.B., (1994)
Prog. Energy Combust. Sci. 20, 431.
CH4-C2H6-O2-Ar (3.29%-0.21%-7%-89.5%)
|
|
IG.St1b | ||
IG.St3a | Seery, D.J., and Bowman, C.T., (1970) Combust.
Flame 14, 37.
CH4-O2-Ar (4.8%-19.1%-76.1%)
|
|
IG.St3b | ||
IG.St4a | Peterson, E.L., Davidson, D.F., Rohrig, M.,
Hanson, R.K., and Bowman, C.T. (1995), to be published.
CH4-O2-Ar, phi = 1.0
|
|
IG.St4b | ||
SHOCK-TUBE SPECIES PROFILE MEASUREMENTS | ||
CH3.C1a | Maximum CH3 concentration | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K.,
and Bowman, C.T. (1994) 25th Symposium (International) on Combustion,
Poster 3 - 23.
CH4-O2-Ar (994 ppm-0.2021%-99.7%)
|
CH3.C1b | ||
CH3.T1a | Time of CH3 maximum | |
CH3.T1b | ||
CH3.C2 | Maximum CH3 concentration | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K.,
and Bowman, C.T. (1994) 25th Symposium (International) on Combustion,
Poster 3 - 23.
C2H6-O2-Ar (295 ppm-0.1055%-99.865%)
|
CH3.T2 | Time of CH3 maximum | |
CH3.C3 | Maximum CH3 concentration | Davidson, D.F., DiRosa, M., Chang, A.Y., Hanson, R.K.,
and Bowman, C.T. (1992) 24th Symposium (International) on Combustion,
p. 589.
C2H6-Ar (206 ppm)
|
CH3.T3 | Time of CH3 maximum | |
CH3.C4 | Maximum CH3 concentration | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K.,
and Bowman, C.T. (1994) 25th Symposium (International) on Combustion,
Poster 3 - 23.
CH4-O2-Ar (0.2021%-0.1%-99.7%)
|
CH3.T4 | Time of CH3 maximum | |
CH3.StC6 | Time of CH3 maximum | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K.,
and Bowman, C.T. (1994) 25th Symposium (International) on Combustion,
Poster 3 - 23.
CH4-O2-Ar (0.1%-0.4%-99.5%)
|
CH3.StC7 | ||
OH.1a | Time to half of OH maximum | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K.,
and Bowman, C.T. (1994) 25th Symposium (International) on Combustion,
Poster 3 - 23.
CH4-O2-Ar (0.1%-0.2%-99.7%)
|
OH.1b | ||
OH.2 | Chang, A.Y., Davidson, D.F., DiRosa, M., Hanson, R.K., and Bowman,
C.T. (1994) 25th Symposium (International) on Combustion, Poster
3 - 23.
C2H6-O2-Ar (300ppm-0.105%-99.89%)
|
|
OH.3a | Time to half of OH maximum | Yu, C.L., Wang, C., and Frenklach, M. (1995) J.Phys.Chem.
99, 14377.
CH4-O2-Ar, phi = 0.04, P = 1.51 atm, T = 1750 K |
CO.C1a | Maximum CO concentration | |
CO.T1a | Time to half of CO maximum | |
OH.3b | Time to half of OH maximum | Yu, C.L., Wang, C., and Frenklach, M. (1995) J.Phys.Chem.
99, 14377.
CH4-O2-Ar, phi = 0.04, P = 1.64 atm, T = 1900 K |
CO.C1b | Maximum CO concentration | |
CO.T1b | Time to half of CO maximum | |
OH.3c | Time to half of OH maximum | Yu, C.L., Wang, C., and Frenklach, M. (1995) J.Phys.Chem.
99, 14377.
CH4-O2-Ar, phi = 0.667, P = 1.51 atm, T = 1750 K |
CO.C1c | Maximum CO concentration | |
CO.T1c | Time to half of CO maximum | |
OH.3d | Time to half of OH maximum | Yu, C.L., Wang, C., and Frenklach, M. (1995) J.Phys.Chem.
99, 14377.
CH4-O2-Ar, phi = 0.667, P = 1.64 atm, T = 1900 K |
CO.C1d | Maximum CO concentration | |
CO.T1d | Time to half of CO maximum | |
OH.ST8 | Time to half of OH maximum | Yu, C.L., Wang, C., and Frenklach, M.,
'Chemical Kinetics of Methyl Oxidation by Molecular Oxygen,'
(1995) J.Phys.Chem. 99, 14377.
CH4-O2-Ar (1%-1.5%-97.5%)
|
CO.ST8 | Time to half of CO maximum | |
CO.SC8 | Maximum CO concentration | |
BCO.T1 | Time to half of CO maximum | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-Ar P = 1.17 atm, T = 2124 K |
BCO.T2 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-Ar P = 1.51 atm, T = 1724 K |
|
BCO.T3 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-O2-Ar, phi = 5.88 P = 2.32 atm, T = 1784 K |
|
BCO.T4 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-O2-Ar, phi = 0.25 P = 1.89 atm, T = 1442 K |
|
BCO.T5 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-O2-Ar, phi = 1.68 P = 0.911 atm, T = 1768 |
|
BCO.T6 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-O2-Ar, phi = 0.17 P = 2 atm, T = 1515 K |
|
BCO.T7 | Eiteneer, B., Yu, C.-L., Goldenberg, M., and Frenklach,
M., J. Phys. Chem. 102, 5196, (1998).
CH2O-O2-Ar, phi = 1 P = 1.51 atm, T = 1720 K |
|
BCH2O.T1 | Time to half of CO maximum | Hidaka, Y., Taniguchi, T., Tanaka, H., Kamesawa, T., Inami,
K., and Kawano, H. (1993) Combust. Flame 92, 365.
CH2O-O2-Ar, phi = 2.0 - 4.0
|
BCH2O.T2 | ||
BCH2O.T3 | ||
REACTORS | ||
SR.10c | Temperature where [CO2] = 500 ppm | Kristensen, P.G., Glarborg, P., and Dam-Johansen, K. (1995), unpublished data. |
LAMINAR FLAME SPEED | ||
F1 | Vagelopoulos, C.N., and Egolfopoulos,F.N., 'Direct Experimental
Determination of Laminar Flame Speeds,' Paper WSS/CI 97S-022, Western States
Section/Combustion Institute Meeting, Livermore, CA, April 1997.
CH4-air, T(0) = 300 K
|
|
F2 | Vagelopoulos, C.M., Egolfopoulos, F.N., and
Law, C.K., 'Further Considerations on the Determination of Laminar Flame
Speeds with the Counterflow Twin-Flame Technique,' (1994) 25th Symposium
(International) on Combustion, p. 1341.
CH4-air, T(0) = 300 K
|
|
F3 | ||
F4 | Egolfopolous, F.N., Cho, P., and Law, C.K., 'Laminar Flame
Speeds of Methane-Air Mixtures under Reduced and Elevated Pressures,' (1989)
Combust. Flame 76, 375.
CH4-air, T(0) = 300 K
|
|
F5 | Just, Th. (1994) Private communication.
CH4-air, T(0) = 400 K
|
|
F6 | ||
SF7 | McLean, I.C., Smith, D.B., and Taylor, S.B., (1994) 25th
Symposium (International) on Combustion, p. 749.
CO-H2-air, (20.8%, 20.8%, 58.4%) T(0) = 300 K
|
|
StF8 | References: Vagelopoulos, C.N., and Egolfopoulos,F.N.,
'Direct Experimental Determination of Laminar Flame Speeds,' Paper WSS/CI
97S-022, Western States Section/Combustion Institute Meeting, Livermore,
CA, April 1997.
C2H6-air, T(0) = 300 K
|
|
PROMPT NO | ||
SNO.C11 | NO (2 cm) concentration | Luque, J., Smith, G.P., and Crosley, D.R. (1996) 26th Symposium
(International) on Combustion, 959.
CH4-O2-N2 (13.8%-25.9%-60.3%)
|
SCH.C11 | Maximum CH concentration | Luque, J., Smith, G.P., and Crosley, D.R. (1996) 26th Symposium
(International) on Combustion, 959
CH4-O2-N2 (13.8%-25.9%-60.3%)
|
SCH.C12 | Maximum CH concentration | Berg, P.A., Hill, D.A., Noble, A.R., Smith, G.P., Jeffries,
J.B., and Crosley, D.R., 'Absolute CH Concentration Measurements in Low
Pressure Hydrocarbon Flames: Comparison with Model Predictions,' (1997)
35th Aerospace Sciences Meeting, Paper 97-0905, Reno.
CH4-O2-N2
|
SCH.C13 | ||
CH.St | Maximum of CH concentration without NO doping | Woiki, D., Votsmeier, M., Davidson, D.F., Hanson, R.K.,and Bowman,
C.T., 'CH-Radical Concentration Measurements in Fuel Rich CH4/O2/Ar and
CH4/O2/NO/Ar Mixtures Behind Shock Waves', (1998) Combust. Flame
113, 624.
CH4-O2-Ar
|
HCN OXIDATION | ||
NFR1 | HCN relative concentration | Flow Reactor: Glarborg, P., and Miller, J.A. (1994) Combust.
Flame 99, 475.
HCN-CO-O2-H2O-N2
|
NFR2 | NO relative concentration | |
NFR3 | N2O relative concentration | |
NF6 | Maximum NO mole fraction | Sandia Flame: Miller, J.A., et al. (1984) 20th Symposium
(International) on Combustion, p. 673.
H2-O2-HCN-Ar
|
NF7 | Maximum CN mole fraction | |
REBURNING | ||
NF11 | Ratio of CH maximum concentrations at two levels of NO doping | NRL Flame: Williams, B.A., and Fleming, J.W. (1994) Combust.
Flame 98, 93.
CH4-O2-Ar-NO or N2O (1.2%)
|
NF12/13 | Ratio of CN maximum concentrations for NO and N2O doping | |
NFR4 | NO mole fraction at reactor exit | Flow Reactor: Alzueta, M.U., Glarborg, P., and Dam-Johansen,
K., 'Low Temeprature Interactions between Hydrocarbons and Nitric Oxide.
An Experimental Study,' Combust. Flame 109, 25.
CH4-C2H6-O2-NO-H2O-N2 (2864 ppm - 298 ppm - 5090 ppm - 947 ppm
- 2.16% - 96.92%)
|
NFR5 | HCN mole fraction at reactor exit | |
CHNO.St | Maximum of CH concentration with NO doping | Woiki, D., Votsmeier, M., Davidson, D.F., Hanson, R.K.,and Bowman,
C.T., 'CH-Radical Concentration Measurements in Fuel Rich CH4/O2/Ar and
CH4/O2/NO/Ar Mixtures Behind Shock Waves', (1998) Combust. Flame
113, 624.
CH4-O2-NO-Ar
|
% CH4 | % O2 | % C2H6 | Density
(10-5 mol/cm3) |
Temperature (K) | Ignition Delay |
3.29 | 7.0 | 0.21 | 5.5 | 1354 - 1698 | Graph
(Targets IG.St1a and IG.St1b) |
3.40 | 7.0 | 0.10 | 5.5 | 1418 - 1742 | Graph |
3.46 | 15.25 | 3.3, 6.3 | 1336 - 1679 | Graph | |
3.50 | 5.6 | 2.2, 4.2 | 1461 - 2025 | Graph | |
3.50 | 7.0 | 2.7, 4.9 | 1516 - 1918 | Graph | |
3.50 | 7.0 | 0.04
0.35 |
5.5 | 1682 - 1723
1445 - 1669 |
Graph |
3.50 | 9.3 | 3.1, 5.8 | 1479 - 1759 | Graph | |
6.00 | 12.0 | 3.9, 8.0 | 1464 - 1708 | Graph |
% CH4 | % O2 | Density (10-5 mol/cm3) | Temperature (K) | Profile Graphs |
0.4 | 5.0 | 1.04 | 1821 | CO,OH |
0.4 | 5.0 | 1.58 | 1941 | CO,OH |
0.4 | 20 | 1.04 | 1711 | CO,OH |
0.4 | 20 | 1.07 | 1944 | CO,OH |
0.5 | 10 | 1.05 | 1752 | CO,OH |
0.5 | 10 | 1.06 | 1843 | CO,OH |
1.0 | 1.5 | 1.56 | 2025 | CO,OH |
1.0 | 1.5 | 1.60 | 1865 | CO,OH
(Targets CO.SC8, C0.SC8, and OH.ST8) |
1.0 | 2.0 | 1.06 | 2040 | CO,OH |
1.0 | 2.0 | 1.57 | 1800 | CO,OH |
1.0 | 3.0 | 1.07 | 1826 | CO,OH |
1.0 | 3.0 | 1.58 | 1856 | CO,OH |
Comments on GRI-Mech 1.2 flame speed computations.
An NO concentration of 1.8+/-0.6E12 cm-3 was measured at 2 cm, while GRI-Mech 3.0 predicts 1.3E12 cm-3.