DIRECT NUMERICAL SIMULATION OF TURBINE CASCADE FLOW WITH HEAT TRANSFER EFFECTS

低雷诺数流动对高空动力装置, 特别是涡轮部件的性能产生重要的影响. 本文采用具有7阶精度的差分格式, 通过直接求解二维瞬态可压缩Navier-Stokes方程组, 对雷诺数为241 800 (基于叶片弦长)时的叶片表面带有热传导效应的平面涡轮叶栅流动进行了二维直接数值模拟, 对低雷诺数平面涡轮叶栅流动的非定常流动现象作了初步的探索.数值结果表明:在叶栅通道入口处, 流场的非定常性很弱;在叶栅尾缘处, 具有正负涡量的尾涡交替地从压力面和吸力面上脱落;周期性的涡脱落使得叶栅通道内和尾迹区的总压发生(准)周期的变化, 并且, 尾迹区总压变化主频率是通道内总压变化主频率的2倍;在时均流场中, 叶片表面压力的分布与实验值吻合良好, 表征热传导效应的斯坦顿数除湍流区外与实验值基本吻合;尾迹区速度脉动的2阶统计量与圆柱绕流尾迹区速度脉动2阶统计量具有基本相似的分布特征.     

Turbulence measurements in the inlet plane of a centrifugal compressor vaneless diffuser

Detailed flow measurements at the inlet of a centrifugal compressor vaneless diffuser are presented. The mean 3-d velocities and six Reynolds stress components tensor are used to determine the turbulence production terms which lead to total pressure loss. High levels of turbulence kinetic energy were observed in both the blade and passage wakes, but these were only associated with high Reynolds stresses in the blade wakes. For this reason the blade wakes mixed out rapidly, whereas the passage wake maintained its size, but was redistributed across the full length of the shroud wall. Peak levels of Reynolds stress occurred in regions of high velocity shear and streamline curvature which would tend to destabilize the shear gradient. Four regions in the flow are identified as potential sources of loss - the blade wake, the shear layers between passage wake and jet, the thickened hub boundary layer and the interaction region between the secondary flow within the blade wake and the passage vortex. The blade wakes generate most turbulence, with smaller contributions from the hub boundary layer and secondary flows, but no significant contribution is apparent from the passage wake shear layers.     

Experimental and numerical investigations of tip injection on tip clearance flow in an axial turbine cascade

This investigation presents detailed experimental measurements of an active tip-clearance control method based on tip injection in a high-turning axial turbine cascade. Besides that, numerical investigations are also conducted to study phenomena which are not easily measured in the experiments. It aims to study the influence of tip injection on tip clearance flow, with emphasis on the effects of injection locations. Detailed flow field measurements were made downstream of the cascade using a three-hole probe. Static pressure distributions were also measured on the blade surface at 50% and 97.5% span, respectively. The results suggest that tip injection can weaken tip clearance flow, reducing the tip clearance mass flow and its associated losses. Meanwhile, the heat transfer condition on the blade tip surface can be also improved significantly. It also can be found that injection chordwise location plays an important role in the redistribution of secondary flow within the cascade passage. When the same number of injection holes and injection mass flow are applied, holes located in the aft part of blade can perform much better than that in the front part.     

Two-equation turbulence models for prediction of heat transfer on a transonic turbine blade

Two versions of the two-equation k–ω model and a shear stress transport (SST) model are used in a three-dimensional, multi-block, Navier–Stokes code to compare the detailed heat transfer measurements on a transonic turbine blade. It is found that the SST model resolves the passage vortex better on the suction side of the blade, thus yielding a better comparison with the experimental data than either of the k–ω models. However, the comparison is still deficient on the suction side of the blade. Use of the SST model does require the computation of distance from a wall, which for a multi-block grid, such as in the present case, can be complicated. However, a relatively easy fix for this problem was devised. Also addressed are issues such as (1) computation of the production term in the turbulence equations for aerodynamic applications, and (2) the relation between the computational and experimental values for the turbulence length scale, and its influence on the passage vortex on the suction side of the turbine blade.     

