Direct Numerical Simulation of Head-On Quenching of Statistically Planar Turbulent Premixed Methane-Air Flames Using a Detailed Chemical Mechanism

A three-dimensional compressible Direct Numerical Simulation (DNS) analysis has been carried out for head-on quenching of a statistically planar stoichiometric methane-air flame by an isothermal inert wall. A multi-step chemical mechanism for methane-air combustion is used for the purpose of detailed chemistry DNS. For head-on quenching of stoichiometric methane-air flames, the mass fractions of major reactant species such as methane and oxygen tend to vanish at the wall during flame quenching. The absence of \(\text {OH}\) at the wall gives rise to accumulation of carbon monoxide during flame quenching because \(\text {CO}\) cannot be oxidised anymore. Furthermore, it has been found that low-temperature reactions give rise to accumulation of \(\text {HO}_{2}\) and \(\mathrm {H}_{2}\mathrm {O}_{2}\) at the wall during flame quenching. Moreover, these low temperature reactions are responsible for non-zero heat release rate at the wall during flame-wall interaction. In order to perform an in-depth comparison between simple and detailed chemistry DNS results, a corresponding simulation has been carried out for the same turbulence parameters for a representative single-step Arrhenius type irreversible chemical mechanism. In the corresponding simple chemistry simulation, heat release rate vanishes once the flame reaches a threshold distance from the wall. The distributions of reaction progress variable c and non-dimensional temperature T are found to be identical to each other away from the wall for the simple chemistry simulation but this equality does not hold during head-on quenching. The inequality between c (defined based on \(\text {CH}_{4}\) mass fraction) and T holds both away from and close to the wall for the detailed chemistry simulation but it becomes particularly prominent in the near-wall region. The temporal evolutions of wall heat flux and wall Peclet number (i.e. normalised wall-normal distance of \(T = 0.9\) isosurface) for both simple and detailed chemistry laminar and turbulent cases have been found to be qualitatively similar. However, small differences have been observed in the numerical values of the maximum normalised wall heat flux magnitude \(\left ({\Phi }_{\max } \right )_{\mathrm {L}}\) and the minimum Peclet number \((Pe_{\min })_{\mathrm {L}}\) obtained from simple and detailed chemistry based laminar head-on quenching calculations. Detailed explanations have been provided for the observed differences in behaviours of \(\left ({\Phi }_{\max }\right )_{\mathrm {L}}\) and \((Pe_{\min })_{\mathrm {L}}\). The usual Flame Surface Density (FSD) and scalar dissipation rate (SDR) based reaction rate closures do not adequately predict the mean reaction rate of reaction progress variable in the near-wall region for both simple and detailed chemistry simulations. It has been found that recently proposed FSD and SDR based reaction rate closures based on a-priori DNS analysis of simple chemistry data perform satisfactorily also for the detailed chemistry case both away from and close to the wall without any adjustment to the model parameters.