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.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

Obtaining bifurcation diagrams with a thermoacoustic network model

Linear techniques can predict whether the non-oscillating (steady) state of a thermoacoustic system is stable or unstable. With a sufficiently large impulse, however, a thermoacoustic system can reach a stable oscillating state even when the steady state is also stable. A nonlinear analysis is required to predict the existence of this oscillating state. Continuation methods are often used for this but they are computationally expensive. In this paper, an acoustic network code called LOTAN is used to obtain the steady and the oscillating solutions for a horizontal Rijke tube. The heat release is modelled as a nonlinear function of the mass flow rate. Several test cases from the literature are analysed in order to investigate the effect of various nonlinear terms in the flame model. The results agree well with the literature, showing that LOTAN can be used to map the steady and oscillating solutions as a function of the control parameters. Furthermore, the nature of the bifurcation between steady and oscillating states can be predicted directly from the nonlinear terms inside the flame model. 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.

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.

Computational engineering design for micro-scale combustors

At present, optimisation is an enabling technology in innovation. Multi-objective and multi-disciplinary design tools are essential in the engineering design process, and have been applied successfully in aerospace and turbomachinery applications extensively. These approaches give insight into the design space and identify the trade-offs between the competing performance measures satisfying a number of constraints at the same time. It is anticipated here that the same benefits can be obtained for the design of micro-scale combustors. In this paper, a multi-disciplinary automated design optimisation system was developed for this purpose, which comprises a commercial computational fluid dynamics package and a multi-objective variant of the Tabu Search optimisation algorithm. The main objectives that are considered in this study are to optimise the main micro-scale combustor design characteristics and to satisfy manufacturability considerations from the very beginning of the whole design operation. Hydrogen-air combustion as well as 14 geometrical and 2 operational parameters are used to describe and model the design problem. Two illustrative test cases will be presented, in which the most important device operational requirements are optimised, and the efficiency of the developed optimisation system is demonstrated. The identification, assessment and suitability of the optimum design configurations are discussed in detail. Copyright © 2012 by ASME.

Angle-of-attack effects on counter-rotating propellers at take-off

Due to their potential for significant fuel consumption savings, Counter-Rotating Open Rotors (CRORs) are currently being considered as an alternative to high-bypass turbofans. When CRORs are mounted on an aircraft, several 'installation effects' arise which are not present when the engine is operated in isolation. This paper investigates how flow features arising from one such effect - The angle-of-attack of the engine centre-line relative to the oncoming flow - can influence the design of CROR engines. Three-dimensional full-annulus unsteady CFD simulations are used to predict the time-varying flow field experienced by each rotor and emphasis is put on the interaction of the frontrotor wake and tip vortex with the rear-rotor. A parametric study is presented that quantifies the rotorrotor interaction as a function of the angle-of-attack. It is shown that angle-of-attack operation significantly changes the flow field and the unsteady lift on both rotors. In particular, a frequency analysis shows that the unsteady lift exhibits sidebands around the rotor-rotor interaction frequencies. Further, a non-linear increase in the total rear-rotor tip unsteadiness is observed for moderate and high angles-of-attack. The results presented in this paper demonstrate that common techniques used to mitigate CROR noise, such as modifying the rotor-rotor axial spacing and rear-rotor crop, can not be applied correctly unless angle-of-attack effects are taken into account. Copyright © 2012 by ASME.

An experimental study of loss sources in a fan operating with continuous inlet stagnation pressure distortion

The viability of Boundary Layer Ingesting (BLI) engines for future aircraft propulsion is dependent on the ability to design robust, efficient engine fan systems for operation with continuously distorted inlet flow. A key step in this process is to develop an understanding of the specific mechanisms by which an inlet distortion affects the performance of a fan stage. In this paper, detailed full-annulus experimental measurements of the flow field within a low-speed fan stage operating with a continuous 60-degree inlet stagnation pressure distortion are presented. These results are used to describe the three-dimensional fluid mechanics governing the interaction between the fan and the distortion and to make a quantitative assessment of the impact on loss generation within the fan. A 5.3 percentage point reduction in stage total-to-total efficiency is observed as a result of the inlet distortion. The reduction in performance is shown to be dominated by increased loss generation in the rotor due to off-design incidence values at its leading edge, an effect which occurs throughout the annulus despite the localised nature of the inlet distortion. Increased loss generation in the stator row is also observed due to flow separations that are shown to be caused by whirl angle distortion at rotor exit. By addressing these losses, it should be possible to achieve improved efficiency in BLI fan systems. 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.

A numerical study into the effects of Bio-synthetic Paraffinic Kerosine blends with Jet-A fuel for civil aircraft engine

Growing concerns regarding fluctuating fuel costs and pollution targets for gas emissions, have led the aviation industry to seek alternative technologies to reduce its dependency on crude oil, and its net emissions. Recently blends of bio-fuel with kerosine, have become an alternative solution as they offer "greener" aircraft and reduce demand on crude oil. Interestingly, this technique is able to be implemented in current aircraft as it does not require any modification to the engine. Therefore, the present study investigates the effect of blends of bio-synthetic paraffinic kerosine with Jet-A in a civil aircraft engine, focusing on its performance and exhaust emissions. Two bio-fuels are considered: Jatropha Bio-synthetic Paraffinic Kerosine (JSPK) and Camelina Bio-synthetic Paraffinic Kerosine (CSPK); there are evaluated as pure fuels, and as 10% and 50% blend with Jet-A. Results obtained show improvement in thrust, fuel flow and SFC as composition of bio-fuel in the blend increases. At design point condition, results on engine emissions show reduction in NO x, and CO, but increases of CO is observed at fixed fuel condition, as the composition of bio-fuel in the mixture increases. Copyright © 2012 by ASME.