Developing Buoyancy-Modified Turbulent Flow in Ducts Rotating in Orthogonal Mode

1995 ◽  
Vol 117 (3) ◽  
pp. 474-484 ◽  
Author(s):  
T. Bo ◽  
H. Iacovides ◽  
B. E. Launder

A numerical study of developing flow through a heated duct of square cross section rotating in orthogonal mode is reported. The two main aims are to explore the effects of rotational buoyancy on the flow development and to assess the ability of available turbulence models to predict such flows. Two test cases have been computed corresponding to values of the rotation number, Ro, of 0.12 and 0.24, which are typical of operating conditions in internal cooling passages of gas turbine blades. Computations from three turbulence models are presented: a k–ε eddy viscosity (EVM) model matched to a low-Reynolds-number one-equation EVM in the near-wall region; a low-Re k–ε EVM and a low-Re algebraic stress model (ASM). Additional computations in which the fluid density is assumed to remain constant allow the distinct contributions from buoyancy and Coriolis forces to be separated. It is thus shown that rotational buoyancy can have a substantial influence on the flow development and that, in the case of outward flow, it leads to a considerable increase of the side-averaged heat transfer coefficient. The Coriolis-induced secondary motion leads to an augmentation of the mean heat transfer coefficient on the pressure surface and a reduction on the suction side. The k–ε/one-equation EVM produces a mostly reasonable set of heat transfer predictions, but some deficiencies do emerge at the higher rotation number. In contrast, predictions with the low-Re k–ε EVM return a spectacularly unrealistic behavior while the low-Re ASM thermal predictions are in encouragingly close agreement with available measurements.

Author(s):  
Vijay K. Garg ◽  
Ali A. Ameri

A three-dimensional Navier-Stokes code has been used to compute the heat transfer coefficient on two film-cooled turbine blades, namely the VKI rotor with six rows of cooling holes including three rows on the shower head, and the C3X vane with nine rows of holes including five rows on the shower head. Predictions of heat transfer coefficient at the blade surface using three two-equation turbulence models, specifically, Coakley’s q-ω model, Chien’s k-ε model and Wilcox’s k-ω model with Menter’s modifications, have been compared with the experimental data of Camci and Arts (1990) for the VKI rotor, and of Hylton et al. (1988) for the C3X vane along with predictions using the Baldwin-Lomax (B-L) model taken from Garg and Gaugler (1995). It is found that for the cases considered here the two-equation models predict the blade heat transfer somewhat better than the B-L model except immediately downstream of the film-cooling holes on the suction surface of the VKI rotor, and over most of the suction surface of the C3X vane. However, all two-equation models require 40% more computer core than the B-L model for solution, and while the q-ω and k-ε models need 40% more computer time than the B-L model, the k-ω model requires at least 65% more time due to slower rate of convergence. It is found that the heat transfer coefficient exhibits a strong spanwise as well as streamwise variation for both blades and all turbulence models.


2014 ◽  
Vol 137 (3) ◽  
Author(s):  
Ignacio Mayo ◽  
Tony Arts ◽  
Ahmed El-Habib ◽  
Benjamin Parres

The convective heat transfer distribution in a rib-roughened rotating internal cooling channel was measured for different rotation and Reynolds numbers, representative of engine operating conditions. The test section consisted of a channel of aspect ratio equal to 0.9 with one wall equipped with eight ribs perpendicular to the main flow direction. The pitch to rib height ratio was 10 and the rib blockage was 10%. The test rig was designed to provide a uniform heat flux boundary condition over the ribbed wall, minimizing the heat transfer losses and allowing temperature measurements at significant rotation rates. Steady-state liquid crystal thermography (LCT) was employed to quantify a detailed 2D distribution of the wall temperature, allowing the determination of the convective heat transfer coefficient along the area between the sixth and eighth rib. The channel and all the required instrumentation were mounted on a large rotating disk, providing the same spatial resolution and measurement accuracy as in a stationary rig. The assembly was able to rotate both in clockwise and counterclockwise directions, so that the investigated wall was acting either as leading or trailing side, respectively. The tested Reynolds number values (based on the hydraulic diameter of the channel) were 15,000, 20,000, 30,000, and 40,000. The maximum rotation number values were ranging between 0.12 (Re = 40,000) and 0.30 (Re = 15,000). Turbulence profiles and secondary flows modified by rotation have shown their impact not only on the average value of the heat transfer coefficient but also on its distribution. On the trailing side, the heat transfer distribution flattens as the rotation number increases, while its averaged value increases due to the turbulence enhancement and secondary flows induced by the rotation. On the leading side, the secondary flows counteract the turbulence reduction and the overall heat transfer coefficient exhibits a limited decrease. In the latter case, the secondary flows are responsible for high heat transfer gradients on the investigated area.


