Effects of turbulent reynolds number on the performance of algebraic flame surface density models for large eddy simulation in the thin reaction zones regime: A direct numerical simulation analysis

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Ret has been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Ret and Ka before reaching an asymptotic value for large values of Ret and Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel'dovich flame thickness ηi/δz does not exhibit any significant dependence on Ret for the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width. © 2012 Mohit Katragadda et al.

DNS of EGR-type turbulent flame in MILD condition

Three-dimensional direct numerical simulation (DNS) of exhaust gas recirculation (EGR)-type turbulent combustion operated in moderate and intense low-oxygen dilution (MILD) condition has been carried out to study the flame structure and flame interaction. In order to achieve adequate EGR-type initial/inlet mixture fields, partially premixed mixture fields which are correlated with the turbulence are carefully preprocessed. The chemical kinetics is modelled using a skeletal mechanism for methane-air combustion. The results suggest that the flame fronts have thin flame structure and the direct link between the mean reaction rate and scalar dissipation rate remains valid in the EGR-type combustion with MILD condition. However, the commonly used canonical flamelet is not fully representative for MILD combustion. During the flame-flame interactions, the heat release rate increases higher than the maximum laminar flame value, while the gradient of progress variable becomes smaller than laminar value. It is also proposed that the reaction rate and the scalar gradient can be used as a marker for the flame interaction. © 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Self-excited circumferential instabilities in a model annular gas turbine combustor: Global flame dynamics

In this paper the global flame dynamics of a model annular gas turbine combustor undergoing strong self-excited circumferential instabilities is presented. The combustor consisted of either 12, 15 or 18 turbulent premixed bluff-body flames arranged around an annulus of fixed circumference so that the effect of flame separation distance, S, on the global heat release dynamics could be investigated. Reducing S was found to produce both an increase in the resonant frequency and the limit-cycle amplitudes of pressure and heat release for the same equivalence ratio. The phase-averaged global heat release, obtained from high-speed OH- chemiluminescence imaging from above, showed that these changes are caused by large-scale modifications to the flame structure around the annulus. For the largest S studied (12 flame configuration) the azimuthal instability produced a helical-like global heat release structure for each flame. When S was decreased, large-scale merging or linking between adjacent flames occurred spanning approximately half of the annulus with the peak heat release concentrated at the outer annular wall. The circumferential nature of the instability was evident from both the pressure measurements and the phase-averaged OH- chemiluminescence showing the phase of the heat release on either side of the annulus to be ≈180°apart and spinning in the counter clockwise direction. Both spinning and standing modes were found but only spinning modes are considered in this paper. To the best of the authors knowledge, these are the first experiments to provide a phase-averaged picture of self-excited azimuthal instabilities in a laboratory-scale annular combustor relevant to gas turbines. © 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Flame interactions in turbulent premixed twin V-flames

Multiple flame-flame interactions in premixed combustion are investigated using direct numerical simulations of twin turbulent V-flames for a range of turbulence intensities and length scales. Interactions are identified using a novel automatic feature extraction (AFE) technique, based on data registration using the dual-tree complex wavelet transform. Information on the time, position, and type of interactions, and their influence on the flame area is extracted using AFE. Characteristic length and time scales for the interactions are identified. The effect of interactions on the flame brush is quantified through a global stretch rate, defined as the sum of flamelet stretch and interaction stretch contributions. The effects of each interaction type are discussed. It is found that the magnitude of the fluctuations in flamelet and interaction stretch are comparable, and a qualitative sensitivity to turbulence length scale is found for one interaction type. Implications for modeling are discussed. © 2013 Copyright Taylor and Francis Group, LLC.

Spray flame study using a model gas turbine swirl burner

A model gas turbine burner was employed to investigate spray flames established under globally lean, continuous, swirling conditions. Two types of fuel were used to generate liquid spray flames: palm biodiesel and Jet-A1. The main swirling air flow was preheated to 350°C prior to mixing with airblast-atomized fuel droplets at atmospheric pressure. The global flame structure of flame and flow field were investigated at the fixed power output of 6 kW. Flame chemiluminescence imaging technique was employed to investigate the flame reaction zones, while particle imaging velocimetry (PIV) was utilized to measure the flow field within the combustor. The flow fields of both flames are almost identical despite some differences in the flame reaction zones. © (2013) Trans Tech Publications, Switzerland.

When will a flame describing function approach to thermoacoustics work well?

In any thermoacoustic analysis, it is important not only to predict linear frequencies and growth rates, but also the amplitude and frequencies of any limit cycles. The Flame Describing Function (FDF) approach is a quasi-linear analysis which allows the prediction of both the linear and nonlinear behaviour of a thermoacoustic system. This means that one can predict linear growth rates and frequencies, and also the amplitudes and frequencies of any limit cycles. The FDF achieves this by assuming that the acoustics are linear and that the flame, which is the only nonlinear element in the thermoacoustic system, can be adequately described by considering only its response at the frequency at which it is forced. Therefore any harmonics generated by the flame's nonlinear response are not considered. This implies that these nonlinear harmonics are small or that they are sufficiently filtered out by the linear dynamics of the system (the low-pass filter assumption). In this paper, a flame model with a simple saturation nonlinearity is coupled to simple duct acoustics, and the success of the FDF in predicting limit cycles is studied over a range of flame positions and acoustic damping parameters. Although these two parameters affect only the linear acoustics and not the nonlinear flame dynamics, they determine the validity of the low-pass filter assumption made in applying the flame describing function approach. Their importance is highlighted by studying the level of success of an FDF-based analysis as they are varied. This is achieved by comparing the FDF's prediction of limit-cycle amplitudes to the amplitudes seen in time domain simulations.

