Multi-objective design and optimisation of bypass outlet-guide vanes

This paper describes the development of an automated design optimization system that makes use of a high fidelity Reynolds-Averaged CFD analysis procedure to minimize the fan forcing and fan BOGV (bypass outlet guide vane) losses simultaneously taking into the account the down-stream pylon and RDF (radial drive fairing) distortions. The design space consists of the OGV's stagger angle, trailing-edge recambering, axial and circumferential positions leading to a variable pitch optimum design. An advanced optimization system called SOFT (Smart Optimisation for Turbomachinery) was used to integrate a number of pre-processor, simulation and in-house grid generation codes and postprocessor programs. A number of multi-objective, multi-point optimiztion were carried out by SOFT on a cluster of workstations and are reported herein.

The effects of a trip wire and unsteadiness on a high speed highly loaded low-pressure turbine blade

This paper presents the effect of a single spanwise 2D wire upon the downstream position of boundary layer transition under steady and unsteady inflow conditions. The study is carried out on a high turning, high-speed, low pressure turbine (LPT) profile designed to take account of the unsteady flow conditions. The experiments were carried out in a transonic cascade wind tunnel to which a rotating bar system had been added. The range of Reynolds and Mach numbers studied includes realistic LPT engine conditions and extends up to the transonic regime. Losses are measured to quantify the influence of the roughness with and without wake passing. Time resolved measurements such as hot wire boundary layer surveys and surface unsteady pressure are used to explain the state of the boundary layer. The results suggest that the effect of roughness on boundary layer transition is a stability governed phenomena, even at high Mach numbers. The combination of the effect of the roughness elements with the inviscid Kelvin-Helmholtz instability responsible for the rolling up of the separated shear layer (Stieger [1]) is also examined. Wake traverses using pneumatic probes downstream of the cascade reveal that the use of roughness elements reduces the profile losses up to exit Mach numbers of 0.8. This occurs with both steady and unsteady inflow conditions.

Heat transfer and aerodynamics of turbine blade tips in a linear cascade

Local measurements of the heat transfer coefficient and pressure coefficient were conducted on the tip and near tip region of a generic turbine blade in a five-blade linear cascade. Two tip clearance gaps were used: 1.6% and 2.8% chord. Data was obtained at a Reynolds number of 2.3 × 10 5 based on exit velocity and chord. Three different tip geometries were investigated: a flat (plain) tip, a suction-side squealer, and a cavity squealer. The experiments reveal that the flow through the plain gap is dominated by flow separation at the pressure-side edge and that the highest levels of heat transfer are located where the flow reattaches on the tip surface. High heat transfer is also measured at locations where the tip-leakage vortex has impinged onto the suction surface of the aerofoil. The experiments are supported by flow visualisation computed using the CFX CFD code which has provided insight into the fluid dynamics within the gap. The suction-side and cavity squealers are shown to reduce the heat transfer in the gap but high levels of heat transfer are associated with locations of impingement, identified using the flow visualisation and aerodynamic data. Film cooling is introduced on the plain tip at locations near the pressure-side edge within the separated region and a net heat flux reduction analysis is used to quantify the performance of the successful cooling design. copyright © 2005 by ASME.

The effect of work processes on the casing heat transfer of a transonic turbine

This paper considers the effect of the rotor tip on the casing heat load of a transonic axial flow turbine. The aim of the research is to understand the dominant causes of casing heat-transfer. Experimental measurements were conducted at engine-representative Mach number, Reynolds number and stage inlet to casing wall temperature ratio. Time-resolved heat-transfer coefficient and gas recovery temperature on the casing were measured using an array of heat-transfer gauges. Time-resolved static pressure on the casing wall was measured using Kulite pressure transducers. Time-resolved numerical simulations were undertaken to aid understanding of the mechanism responsible for casing heat load. The results show that between 35% and 60% axial chord the rotor tip-leakage flow is responsible for more than 50% of casing heat transfer. The effects of both gas recovery temperature and heat transfer coefficient were investigated separately and it is shown that an increased stagnation temperature in the rotor tip gap dominates casing heat-transfer. In the tip gap the stagnation temperature is shown to rise above that found at stage inlet (combustor exit) by as much as 35% of stage total temperature drop. The rise in stagnation temperature is caused by an isentropic work input to the tip-leakage fluid by the rotor. The size of this mechanism is investigated by computationally tracking fluid path-lines through the rotor tip gap to understand the unsteady work processes that occur. Copyright © 2005 by ASME.