Author(s):  
C. Bianchini ◽  
B. Facchini ◽  
L. Mangani ◽  
M. Maritano

Fan-shaped holes are widely used to provide better cooling performances than cylindrical holes over a large range of different operating conditions. Main advantages of such solution include a reduced amount of cooling air for the same performance, increased part lifetime and fewer required holes. As the overall cooling performance of such holes is strictly related to the adiabatic effectiveness and heat transfer coefficient (HTC) increase due to the coolant injection, both issues should be investigated. A numerical analysis has been conducted on a laidback fan-shaped film cooling hole onto a flat plate with the aim of investigating the increase of heat transfer. A steady-state RANS analysis was performed at two different blowing ratios (1.25 and 2.5) with imposed heat flux on the heated wall reproducing the same conditions as in the experimental tests presented in the companion paper. Despite no temperature difference was imposed between main gas and coolant flow, adiabatic effectiveness maps were extracted from tracing distribution over the plate. Performances of four different eddy viscosity turbulence models have been tested: the Two-Layer model by Rodi both in the isotropic original formulation and with an anisotropic algebraic correction based on DNS data fitting as firstly proposed by Lakheal, the k–ω SST by Menter and the ν2–f by Durbin. All calculations were conducted with a 3D unstructured pressure-based compressible solver based on the open-source OpenFOAM® CFD platform. A detailed analysis of both the predicted flow field and thermal distribution in the domain was presented. The obtained results were compared with the experimental measurements showed in the companion paper both in terms of wall heat transfer coefficient and adiabatic effectiveness.


1991 ◽  
Vol 113 (3) ◽  
pp. 346-353 ◽  
Author(s):  
M. E. Taslim ◽  
L. A. Bondi ◽  
D. M. Kercher

Turbine blade cooling is imperative in advanced aircraft engines. The extremely hot gases that operate within the turbine section require turbine blades to be cooled by a complex cooling circuit. This cooling arrangement increases engine efficiency and ensures blade materials a longer creep life. One principle aspect of the circuit involves serpentine internal cooling passes throughout the core of the blade. Roughening the inside surfaces of these cooling passages with turbulence promoters provides enhanced heat transfer rates from the surface. The purpose of this investigation was to study the effect of rotation, aspect ratio, and turbulator roughness on heat transfer in these rib-roughened passages. The investigation was performed in an orthogonally rotating setup to simulate the actual rotation of the cooling passages. Single-pass channels, roughened on two opposite walls, with turbulators positioned at 45 deg angle to the flow, in a criss-cross arrangement, were studied throughout this experiment. The ribs were arranged such that their pitch-to-height ratio remained at a constant value of 10. An aspect ratio of unity was investigated under three different rib blockage ratios (turbulator height/channel hydraulic diameter) of 0.1333, 0.25, and 0.3333. A channel with an aspect ratio of 2 was also investigated for a blockage ratio of 0.25. Air was flown radially outward over a Reynolds number range of 15,000 to 50,000. The rotation number was varied from 0 to 0.3. Stationary and rotating cases of identical geometries were compared. Results indicated that rotational effects are more pronounced in turbulated passages of high aspect and low blockage ratios for which a steady increase in heat transfer coefficient is observed on the trailing side as rotation number increases while the heat transfer coefficient on the leading side shows a steady decrease with rotation number. However, the all-smooth-wall classical pattern of heat transfer coefficient variation on the leading and trailing sides is not followed for smaller aspect ratios and high blockage ratios when the relative artificial roughness is high.


Author(s):  
Ignacio Mayo ◽  
Ahmed El-Habib ◽  
Tony Arts ◽  
Benjamin Parres

