Three-dimensional direct numerical simulations of V-flames interacting with chemically inert walls in a fully developed turbulent channel flow have been performed under adiabatic and isothermal wall boundary conditions using single-step chemistry. These simulations are representative of stoichiometric methane-air mixture at unity Lewis number under atmospheric conditions. The turbulence in the non-reacting channel is representative of the friction velocity based Reynolds number $$Re_{\tau }=110$$ . Differences in the statistical behaviours of the mean values of progress variable, temperature, and density have been reported for different wall boundary conditions. It is found that the mean location of the oblique flame interacting with the wall is affected by the choice of the wall boundary condition used. The influence of these differences on the flame dynamics is investigated by analysing the statistical behaviours of the surface density function (SDF) and the strain rates, which govern the evolution of the SDF. The mean variation of the SDF and the flame displacement speed are strongly affected by the wall boundary condition within the viscous sub-layer region of the boundary layer. The behaviours of the normal and tangential strain rates are found to be influenced by not only the differences in the wall boundary conditions, but also by the distance from the wall. The differences in the displacement speed statistics for different wall boundary conditions and wall distance affect the behaviours of the normal strain rate arising due to flame propagation and curvature stretch. The changes in the SDF behaviour in the near wall region have been explained in terms of the statistics of effective normal strain rate experienced by the progress variable iso-surfaces under different wall boundary conditions and wall normal distances.
Multiscale analysis of wall-bounded turbulent premixed flames is performed using three-dimensional direct numerical simulation data of flame-wall interaction (FWI). The chosen configuration represents head-on quenching of a turbulent statistically planar stoichiometric methane-air flame by an isothermal inert wall. Different turbulence intensities and chemical mechanisms have been analyzed. A bandpass filtering technique is utilised to analyze the influence of turbulent eddies of varying size and the statistics of vorticity and strain rate fields associated with them. It is found that the presence of the flame does not alter the mechanism of vortex stretching in turbulent flows when the flame is away from the wall, but in the case of FWI, the mechanism of vortex stretching is altered due to a reduction in the contribution from non-local strain, and the small scales of turbulence start to contribute to the flame straining process. The results indicate that small scale eddies do not contribute to the tangential strain rate when the flames are away from the walls, whereas the contribution from the small scales to the tangential strain rate increases when the flame is in the vicinity of the wall. It is also found that the choice of chemical mechanism does not influence the underlying fluid mechanical processes involved in FWI.
17 March 2019 in Combustion Science and Technology
The bending effect of turbulent flame speed variation (i.e., the deviation from the linear increase of flame speed with increasing root-mean-square turbulent velocity fluctuation) has been investigated based on a Direct Numerical Simulation database of statistically planar turbulent premixed flames propagating into forced unburned gas turbulence. The validity of Damköhler’s first hypothesis has been utilized to analyze the bending effect in terms of generalized Flame Surface Density (FSD) evolution. The volume-integrated value of the tangential strain rate term of the FSD transport equation remains positive, whereas the volume-integrated value of the curvature term assumes negative values. Under statistically stationary state, the positive value of the volume-integrated tangential strain rate term remains in equilibrium with the negative value of the volume-integrated curvature term. It has been found that the contribution of the normal strain rate to the flame surface area remains negative for small turbulence intensities, which eventually become positive for large turbulence intensities. This is a consequence of the change of collinear alignment of the reaction progress variable gradient from the most extensive principal strain rate direction to the direction of the eigenvector associated with the most compressive principal strain rate with increasing turbulence intensity. An increase in turbulence intensity increases the width of the probability density functions of flame curvature, and thereby increases the surface-averaged curvature squared values. This eventually makes the FSD curvature term due to the tangential diffusion component of displacement speed as the major contributor to the negative contribution of the volume-integrated curvature term in the FSD transport equation for large turbulence intensities. However, the negative contribution of the volume-integrated FSD curvature term does not increase indefinitely with increasing turbulence intensity and the inner cut-off scale, which also limits the maximum possible value of the volume-integrated FSD strain rate term under statistically stationary state, governs the maximum possible destruction of flame surface area. It has been argued that the upper limits of the flame surface area generation and destruction are responsible for the bending effects in the variations of turbulent flame speed and flame surface area.
Turbulence in open channel flows is ubiquitous to hydro-environmental applications and has recently increased in importance with the deployment of tidal stream turbines, as turbulence impacts both the turbine performance and the blade fatigue life. Tidal turbine analysis requires fully developed turbulence characteristics at the inlet of numerical simulations where generally the length scale information is limited. In this study, fully resolved large eddy simulations (LES) with flat beds were undertaken using an open source code at friction Reynolds numbers ( ) of 150, 400 and 1020. It was found that the effects of the free surface on turbulence length scales were felt in approximately the uppermost 10% of the channel only, although the influence on Reynolds stresses extended further downwards. Furthermore, the cross-correlation length scales of both streamwise and spanwise velocities were found to be significantly affected by the free surface where turbulent eddies were flattened to the two-component limit.
A direct numerical simulation (DNS) database for boundary layer flashback of a premixed hydrogen-air flame with an equivalence ratio of 1.5 in a fully developed turbulent channel flow has been considered for this analysis. The nonreacting part of the channel flow is representative of the friction velocity based Reynolds number . A skeletal chemical mechanism with 9 chemical species and 20 reactions is employed for representing hydrogen-air combustion. In this work the flow configuration and the turbulence and flame characteristics are similar to those of Gruber et al. [ J. Fluid Mech. 709, 516 (2012)]. The interaction between the flame structure and the turbulent flow has been investigated for boundary layer flashback for a comparison with the earlier work of Gruber et al. [ J. Fluid Mech. 709, 516 (2012)]. The statistics of wall shear stress, turbulent kinetic energy, and its dissipation have been analyzed to probe the influence of the flame on the underlying turbulence in the channel flow configuration. Furthermore, the budgets for the individual terms in the turbulent kinetic energy transport equation have also been investigated at a given plane in the channel. It is found that the propagation of the flame into the upstream part of the fully developed turbulent boundary layer introduces a flow reversal in some regions upstream of the flame and these regions lead to negative wall shear stress. Interrogation of the DNS data for the budgets of the turbulent kinetic energy transport has revealed that the aforementioned local flow reversal regions have significant influences on the turbulent kinetic energy production, pressure dilatation, and pressure transport terms. It has been found that the flame propagation into the upstream reactants leads to some weak local compressibility effects as demonstrated by the changes in the pressure related terms in the turbulent kinetic energy transport equation. These results indicate that the pressure dilatation and turbulent transport due to pressure are the two dominant terms in the turbulent kinetic energy equation in the case of wall bounded flashback flames.