Experimental study on the effect of unsteadiness on boundary layer development on a linear turbine cascade

 The results from an experimental investigation of unsteady boundary layer behavior on a linear turbine cascade are presented in this paper. To perform a detailed study on unsteady cascade aerodynamics and heat transfer, a new large-scale, high-subsonic research facility for simulating the periodic unsteady flow has been developed. It is capable of sequentially generating up to four different unsteady inlet flow conditions that lead to four different passing frequencies, wake structures, and freestream turbulence intensities. For a given Reynolds number, two different unsteady wake formations are utilized. Detailed unsteady boundary layer velocity. turbulence intensity, and pressure measurements are performed along the suction and pressure surfaces of one blade. The results display the transition and development of the boundary layer, ensemble-averaged velocity, and turbulence intensity. Received: 23 September 1996/Accepted: 19 February 1997     

Heat transfer coefficient measurements on the pressure surface of a transonic airfoil

This paper presents steady-state recovery temperature and heat transfer coefficient measurements on the pressure surface of a modern, highly cambered transonic airfoil. These measurements were collected with a peak Mach number of 1.5 and a maximum turbulence intensity of 30%. We used a single passage model to simulate the idealized two-dimensional flow path between rotor blades in a modern transonic turbine. This set up offered a simpler construction than a linear cascade, yet produced an equivalent flow condition. We performed validated high accuracy (±0.2°C) surface temperature measurements using wide-band thermochromic liquid crystals allowing separate measurements of the previously listed parameters with the same heat transfer surface. We achieved maximum heat transfer coefficient uncertainties that were equivalent to similar investigations (±10%). Two key observations are the heat transfer coefficient along the aft portion of the airfoil is sensitive to the surface heat flux and is highly insensitive to the level of freestream turbulence. Possible explanations for these observations are discussed.     

平面叶栅气动试验研究进展与展望

平面叶栅气动试验传统上是验证压气机、涡轮的基元性能的主要手段, 近年来国内外研究人员利用平面叶栅开展了大量的流动测量试验, 以揭示叶栅内部复杂流动现象的本质和规律、探索减小叶栅内流动损失的方法. 本文从试验装置、测试技术和研究内容三个方面, 综述了近年来平面叶栅气动试验研究的进展情况. 首先介绍了平面叶栅试验装置的发展及提高平面叶栅试验段流场品质的措施; 其次介绍了叶栅气动试验采用的部分流场测试技术, 包括叶片表面压力场、叶片表面温度场、内流速度场及流场可视化等测试技术, 分析了这些测试技术的进展和存在的问题; 然后梳理了近年来平面叶栅试验研究的相关科学问题及进展, 包括跨音速叶栅中的激波研究, 叶顶间隙泄漏流动研究, 叶型优化研究, 多尺度非定常旋涡结构研究, 振动环境下叶栅流场研究等; 最后对平面叶栅气动试验研究方向进行了展望. 通过了解叶栅内复杂流动现象及本质, 为进一步探索和提高压气机、涡轮的气动性能提供技术支撑.      

Measurements of the tip leakage vortex structures and turbulence in the meridional plane of an axial water-jet pump

Particle image velocimetry (PIV) measurements at varying resolutions focus on the flow structures in the tip region of a water-jet pump rotor, including the tip-clearance flow and the rollup process of a tip leakage vortex (TLV). Unobstructed views of these regions are facilitated by matching the optical refractive index of the transparent pump with that of the fluid. High-magnification data reveal the flow non-uniformities and associated turbulence within the tip gap. Instantaneous data and statistics of spatial distributions and strength of vortices in the rotor passage reveal that the leakage flow emerges as a wall jet with a shear layer containing a train of vortex filaments extending from the tip of the blade. These vortices are entrained into the TLV, but do not have time to merge. TLV breakdown in the aft part of the blade passage further fragments these structures, increasing their number and reducing their size. Analogy is made between the circumferential development of the TLV in the blade passage and that of the starting jet vortex ring rollup. Subject to several assumptions, these flows display similar trends, including conditions for TLV separation from the shear layer feeding vorticity into it.     

Large chord turbine cascade testing at engine Mach and Reynolds numbers

