Flame self-interactions in MILD combustion of homogeneous and inhomogeneous mixtures

Flame self-interaction (FSI) events in Moderate or Intense Low-Oxygen Dilution (MILD) combustion of homogeneous and inhomogeneous mixtures of methane and oxidiser have been analysed using three-dimensional Direct Numerical Simulations (DNS). The simulations have been conducted at the same global equivalence ratio ( ⟨ ϕ ⟩ = 0.8 ) for different levels of O 2 concentration (dilution) and initial turbulence intensities. It has been reported that both homogeneous and inhomogeneous mixture MILD combustion cases exhibit significant occurrences of FSI events, with the peak frequency of FSI events occurring towards the burned gas side in all cases. Moreover, the frequency of FSI events increases with increasing dilution level and turbulence intensity, but the presence of mixture inhomogeneity leads to a reduction in total FSI events. In all cases, the cylindrical FSI topologies (i.e. tunnel formation and tunnel closure) were found to have a higher likelihood of occurrence compared to spherical FSI topologies (i.e. unburned and burned gas pockets). The geometries of FSI topologies were also analysed using the mean and Gaussian curvatures. It has been shown that the inward propagating spherical FSI topologies (i.e. unburned gas pockets) are associated with negative mean curvature, while outward propagating spherical FSI topologies (i.e. burned gas pockets) are associated with positive mean curvature. Moreover, tunnel formation (tunnel closure) FSI topologies predominantly exhibit either elliptic geometries with positive (negative) mean curvature or hyperbolic saddle geometries with negative (positive) mean curvature. It has been shown for the first time in MILD combustion that the mean values of kinematic restoration and dissipation terms in the transport equation of the magnitude of the reaction progress variable conditional upon the reaction progress variable tend to cancel each other in the vicinity of the critical points associated with cylindrical topologies. Thus, the singular contributions in these terms, which are obtained from analytical descriptions in the vicinity of tunnel formation and tunnel closure topologies, do not affect the balance equation of the magnitude of the gradient of the reaction progress variable. Consequently, there is no need for a separate model treatment for singularities in modelling approaches based on the magnitude of the gradient of the reaction progress variable. The FSI events in the reaction dominated and propagating flame regi