Passive control of combustion instability in a low emissions aeroderivative gas turbine

This paper describes the conceptual ideas, the theoretical validation, the laboratory testing and the field trials of a recently patented fuel-air mixing device for use in high-pressure ratio, low emissions, gaseous-fueled gas turbines. By making the fuel-air mixing process insensitive to pressure fluctuations in the combustion chamber, it is possible to avoid the common problem of positive feedback between mixture strength and the unsteady combustion process. More specifically, a mixing duct has been designed such that fuel-air ratio fluctuations over a wide range of frequencies can be damped out by passive design means. By scaling the design in such a way that the range of damped frequencies covers the frequency spectrum of the acoustic modes in the combustor, the instability mechanism can be removed. After systematic development, this design philosophy was successfully applied to a 35:1 pressure ratio aeroderivative gas turbine yielding very low noise levels and very competitive NOx and CO measurements. The development of the new premixer is described from conceptual origins through analytic and CFD evaluation to laboratory testing and final field trials. Also included in this paper are comments about the practical issues of mixing, flashback resistance and autoignition.

Spray combustion characteristics of palm biodiesel

The potential of palm methyl esters (PME) as an alternative fuel for gas turbines is investigated using a swirl burner. The main air flow is preheated to 623 K, and a swirling spray flame is established at atmospheric pressure. The spray combustion characteristics of PME are compared to diesel and Jet-A1 fuel under the same burner power output of 6 kW. Investigation of the fuel atomizing characteristics using phase Doppler anemometry (PDA) shows that most droplets are distributed within the flame reaction zone region. PME droplets exhibit higher Sautermean diameter (SMD) values than baseline fuels, and thus higher droplet penetration length and longer evaporation timescales. The PME swirl flame presents a different visible flame reaction zone while combusting with low luminosity and produces no soot. NO x emissions per unit mass of fuel and per unit energy are reduced by using PME relative to those of conventional fuels. © 2012 Copyright Taylor and Francis Group, LLC.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

Combustion rumble prediction with integrated computational-fluid-dynamics/ low-order-model methods

The pressure oscillation within combustion chambers of aeroengines and industrial gas turbines is a major technical challenge to the development of high-performance and low-emission propulsion systems. In this paper, an approach integrating computational fluid dynamics and one-dimensional linear stability analysis is developed to predict the modes of oscillation in a combustor and their frequencies and growth rates. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone, which enables the computational fluid dynamics calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in agreement with experimental measurements. Furthermore, results from the computational fluid dynamics calculation provide the flame transfer function to describe unsteady heat release rate. Departures from ideal one-dimensional flows are described by shape factors. Combined with this information, low-order models can work out the possible oscillation modes and their initial growth rates. The approach developed here can be used in more general situations for the analysis of combustion oscillations. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Design of brushless doubly-fed machines based: On magnetic circuit modeling

This paper proposes a magnetic circuit model (MCM) for the design of a brushless doubly-fed machine (BDFM). The BDFM possesses advantages in terms of high reliability and reduced gearbox stages, and it requires a fractionally-rated power converter. This makes it suitable for utilization in offshore wind turbines. It is difficult for conventional design methods to calculate the flux in the stator because the two sets of stator windings, which have different pole number, form a complex flux pattern which is not easily determined using common analytical approaches. However, it is advantageous to predict the flux density in the teeth and air-gap at the initial design stage for sizing purposes without recourse finite element analysis. Therefore a magnetic circuit model is developed in this paper to calculate the flux density. A BDFM is used as a case study with FEA validation. © 1965-2012 IEEE.

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.

Experimental investigation of dynamic response of acoustically forced turbulent premixed CH4/CO2/air flames

Increasing demand for energy and continuing increase in environmental as well as financial cost of use of fossil fuels drive the need for utilization of fuels from sustainable sources for power generation. Development of fuel-flexible combustion systems is vital in enabling the use of sustainable fuels. It is also important that these sustainable combustion systems meet the strict governmental emission legislations. Biogas is considered as one of the viable sustainable fuels that can be used to power modern gas turbines: However, the change in chemical, thermal and transport properties as well as change in Wobbe index due to the variation of the fuel constituents can have a significant effect on the performance of the combustor. It is known that the fuel properties have strong influence on the dynamic flame response; however there is a lack of detailed information regarding the effect of fuel compositions on the sensitivity of the flames subjected to flow perturbations. In this study, we describe an experimental effort investigating the response of premixed biogas-air turbulent flames with varying proportions of CH4 and CO2 to velocity perturbations. The flame was stabilized using a centrally placed conical bluff body. Acoustic perturbations were imposed to the flow using loud speakers. The flame dynamics and the local heat release rate of these acoustically excited biogas flames were studied using simultaneous measurements of OH and H2CO planar laser induced fluorescence. OH* chemiluminescence along with acoustic pressure measurements were also recorded to estimate the total flame heat release modulation and the velocity fluctuations. The measurements were carried out by keeping the theoretical laminar flame speed constant while varying the bulk velocity and the fuel composition. The results indicate that the flame sensitivity to perturbations increased with increased dilution of CH4 by CO2 at low amplitude forcing, while at high amplitude forcing conditions the magnitude of the flame response was independent of dilution.

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.

