Large eddy simulation is used to investigate the flashback mechanism caused by the combustion-induced vortex breakdown (CIVB) in a high-pressure lean-burn annular combustor with lean direct injection of kerosene. A single sector of the geometry, including a central pilot flame surrounded by a main flame, is simulated at take-off conditions. A previously-developed flamelet-based approach is used to model turbulence-combustion interactions due to its relatively low cost, allowing to simulate a sufficiently long time window. In stable operations, the flame stabilises in an M-shape configuration and a periodic movement of the pilot jet, with the corresponding formation of a small recirculation bubble, is observed. Flashback is then observed, with the flame accelerating upstream towards the injector as already described in other studies. This LES, however, reveals a precursor partial blow-out of the main flame induced by a cluster of vortices appearing in the outer recirculation region. The combined effect of vortices and sudden quenching alters the mixing level close to the injector, causing first the main, then the pilot flame, to accelerate upstream and initiate the CIVB cycle before the quenched region can re-ignite. Main and pilot flames partly extinguish as they cross their respective fuel injection point, and re-ignition follows due to the remnants of the reaction in the pilot stream. The process is investigated in detail, discussing the causes of CIVB-driven flashback in realistic lean-burn systems.

In gas turbines, combustor inlets are characterised by significant levels of unsteady circumferential distortion due to compressor wakes and secondary flows, together with additional radial non-uniformity induced by the adverse pressure gradients in the pre-diffuser. This can cause non-uniform velocity distributions across the fuel injector, although the exact interaction mechanism, and the effects it has on the downstream air-fuel mixing, is not fully understood. This paper investigates the flow in an a single sector of a fully featured isothermal rig comprising of compression and combustion systems, exploiting the synchronous coupling of a compressible unsteady RANS simulation with a low-Mach LES. Validation against five-hole probe measurements shows that the coupled approach can correctly predict distortion onset and development, with no solution discontinuity at the coupling interface, and is able to preserve unsteady information. The coupled prediction is then compared against a standalone combustor simulation carried out using a circumferentially uniform inlet profile, showing that the additional turbulence from the wakes interacts with the injector, reducing the coherence of the precessing vortex core and potentially affecting the air-fuel mixing characteristics.

The physical mechanism leading to flame local extinction remains a key issue to be further understood. An analysis of large eddy simulation (LES) data with presumed probability density function (PDF) based closure (Chen et al., 2020, Combust. Flame, vol. 212, pp. 415) indicated the presence of localised breaks of the flame front along the stoichiometric line. These observations and their relation to local quenching of burning fluid particles, together with the possible physical mechanisms and conditions allowing their appearance in LES with a simple flamelet model, are investigated in this work using a combined Lagrangian-Eulerian analysis. The Sidney/Sandia piloted jet flames with compositionally inhomogeneous inlet and increasing bulk speeds, amounting to respectively 70 and 90% of the experimental blow-off velocity, are used for this analysis. Passive flow tracers are first seeded in the inlet streams and tracked for their lifetime. The critical scenario observed in the Lagrangian analysis, i.e., burning particles crossing extinction holes on the stoichiometric iso-surface, is then investigated using the Eulerian control-volume approach. For the 70% blow-off case the observed flame front breaks/extinction holes are due to cold and inhomogeneous reactants that are cast onto the stoichiometric iso-surface by large vortices initiated in the jet/pilot shear layer. In this case an extinction hole forms only when the strain effect is accompanied by strong subgrid mixing. This mechanism is captured by the unstrained flamelets model due to the ability of the LES to resolve large-scale strain and considers the SGS mixture fraction variance weakening effect on the reaction rate through the flamelet manifold. Only at 90% blow-off speed the expected limitation of the underlying combustion model assumption become apparent, where the amount of local extinctions predicted by the LES is underestimated compared to the experiment. In this case flame front breaks are still observed in the LES and are caused by a stronger vortex/strain interaction yet without the aid of mixture fraction variance. The reasons for these different behaviours and their implications from a physical and modelling point of view are discussed in this study.