A technique is presented for producing a flow through a linear cascade of turbine blades of large chord which gives the pressure distribution around a blade the same as that obtained in an infinite cascade for Mach and Reynolds numbers typical of gas turbine operating conditions. Results of experiments with a cascade of three blades of large chord are compared with results from a cascade of nine blades of smaller chord to confirm the validity of the technique. Experiments are performed on the large-chord cascade to examine surface phenomena with high spatial resolution. Boundary layer scales are also increased and profiles on both suction and pressure surfaces of the blade are obtained.     

Turbine vane film cooling: Heat transfer evaluation using high-order discontinuous Galerkin RANS computations

A high-order accurate CFD solver, based on the Discontinuous Galerkin (DG) finite element method, is here employed to compute the heat transfer, with and without film coolant injection, around a turbine vane extensively tested in a wind tunnel. The numerical solution makes also use of a high-order polynomial representation of the airfoil curved boundary in order to minimize the numerical sources of error, leaving possibly only those related to the physical model adopted. The objective of the work is therefore twofold: on the one hand to provide a detailed investigation, often beyond the reach of the experiments, of the complex flow field arising in a film-cooled gas turbine cascade, on the other hand to ascertain the limits of the Reynolds-averaged Navier-Stokes (RANS) approach and its associated turbulence model when using high-order accurate methods. The DG formulation is briefly reviewed, as well as the experimental apparatus and the measuring technique, and then the code is applied to the computation of various test cases characterized by different reference Reynolds and Mach numbers. Two-dimensional results (up to seventh-order accurate) obtained both with the high- and low-Reynolds version of the k-ω model employed are presented. Reasonably good agreement between experimental and numerical results is obtained, even though the outcomes are far from being completely satisfactory especially for flow regimes in the low Reynolds number range. This is due to the lack of suitable modeling of the laminar-turbulent transition process taking place around the blade leading edge. Such a complex phenomenon is out of reach of the modeling capabilities of the high-Re k-ω model, while can be roughly mimicked by the low-Re version of the model, which is able to provide a delayed onset of the turbulence quantities along the blade surface.Third-order accurate computation of the three-dimensional turbine vane are also presented in this work and compared with available measurements to investigate the relevant fluid flow phenomena occurring and to discuss significant issues related to an accurate prediction of the turbine wall heat transfer.     

Calibration of a four-hole pyramid probe and area traverse measurements in a short-duration transonic turbine cascade tunnel

A four-hole pyramid probe has been calibrated for use in a short-duration transonic turbine cascade tunnel. The probe is used to create area traverse maps of total and static pressure, and pitch and yaw angles of the flow downstream of a transonic annular cascade. This data is unusual in that it was acquired in a short-duration (5 s of run time) annular cascade blowdown tunnel. A four-hole pyramid probe was used which has a 2.5 mm section head, and has the side faces inclined at 60° to the flow to improve transonic performance. The probe was calibrated in an ejector driven, perforated wall transonic tunnel over the Mach number range 0.5–1.2, with pitch angles from -20° to + 20° and yaw angles from-23° to +23°. A computer driven automatic traversing mechanism and data collection system was used to acquire a large probe calibration matrix (～ 10,000 readings) of non dimensional pitch, yaw, Mach number, and total pressure calibration coefficients. A novel method was used to transform the probe calibration matrix of the raw coefficients into a probe application matrix of the physical flow variables (pitch, yaw, Mach number etc.). The probe application matrix is then used as a fast look-up table to process probe results. With negligible loss of accuracy, this method is faster by two orders of magnitude than the alternative of global interpolation on the raw probe calibration matrix. The blowdown tunnel (mean nozzle guide vane blade ring diameter 1.1 m) creates engine representative Reynolds numbers, transonic Mach numbers and high levels (≈ 13%) of inlet turbulence intensity. Contours of experimental measurements at three different engine relevant conditions and two axial positions have been obtained. An analysis of the data is presented which includes a necessary correction for the finite velocity of the probe. Such a correction is non trivial for the case of fast moving probes in compressible flow.     

An investigation of the heat transfer and static pressure on the over-tip casing wall of an axial turbine operating at engine representative flow conditions. (II). Time-resolved results

