The behaviour of laminar stratified methane/air flames in counterflow

Flames propagating through a mixture with a gradient of equivalence ratio have been previously demonstrated to travel faster or slower than their equivalent premixed flames. The present study aims to numerically investigate the response of strained laminar methane-air flames to such gradients. The flames are simulated in a counterflow configuration where a premixed reactant stream at equivalence ratio φR opposes a hot equilibrium stream at equivalence ratio φP. Premixed and stratified flames are compared with respect to the equivalence ratio φ* and the corresponding gradient ∇φ* at the point of peak heat release rate, for three strain rates, a=50, 300 and 500s-1 and a range of φ*. The effect of different stratification levels is also investigated by varying the ratio of φP to φR, Θ. Results indicate that, as long as flames stabilize within the diffusion layer and Θ>1, increased heat release rate Q is seen throughout the progress variable space in comparison to the premixed state. In contrast, an attenuation of heat release rate is seen for Θ<1. The enhancement (or attenuation) of heat release varies monotonically with Θ. The effect of stratification on flame behavior becomes more pronounced as the strain rate increases. The present study reveals the mechanisms for the propagation of quasi-steady stratified flames under lean and rich conditions: stratified flames are primarily dominated by the diffusion of heat under lean conditions, and diffusion of H2 under rich conditions. Thanks to species and thermal support, stratified flames continue to burn beyond the premixed lean and rich flammability limits. Further investigation on the unsteady response of flames to the fluctuating equivalence ratio implies that the steady results represent the unsteady response well, as long as φ* and ∇φ* are similar in both steady and unsteady cases. © 2013 The Combustion Institute.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

Absolute and convective instability in gas turbine fuel injectors

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier-Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re - 106, and a lean premixed gas turbine injector with five swirling streams at Re - 106. Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes. Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable. Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability. Copyright © 2012 by ASME.

On Mixing and Shaking : Structure and Dynamics of Turbulent Stratified Flames

The drive for low emission combustion systems encourages applications using premixed flames. Yet in many applications, considerations of flame stability or mixing times lead to systems with neither premixed nor diffusion flames, which are often called technically premixed or stratified flames. In this talk we discuss the current state of understanding of the effect of mixing and extent of stratification on the structure, microstructure and dynamics of selected turbulent stratified flames. Over the past few years, a significant database of scalar and velocity data has been built to analyze the effects of unmixedness on local and global flame structure. Microscale studies of the flame structures show in detail how the effect of local stratification affects (or not!) the flame structure, flame surface density and scalar dissipation rates, and production of selected species. The experiments place exacting demands on current spectroscopic diagnostics, and reveal the progress and limits to our understanding of turbulent flames in general. The dynamics of stratified flames with respect to instabilities is also shown to be very rich, as the particular shape of the flames and the stabilization points are is significantly affected by the fuel distribution, modifying the rate and location of heat release, and thus the coupling with the surrounding acoustics and determining the onset of self-excitations.

Application of scalar dissipation rate model to siemens DLE combustors

The standard design process for the Siemens Industrial Turbomachinery, Lincoln, Dry Low Emissions combustion systems has adopted the Eddy Dissipation Model with Finite Rate Chemistry for reacting computational fluid dynamics simulations. The major drawbacks of this model have been the over-prediction of temperature and lack of species data limiting the applicability of the model. A novel combustion model referred to as the Scalar Dissipation Rate Model has been developed recently based on a flamelet type assumption. Previous attempts to adopt the flamelet philosophy with alternative closure models have failed, with the prediction of unphysical phenomenon. The Scalar Dissipation Rate Model (SDRM) was developed from a physical understanding of scalar dissipation rate, signifying the rate of mixing of hot and cold fluids at scales relevant to sustain combustion, in flames and was validated using direct numerical simulations data and experimental measurements. This paper reports on the first industrial application of the SDRM to SITL DLE combustion system. Previous applications have considered ideally premixed laboratory scale flames. The industrial application differs significantly in the complexity of the geometry, unmixedness and operating pressures. The model was implemented into ANSYS-CFX using their inbuilt command language. Simulations were run transiently using Scale Adaptive Simulation turbulence model, which switches between Large Eddy Simulation and Unsteady Reynolds Averaged Navier Stokes using a blending function. The model was validated in a research SITL DLE combustion system prior to being applied to the actual industrial geometry at real operating conditions. This system consists of the SGT-100 burner with a glass square-sectioned combustor allowing for detailed diagnostics. This paper shows the successful validation of the SDRM against time averaged temperature and velocity within measurement errors. The successful validation allowed application of the SDRM to the SGT-100 twin shaft at the relevant full load conditions. Limited validation data was available due to the complexity of measurement in the real geometry. Comparison of surface temperatures and combustor exit temperature profiles showed an improvement compared to EDM/FRC model. Furthermore, no unphysical phenomena were predicted. This paper presents the successful application of the SDRM to the industrial combustion system. The model shows a marked improvement in the prediction of temperature over the EDM/FRC model previously used. This is of significant importance in the future applications of combustion CFD for understanding of hardware mechanical integrity, combustion emissions and dynamics of the flame. Copyright © 2012 by ASME.

