Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

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.

A novel system for reducing turbo-lag by injection of compressed gas into the exhaust manifold

A key challenge in achieving good transient performance of highly boosted engines is the difficulty of accelerating the turbocharger from low air flow conditions (turbo lag). Multi-stage turbocharging, electric turbocharger assistance, electric compressors and hybrid powertrains are helpful in the mitigation of this deficit, but these technologies add significant cost and integration effort. Air-assist systems have the potential to be more cost-effective. Injecting compressed air into the intake manifold has received considerable attention, but the performance improvement offered by this concept is severely constrained by the compressor surge limit. The literature describes many schemes for generating the compressed gas, often involving significant mechanical complexity and/or cost. In this paper we demonstrate a novel exhaust assist system in which a reservoir is charged during braking. Experiments have been conducted using a 2.0 litre light-duty Diesel engine equipped with exhaust gas recirculation (EGR) and variable geometry turbine (VGT) coupled to an AC transient dynamometer, which was controlled to mimic engine load during in-gear braking and acceleration. The experimental results confirm that the proposed system reduces the time to torque during the 3rd gear tip-in by around 60%. Such a significant improvement was possible due to the increased acceleration of turbocharger immediately after the tip-in. Injecting the compressed gas into the exhaust manifold circumvents the problem of compressor surge and is the key enabler of the superior performance of the proposed concept. Copyright © 2013 SAE International.

Assessing boost-assist options for turbocharged engines using 1-D engine simulation and model predictive control

Delivering acceptable low end torque and good transient response is a significant challenge for all turbocharged engines. As downsized gasoline engines and Diesel engines make up a larger and larger proportion of the light-duty engines entering the market, the issue takes on greater significance. Several schemes have been proposed to improve torque response in highly boosted engines, including the use of electrical assist turbochargers and compressed air assist. In this paper we examine these methods with respect to their effectiveness in improving transient response and their relative performance along with some of the practical considerations for real world application. Results shown in this paper are from 1-D simulations using the Ricardo WAVE software package. The simulation model is based on a production light-duty Diesel engine modified to allow the introduction of compressed air at various points in the air-path as well as direct torque application to the turbocharger shaft (such as might be available from an electrical assist turbocharger). Whilst the 1-D simulation software provides a suitable environment for investigating the various boost assistance options, the overall air path performance also depends upon the control system. The introduction of boost assistance complicates the control in two significant ways: the system may run into constraints (such as compressor surge) that are not encountered in normal operation and the assistance introduces an additional control input. Production engine controllers are usually based on gain-scheduled PID control and extensive calibration. For this study, the non-linear nature of the engine together with the multiple configurations considered and the slower than real-time execution of 1-D models makes such an approach time consuming. Moreover, an ad-hoc approach would leave some doubt as to the fairness of comparisons between the different boost-assist options. Model Predictive Control has been shown to offer a convenient approach to controlling the 1-D simulations in a close to optimal manner for a typical Diesel VGT-EGR air path configuration. We show that the same technique can be applied to all the considered assistance methods with only modest calibration effort required. Copyright © 2012 SAE International.

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.

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.

Multi-objective optimisation of horizontal axis wind turbine structure and energy production using aerofoil and blade properties as design variables

The design of wind turbine blades is a true multi-objective engineering task. The aerodynamic effectiveness of the turbine needs to be balanced with the system loads introduced by the rotor. Moreover the problem is not dependent on a single geometric property, but besides other parameters on a combination of aerofoil family and various blade functions. The aim of this paper is therefore to present a tool which can help designers to get a deeper insight into the complexity of the design space and to find a blade design which is likely to have a low cost of energy. For the research we use a Computational Blade Optimisation and Load Deflation Tool (CoBOLDT) to investigate the three extreme point designs obtained from a multi-objective optimisation of turbine thrust, annual energy production as well as mass for a horizontal axis wind turbine blade. The optimisation algorithm utilised is based on Multi-Objective Tabu Search which constitutes the core of CoBOLDT. The methodology is capable to parametrise the spanning aerofoils with two-dimensional Free Form Deformation and blade functions with two tangentially connected cubic splines. After geometry generation we use a panel code to create aerofoil polars and a stationary Blade Element Momentum code to evaluate turbine performance. Finally, the obtained loads are fed into a structural layout module to estimate the mass and stiffness of the current blade by means of a fully stressed design. For the presented test case we chose post optimisation analysis with parallel coordinates to reveal geometrical features of the extreme point designs and to select a compromise design from the Pareto set. The research revealed that a blade with a feasible laminate layout can be obtained, that can increase the energy capture and lower steady state systems loads. The reduced aerofoil camber and an increased L/. D-ratio could be identified as the main drivers. This statement could not be made with other tools of the research community before. © 2013 Elsevier Ltd.