The convective heat transfer distribution in a rib-roughened rotating internal cooling channel was measured for different Rotation and Reynolds numbers, representative of engine operating conditions. The test section consisted of a channel of aspect ratio equal to 0.9 with one wall equipped with 8 ribs perpendicular to the main flow direction. The pitch to rib height ratio was 10 and the rib blockage was 10 per cent. The test rig was designed to provide a uniform heat flux boundary condition over the ribbed wall, minimizing the heat transfer losses and allowing temperature measurements at significant rotation rates. Steady-state Liquid Crystal Thermography was employed to quantify a detailed two dimensional distribution of the wall temperature, allowing the determination of the convective heat transfer coefficient along the area between the 6th and 8th rib. The channel and all the required instrumentation were mounted on a large rotating disk, providing the same spatial resolution and measurement accuracy as in a stationary rig. The assembly was able to rotate both in clockwise and counterclockwise directions, so that the investigated wall was acting either as leading or trailing side, respectively. The tested Reynolds number values (based on the hydraulic diameter of the channel) were 15000, 20000, 30000 and 40000. The maximum Rotation number values were ranging between 0.12 (Re = 40000) and 0.30 (Re = 15000). Turbulence profiles and secondary flows modified by rotation have shown their impact not only on the average value of the heat transfer coefficient but also on its distribution. On the trailing side, the heat transfer distribution flattens as the Rotation number increases, while its averaged value increases due to the turbulence enhancement and secondary flows induced by the rotation. On the leading side, the secondary flows counteract the turbulence reduction and the overall heat transfer coefficient exhibits a limited decrease. In the latter case the secondary flows are responsible for high heat transfer gradients on the investigated area.


Author(s):  
Antonio Andreini ◽  
Alberto Ceccherini ◽  
Bruno Facchini ◽  
Daniele Coutandin

Due to the higher cooling requirements of novel combustor liners a comprehensive understanding of the phenomena concerning the interaction of hot gases with different coolant flows plays a major role in the definition of a well performing liner. A numerical study of a real engine cooling scheme was performed on a test article replicating a slot injection and an effusion array with a central large dilution hole. Geometry consists of a rectangular cross-section duct with a flat plate comprised 272 holes arranged in 29 staggered rows (d = 1.65 mm, Sx/d = 7.6, Sy/d = 6, L/d = 5.5, α = 30 deg); a dilution hole (D = 18.75 mm) is located at the 14th row. A detailed experimental survey has been performed on this test article making possible to compare both predicted adiabatic effectiveness and heat transfer coefficient. The study has a twofold objective. On one hand it aims to assess the accuracy of standard industrial CFD analysis in the prediction of heat transfer on the hot side of realistic effusion cooled plates, and, on the other hand, it allows to better understand the structure of flow field, not investigated with experiments. Steady state RANS calculations have been performed on 3D computational domain with a full explicit discretization of effusion holes, with a sensitivity to standard two-equation turbulence models. Numerical results have pointed out a large dependence on effusion velocity ratio and, despite the well known deficiency of eddy viscosity models in the prediction of film effectiveness, CFD results have shown an excellent agreement with experiments in the prediction of hot side heat transfer coefficient. The entity of local heat transfer augmentation due to gas-jets interaction and its dependence on jets velocity ratio were predicted with very satisfactory agreement. The increase of heat transfer is usually located very close to jet exits and it is mainly due to local flow acceleration and vortices whose calculation is not affected by the inaccurate jet mixing prediction of first order turbulence models. Besides the comparison with experimental data of the companion paper, an additional numerical investigation was performed to assess the effect of a variable density ratio. Obtained results point out the opportunity to scale the increase in heat transfer coefficient with effusion jets velocity ratio.


Author(s):  
M. E. Taslim ◽  
L. A. Bondi ◽  
D. M. Kercher

