Choice of reaction progress variable under preferential diffusion effects in turbulent syngas combustion based on detailed chemistry direct numerical simulations

The combustion of hydrogen and carbon-monoxide mixtures, so-called syngas, plays an increasingly important role in the safety context of non-fossil energy generation, more specifically in the risk management of incidents in process engineering plants for ammonia synthesis and in nuclear power plants. In order to characterize and simulate syngas/air combustion on industrially relevant scales, subgrid modelling is required, which is often based on a reaction progress variable. To understand the influence of different fuel compositions, turbulence intensities and flame topologies on different possible definitions of reaction progress variable, detailed chemistry direct numerical simulations data of premixed, lean hydrogen/air and syngas/air flames has been considered. A reaction progress variable based on normalized molecular oxygen mass fraction has been found not to capture the augmentation of the normalized burning rate per unit flame surface area in comparison to the corresponding 1D unstretched premixed flame due to preferential diffusion effects. By contrast, reaction progress variables based on other individual species, such as hydrogen, can capture the augmentation of the rate of burning well, but exhibit a pronounced sensitivity to preferential diffusion effects, especially in response to flame curvatures. However, a reaction progress variable based on the linear combination of the main products can accurately represent the temperature evolution of the flame for different mixtures, turbulence intensities and varying local flame topology, while effectively capturing the augmentation of burning rate due to preferential diffusion effects. However, its tendency to assume values larger than 1.0 in the regions of super-adiabatic temperatures poses challenges for future modeling approaches, whereas the reaction progress variable based on hydrogen mass fraction remains bound between 0.0 and 1.0 despite showing deviations in comparison to corresponding variations obtained from the unstretched laminar flame depending on flame curvature variations.