Statistical Analysis of Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Abstract

Statistically planar turbulent premixed and partially premixed flames for different initial turbulence intensities are simulated for global equivalence ratios = 0.7 and = 1.0 using three-dimensional Direct Numerical Simulations (DNS) with simplified chemistry. For the simulations of partially premixed flames, a random distribution of equivalence ratio following a bimodal distribution of equivalence ratio is introduced in the unburned reactants ahead of the flame. The simulation parameters in all of the cases were chosen such that the combustion situation belongs to the thin reaction zones regime. The DNS data has been used to analyse the behaviour of the dissipation rate transports of both active and passive scalars (i.e. the fuel mass fraction Y F and the mixture fraction ξ) in the context of Reynolds Averaged Navier–Stokes (RANS) simulations. The behaviours of the unclosed terms of the Favre averaged scalar dissipation rates of fuel mass fraction and mixture fraction (i.e. $\widetilde {\varepsilon }_Y =\overline {\rho D\nabla Y_F^{\prime \prime } \cdot \nabla Y_F^{\prime \prime } } /\overline{\rho }$ and $\widetilde {\varepsilon }_\xi =\overline {\rho D\nabla \xi ^{\prime \prime }\cdot \nabla \xi ^{\prime \prime }} /\overline {\rho })$ transport equations have been analysed in detail. In the case of the $\widetilde {\varepsilon }_Y $ transport, it has been observed that the turbulent transport term of scalar dissipation rate remains small throughout the flame brush whereas the terms due to density variation, scalar–turbulence interaction, reaction rate and molecular dissipation remain the leading order contributors. The term arising due to density variation remains positive throughout the flame brush and the combined contribution of the reaction and molecular dissipation to the $\widetilde {\varepsilon }_Y $ transport remains negative throughout the flame brush in all cases. However, the behaviour of scalar–turbulence interaction term of the $\widetilde {\varepsilon }_Y $ transport equation is significantly affected by the relative strengths of turbulent straining and the straining due to chemical heat release. In the case of the $\widetilde {\varepsilon }_\xi $ transport, the turbulent transport term remains small throughout the flame brush and the density variation term is found to be negligible in all cases, whilst the reaction rate term is exactly zero. The scalar–turbulence interaction term and molecular dissipation term remain the leading order contributors to the $\widetilde {\varepsilon }_\xi $ transport throughout the flame brush in all cases that have been analysed in the present study. Performances of existing models for the unclosed terms of the transport equations of $\widetilde {\varepsilon }_Y $ and $\widetilde {\varepsilon }_\xi $ are assessed with respect to the corresponding quantities obtained from DNS data. Based on this exercise either suitable models have been identified or new models have been proposed for the accurate closure of the unclosed terms of both $\widetilde {\varepsilon }_Y $ and $\widetilde {\varepsilon }_\xi $ transport equations in the context of Reynolds Averaged Navier–Stokes (RANS) simulations.