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.

Investigation of flame stretch in turbulent lifted jet flame

DNS data of a laboratory-scale turbulent lifted hydrogen jet flame has been analyzed to show that this flame has mixed mode combustion not only at the flame base but also in downstream locations. The mixed mode combustion is observed in instantaneous structures as in earlier studies and in averaged structure, in which the predominant mode is found to be premixed combustion with varying equivalence ratio. The non-premixed combustion in the averaged structure is observed only in a narrow region at the edge of the jet shear layer. The analyzes of flame stretch show large probability for negative flame stretch leading to negative surface averaged flame stretch. The displacement speed-curvature correlation is observed to be negative contributing to the negative flame stretch and partial premixing resulting from jet entrainment acts to reduce the negative correlation. The contribution of turbulent straining to the flame stretch is observed to be negative when the scalar gradient aligns with the most extensive principal strain rate. The physics behind the negative flame stretch resulting from turbulent straining is discussed and elucidated through a simple analysis of the flame surface density transport equation. © 2014 Copyright Taylor and Francis Group, LLC.

A flame surface density approach to large-eddy simulation of premixed turbulent combustion

The flame surface density approach to the modeling of premixed turbulent combustion is well established in the context of Reynolds-averaged simulations. For the future, it is necessary to consider large-eddy simulation (LES), which is likely to offer major advantages in terms of physical accuracy, particularly for unsteady combustion problems. LES relies on spatial filtering for the removal of unresolved phenomena whose characteristic length scales are smaller than the computational grid scale. Thus, there is a need for soundly based physical modeling at the subgrid scales. The aim of this paper is to explore the usefulness of the flame surface density concept as a basis for LES modeling of premixed turbulent combustion. A transport equation for the filtered flame surface density is presented, and models are proposed for unclosed terms. Comparison with Reynolds-averaged modeling is shown to reveal some interesting similarities and differences. These were exploited together with known physics and statistical results from experiment and from direct numerical stimulation in order to gain insight and refine the modeling. The model has been implemented in a combustion LES code together with standard models for scalar and momentum transport. Computational results were obtained for a simple three-dimensional flame propagation test problem, and the relative importance of contributing terms in the modeled equation for flame surface density was assessed. Straining and curvature are shown to have a major influence at both the resolved and subgrid levels.

Interaction of a premixed flame with a liquid fuel film on a wall

© 2004 The Combustion Institute. Published by Elsevier Inc. All rights reserved. In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.

Heat release rate measurement in turbulent flames

The dominant industrial approach for the reduction of NO x emissions in industrial gas turbines is the lean pre-mixed prevaporized concept. The main advantage of this concept is the lean operation of the combustion process; this decreases the heat release rate from the flame and results in a reduction in operating temperature. The direct measurement of heat release rates via simultaneous laser induced fluorescence of OH and CH 2O radicals using planar laser induced fluorescence. The product of the two images correlated with the forward production rate of the HCO radical, which in turn has correlated well with heat release rates from premixed hydrocarbon flames. The experimental methodology of the measurement of heat release rate and applications in different turbulent premixed flames were presented. This is an abstract of a paper presented at the 7th World Congress of Chemical Engineering (Glasgow, Scotland 7/10-14/2005).

Forcing of self-excited round jet diffusion flames

In this experimental and numerical study, two types of round jet are examined under acoustic forcing. The first is a non-reacting low density jet (density ratio 0.14). The second is a buoyant jet diffusion flame at a Reynolds number of 1100 (density ratio of unburnt fluids 0.5). Both jets have regions of strong absolute instability at their base and this causes them to exhibit strong self-excited bulging oscillations at welldefined natural frequencies. This study particularly focuses on the heat release of the jet diffusion flame, which oscillates at the same natural frequency as the bulging mode, due to the absolutely unstable shear layer just outside the flame. The jets are forced at several amplitudes around their natural frequencies. In the non-reacting jet, the frequency of the bulging oscillation locks into the forcing frequency relatively easily. In the jet diffusion flame, however, very large forcing amplitudes are required to make the heat release lock into the forcing frequency. Even at these high forcing amplitudes, the natural mode takes over again from the forced mode in the downstream region of the flow, where the perturbation is beginning to saturate non-linearly and where the heat release is high. This raises the possibility that, in a flame with large regions of absolute instability, the strong natural mode could saturate before the forced mode, weakening the coupling between heat release and incident pressure perturbations, hence weakening the feedback loop that causes combustion instability. © 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Modelling soot formation in a DISI engine

In this work, the formation of soot in a Direct Injection Spark Ignition (DISI) engine is simulated using the Stochastic Reactor Model (SRM) engine code. Volume change, convective heat transfer, turbulent mixing, direct injection and flame propagation are accounted for. In order to simulate flame propagation, the cylinder is divided into an unburned, entrained and burned zone, with the rate of entrainment being governed by empirical equations but combustion modelled with chemical kinetics. The model contains a detailed chemical mechanism as well as a highly detailed soot formation model, however computation times are relatively short. The soot model provides information on the morphology and chemical composition of soot aggregates along with bulk quantities, including soot mass, number density, volume fraction and surface area. The model is first calibrated by simulating experimental data from a Gasoline Direct Injection (GDI) Spark Ignition (SI) engine. The model is then used to simulate experimental data from the literature, where the numbers, sizes and derived mass particulate emissions from a 1.83 L, 4-cylinder, 4 valve production DISI engine were examined. Experimental results from different injection and spark timings are compared with the model and the qualitative trends in aggregate size distribution and emissions match the exhaust gas measurements well. © 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.