An investigation into numerical analysis alternatives for predicting reingestion in turbine disc rim cavities

Reliable means of predicting ingestion in cavities adjacent to the main gas path are increasingly being sought by engineers involved in the design of gas turbines. In this paper, analysis is to be presented that results from an extended research programme, MAGPI, sponsored by the EU and several leading gas turbine manufactures and universities. Extensive use is made of CFD modelling techniques to understand the aerodynamic behaviour of a turbine stator well cavity, focusing on the interaction of cooling air supply with the main annulus gas. The objective of the study has been to benchmark a number of CFD codes and numerical techniques covering RANS and URANS calculations with different turbulence models in order to assess the suitability of the standard settings used in the industry for calculating the mechanics of the flow travelling between cavities in a turbine through the main gas path. The modelling methods employed have been compared making use of experimental data gathered from a dedicated two-stage turbine rig, running at engine representative conditions. Extensive measurements are available for a range of flow conditions and alternative cooling arrangements. The limitations of the numerical methods in calculating the interaction of the cooling flow egress and the main stream gas, and subsequent ingestion into downstream cavities in the engine (i.e. re-ingestion), have been exposed. This has been done without losing sight of the validation of the CFD for its use for predicting heat transfer, which was the main objective of the partners of the MAGPI Work- Package 1 consortium. Copyright © 2012 by ASME.

Multi-objective shape optimisation for horizontal-axis wind turbine blades

We are developing a wind turbine blade optimisation package CoBOLDT (COmputa- tional Blade Optimisation and Load De ation Tool) for the optimisation of large horizontal- axis wind turbines. The core consists of the Multi-Objective Tabu Search (MOTS), which controls a spline parameterisation module, a fast geometry generation and a stationary Blade Element Momentum (BEM) code to optimise an initial wind turbine blade design. The objective functions we investigate are the Annual Energy Production (AEP) and the fl apwise blade root bending moment (MY0) for a stationary wind speed of 50 m/s. For this task we use nine parameters which define the blade chord, the blade twist (4 parameters each) and the blade radius. Throughout the optimisation a number of binary constraints are defined to limit the noise emission, to allow for transportation on land and to control the aerodynamic conditions during all phases of turbine operation. The test case shows that MOTS is capable to find enhanced designs very fast and eficiently and will provide a rich and well explored Pareto front for the designer to chose from. The optimised blade de- sign could improve the AEP of the initial blade by 5% with the same flapwise root bending moment or reduce MY0 by 7.5% with the original energy yield. Due to the fast runtime of order 10 seconds per design, a huge number of optimisation iterations is possible without the need for a large computing cluster. This also allows for increased design flexibility through the introduction of more parameters per blade function or parameterisation of the airfoils in future. © 2012 by Nordex Energy GmbH.

Multi-objective shape optimisation for horizontal-axis wind turbine blades

We are developing a wind turbine blade optimisation package CoBOLDT (COmputa- tional Blade Optimisation and Load Deation Tool) for the optimisation of large horizontal- axis wind turbines. The core consists of the Multi-Objective Tabu Search (MOTS), which controls a spline parameterisation module, a fast geometry generation and a stationary Blade Element Momentum (BEM) code to optimise an initial wind turbine blade design. The objective functions we investigate are the Annual Energy Production (AEP) and the apwise blade root bending moment (MY0) for a stationary wind speed of 50 m/s. For this task we use nine parameters which define the blade chord, the blade twist (4 parameters each) and the blade radius. Throughout the optimisation a number of binary constraints are defined to limit the noise emission, to allow for transportation on land and to control the aerodynamic conditions during all phases of turbine operation. The test case shows that MOTS is capable to find enhanced designs very fast and efficiently and will provide a rich and well explored Pareto front for the designer to chose from. The optimised blade de- sign could improve the AEP of the initial blade by 5% with the same apwise root bending moment or reduce MY0 by 7.5% with the original energy yield. Due to the fast runtime of order 10 seconds per design, a huge number of optimisation iterations is possible without the need for a large computing cluster. This also allows for increased design flexibility through the introduction of more parameters per blade function or parameterisation of the airfoils in future. © 2012 AIAA.