Modal dynamics of self-excited azimuthal instabilities in an annular combustion chamber

In this paper we describe the time-varying amplitude and its relation to the global heat release rate of self-excited azimuthal instabilities in a simple annular combustor operating under atmospheric conditions. The combustor was modular in construction consisting of either 12, 15 or 18 equally spaced premixed bluff-body flames around a fixed circumference, enabling the effect of large-scale interactions between adjacent flames to be investigated. High-speed OH* chemiluminescence imaged from above the annulus and pressure measurements obtained at multiple locations around the annulus revealed that the limit cycles of the modes are degenerate in so much as they undergo continuous transitions between standing and spinning modes in both clockwise (CW) and anti-clockwise (ACW) directions but with the same resonant frequency. Similar behaviour has been observed in LES simulations which suggests that degenerate modes may be a characteristic feature of self-excited azimuthal instabilities in annular combustion chambers. By modelling the instabilities as two acoustic waves of time-varying amplitude travelling in opposite directions we demonstrate that there is a statistical prevalence for either standing m=1 or spinning m=±1 modes depending on flame spacing, equivalence ratio, and swirl configuration. Phase-averaged OH* chemiluminescence revealed a possible mechanism that drives the direction of the spinning modes under limit-cycle conditions for configurations with uniform swirl. By dividing the annulus into inner and outer annular regions it was found that the spin direction coincided with changes in the spatial distribution of the peak heat release rate relative to the direction of the bulk swirl induced along the annular walls. For standing wave modes it is shown that the globally integrated fluctuations in heat release rate vary in magnitude along the acoustic mode shape with negligible contributions at the pressure nodes and maximum contributions at the pressure anti-nodes. © 2013.

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.

Scalar gradient behaviour in MILD combustion

The results of three-dimensional Direct Numerical Simulation (DNS) of Moderate, Intense Low-oxygen Dilution (MILD) and conventional premixed turbulent combustion conducted using a skeletal mechanism including the effects of non-unity Lewis numbers and temperature dependent transport properties are analysed to investigate combustion characteristics using scalar gradient information. The DNS data is also used to synthesise laser induced fluorescence (LIF) signals of OH, CH2O, and CHO. These signals are analysed to verify if they can be used to study turbulent MILD combustion and it has been observed that at least two (OH and CH2O) LIF signals are required since the OH increase across the reaction zone is smaller in MILD combustion compared to premixed combustion. The scalar gradient PDFs conditioned on the reaction rate obtained from the DNS data and synthesised LIF signals suggests a strong gradient in the direction normal to the MILD reaction zone with moderate reaction rate implying flamelet combustion. However, the PDF of the normal gradient is as broad as for the tangential gradient when the reaction rate is high. This suggests a non-flamelet behaviour, which is due to interaction of reaction zones. The analysis of the conditional PDFs for the premixed case confirms the expected behaviour of scalar gradient in flamelet combustion. It has been shown that the LIF signals synthesised using 2D slices of DNS data also provide very similar insights. These results demonstrate that the so-called flameless combustion is not an idealised homogeneous reactive mixture but has common features of conventional combustion while containing distinctive characteristics. © 2013 The Combustion Institute.

Adaptive combustion instability control with saturation

