Investigation of boundary layer development in a multistage LP turbine

This paper describes an investigation of the behavior of suction surface boundary layers in a modern multistage Low-Pressure turbine. An array of 18 surface-mounted hot-film anemometers was mounted on a stator blade of the third stage of a four-stage machine. Data were obtained at Reynolds numbers between 0.9 × 105 and 1.8 × 105. At the majority of the test conditions, wakes from upstream rotors periodically initiated transition at about 40 percent surface length. In between these events, laminar separation occurred at about 75 percent surface length. Because the wake-affected part of the flow appeared to be only intermittently turbulent, laminar separation also occurred at about 75 percent surface length while this flow was instantaneously laminar. At all but the lowest Reynolds numbers, the time-mean boundary layer appeared to have re-attached by the trailing edge even though it was not fully turbulent. It is inferred that the effect of the wakes on the performance of the blade row is limited and that steady flow design methods should provide an adequate assessment of LP turbine performance during design.

A comparison of different failure assessment methodologies applied to the NESC-1 test

Predictions for a 75x205mm surface semi-elliptic defect in the NESC-1 spinning cylinder test have been made using BS PD 6493:1991, the R6 procedure, non-linear cracked body finite element analysis techniques and the local approach to fracture. All the techniques agree in predicting ductile tearing near the inner surface of the cylinder followed by cleavage initiation. However they differ in the amount of ductile tearing, and the exact location and time of any cleavage event. The amount of ductile tearing decreases with increasing sophistication in the analysis, due to the drop in peak crack driving force and more explicit consideration of constraint effects. The local approach predicts a high probability of cleavage in both HAZ and base material after 190s, while the other predictions suggest that cleavage is unlikely in the HAZ due to constraint loss, but likely in the underlying base material. The timing of this event varies from ∼150s for R6 predictions to ∼250-300s using non-linear cracked body analysis.

An experimental and computational study of the formation of a streamwise shed vortex in a turbine stage

Shear layers shed by aircraft wings roll up into vortices. A similar, though far less common, phenomenon can occur in the wake of a turbomachine blade. This paper presents experimental data from a new single stage turbine that has been commissioned at the Whittle Laboratory. Two low aspect ratio stators have been tested with the same rotor row. Surface flow visualisation illustrates the extremely strong secondary flows present in both NGV designs. These secondary flows lead to conventional passage vortices but also to an intense vortex sheet which is shed from the trailing edge of the blades. Pneumatic probe traverse show how this sheet rolls up into a concentrated vortex in the second stator design, but not in the first. A simple numerical experiment is used to model the shear layer instability and the effects of trailing edge shape and exit yaw angle distribution are investigated. It is found that the latter has a strong influence on shear layer rollup: inhibiting the formation of a vortex downstream of NGV 1 but encouraging it behind NGV 2.

The development of turbine exit flow in a swan-necked inter-stage diffuser

This paper describes both the migration and dissipation of flow phenomena downstream of a transonic high-pressure turbine stage. The geometry of the HP stage exit duct considered is a swan-necked diffuser similar to those likely to be used in future engine designs. The paper contains results both from an experimental programme in a turbine test facility and from numerical predictions. Experimental data was acquired using three fast-response aerodynamic probes capable of measuring Mach number, whirl angle, pitch angle, total pressure and static pressure. The probes were used to make time-resolved area traverses at two axial locations downstream of the rotor trailing edge. A 3D time-unsteady viscous Navier-Stokes solver was used for the numerical predictions. The unsteady exit flow from a turbine stage is formed from rotordependent phenomena (such as the rotor wake, the rotor trailing edge recompression shock, the tip-leakage flow and the hub secondary flow) and vane-rotor interaction dependant phenomena. This paper describes the time-resolved behaviour and three-dimensional migration paths of both of these phenomena as they convect downstream. It is shown that the inlet flow to a downstream vane is dominated by two corotating vortices, the first caused by the rotor tip-leakage flow and the second by the rotor hub secondary flow. At the inlet plane of the downstream vane the wake is extremely weak and the radial pressure gradient is shown to have caused the majority of the high loss wake fluid to be located between the mid-height of the passage and the casing wall. The structure of the flow indicates that between a high pressure stage and a downstream vane simple two-dimensional blade row interaction does not occur. The results presented in this paper indicate that the presence of an upstream stage is likely to significantly alter the structure of the secondary flow within a downstream vane. The paper also shows that vane-rotor interaction within the upstream stage causes a 10° circumferential variation in the inlet flow angle of the 2nd stage vane.

