An advanced multi-configuration stator well cooling test facility

Optimisation of cooling systems within gas turbine engines is of great interest to engine manufacturers seeking gains in performance, efficiency and component life. The effectiveness of coolant delivery is governed by complex flows within the stator wells and the interaction of main annulus and cooling air in the vicinity of the rim seals. This paper reports the development of a test facility which allows the interaction of cooling air and main gas paths to be measured at conditions representative of those found in modern gas turbine engines. The test facility features a two stage turbine with an overall pressure ratio of approximately 2.6:1. Hot air is supplied to the main annulus using a Rolls-Royce Dart compressor driven by an aero-derivative engine plant. Cooling air can be delivered to the stator wells at multiple locations and at a range of flow rates which cover bulk ingestion through to bulk egress. The facility has been designed with adaptable geometry to enable rapid changes of cooling air path configuration. The coolant delivery system allows swift and accurate changes to the flow settings such that thermal transients may be performed. Particular attention has been focused on obtaining high accuracy data, using a radio telemetry system, as well as thorough through-calibration practices. Temperature measurements can now be made on both rotating and stationary discs with a long term uncertainty in the region of 0.3 K. A gas concentration measurement system has also been developed to obtain direct measurement of re-ingestion and rim seal exchange flows. High resolution displacement sensors have been installed in order to measure hot running geometry. This paper documents the commissioning of a test facility which is unique in terms of rapid configuration changes, non-dimensional engine matching and the instrumentation density and resolution. Example data for each of the measurement systems is presented. This includes the effect of coolant flow rate on the metal temperatures within the upstream cavity of the turbine stator well, the axial displacement of the rotor assembly during a commissioning test, and the effect of coolant flow rate on mixing in the downstream cavity of the stator well. Copyright © 2010 by ASME.

Electricity Storage Using a Thermal Storage Scheme

The increasing use of renewable energy technologies for electricity generation, many of which have an unpredictably intermittent nature, will inevitably lead to a greater demand for large-scale electricity storage schemes. For example, the expanding fraction of electricity produced by wind turbines will require either backup or storage capacity to cover extended periods of wind lull. This paper describes a recently proposed storage scheme, referred to here as Pumped Thermal Storage (PTS), and which is based on "sensible heat" storage in large thermal reservoirs. During the charging phase, the system effectively operates as a high temperature-ratio heat pump, extracting heat from a cold reservoir and delivering heat to a hot one. In the discharge phase the processes are reversed and it operates as a heat engine. The round- trip efficiency is limited only by process irreversibilities (as opposed to Second Law limitations on the coefficient of performance and the thermal efficiency of the heat pump and heat engine respectively). PTS is currently being developed in both France and England. In both cases, the schemes operate on the Joule-Brayton (gas turbine) cycle, using argon as the working fluid. However, the French scheme proposes the use of turbomachinery for compression and expansion, whereas for that being developed in England reciprocating devices are proposed. The current paper focuses on the impact of the various process irreversibilities on the thermodynamic round-trip efficiency of the scheme. Consideration is given to compression and expansion losses and pressure losses (in pipe-work, valves and thermal reservoirs); heat transfer related irreversibility in the thermal reservoirs is discussed but not included in the analysis. Results are presented demonstrating how the various loss parameters and operating conditions influence the overall performance.

A 5-step reduced mechanism for combustion of CO/H2/H2O/CH4/CO2 mixtures with low hydrogen/methane and high H2O content

In this study a 5-step reduced chemical kinetic mechanism involving nine species is developed for combustion of Blast Furnace Gas (BFG), a multi-component fuel containing CO/H2/CH4/CO2, typically with low hydrogen, methane and high water fractions, for conditions relevant for stationary gas-turbine combustion. This reduced mechanism is obtained from a 49-reaction skeletal mechanism which is a modified subset of GRI Mech 3.0. The skeletal and reduced mechanisms are validated for laminar flame speeds, ignition delay times and flame structure with available experimental data, and using computational results with a comprehensive set of elementary reactions. Overall, both the skeletal and reduced mechanisms show a very good agreement over a wide range of pressure, reactant temperature and fuel mixture composition. © 2012 The Combustion Institute..

Zonal RANS-LES modeling for turbines in aeroengines

The cost of large-eddy simulation (LES) modeling in various zones of gas turbine aeroengines is outlined. This high cost clearly demonstrates the need to perform hybrid Reynolds-averaged Navier-Stokes-LES (RANS-LES) over the majority of engine zones because the Reynolds number is too high for pure LES. The RANS layer is used to cover over the fine streaks found in the inner part of the boundary layer. The hybrid strategy is applied to various engine zones, which is shown to typically give much greater predictive accuracy than pure RANS simulations. However, the cost estimates show that the RANS layer should be disposed within the low-pressure turbine zone. Also, the nature of the flow physics in this zone makes LES most sensible. © 2012 by Begell House, Inc.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

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.

Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

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.

Hybrid-electric propulsion for aircraft

Against a background of increasing energy demand and rising fuel prices, hybrid-electric propulsion systems (HEPS) have the potential to significantly reduce fuel consumption in the aviation industry, particularly in the lighter sectors. By taking advantage of both Electric Motor (EM) and Internal Combustion Engine (ICE), HEPS provide not only a benefit in fuel saving but also a reduction in take-off noise and the emission levels. This research considers the design and sizing process of a hybrid-electric propulsion system for a single-seat demonstrator aircraft, the experimental derivation of the ICE map and the EM parameters. In addition to the experimental data, a novel modeling approach including several linked desktop PC software packages is presented to analyze and optimize hybrid-electric technology for aircraft. Further to the analysis of a parallel hybrid-electric, mid-scale aircraft, this paper also presents a scaling approach for a 20 kg UAV and a 50 tonne inter-city airliner. At the smaller scale, two different mission profiles are analyzed: an ISR mission profile, where the simulation routine optimizes the component size of the hybrid-electric propulsion system with respect to fuel saving, and a maximum duration profile; where the flight endurance is determined as a function of payload weight. At the larger scale, the performance of a 50 tonne inter-city airliner is modeled, based on a hybrid-electric gas-turbine, assuming a range of electric boost powers and battery masses.

Prediction of combustion noise for an aeroengine combustor

Combustion noise may become an important noise source for lean-burn gas turbine engines, and this noise is usually associated with highly unsteady flames. This work aims to compute the broadband combustion noise spectrum for a realistic aeroengine combustor and to compare with available measured noise data on a demonstrator aeroengine. A low-order linear network model is applied to a demonstrator engine combustor to obtain the transfer function that relates to unsteadiness in the rate of heat release, acoustic, entropic, and vortical fluctuations. A spectral model is used for the heat release rate fluctuation, which is the source of the noise. The mean flow of the aeroengine combustor required as input data to this spectral model is obtained from Reynolds-averaged Navier-Stokes simulations. The computed acoustic field for a low-to-medium power setting indicates that the models used in this study capture the main characteristics of the broadband spectral shape of combustion noise. Reasonable agreement with the measured spectral level is achieved. © 2012 AIAA.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.