Aerodynamic design of high end wall angle turbine stages-part I: Methodology development

A design methodology is presented for turbines in an annulus with high end wall angles. Such stages occur where large radial offsets between the stage inlet and stage outlet are required, for example in the first stage of modern low pressure turbines, and are becoming more prevalent as bypass ratios increase. The turbine vanes operate within s-shaped ducts which result in meridional curvature being of a similar magnitude to the bladeto-blade curvature. Through a systematic series of idealized computational cases, the importance of two aspects of vane design are shown. First, the region of peak end wall meridional curvature is best located within the vane row. Second, the vane should be leant so as to minimize spanwise variations in surface pressure-this condition is termed "ideal lean." This design philosophy is applied to the first stage of a low pressure turbine with high end wall angles. © 2014 by ASME.

Improving intermediate pressure turbine performance by using a nonorthogonal stator

Intermediate pressure (IP) turbines in high bypass ratio civil aeroengines are characterized by a significant increase in radius and a low aspect ratio stator. Conventional aerodynamic designs for the IP turbine stator have had leading and trailing edges orthogonal to the hub and casing end walls. The IP turbine rotor, however, is stacked radially due to stress limits. These choices inevitably lead to a substantial gap between the IP stator and rotor at the outer diameter in a duct that is generally diffusing the flow due to the increasing radius. In this low Mach number study, the IP stator is redesigned, incorporating compound sweep so that the leading and trailing edges are no longer orthogonal to the end walls. Computational investigations showed that the nonorthogonal stator reduces the flow diffusion between the stator and rotor, which yields two benefits: The stator trailing edge velocity was reduced, as was the boundary layer growth on the casing end wall within the gap. Experimental measurements confirm that the turbine with the nonorthogonal stator has an increased efficiency by 0.49%, while also increasing the work output by 4.6%, at the design point. © 2014 by ASME.

Aerodynamic design of high end wall angle turbine stages- part II: Experimental verification

An experimental investigation of a turbine stage featuring very high end wall angles is presented. The initial turbine design did not achieve a satisfactory performance and the difference between the design predictions and the test results was traced to a large separated region on the rear suction-surface. To improve the agreement between computational fluid dynamics (CFD) and experiment, it was found necessary to modify the turbulence modeling employed. The modified CFD code was then used to redesign the vane, and the changes made are described. When tested, the performance of the redesigned vane was found to have much closer agreement with the predictions than the initial vane. Finally, the flowfield and performance of the redesigned stage are compared to a similar turbine, designed to perform the same duty, which lies in an annulus of moderate end wall angles. A reduction in stage efficiency of at least 2.4% was estimated for the very high end wall angle design. © 2014 by ASME.

Unsteady numerical study of wet steam flow in a low pressure steam turbine

In steam power plants condensation already starts in the flow path of the low pressure part of the steam turbine, which leads to a complex three-dimensional two-phase flow. Wetness losses are caused due to thermodynamic and mechanical relaxation processes during condensation and droplet transport. The present investigation focuses on the unsteady effects due to rotor-stator interaction on the droplet formation process. Results of unsteady three dimensional flow simulations of a two-stage steam turbine are presented, whereby this is the first time that non-equilibrium condensation is considered in such simulations. The numerical approach is based on RANS equations, which are extended by a wet steam specific nucleation and droplet growth model. Despite the use of a high performance cluster the unsteady simulation has a considerably high simulation time of approximately 60 days by use of 48 CPUs. © Springer-Verlag Berlin Heidelberg 2012.