Large Eddy Simulations with flamelet-based thermochemistry are used to investigate the behavior of a premixed hydrogen-air flame stabilized by a bluff-body. Validation against experimental data is carried out first to demonstrate the model’s ability to predict both velocity field and flame structure. The capability of the model in predicting differential diffusion effects is then assessed, in particular regarding the coupling between differential diffusion, tangential strain and curvature, and their effect on mixture fraction redistribution and reaction rate variation. Results indicate that unstretched flamelet thermochemistry is capable of capturing the increase in mixture fraction caused by positive resolved strain, as well as negative variations of mixture fraction due to negative curvature. Furthermore, the model is observed to mimic the effects of negative Markstein length to a certain extent, so that positive tangential strain causes reaction rate increase. The interplay between resolved stretch and preferential diffusion is also shown to lead to a shorter flame length which is in better agreement with experimental observations as compared to simulations under unity Lewis number assumption. These findings highlight that the macroscopic effects of differential diffusion and stretch on the premixed hydrogen flame, characterized by significant strain levels, can be predicted using a flamelet-based approach and without recurring to strained flamelets database, which implies important simplifications in the combustion modeling of turbulent hydrogen-premixed flames and offers valuable insights for the design of novel combustors.

Large eddy simulations (LES) with flamelet and presumed filtered density function closure are used to simulate turbulent premixed and partially premixed hydrogen flames. Different approaches to model differential diffusion are investigated and compared. In particular, two existing models are extended to the LES framework to correct the resolved diffusive flux of the controlling variables due to differential diffusion. A lean premixed turbulent hydrogen flame in a slotted burner configuration is simulated first to compare the capability of the considered models in capturing local mixture fraction redistribution, super-adiabatic temperatures and thermo-diffusive instabilities. Results show that both models describe the formation of cellular burning structures. Next, a partially premixed lifted hydrogen flame in vitiated hot coflow is simulated to gain insight on the relevance of differential diffusion modelling at a higher turbulence level, a different combustion mode and in the presence of a complex stabilisation mechanism. Good predictions of the turbulent mixing and temperature fields are observed. Moreover, results show that the flame lift-off height has an appreciable sensitivity to the differential diffusion model. When differential diffusion is included only in the thermochemistry database, only mild effects on the predicted temperature fields, mixing and flame height are observed. On the contrary, a considerable shift of the flame base is observed when corrections are applied in the LES at the resolved level, depending on what controlling variables are considered. Further analyses reveal how the corrections of diffusive fluxes in the thermochemistry and at the LES level affect differently the flame burning mode, whose details are given throughout the paper.

Large eddy simulation (LES) paradigms are used in the present work to predict premixed and partially premixed turbulent flames with flamelets based thermochemistry and presumed filtered density function approach for turbulence-chemistry interaction modelling. The combustion model requires a closure for the scalar dissipation rate of a progress variable, in which a modelling constant must be chosen. The present work focuses on the computation of the model constant through dynamic procedures based on the scale-similarity assumption, which requires the application of test-filters. In particular, two test-filtering approaches for LES, based respectively on an algebraic formulation and a newly proposed differential equation, are tested for flame configurations at different levels of turbulence, and using block-structured and unstructured meshes. The analysis shows that the differential filter, unlike the algebraic one, is handled well in situations of weak turbulence at comparable computational costs. At higher turbulence conditions the outcome looks less dependent on the test-filter and mesh topology used, although quantitative differences in the behaviour of the dynamically-computed model constant are still observable and discussed. Further analyses to understand the behaviour of the two filters are presented in the paper.

Large eddy simulation (LES) paradigms are employed to analyse the internal flow field of a lean premixed swirl-stabilized combustor with axial air injection at both non-reacting and reacting conditions, for a methane and a methane-hydrogen fuel mixture. The thickened flame combustion model (TFM) with detailed chemical kinetic mechanism is employed to simulate the flow. An adaptive mesh strategy is used to maximise the mesh resolution in the flame and boundary layer regions. The numerical results for the methane flame are firstly validated against experimental velocity measurements obtained via particle image velocimetry (PIV). Subsequently the LES is employed to simulate hydrogen-enriched methane flames by keeping the same output power in the combustor, in order to obtain insights on the flow behaviour when hydrogen is added, in terms of flame stability and emissions. A POD analysis reveals the presence of a precessing vortex core (PVC) in both reacting and non-reacting conditions, and how this PVC is affected by the reactants mixture is discussed in the paper. Moreover, the flame is observed to propagate upstream in the jet core despite the use of axial air injection, although flashback is not observed. In terms of emissions, significant reduction in CO and NO_x is observed when adding the hydrogen to the reactants mixture despite the higher flame speed, the reason for are discussed in the paper.