Numerical simulation of post-flame oxidation of hydrocarbons in spark ignition engines

About 50-90 percent of the hydrocarbons that escape combustion during flame passage in spark-ignition engine operation are oxidized in the cylinder before leaving the system. The process involves the transport of unreacted fuel from cold walls towards the hotter burned gas regions and subsequent reaction. In order to understand controlling factors in the process, a transient one-dimensional reactive-diffusive model has been formulated for simulating the oxidation processes taking place in the reactive layer between hot burned gases and cold unreacted air/fuel mixture, with initial and boundary conditions provided by the emergence of hydrocarbons from the piston top land crevice. Energy and species conservation equations are solved for the entire process, using a detailed chemical kinetic mechanism for propane. Simulation results show that the post-flame oxidation process takes place within a reactive layer where intermediate hydrocarbon products are formed at temperatures above 1100-1200 K, followed by a carbon monoxide conversion region closer to the hot burned gases. Model results show that most of hydrocarbons leaving the crevice are completely oxidized inside the cylinder. The largest contribution of remaining hydrocarbons are those leaving the crevice at temperatures below 1400 K. The largest fraction of non-fuel (intermediate) hydrocarbons results from hydrocarbons leaving the crevice when core temperatures are around 1400 K Copyright © 1997 Society of Automotive Engineers, Inc.

The effects of non-unity lewis numbers on turbulent premixed flame interactions in a twin V-flame configuration

The influence of Lewis number on turbulent premixed flame interactions is investigated using automatic feature extraction (AFE) applied to high-resolution flame simulation data. Premixed turbulent twin V-flames under identical turbulence conditions are simulated at global Lewis numbers of 0.4, 0.8, 1.0, and 1.2. Information on the position, frequency, and magnitude of the interactions is compared, and the sensitivity of the results to sample interval is discussed. It is found that both the frequency and magnitude of normal type interactions increases with decreasing Lewis number. Counternormal type interactions become more likely as the Lewis number increases. The variation in both the frequency and the magnitude of the interactions is found to be caused by large-scale changes in flame wrinkling resulting from differences in the thermo-diffusive stability of the flames. During flame interactions, thermo-diffusive effects are found to be insignificant due to the separation of time scales. © 2013 Copyright Taylor and Francis Group, LLC.

Spray flame structure of rapeseed biodiesel and Jet-A1 fuel

The spray combustion characteristics of rapeseed methyl esters (RME) were compared to Jet-A1 fuel using a gas turbine type combustor. The swirling spray flames for both fuels were established at a constant power output of 6 kW. The main swirling air flow was preheated to 350 C prior to coaxially enveloping the airblast-atomized liquid fuel spray at atmospheric pressure. Investigation of the fundamental spray combustion was performed via measurements of the fuel droplet sizes and velocities, gas phase flow fields and flame reaction zones. The spray flame droplets and flow fields in the combustors were characterised using phase Doppler anemometry (PDA) and particle imaging velocimetry (PIV) respectively. Flame chemiluminescence imaging was employed to identify the flame reaction zones. The highest droplet concentration zone extends along a 30 angle from the symmetry axis, inside the flame zone. Only small droplets(<17 μ) (<17 μm)are found around the centreline region, while larger droplets are found at the edge of the spray outside the flame reaction zone. RME exhibits spray characteristics similar to Jet-A1 but with droplet concentration and volume fluxes four times higher, consistent with the expected longer droplet evaporation timescale. The flow field characteristics for both RME and Jet-A1 spray flames are very similar despite the significantly different visible characteristics of the flame reaction zones. © 2013 Elsevier Ltd. All rights reserved.

CFD investigation of turbulent premixed flame response to transverse forcing

Screech is a high frequency oscillation that is usually characterized by instabilities caused by large-scale coherent flow structures in the wake of bluff-body flameholders and shear layers. Such oscillations can lead to changes in flame surface area which can cause the flame to burn unsteadily, but also couple with the acoustic modes and inherent fluid-mechanical instabilities that are present in the system. In this study, the flame response to hydrodynamic oscillations is analyzed in a controlled manner using high-fidelity Computational Fluid Dynamics (CFD) with an unsteady Reynolds-averaged Navier-Stokes approach. The response of a premixed flame with and without transverse velocity forcing is analyzed. When unforced, the flame is shown to exhibit a self-excitation that is attributed to the anti-symmetric shedding of vortices in the wake of the flameholder. The flame is also forced using two different kinds of low-amplitude out-of-phase inlet velocity forcing signals. The first forcing method is harmonic forcing with a single characteristic frequency, while the second forcing method involves a broadband forcing signal with frequencies in the range of 500 - 1000 Hz. For the harmonic forcing method, the flame is perturbed only lightly about its mean position and exhibits a limit cycle oscillation that is characteristic of the forcing frequency. For the broadband forcing method, larger changes in the flame surface area and detachment of the flame sheet can be seen. Transition to a complicated trajectory in the phase space is observed. When analyzed systematically with system identification methods, the CFD results, expressed in the form of the Flame Transfer Function (FTF) are capable of elucidating the flame response to the imposed perturbation. The FTF also serves to identify, both spatially and temporally, regions where the flame responds linearly and nonlinearly. Locking-in between the flame's natural self-excited frequency and the subharmonic frequencies of the broadband forcing signal is found to alter the dynamical behaviour of the flame. Copyright © 2013 by ASME.