Secondary flows and loss caused by blade row interaction in a turbine stage

A study of the three-dimensional stator-rotor interaction in a turbine stage is presented. Experimental data reveal vortices downstream of the rotor which are stationary in the absolute frame - indicating that they are caused by the stator exit flowfield. Evidence of the rotor hub passage vortices is seen, but additional vortical structures away from the endwalls, which would not be present if the rotor were tested in isolation, are also identified. An unsteady computation of the rotor row is performed using the measured stator exit flowfield as the inlet boundary condition. The strength and location of the vortices at rotor exit are predicted. A formation mechanism is proposed whereby stator wake fluid with steep spanwise gradients of absolute total pressure is responsible for all but one of the rotor exit vortices. This mechanism is then verified computationally using a passive-scalar tracking technique. The predicted loss generation through the rotor row is then presented and a comparison made with a steady calculation where the inlet flow has been mixed out to pitchwise uniformity. The loss produced in the steady simulation, even allowing for the mixing loss at inlet, is 10% less than that produced in the unsteady simulation. This difference highlights the importance of the time-accurate calculation as a tool of the turbomachine designer.

The transition mechanism of highly-loaded LP turbine blades

A detailed experimental investigation was conducted into the interaction of a converted wake and a separation bubble on the rear suction surface of a highly loaded low-pressure (LP) turbine blade. Boundary layer measurements, made with 2D LDA, revealed a new transition mechanism resulting from this interaction. Prior to the arrival of the wake, the boundary layer profiles in the separation region are inflexional. The perturbation of the separated shear layer caused by the converting wake causes an inviscid Kelvin-Helmholtz rollup of the shear layer. This results in the breakdown of the laminar shear layer and a rapid wake-induced transition in the separated shear layer.

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.

Vortex filament simulation of trailing vortex merger

Aircraft in high-lift configuration shed multiple vortices. These generally merge to form a downstream wake consisting of two counter-rotating vortices of equal strength. The understanding of the merger of two co-rotating trailing vortices is important in evaluating the separation criteria for different aircraft to prevent wake vortex hazards during landing and take-off. There is no existing theoretical method on the basis of which such norms can be set. The present study is aimed at gaining a better understanding of the behaviour of wake vortices behind the aircraft. Two dimensional studies are carried out using the vortex blob method and compared with Bertenyi's experiment. It is shown that inviscid two dimensional effects are insufficient to explain the observations. Three dimensional studies, using the vortex filament method, are applied to the same test case. Two Lamb-Oseen profile vortices of the same dimensions and initial separation as the experiment are allowed to evolve from a straight starting condition until a converged steady flow is achieved. The results obtained show good agreement with the experimental distance to merger. Core radius and separation behaviour is qualitatively similar to experiment, with the exception of rapid increases at first. This may be partially attributable to the choice of filament-element length, and recommended further work includes a convergence study for this parameter. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

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.

遥感研究南海台风引起的冷迹特征

海表面温度（Sea Surface Temperature, SST）是海洋生态环境变化的重要因子，也是海洋观测研究的最重要的因子之一。海表温度场表征了海洋热力、动力过程和海洋与大气相互作用的综合结果，是反映海气热量、动力和水汽交换的重要参量，也是气候变化及海洋营养盐浓度和初级生产力的重要影响因子。台风过后通常会在上层海洋引起冷迹，即台风路径附近的海表温度降低。海表温度的变化是台风与海洋之间能量交换的关键，对于台风的产生，演化和强度变化有非常重要的作用。因此，研究上层海洋对台风的响应，特别是海表温度对台风的响应过程有非常重要的科学意义。本研究主要利用微波SST逐日数据资料统计研究了海表面温度对南海台风的响应过程，并分析了中尺度涡（特别是气旋性涡）对台风所引起的降温幅度和降温位置的影响作用。在此研究基础上初步尝试提出了一种新的量化台风降温效应的方法。 首先，本研究以1998-2009年总共12年内经过中国南海的92个台风为研究对象，利用逐日连续的微波TMI及TMI-AMSER海表温度数据资料分析了92个台风所引起的最大降温幅度以及冷迹位置的分布特征。统计结果表明有64个（69.6%）台风引起了海表温度明显的降低（降温≥2oC），而对最大降温出现的位置进行分析发现：有超过半数（56个）的台风引起的最大降温位于台风路径的附近或路径右侧，同时还有8个（8.7%）台风引起的最大降温明显位于路径的左侧。文中选取了冷迹相对于台风路径出现在不同位置（左、右、附近）的4个台风为例来说明台风引起的降温特征。通过量化降温中心到台风路径的距离发现台风所引起的最大降温出现的位置主要集中分布在路径左右两侧100km范围内。另外，通过对比台风来之前的海平面高度距平发现先于台风之前处在于上层海洋环境的冷涡，特别是强冷涡，在台风引起的海表面降温幅度和位置分布中可能起着重要的作用，其中冷涡的位置与台风所引起的最大降温中心出现的位置有很好的相关性。 其次，本研究提出一个新的描述降温效应的指数。台风所经过的海域往往会出现海表温度的降低，而海表温度的变化对海洋生态系统有很重要的影响。台风的降温效应因台风、背景温度场和台风前的海洋环境条件而异。目前，海表温度对台风的响应效应主要是通过三个物理量来描述，降温幅度、冷迹范围、冷迹持续时间。但是还没有一个综合的量来对台风的降温效应进行定量的描述。本研究利用遥感微波SST数据尝试初步