This article reports the measurements of time-resolved heat transfer rate and time-resolved static pressure that have been made on the over-tip casing of a transonic axial-flow turbine operating at flow conditions that are representative of those found in modern gas turbine engines. This data is discussed and analysed in the context of explaining the physical mechanisms that influence the casing heat flux. The physical size of the measurement domain was one nozzle guide vane-pitch and from −20% to +80% rotor axial chord. Additionally, measurements of the time-resolved adiabatic wall temperature are presented. The time-mean data from the same set of experiments is presented and discussed in Part I of this article. The nozzle guide vane exit flow conditions in these experiments were a Mach number of 0.93 and a Reynolds number of 2.7 × 10⁶ based on nozzle guide vane mid-height axial chord. The data reveal large temporal variations in heat transfer characteristics to the casing wall that are associated with blade-tip passing events and in particular the blade over-tip leakage flow. The highest instantaneous heat flux to the casing wall occurs within the blade-tip gap, and this has been found to be caused by a combination of increasing flow temperature and heat transfer coefficient. The time-resolved static pressure measurements have enabled a detailed understanding of the tip-leakage aerodynamics to be established, and the physical mechanisms influencing the casing heat load have been determined. In particular, this has focused on the role of the unsteady blade lift distribution that is produced by upstream vane effects. This has been seen to modulate the tip-leakage flow and cause subsequent variations in casing heat flux. The novel experimental techniques employed in these experiments have allowed the measurement of the time-resolved adiabatic wall temperature on the casing wall. These data clearly show the falling flow temperatures as work is extracted from the gas by the turbine. Additionally, these temperature measurements have revealed that the absolute stagnation temperature within the tip-gap flow can be above the turbine inlet total temperature, and indicates the presence of a work process that leads to high adiabatic wall temperatures as a blade-tip passes a point on the casing wall. It is shown that this phenomena can be explained by consideration of the flow vectors within the tip-gap, and that these in turn are related to the local blade loading distribution. The paper also assesses the relative importance of different time-varying phenomena to the casing heat load distribution. This analysis has indicated that up to half of the casing heat load is associated with the over-tip leakage flow. Finally, the implications of the experimental findings are discussed in relation to future shroudless turbine design, and in particular the importance of accounting for the high heat fluxes found within the tip-gap.     

Calculation of secondary flow in cascades including effects of Bernoulli surface distortion

A method is proposed for calculating the secondary flow in a cascade, including effects due to the distortion of the Bernoulli surfaces. An approximation to the secondary flow is first used to determine the displacement of the streamlines and Bernoulli surfaces in the blade passage. From this, improved values for the secondary vorticity distribution can be obtained, and the secondary flow recalculated. An example is presented in which the Bernoulli surfaces are found to be considerably distorted, and a vortex core is shown to form near the corner between the blade suction surface and the end wall     

Three-dimensional flow structures and associated turbulence in the tip region of a waterjet pump rotor blade

Stereo particle image velocimetry measurements focus on the flow structure and turbulence within the tip leakage vortex (TLV) of an axial waterjet pump rotor. Unobstructed optical access to the sample area is achieved by matching the optical refractive index of the transparent pump with that of the fluid. Data obtained in closely spaced planes enable us to reconstruct the 3D TLV structure, including all components of the mean vorticity and strain-rate tensor along with the Reynolds stresses and associated turbulence production rates. The flow in the tip region is highly three-dimensional, and the characteristics of the TLV and leakage flow vary significantly along the blade tip chordwise direction. The TLV starts to roll up along the suction side tip corner of the blade, and it propagates within the passage toward the pressure side of the neighboring blade. A shear layer with increasing length connects the TLV to the blade tip and initially feeds vorticity into it. During initial rollup, the TLV involves entrainment of a few vortex filaments with predominantly circumferential vorticity from the blade tip. Being shed from the blade, these filaments also have high circumferential velocity and appear as swirling jets. The circumferential velocity in the TLV core is also substantially higher than that in the surrounding passage flow, but the velocity peak does not coincide with the point of maximum vorticity. When entrainment of filaments stops in the aft part of the passage, newly forming filaments wrap around the core in helical trajectories. In ensemble-averaged data, these filaments generate a vortical region that surrounds the TLV with vorticity that is perpendicular to that in the vortex core. Turbulence within the TLV is highly anisotropic and spatially non-uniform. Trends can be traced to high turbulent kinetic energy and turbulent shear stresses, e.g., in the shear layer containing the vortex filaments and the contraction region situated along the line where the leakage backflow meets the throughflow, causing separation of the boundary layer at the pump casing. Upon exposure to adverse pressure gradients in the aft part of the passage, at 0.65–0.7 chord fraction in the present conditions, the TLV bursts into a broad turbulent array of widely distributed vortex filaments.     

