Estimating the joint distribution of rate parameters across multiple reactions in the absence of experimental data

Abstract A procedure for determining the joint uncertainty of Arrhenius parameters across multiple combustion reactions of interest is demonstrated. This approach is capable of constructing the joint distribution of the Arrhenius parameters arising from the uncertain measurements performed in specific target experiments without having direct access to the underlying experimental data. The method involves constructing an ensemble of hypothetical data sets with summary statistics consistent with the available information reported by the experimentalists, followed by a fitting procedure that learns the structure of the joint parameter density across reactions using this consistent hypothetical data as evidence. The procedure is formalized in a Bayesian statistical framework, employing maximum-entropy and approximate Bayesian computation methods and utilizing efficient Markov chain Monte Carlo techniques to explore data and parameter spaces in a nested algorithm. We demonstrate the application of the method in the context of experiments designed to measure the rates of selected chain reactions in the H2-O2 system and highlight the utility of this approach for revealing the critical correlations between the parameters within a single reaction and across reactions, as well as for maximizing consistency when utilizing rate parameter information in predictive combustion modeling of systems of interest.

[1]  David A. Sheen,et al.  The method of uncertainty quantification and minimization using polynomial chaos expansions , 2011 .

[2]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[3]  Habib N. Najm,et al.  Data-free inference of the joint distribution of uncertain model parameters , 2010, J. Comput. Phys..

[4]  Dirk P. Kroese,et al.  Kernel density estimation via diffusion , 2010, 1011.2602.

[5]  Raymond W. Walker,et al.  Evaluated kinetic data for combustion modelling supplement I , 1994 .

[6]  Habib N. Najm,et al.  Data free inference with processed data products , 2016, Stat. Comput..

[7]  H. Najm,et al.  Inference of reaction rate parameters based on summary statistics from experiments , 2017 .

[8]  Yelon,et al.  Microscopic explanation of the compensation (Meyer-Neldel) rule. , 1990, Physical review letters.

[9]  Y. Marzouk,et al.  Uncertainty quantification in chemical systems , 2009 .

[10]  Tamás Varga,et al.  Determination of rate parameters based on both direct and indirect measurements , 2012 .

[11]  Yelon,et al.  Origin and consequences of the compensation (Meyer-Neldel) law. , 1992, Physical review. B, Condensed matter.

[12]  Zhenwei Zhao,et al.  An updated comprehensive kinetic model of hydrogen combustion , 2004 .

[13]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[14]  Tamás Turányi,et al.  Uncertainty of Arrhenius parameters , 2011 .

[15]  H. Haario,et al.  An adaptive Metropolis algorithm , 2001 .

[16]  Tamás Turányi,et al.  Uncertainty of the rate parameters of several important elementary reactions of the H2 and syngas combustion systems , 2015 .