Determination of rate parameters based on both direct and indirect measurements

Abstract

The determination of rate parameters of gas-phase elementary reactions is usually based on direct measurements. The rate parameters obtained in many independent direct measurements are then used in reaction mechanisms, which are tested against the results of indirect experiments, like time-to-ignition or laminar flame velocity measurements. We suggest a new approach that takes into account both direct and indirect measurements and optimizes all influential rate parameters. First, the domain of feasibility of the Arrhenius parameters is determined from all of the available direct measurements. Thereafter, the optimal Arrhenius parameters are sought within this domain to reproduce the selected direct and indirect measurements. Other parameters of a complex mechanism (third-body efficiencies, enthalpies of formation, parameters of pressure dependence, etc.) can also be taken into account in a similar way. A new fitting algorithm and a new method for error calculation were developed to determine the optimal mean values and the covariance matrix of all parameters. The approach is demonstrated on the calculation of Arrhenius parameters of reactions (R1): H + O 2 = OH + O and (R2): H + O 2 + M = HO 2 + M (low-pressure limit, M = N 2 or Ar). In total, 9 direct measurements for reaction (R1) (745 data points), 10 direct measurements for reaction (R2) (258 data points), and 11 ignition time measurements (79 data points) were taken into account. The application of the method resulted in the following rate parameters for the investigated reactions-(R1): A = 3.003 × 10 10 cm 3 mol −1 s −1 , n = 0.965, E/R = 6158 K (T = 950-3550 K) DETERMINATION OF RATE PARAMETERS BASED ON DIRECT AND INDIRECT MEASUREMENTS 285and (R2): A = 7.856 × 10 18 cm 6 mol −2 s −1 , n = −1.100, E/R = 0 K (low-pressure limit, M = N 2 , T = 300-1850 K). The optimized third-body efficiency of Ar relative to N 2 is m = 0.494 (standard deviation σ = 0.010). The uncertainty parameter f as a function of temperature was also calculated. Average uncertainty parameter values are f = 0.025 and 0.049 for reactions (R1) and (R2) (corresponding to 6% and 12%), respectively, which are much lower than those of the previous evaluations. C 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 284-302, 2012

The effect of temperature on the adiabatic burning velocities of diluted hydrogen flames: A kinetic study using an updated mechanism

“…Improved presentation of the uncertainty information has been discussed by Nagy and Turanyi [60] and others. It was demonstrated [61] that dedicated optimization procedures may bring significantly reduced uncertainties varying with temperature and lead to an improvement in predicting ability [62].…”

Comparison of the performance of several recent syngas combustion mechanisms

“…Error function values E i and E are expected to be near unity if the chemical kinetic model is accurate, and deviations of the measured and simulated results are caused by the scatter of the experimental data only. Note that due to the squaring in the definition of E, a twice as high deviation of the simulated and experimental values of one mechanism in comparison to another leads to a four times higher value of E. This objective function has been used in our previous studies on the comparison of reaction mechanisms [2] and the estimation of rate parameters from experimental data [3][4][5][6].…”

Global Uncertainty Propagation and Sensitivity Analysis in the CH3OCH2 + O2 System: Combining Experiment and Theory To Constrain Key Rate Coefficients in DME Combustion

“…The analysis supports previous findings by Nagy and Turányi concerning the treatment of correlations in the uncertainty parameters for Arrhenius expressions. [41][42][43] The correlations will also affect the interpretation of the sensitivity information as highlighted above. A high total sensitivity to E(RO 2 ) or E(TS2) will lead to a high sensitivity in the other because of their high correlation.…”