Effects of Mean Inflow Velocity and Droplet Diameter on the Propagation of Turbulent V-Shaped Flames in Droplet-Laden Mixtures

Three-dimensional carrier phase Direct Numerical Simulations of V-shaped n-heptane spray flames have been performed for different initially mono-sized droplet diameters to investigate the influence of mean flow velocity on the burning rate and flame structure at different axial locations from the flame holder. The fuel is supplied as liquid droplets through the inlet and an overall (i.e., liquid + gaseous) equivalence ratio of unity is retained in the unburned gas. Additionally, turbulent premixed stoichiometric V-shaped n-heptane flames under the same turbulent flow conditions have been simulated to distinguish the differences in combustion behaviour of the pure gaseous phase premixed combustion in comparison to the corresponding behaviour in the presence of liquid n-heptane droplets. It has been found that reacting gaseous mixture burns predominantly under fuel-lean mode and the availability of having fuel-lean mixture increases with increasing mean flow velocity. The extent of flame wrinkling for droplet cases has been found to be greater than the corresponding gaseous premixed flames due to flame-droplet-interaction, which is manifested by dimples on the flame surface, and this trend strengthens with increasing droplet diameter. As the residence time of the droplets within the flame decreases with increasing mean inflow velocity, the droplets can survive for larger axial distances before the completion of their evaporation for the cases with higher mean inflow velocity and this leads to greater extents of flame-droplet interaction and droplet-induced flame wrinkling. Mean inflow velocity, droplet diameter and the axial distance affect the flame brush thickness. The flame brush thickens with increasing droplet diameter for the cases with higher mean inflow velocity due to the predominance of fuel-lean gaseous mixture within the flame. However, an opposite behaviour has been observed for the cases with lower mean inflow velocity where the smaller extent of flame wrinkling due to smaller values of integral length scale to flame thickness ratio arising from higher likelihood of fuel-lean combustion for larger droplets dominates over the thickening of the flame front. It has been found that the major part of the heat release arises due to premixed mode of combustion for all cases but the contribution of non-premixed mode of combustion to the total heat release has been found to increase with increasing mean inflow velocity and droplet diameter. The increase in the mean inflow velocity yields an increase in the mean values of consumption and density-weighted displacement speed for the droplet cases but leads to a decrease in turbulent burning velocity. By contrast, an increase in droplet diameter gives rise to decreases in turbulent burning velocity, and the mean values of consumption and density-weighted displacement speeds. Detailed physical explanations have been provided to explain the observed mean inflow velocity and droplet diameter dependences of the flame propagation behaviour.

Scaling of Second-Order Structure Functions in Turbulent Premixed Flames in the Flamelet Combustion Regime

The second-order velocity structure function statistics have been analysed using a DNS database of statistically planar turbulent premixed flames subjected to unburned gas forcing. The flames considered here represent combustion for moderate values of Karlovitz number from the wrinkled flamelets to the thin reaction zones regimes of turbulent premixed combustion. It has been found that the second-order structure functions exhibit the theoretical asymptotic scalings in the dissipative and (relatively short) inertial ranges. However, the constant of proportionality for the theoretical asymptotic variation for the inertial range changes from one case to another, and this value also changes with structure function orientation. The variation of the structure functions for small length scale separation remains proportional to the square of the separation distance. However, the constant of proportionality for the limiting behaviour according to the separation distance square remains significantly different from the theoretical value obtained in isotropic turbulence. The disagreement increases with increasing turbulence intensity. It has been found that turbulent velocity fluctuations within the flame brush remain anisotropic for all cases considered here and this tendency strengthens towards the trailing edge of the flame brush. It indicates that the turbulence models derived based on the assumptions of homogeneous isotropic turbulence may not be fully valid for turbulent premixed flames.

Closure Relations for Fluxes of Flame Surface Density and Scalar Dissipation Rate in Turbulent Premixed Flames

