Bayesian Error Propagation for a Kinetic Model of n‐Propylbenzene Oxidation in a Shock Tube |
| |
Authors: | Sebastian Mosbach Je Hyeong Hong George P E Brownbridge Markus Kraft Soumya Gudiyella Kenneth Brezinsky |
| |
Institution: | 1. Department of Chemical Engineering and Biotechnology, University of Cambridge, , Cambridge, CB2 3RA United Kingdom;2. Department of Engineering, University of Cambridge, , Cambridge, CB2 1PZ United Kingdom;3. Department of Chemical Engineering, University of Illinois at Chicago, , Chicago, 60607;4. Department of Mechanical and Industrial Engineering, University of Illinois at Chicago, , Chicago, 60607 |
| |
Abstract: | We apply a Bayesian parameter estimation technique to a chemical kinetic mechanism for n‐propylbenzene oxidation in a shock tube to propagate errors in experimental data to errors in Arrhenius parameters and predicted species concentrations. We find that, to apply the methodology successfully, conventional optimization is required as a preliminary step. This is carried out in two stages: First, a quasi‐random global search using a Sobol low‐discrepancy sequence is conducted, followed by a local optimization by means of a hybrid gradient‐descent/Newton iteration method. The concentrations of 37 species at a variety of temperatures, pressures, and equivalence ratios are optimized against a total of 2378 experimental observations. We then apply the Bayesian methodology to study the influence of uncertainties in the experimental measurements on some of the Arrhenius parameters in the model as well as some of the predicted species concentrations. Markov chain Monte Carlo algorithms are employed to sample from the posterior probability densities, making use of polynomial surrogates of higher order fitted to the model responses. We conclude that the methodology provides a useful tool for the analysis of distributions of model parameters and responses, in particular their uncertainties and correlations. Limitations of the method are discussed. For example, we find that using second‐order response surfaces and assuming normal distributions for propagated errors is largely adequate, but not always. |
| |
Keywords: | |
|
|