Scalar Forcing Methodology for Direct Numerical Simulations of Turbulent Stratified Mixture Combustion

Scalar forcing in the context of turbulent stratified flame simulations aims to maintain the fuel-air inhomogeneity in the unburned gas. With scalar forcing, stratified flame simulations have the potential to reach a statistically stationary state with a prescribed mixture fraction distribution and root-mean-square value in the unburned gas, irrespective of the turbulence intensity. The applicability of scalar forcing for Direct Numerical Simulations of stratified mixture combustion is assessed by considering a recently developed scalar forcing scheme, known as the reaction analogy method, applied to both passive scalar mixing and the imperfectly mixed unburned reactants of statistically planar stratified flames under low Mach number conditions. The newly developed method enables statistically symmetric scalar distributions between bell-shaped and bimodal to be maintained without any significant departure from the specified bounds of the scalar. Moreover, the performance of the newly proposed scalar forcing methodology has been assessed for a range of different velocity forcing schemes (Lundgren forcing and modified bandwidth forcing) and also without any velocity forcing. It has been found that the scalar forcing scheme has no adverse impact on flame-turbulence interaction and it only maintains the prescribed root-mean-square value of the scalar fluctuation, and its distribution. The scalar integral length scale evolution is shown to be unaffected by the scalar forcing scheme studied in this paper. Thus, the scalar forcing scheme has a high potential to provide a valuable computational tool to enable analysis of the effects of unburned mixture stratification on turbulent flame dynamics.

Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation

The effects of varying turbulence intensity and turbulence length scale on premixed turbulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation.

Analysis of Flame Front Breaks Appearing in LES of Inhomogeneous Jet Flames Using Flamelets

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.

Flow Separation Dynamics in Three-Dimensional Asymmetric Diffusers

The mean and instantaneous flow separation of two different three-dimensional asymmetric diffusers is analysed using the data of large-eddy simulations. The geometry of both diffusers under investigation is based on the experimental configuration of Cherry et al. (Int J Heat Fluid Flow 29(3):803–811, 2008). The two diffusers feature similar area ratios of $$\mathrm{AR}=4.8$$ AR = 4.8 and $$\mathrm{AR}=4.5$$ AR = 4.5 while exhibiting differing asymmetric expansion ratios of $$\mathrm{AER}=4.5$$ AER = 4.5 or $$\mathrm{AER}=2.0$$ AER = 2.0 , respectively. The Reynolds number based on the averaged inlet velocity and height of the inlet duct is approximately $${\textit{Re}}=10{,}000$$ Re = 10 , 000 . The time-averaged flow in both diffusers in terms of streamwise velocity profiles or the size and location of the mean backflow region are validated using experimental data. In general good agreement of simulated results with the experimental data is found. Further quantification of the flow separation behaviour and unsteadiness using the backflow coefficient reveals the volume portion in which the instantaneous reversal flow evolves. This new approach investigates the cumulative fractional volume occupied by the instantaneous backflow throughout the simulation, a power density spectra analysis of their time series reveals the periodicity of the growth and reduction phases of the flow separation within the diffusers. The dominating turbulent events responsible for the formation of the energy-containing motions including ejection and sweep are examined using the quadrant analysis at various locations. Finally, isourfaces of the Q-criterion visualise the instantaneous flow and the origin and fate of coherent structures in both diffusers.