LES and hybrid LES/RANS simulations for conjugate heat transfer over a matrix of cubes

Large Eddy Simulation (LES) and a novel k -l based hybrid LES/RANS approach have been applied to simulate a conjugate heat transfer problem involving flow over a matrix of surface mounted cubes. In order to assess the capability and reliability of the newly developed k -l based hybrid LES/RANS, numerical results are compared with new LES and existing RANS results. Comparisons include mean velocity profiles, Reynolds stresses and conjugate heat transfer. As well as for hybrid LES/RANS validation purposes, the LES results are used to gain insights into the complex flow physics and heat transfer mechanisms. Numerical simulations show that the hybrid LES/RANS approach is effective. Mean and instantaneous fluid temperatures adjacent to the cube surface are found to strongly correlate with flow structure. Although the LES captures more mean velocity field complexities, broadly time averaged wake temperature fields are found similar for the LES and hybrid LES/RANS. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Vortex generators near shock/ boundary layer interactions

Experiments have been performed in a blowdown supersonic wind tunnel to investigate the effect of arrays of sub-boundary layer vortex generators placed upstream of a normal shock/ boundary layer interaction. The investigation makes use of a recovery shock wave and the naturally grown turbulent boundary layer on the wind tunnel floor. Experiments were performed at Mach numbers of 1.5 and 1.3 and a freestream Reynolds number of 28 × 106. Two types of vortex generators were investigated - wedge-shaped and arrays of counter-rotating vanes. It was found that at Mach 1.5 the vane-type VGs eliminated and the wedge-type VGs greatly reduced the separation bubble under the shock. When placed in the supersonic part of the flow both VGs caused a wave pattern consisting of a shock, re-expansion and shock. The re-expansion and double shocks are undesirable features since they equate to increased total pressure losses and hence increased -wave drag. Furthermore there are indications that the vortex intensity is reduced by the normal shock/ boundary layer interaction. When the shock was located directly over the VGs there was no re-expansion present, but the 'damping' effect of the shock on the vortex persisted. It appears that the vortices produced by the wedge-shaped VGs lift off the surface more rapidly. Similar results were observed at Mach 1.3, where the flow was unseparated.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Adjoint linearised Euler solver in the frequency domain for jet noise modelling

The paper is devoted to extending the new efficient frequency-domain method of adjoint Green's function calculation to curvilinear multi-block RANS domains for middle and farfield sound computations. Numerical details of the method such as grids, boundary conditions and convergence acceleration are discussed. Two acoustic source models are considered in conjunction with the method and acoustic modelling results are presented for a benchmark low-Reynolds-number jet case.

The prediction of hot streak migration in a high-pressure turbine

Accurate predictions of combustor hot streak migration enable the turbine designer to identify high-temperature regions that can limit component life. It is therefore important that these predictions are achieved within the short time scales of a design process. This article compares temperature measurements of a circular hot streak through a turning duct and a research turbine with predictions using a three-dimensional Reynolds-averaged Navier-Stokes solver. It was found that the mixing length turbulence model did not predict the hot streak dissipation accurately. However, implementation of a very simple model of the free stream turbulence (FST) significantly improved the exit temperature predictions on both the duct and research turbine. One advantage of the simple FST model described over more complex alternatives is that no additional equations are solved. This makes the method attractive for design purposes, as it is not associated with any increase in computational time.

Velocity distributions for low pressure turbines

A parametric set of velocity distributions has been investigated using a flat plate experiment. Three different diffusion factors and peak velocity locations were tested. These were designed to mimic the suction surfaces of Low Pressure (LP) turbine blades. Unsteady wakes, inherent in real turbomachinery flows, were generated using a moving bar mechanism. A turbulence grid generated a freestream turbulence level that is believed to be typical of LP turbines. Measurements were taken across a Reynolds number range of 50,000-220,000 at three reduced frequencies (0.314, 0.628, 0.942). Boundary layer traverses were performed at the nominal trailing edge using a Laser Doppler Anemometry system and hot-films were used to examine the boundary layer behaviour along the surface. For every velocity distribution tested, the boundary layer separated in the diffusing flow downstream of the peak velocity. The loss production is dominated by the mixing in the reattachment process, mixing in the turbulent boundary layer downstream of reattachment and the effects of the unsteady interaction between the wakes and the boundary layer. A sensitive balance governs the optimal location of peak velocity on the surface. Moving the velocity peak forwards on the blade was found to be increasingly beneficial when bubblegenerated losses are high, i.e. at low Reynolds number, at low reduced frequency and at high levels of diffusion. Copyright © 2008 by ASME.

