This study investigates the effect of a rotating axial swirler on the flame shape and combustion dynamics of premixed methane-air flames in a dual-swirler burner by using simultaneous hydroxyl (OH) planar laser- induced fluorescence (PLIF) and stereoscopic particle image velocimetry (sPIV). The rapid expansion from the nozzle to the combustor induces shear layer vortices that fold the flame, resulting in subsequent flame-flame interaction. This flame-flame interaction results in flame surface annihilation, causing heat release fluctuations that generate coherent pressure fluctuations, which is explored through a phase-averaged analysis. Increasing the rotational speed of the axial swirler diminishes vortex strength at the shear layer and shifts the vortex trajectory away from the central axis. This is caused by an increase in radial velocity and a decrease in the difference in momentum flux between the inside and outside shear layer. As the momentum flux difference decreases, the shear force weakens, leading to a reduction in the strength of the shear layer vortex. The shear layer vortex structure propagates outward and eventually dissipates due to the increased radial velocity with an increase in swirl number. This, in turn, substantially suppresses flame-flame interactions, resulting in the of the fluctuations.

Local flame displacement speed Sd of a turbulent premixed flame is of fundamental and practical interest. For H2-air flames, the interest is further accentuated given the recent drive towards the development of zero-carbon combustors for both power and aircraft engine applications. The present study investigates several three-dimensional Direct Numerical Simulation (3D DNS) cases of premixed H2-air turbulent flames to theoretically model the Sd at negative curvatures, building upon recent works. Two of the four DNS cases presented are simulated at atmospheric pressure and two at elevated pressure. The DNS cases at different turbulence Reynolds numbers (Ret) and Karlovitz numbers (Ka) are generated using detailed chemistry. It has been shown in the previous studies that at atmospheric pressure, the density-weighted flame displacement speed Sd˜ is enhanced significantly over its laminar value (SL) at large negative curvature κ due to flame-flame interactions. The current work justifiably employs an imploding cylindrical laminar flame configuration to represent the local flame surfaces undergoing flame-flame interaction in a 3D turbulent flame. Therefore, to acquire a deep understanding of the interacting flame dynamics at large negative curvatures, one-dimensional (1D) simulations of an inwardly propagating cylindrical H2-air laminar premixed flame, with detailed chemistry at the corresponding atmospheric and elevated pressure conditions are performed. In particular, the 1D simulations emphasized the transient nature of the flame structure during these interactions. Based on the insights from the 1D simulations, we utilize an analytical approach to model the Sd˜ at these regions of extreme negative κ of the 3D DNS. The analytical approach is formulated to include the effect of variable density, convection and the inner reaction zone motion. The joint probability density function (JPDF) of Sd˜ and κ and the corresponding conditional averages obtained from 3D DNS showed clear negative correlation between Sd˜ and κ at all pressures. The obtained model successfully predicts the variation of 〈Sd˜|κ〉 with κ for the regions on the flame surface with large negative curvature (κδL≪−1) at atmospheric as well as at elevated pressure, with good accuracy. This showed that the 1D cylindrical, interacting flame model is a fruitful representation of a local flame-flame interaction that persists in a 3D turbulent flame, and is able to capture the intrinsically transient dynamics of the local flame-flame interaction. The 3D DNS cases further showed that even in the non-interacting state at κ=0, on average Sd˜ can deviate from SL. Sd˜ at κ=0 is a manifestation of the internal flame structure, controlled by turbulence transport in the large Ka regime. Therefore, the correlation of 〈Sd˜〉/SL with the the normalized gradient of the progress variable, 〈|∇c^|c0〉 at κ=0 is explored.