A Bayesian approach to calibrating hydrogen flame kinetics using many experiments and parameters

TL;DR

This paper employs Bayesian Markov Chain Monte Carlo sampling to quantify and propagate uncertainties in hydrogen flame kinetics parameters, linking experimental data to predictive combustion modeling.

Contribution

It introduces a Bayesian calibration framework for hydrogen kinetics parameters using extensive experimental data and demonstrates its impact on complex flame predictions.

Findings

Posterior distribution reveals strong parameter correlations.

Uncertainty in ignition time exceeds 10% for hydrogen jet ignition.

Calibration constrains predictions in different combustion regimes.

Abstract

First-principles Markov Chain Monte Carlo sampling is used to investigate uncertainty quantification and uncertainty propagation in parameters describing hydrogen kinetics. Specifically, we sample the posterior distribution of thirty-one parameters focusing on the H2O2 and HO2 reactions resulting from conditioning on ninety-one experiments. Established literature values are used for the remaining parameters in the mechanism. The samples are computed using an affine invariant sampler starting with broad, noninformative priors. Autocorrelation analysis shows that O(1M) samples are sufficient to obtain a reasonable sampling of the posterior. The resulting distribution identifies strong positive and negative correlations and several non-Gaussian characteristics. Using samples drawn from the posterior, we investigate the impact of parameter uncertainty on the prediction of two more complex flames: a 2D premixed flame kernel and the ignition of a hydrogen jet issuing into a heated chamber. The former represents a combustion regime similar to the target experiments used to calibrate the mechanism and the latter represents a different combustion regime. For the premixed flame, the net amount of product after a given time interval has a standard deviation of less than 2% whereas the standard deviation of the ignition time for the jet is more than 10%. The samples used for these studies are posted online. These results indicate the degree to which parameters consistent with the target experiments constrain predicted behavior in different combustion regimes. This process provides a framework for both identifying reactions for further study from candidate mechanisms as well as combining uncertainty quantification and propagation to, ultimately, tie uncertainty in laboratory flame experiments to uncertainty in end-use numerical predictions of more complicated scenarios.