In this study, closure relations for total and turbulent convection fluxes of flame surface density and scalar dissipation rate were developed (i) by placing the focus of consideration on the flow velocity conditioned to the instantaneous flame within the mean flame brush and (ii) by considering the limiting behavior of this velocity at the leading and trailing edges of the flame brush. The model was tested against direct numerical simulation (DNS) data obtained from three statistically stationary, one-dimensional, planar, premixed turbulent flames associated with the flamelet regime of turbulent burning. While turbulent fluxes of flame surface density and scalar dissipation rate, obtained in the DNSs, showed the countergradient behavior, the model predicted the total fluxes reasonably well without using any tuning parameter. The model predictions were also compared with results computed using an alternative closure relation for the flame-conditioned velocity.

Investigation of Extinction and Reignition Events Using the Flamelet Generated Manifold Model

Simulations for the Cambridge swirl bluff-body spray burner are performed near blow-out conditions. A hybrid stress blended eddy simulation (SBES) model is utilized for sub-grid turbulence closure. SBES blends the RANS-SST model at the boundary layer with large eddy simulation dynamic Smagorinsky model outside the boundary layer. The injected N-heptane spray droplets are tracked using a typical Eulerian-Lagrangian approach. Heat transfer coupling between the bluff-body walls and the near-walls fluid is accounted for by coupling the solid and fluid energy equations at the bluff-body surface. Mixing and chemistry are modeled using the Flamelet Generated Manifold (FGM) model. The study investigates how successful the FGM model is in predicting finite rate effects like local extinction and flame lift-off height. To this end, two near blow-out spray flames, the H1S1 (75% to blow-out) and H1S2 (88% to blow-out) are simulated. Good results are shown matching the spray Sauter mean diameter (SMD) and axial velocity mean and rms experimental data. The results also show that the FGM model captured reasonably well the flame structure and lift-off height as well as the spray pattern. Overall the spray droplets mean D32 and mean axial velocity were under-predicted, while the rms distribution matched reasonably well for the H1S1 flame. The mean flame brush lift-off height is estimated based on the statistically stationary mean flame brush and is estimated to be around 6 mm from the bluff-body base. Instantaneous local flame extinction is observed. The H1S2 flame, however, showed similar but slightly better match with the measurements for the mean spray data compared to the H1S1 flame, with slight under-prediction for D32 at Z = 10 mm and Z = 20 mm. Future work will investigate the sensitivity of the simulation to the spray boundary conditions and grid resolution.

Spatially localized multi-scale energy transfer in turbulent premixed combustion

A three-dimensional wavelet multi-resolution analysis of direct numerical simulations of a turbulent premixed flame is performed in order to investigate the spatially localized spectral transfer of kinetic energy across scales in the vicinity of the flame front. A formulation is developed that addresses the compressible spectral dynamics of the kinetic energy in wavelet space. The wavelet basis enables the examination of local energy spectra, along with inter-scale and subfilter-scale (SFS) cumulative energy fluxes across a scale cutoff, all quantities being available either unconditioned or conditioned on the local instantaneous value of the progress variable across the flame brush. The results include the quantification of mean spectral values and associated spatial variabilities. The energy spectra undergo, in most locations in the flame brush, a precipitous drop that starts at scales of the same order as the characteristic flame scale and continues to smaller scales, even though the corresponding decrease of the mean spectra is much more gradual. The mean convective inter-scale flux indicates that convection increases the energy of small scales, although it does so in a non-conservative manner due to the high aspect ratio of the grid, which limits the maximum scale level that can be used in the wavelet transform, and to the non-periodic boundary conditions, which exchange energy through surface forces, as explicitly elucidated by the formulation. The mean pressure-gradient inter-scale flux extracts energy from intermediate scales of the same order as the characteristic flame scale, and injects energy in the smaller and larger scales. The local SFS-cumulative contribution of the convective and pressure-gradient mechanisms of energy transfer across a given cutoff scale imposed by a wavelet filter is analysed. The local SFS-cumulative energy flux is such that the subfilter scales upstream from the flame always receive energy on average. Conversely, within the flame brush, energy is drained on average from the subfilter scales by convective and pressure-gradient effects most intensely when the filter cutoff is larger than the characteristic flame scale.