Global sensitivity analysis of n-butanol reaction kinetics using rate rules

We investigate the sensitivity of the ignition delay time, τ_ign, of n-butanol to uncertainties in the rate-rule kinetic rate parameters, at various initial temperatures (600–1000 K), pressures (10–80 bar) and equivalence ratios (0.5–2.0). We start by considering a 30-dimensional stochastic germ in which each random variable is associated with one reaction class, and build a surrogate model for the ignition delay time using polynomial chaos expansions. The adaptive pseudo-spectral projection technique is used for this purpose. The surrogate model is used to estimate first-order and total-order sensitivity indices characterizing the dependence of the ignition delay time on the uncertain inputs. Results indicate that τ_ign is mostly sensitive to variations in four dominant reaction classes, namely, H–atom abstraction from the fuel (reaction class 2), addition of O₂ to the fuel radicals (reaction class 11), fuel radical isomerization including Waddington type reactions (reaction class 15), and concerted elimination reactions (reaction class 16). We consequently focus our attention on these four reaction classes, and consider variations within corresponding subrules. We explore two approaches to define the subrules of reaction class 2, one based on the radical abstracting from the fuel resulting in eleven subrules and another based on the abstraction site resulting in five subrules. Hence, we investigate the sensitivity of τ_ign due to variability in the rate parameters of 26 and 20 subrules of the resulting models. In particular, the simulations indicate that in reaction class 2 H–atom abstraction by HO₂ dominates the variability in τ_ign at all initial conditions considered. Analysis of this finding reveals that correlations inherent in the rate rule construction plays an important role in the resulting sensitivity predictions, and suggests a hierarchical approach to the calibration of elementary reaction rates.