Shock / boundary-layer interaction control using three-dimensional bumps for transonic wings

Three-dimensional bumps have been developed and investigated, aiming at the two major objectives of shock-wave / boundary-layer interaction control, i.e. drag reduction and suppression of separation, simultaneously. An experimental investigation has been conducted for a default rounded bump in channel now at University of Cambridge and a computational study has been performed for a spanwise series of rounded bumps mounted on a transonic aerofoil at University of Stuttgart. Observed in both cases are wave drag reduction owing to A-shock structures produced by three-dimensional surface bumps and mild control effects on the boundary layer. The effects of rough surface and tall extension have been investigated as well as several geometric variations and multiple bump configurations. A double configuration of narrow rounded bumps has been found to best perform amongst the tested, considerably reducing wave drag through a well-established A-shock structure with little viscous penalty and thus achieving substantial overall drag reduction. Counter-rotating streamwise vortex pairs have been produced by some configurations as a result of local flow separation, but they have been observed to be confined in relatively narrow wake regions, expected to be beneficial in suppressing large-scale separation under off-design condition despite increase of viscous drag. On the whole a large potential of three-dimensional control with discrete rounded bumps has been demonstrated both experimentally and numerically, and experimental investigation of bumps fitted on a transonic aerofoil or wing is suggested toward practical application.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Detached-eddy simulation of transonic flow past a fan-blade section

Detached-eddy simulation of transonic flow past a thin section of a fan blade has been carried out. The inflow Mach number is 1.03, and a bow shock forms upstream of the blade. The shock (corresponding to an adjacent blade) impinges on the suction-side boundary layer which causes separation and rapid transition to turbulence. The boundary layer later re-attaches near the trailing edge. The pressure-side boundary layer transitions near the leading edge and remains attached. Mean surface pressure shows basic agreement with a steady RANS calculation; strong shock motion in the DES is the major cause of discrepancy. Surface pressure spectra are investigated, and low-frequency two-dimensional disturbances associated with the shock motion are dominant. Removing the two-dimensional component from the spectra, the pressure-side three-dimensional spectra reproduce the spectral shape given by a correlation for flat-plate boundary layer wall-pressure spectra developed by Goody. 1 The suction-side disturbances produce similar high- and intermediate-frequency scalings despite substantially different boundary layer development. Near-wake results show that disturbance kinetic energy peaks at the suction-side inflection point of the mean profile, and that the energy is concentrated at low frequencies relative to the near-trailing edge surface pressure. Copyright © 2009 by the authors.

UEGO based measurement of EGR rate and residual gas fraction

A technique is presented for measuring the exhaust gas recirculation (EGR) and residual gas fraction (RGF) using a fast UEGO based O2 measurement of the manifold or in-cylinder gases, and of the exhaust gases. The technique has some advantages over the more common CO2-based method. In the case of an RGF measurement, fuel interference must be eliminated and special fuelling arrangements are is required. It is shown how a UEGO-based measurement, though sensitive to reactive species in the exhaust (such as H 2), as a system reports EGR/ RGF rates faithfully. Preliminary tests showed that EGR and RGF measurements using the O2 approach agreed well with CO2-based measurements. © 2011 SAE International.

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.