Turbine blade cooling is imperative in advanced aircraft engines. The extremely hot gases that operate within the turbine section require turbine blades to be cooled by a complex cooling circuit. This cooling arrangement increases engine efficiency and ensures blade materials a longer creep life. One principle aspect of the circuit involves serpentine internal cooling passes throughout the core of the blade. Roughening the inside surfaces of these cooling passages with turbulence promoters provides enhanced heat transfer rates from the surface. The purpose of this investigation was to study the effects of rotation, aspect ratio, and turbulator roughness on heat transfer in these rib-roughened passage. The investigation was performed in an orthogonally rotating setup to simulate the actual rotation of the cooling passages. Single pass channels, roughened on two opposite walls, with turbulators positioned at a 45° angle to the flow, in a criss-cross arrangement, were studied throughout this experiment. The ribs were arranged such that their pitch-to-height ratio remained at a constant value of 10. An aspect ratio of unity was investigated under three different rib blockage ratios (turbulator height/channel hydraulic diameter) of 0.1333, 0.25, and 0.3333. A channel with an aspect ratio of 2 was also investigated for a blockage ratio of 0.25. Air was flown radially outward over a Reynolds number range of 15000 to 50000. Rotation number was varied from 0 to 0.3. Stationary and rotating cases of identical geometries were compared. Results indicate that rotational effects are more pronounced in turbulated passages of high aspect and low blockage ratios for which a steady increase in heat transfer coefficient is observed on the trailing side as rotation number increases while the heat transfer coefficient on the leading side shows a steady decrease with rotation number. However, the all-smooth-wall classical pattern of heat transfer coefficient variation on the leading and trailing sides is not followed for smaller aspect ratios and high blockage ratios when the relative artificial roughness is high.


Author(s):  
Bo-lun Zhang ◽  
Li Zhang ◽  
Hui-ren Zhu ◽  
Jian-sheng Wei ◽  
Zhong-yi Fu

Film cooling performance of the double-wave trench was numerically studied to improve the film cooling characteristics. Double-wave trench was formed by changing the leading edge and trailing edge of transverse trench into cosine wave. The film cooling characteristics of transverse trench and double-wave trench were numerically studied using Reynolds Averaged Navier Stokes (RANS) simulations with realizable k-ε turbulence model and enhanced wall treatment. The film cooling effectiveness and heat transfer coefficient of double-wave trench at different trench width (W = 0.8D, 1.4D, 2.1D) conditions are investigated, and the distribution of temperature field and flow field were analyzed. The results show that double-wave trench effectively improves the film cooling effectiveness and the uniformity of jet at the downstream wall of the trench. The span-wise averaged film cooling effectiveness of the double-wave trench model increases 20–63% comparing with that of the transverse trench at high blowing ratio. The anti-counter-rotating vortices which can press the film on near-wall are formed at the downstream wall of the double-wave trench. With the double-wave trench width decreasing, the film cooling effectiveness gradually reduces at the hole center-line region of the downstream trench. With the increase of the blowing ratio, the span-wise averaged heat transfer coefficient increases. The span-wise averaged heat transfer coefficient of the double-wave trench with 0.8D and 2.1D trench width is higher than that of the double-wave trench with 1.4D trench width at the high blowing ratio conditions.


Author(s):  
Godwin Ita Ekong ◽  
Christopher A. Long ◽  
Peter R. N. Childs

Compressor tip clearance for a gas turbine engine application is the radial gap between the stationary compressor casing and the rotating blades. The gap varies significantly during different operating conditions of the engine due to centrifugal forces on the rotor and differential thermal expansions in the discs and casing. The tip clearance in the axial flow compressor of modern commercial civil aero-engines is of significance in terms of both mechanical integrity and performance. In general, the clearance is of critical importance to civil airline operators and their customers alike because as the clearance between the compressor blade tips and the casing increases, the aerodynamic efficiency will decrease and therefore the specific fuel consumption and operating costs will increase. This paper reports on the development of a range of concepts and their evaluation for the reduction and control of tip clearance in H.P. compressors using an enhanced heat transfer coefficient approach. This would lead to improvement in cruise tip clearances. A test facility has been developed for the study at the University of Sussex, incorporating a rotor and an inner shaft scaled down from a Rolls-Royce Trent aero-engine to a ratio of 0.7:1 with a rotational speed of up to 10000 rpm. The idle and maximum take-off conditions in the square cycle correspond to in-cavity rotational Reynolds numbers of 3.1×106 ≤ Reφ ≤ 1.0×107. The project involved modelling of the experimental facilities, to demonstrate proof of concept. The analysis shows that increasing the thermal response of the high pressure compressor (HPC) drum of a gas turbine engine assembly will reduce the drum time constant, thereby reducing the re-slam characteristics of the drum causing a reduction in the cold build clearance (CBC), and hence the reduction in cruise clearance. A further reduction can be achieved by introducing radial inflow into the drum cavity to further increase the disc heat transfer coefficient in the cavity; hence a further reduction in disc drum time constant.


Sign in / Sign up

Export Citation Format

Share Document