Flame Surface Density based mean reaction rate closure for Reynolds averaged Navier Stokes methodology in turbulent premixed Bunsen flames with non-unity Lewis number

The Flame Surface Density (FSD) based mean reaction rate closure proposed by Prof. K.N.C. Bray and co-workers in the context of Reynolds Averaged Navier-Stokes simulations has been a-priori analysed using a Direct Numerical Simulation (DNS) database of turbulent premixed Bunsen flames with different characteristic Lewis numbers representing the strict flamelet regime (i.e. high Damkohler number and low Karlovitz number combustion). The statistical behaviours of stretch factor, orientation factor and wrinkling length scale have been assessed for non-unity Lewis number conditions to identify their Lewis number dependencies. The assumption of presumed bimodal distribution has been found to be rendered invalid close to the nozzle exit where the unresolved wrinkling remains relatively small even when the flow parameters at the nozzle exit represent high Damkohler number and low Karlovitz number conditions. Although the PDF of reaction progress variable shows some resemblance to a bimodal distribution away from the nozzle exit, the Bray-Moss-Libby expressions which can be derived for infinitely large values of Damkohler number have been found to show considerable deviations from the Reynolds averaged reaction progress variable and reaction progress variable variance extracted from DNS data even though these cases represent the wrinkled/corrugated flamelets regime combustion based on nozzle exit conditions. This suggests that it might be necessary to solve a modelled scalar variance transport equation even in the strict flamelet regime instead of using the algebraic relation for the scalar variance. It has been found that the characteristic Lewis number has a major influence on the orientation factor, wrinkling length scale, and the stretch factor. Furthermore, these parameters are found to be (sometimes strong) functions of the axial distance from the nozzle exit. Known parameterisations for the wrinkling length scale and the stretch factor have been shown to be unable to capture the correct variation across the flame brush based on a-priori analysis of DNS data. The variability of the orientation factor and the inadequacy of existing relations to approximate other quantities such as the stretch factor, wrinkling factor and the Reynolds averaged reaction progress variable have the potential to severely limit the performance of algebraic FSD based mean reaction rate closures in turbulent premixed flames within the flamelet regime. (C) 2021 The Combustion Institute. Published by Elsevier Inc. All rights reserved.