Intermittency transport modeling of separated flow transition

An intermittency transport model is proposed for modeling separated-flow transition. The model is based on earlier work on prediction of attached flow bypass transition and is applied for the first time to model transition in a separation bubble at various degrees of free-stream turbulence. The model has been developed so that it takes into account the entrainment of the surrounding fluid. Experimental investigations suggest that it is this phenomena which ultimately determines the extent of the separation bubble. Transition onset is determined via a boundary layer correlation based on momentum thickness at the point of separation. The intermittent flow characteristic of the transition process is modeled via an intermittency transport equation. This accounts for both normal and streamwise variation of intermittency and hence models the entrainment of surrounding flow in a more accurate manner than alternative prescribed intermittency models. The model has been validated against the well established T3L semicircular leading edge flat plate test case for three different degrees of free-stream turbulence characteristic of turbomachinery blade applications.

The effect of wake passing on a flow separation in a low pressure turbine cascade

This paper describes an investigation into the effect that passing wakes have on a separation bubble that exists on the pressure surface and near the leading edge of a low pressure turbine blade. Previous experimental studies have shown that the behaviour of this separation is strongly incidence dependent and that it responds to its disturbance environment. The results presented in this paper examine the effect of wake passing in greater detail. Two dimensional, Reynolds averaged, numerical predictions are first used to examine qualitatively the unsteady interaction between the wakes and the separation bubble. The separation is predicted to consist of spanwise vortices whose development is in phase with the wake passing. However, comparison with experiments shows that the numerical predictions exaggerate the coherence of these vortices and also overpredict the time-averaged length of the separation. Nonetheless, experiments strongly suggest that the predicted phase locking of the vortices in the separation onto the wake passing is physical.

Deposition of small particles from turbulent flows

This paper deals with particle deposition onto solid walls from turbulent flows. The aim of the study is to model particle deposition in industrial flows, such as the one in gas turbines. The numerical study has been carried out with a two fluid approach. The possible contribution to the deposition from Brownian diffusion, turbulent diffusion and shear-induced lift force are considered in the study. Three types of turbulent two-phase flows have been studied: turbulent channel flow, turbulent flow in a bent duct and turbulent flow in a turbine blade cascade. In the turbulent channel flow case, the numerical results from a two-dimensional code show good agreement with numerical and experimental results from other resources. Deposition problem in a bent duct flow is introduced to study the effect of curvature. Finally, the deposition of small particles on a cascade of turbine blades is simulated. The results show that the current two fluid models are capable of predicting particle deposition rates in complex industrial flows.

Towards topology-free optimization - An application to turbine internal cooling geometries

The application of automated design optimization to real-world, complex geometry problems is a significant challenge - especially if the topology is not known a priori like in turbine internal cooling. The long term goal of our work is to focus on an end-to-end integration of the whole CFD Process, from solid model through meshing, solving and post-processing to enable this type of design optimization to become viable & practical. In recent papers we have reported the integration of a Level Set based geometry kernel with an octree-based cut- Cartesian mesh generator, RANS flow solver, post-processing & geometry editing all within a single piece of software - and all implemented in parallel with commodity PC clusters as the target. The cut-cells which characterize the approach are eliminated by exporting a body-conformal mesh guided by the underpinning Level Set. This paper extends this work still further with a simple scoping study showing how the basic functionality can be scripted & automated and then used as the basis for automated optimization of a generic gas turbine cooling geometry. Copyright © 2008 by W.N.Dawes.

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.

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.

Effects of multistage coupling and disk flexibility on mistuned bladed disk dynamics

The effects ofdisk flexibility and multistage coupling on the dynamics of bladed disks with and without blade mistuning are investigated. Both free and forced responses are examined using finite element representations of example single and two-stage rotor models. The reported work demonstrates the importance of proper treatment of interstage (stage-to-stage) boundaries in order to yield adequate capture of disk-blade modal interaction in eigenfrequency veering regions. The modified disk-blade modal interactions resulting from interstage-coupling-induced changes in disk flexibility are found to have a significant impact on (a) tuned responses due to excitations passing through eigenfrequency veering regions, and (b) a design's sensitivity to blade mistuning. Hence, the findings in this paper suggest that multistage analyses may be required when excitations are expected to fall in or near eigenfrequency veering regions or when the sensitivity to blade mistuning is to be accounted for Conversely, the observed sensitivity to disk flexibility also indicates that the severity of unfavorable structural interblade coupling may be reduced significantly by redesigning the disk(s) and stage-to-stage connectivity. The relatively drastic effects of such modifications illustrated in this work indicate that the design modifications required to alleviate veering-related response problems may be less comprehensive than what might have been expected.

Fast Response Exhaust Gas Temperature Measurement in IC Engines

The harsh environment presented by engines, particularly in the exhaust systems, often necessitates the use of robust and therefore low bandwidth temperature sensors. Consequently, high frequencies are attenuated in the output. One technique for addressing this problem involves measuring the gas temperature using two sensors with different time-constants and mathematically reconstructing the true gas temperature from the resulting signals. Such a technique has been applied in gas turbine, rocket motor and combustion research. A new reconstruction technique based on difference equations has been developed and its effectiveness proven theoretically. The algorithms have been successfully tested and proven on experimental data from a rig that produces cyclic temperature variations. These tests highlighted that the separation of the thermocouple junctions must be very small to ensure that both sensors are subjected to the same gas temperatures. Exhaust gas temperatures were recorded by an array of thermocouples during transient operation of a high performance two-stroke engine. The results show that the increase in bandwidth arising from the dual sensor technique allowed accurate measurement of exhaust gas temperature with relatively robust thermocouples. Finally, an array of very fine thermocouples (12.5 - 50 microns) was used to measure the in-cycle temperature variation in the exhaust.

Detached-Eddy Simulation of a Louver-Cooling Scheme for Turbine Blades

The performance of a louver-cooling scheme on a flat plate was analyzed using a detached-eddy-simulation turbulence model. It was assumed that the louver-cooling scheme was tested in a wind tunnel with the mainstream flow velocity of 20 m/s, equivalent to a Reynolds number of 16,200, based on the jet diameter. Turbulence closure was achieved by a realizable k-e-based detached-eddy-simulation turbulence model. Solutions of two blowing ratios of 0.5 and 1 were successfully obtained by running parallel on 16 nodes on a computer cluster. The flowfields were found to be highly unsteady and oscillatory in nature, with the maximum fluctuation of the adiabatic effectiveness as high as 15% of the time-averaged value. It is shown that the fluctuations in the adiabatic effectiveness are mainly caused by the spanwise fluctuation of the coolant jet and the unsteady vortical structures created by the interaction of the jet and the mainstream.