The control of a class of combustion systems, suceptible to damage from self-excited combustion oscillations, is considered. An adaptive stable controller, called Self-Tuning Regulator (STR), has recently been developed, which meets the apparently contradictory challenge of relying as little as possible on a particular combustion model while providing some guarantee that the controller will cause no harm. The controller injects some fuel unsteadily into the burning region, thereby altering the heat release, in response to an input signal detecting the oscillation. This paper focuses on an extension of the STR design, when, due to stringent emission requirements and to the danger of flame extension, the amount of fuel used for control is limited in amplitude. A Lyapunov stability analysis is used to prove the stability of the modified STR when the saturation constraint is imposed. The practical implementation of the modified STR remains straightforward, and simulation results, based on the nonlinear premixed flame model developed by Dowling, show that in the presence of a saturation constraint, the self-excited oscillations are damped more rapidly with the modified STR than with the original STR. © 2001 by S. Evesque. Published by the American Institute of Aeronautics and Astronautics, Inc.

Femtosecond laser induced breakdown for combustion diagnostics

The focused beam of a 100 fs, 800 nm laser is used to induce a spark in some laminar premixed air-methane flames operating with variable fuel content (equivalence ratio). The analysis of the light escaping from the plasma revealed that the Balmer hydrogen lines, H α and H β, and some molecular origin emissions were the most prominent spectral features, while the CN (B 2Σ +-X 2Σ +) band intensity was found to depend linearly with methane content, suggesting that femtosecond laser induced breakdown spectroscopy can be a useful tool for the in-situ determination and local mapping of fuel content in hydrocarbon-air combustible mixtures. © 2012 American Institute of Physics.

Modelling of turbulent lifted jet flames using flamelets: A priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z-c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z-c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences. © 2014 Taylor & Francis.

Temperature measurements of the bluff body surface of a Swirl Burner using phosphor thermometry

Flames are often stabilised on bluff-bodies, yet their surface temperatures are rarely measured. This paper presents temperature measurements for the bluff body surface of the Cambridge/Sandia Stratified Swirl Burner. The flame is stabilized by a bluff body, designed to provide a series of turbulent premixed and stratified methane/air flames with a variable degree of swirl and stratification. Recently, modellers have raised concerns about the role of surface temperature on the resulting gas temperatures and the overall heat loss of the burner. Laser-induced phosphorescence is used to measure surface temperatures, with Mg4GeO6F:Mn as the excitation phosphor, creating a spatially resolved temperature map. Results show that the temperature of the bluff body is in the range 550-900 K for different operating conditions. The temperature distribution is strongly correlated with the degree of swirl and local equivalence ratio, reflecting the temperature distribution obtained in the gas phase. The overall heat loss represents only a small fraction (<0.5%) of the total heat load, yet the local surface temperature may affect the local heat transfer and gas temperatures. © 2014 The Combustion Institute.

Thermal performance of a micro-combustor for micro-gas turbine system

Premixed combustion of hydrogen gas and air was performed in a stainless steel based micro-annular combustor for a micro-gas turbine system. Micro-scale combustion has proved to be stable in the micro-combustor with a gap of 2 mm. The operating range of the micro-combustor was measured, and the maximum excess air ratio is up to 4.5. The distribution of the outer wall temperature and the temperature of exhaust gas of the micro-conbustor with excess air ratio were obtained, and the wall temperature of the micro-combustor reaches its maximum value at the excess air ratio of 0.9 instead of 1 (stoichiometric ratio). The heat loss of the micro-combustor to the environment was calculated and even exceeds 70% of the total thermal power computed from the consumed hydrogen mass flow rate. Moreover, radiant hunt transfer covers a large fraction of the total heat loss. Measures used to reduce the heat loss were proposed to improve the thermal performance of the micro-combustor. The optimal operating status of the micro-combustor and micro-gas turbine is analyzed and proposed by analyzing the relationship of the temperature of the exhaust gas of the micro-combustor with thermal power and excess air ratio. The investigation of the thermal performance of the micro-combustor is helpful to design an improved microcombustor.

Spatial resolved temperature measurement based on absorption spectroscopy using a single tunable diode laser

A novel method based on wavelength-multiplexed line-of-sight absorption and profile fitting for non-uniform flow field measurement is reported. A wavelength scanning combing laser temperature and current modulation WMS scheme is used to implement the wavelength-multiplexed-profile fitting method. Second harmonic (2f) signal of eight H2O transitions features near 7,170 cm(-1) are measured in one period using a single tunable diode laser. Spatial resolved temperature distribution upon a CH4/air premixed flat flame burner is obtained. The result validates the feasibility of strategy for non-uniform flow field diagnostics by means of WMS-2f TDLAS.