Assessment of CFD models for shock boundary-layer interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.

Three-dimensional effects on sliding and waving wings

Like large insects, micro air vehicles operate at low Reynolds numbers O(1; 000 - 10; 000) in a regime characterized by separated flow and strong vortices. The leading-edge vortex has been identified as a significant source of high lift on insect wings, but the conditions required for the formation of a stably attached leading-edge vortex are not yet known. The waving wing is designed to model the translational phase of an insect wing stroke by preserving the unsteady starting and stopping motion as well as three-dimensionality in both wing geometry (via a finite-span wing) and kinematics (via wing rotation). The current study examines the effect of the spanwise velocity gradient on the development of the leading-edge vortex along the wing as well as the effects of increasing threedimensionalityby decreasing wing aspect ratio from four to two. Dye flow visualization and particle image velocimetry reveal that the leading-edge vortices that form on a sliding or waving wing have a very high aspect ratio. The structure of the flow is largely two-dimensional on both sliding and waving wings and there is minimal interaction between the leading-edge vortices and the tip vortex. Significant spanwise flow was observed on the waving wing but not on the sliding wing. Despite the increased three-dimensionality on the aspect ratio 2 waving wing, there is no evidence of an attached leading-edge vortex and the structure of the flow is very similar to that on the higher-aspect-ratio wing and sliding wing. © Copyright 2010.

RANS-URANS in axial compressor, a design methodology

Turbomachinery flows are inherently unsteady. Until now during the design process, unsteadiness has been neglected, with resort merely to steady numerical simulations. Despite the assumption involved, the results obtained with steady simulations have been used with success. One of the questions arising in recent years is can unsteady simulations be used to improve the design of turbomachines? In this work the numerical simulation of a multi-stage axial compressor is carried out. Comparison of Reynolds averaged Navier-Stokes (RANS) and unsteady Reynolds averaged Navier-Stokes (URANS) calculation shows that the unsteadiness affects pressure losses and the prediction of stall limit. The unsteady inflow due to the wake passing mainly modifies the losses and whirl angle near the endwalls. The computational cost of the fully unsteady compared with a steady simulation is about four times in terms of mesh dimension and two orders of magnitude as number of iterations. A mixed RANS-URANS solution has been proposed to give the designer the possibility to simulate an unsteady stage embedded in a steady-state simulation. This method has been applied to the simulation of a four-stage axial compressor rig. The mixed RANS-URANS approach has been developed using sliding and mixing planes as interface conditions. The rotor-stator interaction has been captured physically while reducing the computational time and mesh size.

Effects of Lewis number on scalar variance transport in premixed flames

The influences of differential diffusion rates of heat and mass on the transport of the variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been studied using three-dimensional simplified chemistry based Direct Numerical Simulation (DNS) data of statistically planar turbulent premixed flames with global Lewis number ranging from Le = 0.34 to 1.2. The Lewis number effects on the statistical behaviours of the various terms of the transport equations of variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been analysed in the context of Reynolds Averaged Navier Stokes (RANS) simulations. It has been found that the turbulent fluxes of the progress variable and temperature variances exhibit counter-gradient transport for the flames with Lewis number significantly smaller than unity whereas the extent of this counter-gradient transport is found to decrease with increasing Lewis number. The Lewis number is also shown to have significant influences on the magnitudes of the chemical reaction and scalar dissipation rate contributions to the scalar variance transport. The modelling of the unclosed terms in the scalar variance equations for the non-unity Lewis number flames have been discussed in detail. The performances of the existing models for the unclosed terms are assessed based on a-priori analysis of DNS data. Based on the present analysis, new models for the unclosed terms of the active scalar variance transport equations are proposed, whenever necessary, which are shown to satisfactorily capture the behaviours of unclosed terms for all the flames considered in this study. © 2010 Springer Science+Business Media B.V.

Hybrid LES - RANS of complex geometry jets

Large eddy simulation (LES) type studies are made of a realistic geometry coaxial nozzle with a pylon. For the LES, since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving Numerical LES (NLES). To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier Stokes) model is used giving a hybrid NLES-RANS approach.The pylon is shown to influence the flow development, having a significant impact on peak turbulence levels and spreading rates. The results show that real geometry effects are influential and should be taken into account when moving towards real engine simulations. If their effects are ignored then, based on the studies here, key turbulence parameters will have significant error.