Structure of Tip Leakage Flow in a Forward-Swept Axial-Flow Fan

An experimental analysis using three-dimensional laser Dopplervelocimetery (LDV) measurements and computational analysis usingthe Reynolds stress model of the commercial code, FLUENT, wereconducted to give a clear understanding on the structure of thetip leakage flow in a forward-swept axial-flow fan operating atthe peak efficiency condition, and to emphasize the necessity ofusing an anisotropic turbulence model for the accurate predictionof the tip leakage vortex. The rolling-up of the tip leakage flowwas initiated near the position of the maximum static pressuredifference, which was located at approximately 12% axial tipchord downstream from the leading edge of the blade, and developedalong the centerline of the pressure trough on the casing. Areverse flow between the blade tip and the casing due to the tipleakage vortex acted as a blockage on the through-flow. As aresult, high momentum flux was observed below the tip leakagevortex. As the tip leakage vortex proceeded to the aft part of theblade passage, the strength of the tip leakage vortex decreaseddue to the strong interaction with the through-flow and the casingboundary layer, and the diffusion of the tip leakage vortex byhigh turbulence. Through the comparative study of turbulencemodels, it was clearly shown that an anisotropic turbulence model,e.g., Reynolds stress model, should be used to predict reasonablyan anisotropic nature of the turbulent flow fields inside the tipleakage vortex. In comparison with LDV measurement data, thecomputed results predicted the complex viscous flow patternsinside the tip region in a reliable level.     

High resolution PIV measurements around a model turbine blade trailing edge film-cooling breakout

Film cooling downstream of a model turbine blade trailing edge has been studied experimentally. High resolution particle image velocimetry was used to obtain spatially resolved mean velocity and turbulence measurements in the immediate vicinity of the trailing edge breakout. The mean velocity measurements imply the presence of a pair of counter-rotating longitudinal vortices shed from the sides of the breakout lands. The turbulent shear stress measurements above the breakout are significantly intensified as blowing ratio is increased. These results suggest that there is a strong mixing between the film cooling slot jets and the mainstream flow which degrades the film cooling effectiveness.     

Recognition of coherent structures in the boundary layer of a low-pressure-turbine blade for different free-stream turbulence intensity levels

Particle Image Velocimetry (PIV) has been adopted to analyze the instantaneous flow field developing on a high-lift turbine blade profile operating under low and elevated free-stream turbulence conditions (FSTI). Results reported in the paper allow us to analyze the dynamics leading to transition and separation of the suction side boundary layer, looking to generation, propagation and breakdown of coherent structures observed in the two different FSTI cases. To this end, measurements have been performed in two orthogonal planes. Results obtained in the blade-to-blade plane allow the detailed characterization of the propagation of Kelvin–Helmholtz (KH) rolls generating, at low FSTI condition, as a consequence of a non-reattaching separation. Otherwise, data in the wall-parallel plane allow recognizing the presence of three-dimensional disuniformities induced at high FSTI by low and high speed streaks (Klebanoff mode). The sinuous breakdown of boundary layer streaks generates other complex three-dimensional coherent structures such as hairpin or cane-like vortices that induce transition. Proper Orthogonal Decomposition (POD) has been adopted to in depth characterize these structures, thus further explaining the mechanisms through which the free-stream turbulence intensity modify the transition/separation processes of the suction side boundary layer of an highly loaded low pressure turbine blade.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.