TY - JOUR

T1 - Chemical kinetic uncertainty quantification for Large Eddy Simulation of turbulent nonpremixed combustion

AU - Mueller, Michael E.

AU - Iaccarino, Gianluca

AU - Pitsch, Heinz

N1 - Funding Information:
The authors gratefully acknowledge funding from the Predictive Science Academic Alliance Program (PSAAP) . In addition, M.E.M. gratefully acknowledges funding from the National Defense Science and Engineering Graduate (NDSEG) Fellowship and the National Science Foundation (NSF) Graduate Research Fellowship Program .
Copyright:
Copyright 2017 Elsevier B.V., All rights reserved.

PY - 2013

Y1 - 2013

N2 - While the accuracy of chemical kinetic mechanisms continues to improve, these mechanisms are still models with, sometimes considerable, uncertainty. In order to rigorously validate turbulent combustion simulations against experimental data, this uncertainty must be separated from deficiencies in the turbulent combustion model itself. In this work, a method is developed for quantifying the uncertainty in turbulent flame simulations due to input uncertainty in the chemical mechanism. Here the method is developed for Large Eddy Simulation (LES) combined with a steady flamelet model. Rather than a brute force probabilistic approach in which hundreds or thousands of LES runs are required to compute statistics of outputs of interest, the method takes advantage of the actual algorithm employed with the steady flamelet model. First, the high-dimensional uncertainty in the chemical kinetics is propagated through the flamelet equations, and the resulting lower-dimensional joint distribution of the temperature, species mass fractions, and other derived quantities is used as a stochastic equation of state in the LES. Since only a few 'active' quantities are needed to evolve the LES governing equations, efficient non-intrusive stochastic collocation is used to propagate the uncertainty in the density, requiring only a few LES runs. This process captures the uncertainty in the flow field induced by the uncertainty in the chemical kinetic rates. The remaining uncertainty in 'passive' quantities, that is, quantities needed only for post-processing such as the temperature and species mass fractions, is computed with random sampling during the LES runs. The uncertainty quantification algorithm is demonstrated with Sandia flame D, and it is shown that the uncertainty in the simulation results caused by uncertainties in the kinetic rates is sufficiently large to account for the discrepancies with the experimental measurements. The implication is that the turbulent combustion model cannot be fairly assessed with such a large uncertainty.

AB - While the accuracy of chemical kinetic mechanisms continues to improve, these mechanisms are still models with, sometimes considerable, uncertainty. In order to rigorously validate turbulent combustion simulations against experimental data, this uncertainty must be separated from deficiencies in the turbulent combustion model itself. In this work, a method is developed for quantifying the uncertainty in turbulent flame simulations due to input uncertainty in the chemical mechanism. Here the method is developed for Large Eddy Simulation (LES) combined with a steady flamelet model. Rather than a brute force probabilistic approach in which hundreds or thousands of LES runs are required to compute statistics of outputs of interest, the method takes advantage of the actual algorithm employed with the steady flamelet model. First, the high-dimensional uncertainty in the chemical kinetics is propagated through the flamelet equations, and the resulting lower-dimensional joint distribution of the temperature, species mass fractions, and other derived quantities is used as a stochastic equation of state in the LES. Since only a few 'active' quantities are needed to evolve the LES governing equations, efficient non-intrusive stochastic collocation is used to propagate the uncertainty in the density, requiring only a few LES runs. This process captures the uncertainty in the flow field induced by the uncertainty in the chemical kinetic rates. The remaining uncertainty in 'passive' quantities, that is, quantities needed only for post-processing such as the temperature and species mass fractions, is computed with random sampling during the LES runs. The uncertainty quantification algorithm is demonstrated with Sandia flame D, and it is shown that the uncertainty in the simulation results caused by uncertainties in the kinetic rates is sufficiently large to account for the discrepancies with the experimental measurements. The implication is that the turbulent combustion model cannot be fairly assessed with such a large uncertainty.

KW - Large Eddy Simulation

KW - Sandia flame D

KW - Turbulent nonpremixed flames

KW - Uncertainty quantification

UR - http://www.scopus.com/inward/record.url?scp=84874693948&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84874693948&partnerID=8YFLogxK

U2 - 10.1016/j.proci.2012.07.054

DO - 10.1016/j.proci.2012.07.054

M3 - Conference article

AN - SCOPUS:84874693948

SN - 1540-7489

VL - 34

SP - 1299

EP - 1306

JO - Proceedings of the Combustion Institute

JF - Proceedings of the Combustion Institute

IS - 1

ER -