Hazem S. A. M. Awad

A priori Direct Numerical Simulation (DNS) assessment of mean reaction rate closures for reaction progress variable in the context of Reynolds Averaged Navier–Stokes (RANS) simulations has been conducted for MILD combustion of homogeneous (i.e., constant equivalence ratio), methane-air mixtures. The reaction rate predictions according to statistical (e.g., presumed probability density function), phenomenological (e.g., eddy-break up (EBU), eddy dissipation concept (EDC) and the scalar dissipation rate (SDR) based approaches), and flame surface description (e.g., Flame Surface Density) based closures are compared. The performance of the various reaction rate closures has been assessed by comparing the models’ predictions to the corresponding quantities extracted from DNS data. It has been found that the usual presumed probability density function (PDF) approach using the beta-function predicts the PDF of the reaction progress variable in homogenous mixture MILD combustion throughout the flame brush for all cases considered here provided that the scalar variance is accurately predicted. The accurate estimation of scalar variance requires the solution of a modelled transport equation, which depends on the closure of Favre-averaged SDR. A linear relaxation based algebraic closure for the Favre-averaged SDR has been found to capture the behaviour of the Favre-averaged SDR in the current homogenous mixture MILD combustion setup. It has been found that the EBU, SDR and FSD-based mean reaction rate closures do not adequately predict the mean reaction rate of the reaction progress variable for the parameter range considered here. However, a variant of the EDC closure, with model coefficients expressed as functions of micro-scale Damköhler and turbulent Reynolds numbers, has been found to be more successful in predicting the mean reaction rate of reaction progress variable compared to other modelling methodologies for the range of turbulence intensities and dilution levels considered here.

Three-dimensional Direct Numerical Simulations (DNS) data has been utilised to analyse statistical behaviours of the scalar dissipation rate (SDR) and its transport for homogeneous methane-air mixture turbulent Moderate or Intense Low oxygen Dilution (MILD) combustion for different O₂ dilution levels and turbulence intensities for different reaction progress variable definitions. Additional DNS has been conducted for turbulent premixed flames and passive scalar mixing for the purpose of comparison with the SDR statistics of the homogeneous mixture MILD combustion with that in conventional premixed combustion and passive scalar mixing. It has been found that the peak mean value of the scalar dissipation rate decreases with decreasing O₂ concentration for MILD combustion cases. Moreover, SDR magnitudes increase with increasing turbulence intensity for both MILD and conventional premixed combustion cases. The profiles and mean values of the scalar dissipation rate conditioned upon the reaction progress variable are found to be sensitive to the choice of the reaction progress variable definition. This behaviour arises due to the differences in the distributions of the species mass fractions within the flame. The strain rate contribution and the molecular dissipation term are found to be the leading order contributors in the scalar dissipation rate transport for MILD combustion; whereas, in conventional premixed flames, the terms rising from density variation and reaction rate gradient also play leading roles in addition to the strain rate and molecular dissipation contributions. By contrast, the terms due to density gradient and reaction rate gradient remain negligible in comparison to the leading order contributors in MILD combustion cases due to small density variation because of moderate temperature rise and small reaction rate gradient magnitudes. Furthermore, the qualitative behaviour of the strain rate contribution to the SDR transport in premixed flames is significantly different to that in the case of MILD combustion and passive scalar mixing. The findings of the current analysis indicate that the scalar dissipation rate statistics in MILD combustion show several qualitative similarities to the passive scalar mixing despite major differences with the SDR transport in conventional turbulent premixed flames. This further suggests that the scalar dissipation rate models, which were originally proposed in the context of passive scalar mixing, have the potential to be applicable for MILD combustion but the models for the premixed turbulent combustion may not be applicable for MILD combustion of homogeneous mixtures.

Moderate or intense low-oxygen dilution (MILD) combustion is a novel combustion technique that can simultaneously improve thermal efficiency and reduce emissions. This paper focuses on the differences in statistical behaviours of the surface density function (SDF = magnitude of the reaction progress variable gradient) between conventional premixed flames and exhaust gas recirculation (EGR) type homogeneous-mixture combustion under MILD conditions using direct numerical simulations (DNS) data. The mean values of the SDF in the MILD combustion cases were found to be significantly smaller than those in the corresponding premixed flame cases. Moreover, the mean behaviour of the SDF in response to the variations of turbulence intensity were compared between MILD and premixed flame cases, and the differences are explained in terms of the strain rates induced by fluid motion and the ones arising from flame displacement speed. It was found that the effects of dilatation rate were much weaker in the MILD combustion cases than in the premixed flame cases, and the reactive scalar gradient in MILD combustion cases preferentially aligns with the most compressive principal strain-rate eigendirection. By contrast, the reactive scalar gradient preferentially aligned with the most extensive principal strain-rate eigendirection within the flame in the premixed flame cases considered here, but the extent of this alignment weakened with increasing turbulence intensity. This gave rise to a predominantly positive mean value of normal strain rate in the premixed flames, whereas the mean normal strain rate remained negative, and its magnitude increased with increasing turbulence intensity in the MILD combustion cases. The mean value of the reaction component of displacement speed assumed non-negligible values in the MILD combustion cases for a broader range of reaction progress variable, compared with the conventional premixed flames. Moreover, the mean displacement speed increased from the unburned gas side to the burned gas side in the conventional premixed flames, whereas the mean displacement speed in MILD combustion cases decreased from the unburned gas side to the middle of the flame before increasing mildly towards the burned gas side. These differences in the mean displacement speed gave rise to significant differences in the mean behaviour of the normal strain rate induced by the flame propagation and effective strain rate, which explains the differences in the SDF evolution and its response to the variation of turbulence intensity between the conventional premixed flames and MILD combustion cases. The tangential fluid-dynamic strain rate assumed positive mean values, but it was overcome by negative mean values of curvature stretch rate to yield negative mean values of stretch rate for both the premixed flames and MILD combustion cases. This behaviour is explained in terms of the curvature dependence of displacement speed. These findings suggest that the curvature dependence of displacement speed and the scalar gradient alignment with local principal strain rate eigendirections need to be addressed for modelling EGR-type homogeneous-mixture MILD combustion.