Tip gap height effects on flow structure and heat/mass transfer over plane tip of a high-turning turbine rotor blade

The effects of tip gap height-to-chord ratio, h/c, on the flow structure and heat/mass transfer over the plane tip surface of a large-scale high-turning turbine rotor blade have been investigated for h/c = 1.0%, 2.0%, 3.0% and 4.0%. For near-wall tip gap flow visualizations, a high-resolution oil film method is employed, and the naphthalene sublimation technique is used for local heat/mass transfer rate measurements. From the tip surface visualizations, a pair of vortices named “tip gap vortices” is identified in the leading edge region within the tip gap. The overall tip gap flow is characterized not only by the tip gap vortices but also by the flow separation/re-attachment process along the pressure-side tip edge. Within the separation bubble, there exist complicated near-wall flows moving toward a mid-chord flow converging area. With increasing h/c, the tip gap vortices, the flow separation/re-attachment, and the converging flows within the separation bubble tend to be intensified. In general, higher thermal load is found along the loci of the tip gap vortices and along the re-attachment line, while lower thermal load is observed behind the tip gap vortex system and near the mid-chord flow converging area. Heat/mass transfer characteristics with the variation of h/c are discussed in detail in conjunction with the tip gap flow features. Based on the flow visualizations and heat/mass transfer data, new realistic tip gap flow models have been proposed for h/c = 1.0 and 4.0%.     

Simulation of turbulent flow and heat transfer in channels between rod bundles

The turbulent flow and heat transfer in triangular rod bundles are investigated theoretically with CFD code FLUENT. The unsteady Reynolds Stress Model is adopted as turbulence modeling. The wall function is used for near wall boundary layer. The calculation results were in agreement with experimental data. The effects of the Reynolds number and pitch to diameter ratio on the flow and heat transfer in the lattice are significant. The traditional theoretical models could not predict the flow and heat transfer in the lattice. The P/D = 1.03 is a critical point. In this case, the flow and heat transfer in the lattice is the most desirable and most efficient, and the nuclear power could also reach its maximum. The variation of large scale coherent structure with pitch to diameter ratio is consistent with the variation of the Nusselt number with pitch to diameter ratio.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

Turbulent fuel-air mixing study of jet in crossflow at different velocity ratios using LES

In order to study the mixing mechanism of fuel and air in gas turbine, large eddy simulation has been used to investigate the methane jet-in-crossflow with the velocity ratio (R) of 1.5 and 4. This study aims to explore the formation mechanism of vortices such as the hairpin vortices, hovering vortices and horseshoe vortices, the relationship between the fuel–air mixing and flow characteristics at different velocity ratios. The numerical methods in the present work are firstly validated with the experimental data in terms of mean and root mean square values of velocity. For R = 4, the shear layer vortices, horseshoe vortices, counter-rotating vortices pairs (CVP) and wake vortices can be observed, while the jet shear layer cannot be observed for R = 1.5. The hairpin vortices originating from the vortice-ring are lifted and shed from the downstream of the jet-outlet due to Kutta-Joukowski lift. The hairpin vortices are similar to CVP. The horseshoe vortices in R = 1.5 and 4 are formed due to the blockage of the jet (CH₄) and the crossflow (air) respectively, and its evolution is associated with the hovering vortices which only exist for R = 1.5. The uniform index and pr-obability density function are used for quantitative analysis of the mixing performance. The uniform index at X/D = 0 (fuel-inlet) and at X/D = 25 (outlet) are 0.033 and 0.335 for R = 1.5 and 0.130 and 0.047 for R = 4. For R = 4, the jet penetration is higher and the deflection angle of jet is smaller than that in case of R = 1.5. Higher R will provide more region for mixing, therefore uniform index is higher and the mixing is more uniform in the downstream.     

Coherent structures in unsteady swirling jet flow

An LDA technique and phase-averaging analysis were used to study unsteady precessing flow in a model vortex burner. Detailed measurements were made for Re=15,000 and S=1.01. On the basis of the analysis of phase-averaged data and vortex detection by the λ₂-technique of Joeng and Hussain (1995), three precessing spiral vortex structures were identified: primary vortex (PV), inner secondary vortex (ISV), and outer secondary vortex (OSV). The PV is the primary and most powerful structure as it includes primary vorticity generated by the swirler; the ISV and OSV are considered here as secondary vortical structures. The jet breakdown zone is the conjunction of a pair of co-rotating co-winding spiral vortices, PV and ISV. The interesting new feature described is that the secondary vortices form a three-dimensional vortex dipole with a helical geometry. The effect of coupling of secondary vortices was suggested as a mechanism of enhanced stability reflected in their increased axial extent.     

Turbulent film condensation on a horizontal elliptical tube

An investigation has been made of turbulent film condensation on a horizontal elliptical tube. The present study is based on Colburn analogy [1] and potential flow theory to determine the high tangential velocity of vapor flow at the boundary layer and to define the local interfacial shear owing to high velocity vapor flow across the tube surface. The condensate film flow and local/or mean heat transfer characteristics from a horizontal elliptical tube with variable ellipticities, e, under the influence of Froude number, sub-cooling parameter and system pressure have been performed. The present result for dimensionless mean heat transfer coefficient reduces to the same result obtained by Sarma et al.s [2] e=0 (circular tube). Compared with laminar model by Yang and Hsu [3], the present turbulent model shows in better agreement with Michaels experimental data [4] (for e=0). The dependence of mean Nusselt coefficient on the effect of n (power of Reynolds) [1] is also discussed.     

Effects of cross ribs on heat/mass transfer in a two-pass rotating duct

Two-pass internal cooling passage with rib turbulators has been investigated for convective heat/mass transfer under rotating conditions. The flow and heat transfer characteristics in the cooling passage are very complicated so that it is required the detail analysis to design more efficient gas turbine blades. A naphthalene sublimation technique is employed to determine detailed local heat transfer coefficients using the heat and mass transfer analogy. The local heat/mass transfer and flow pattern in the cooling passage are changed significantly according to rib configurations, duct turning geometries and duct rotation speeds. Four different rib configurations are investigated to obtain the combined effects of the angled rib, duct turning and rotation. The results show that the duct rotation generates the heat/mass transfer discrepancy between the leading and trailing walls due to the secondary flows induced by the Coriolis force. The angled ribs generate a single rotating secondary flow with the cross-rib arrangement and the duct turning makes a strong Dean-type vortex. These vortices affect significantly the heat/mass transfer on the duct wall. The overall heat transfer pattern on the leading and trailing surfaces for the first and second passes are dependent on the duct rotation, but the local heat transfer trend is affected mainly by the rib arrangements. In addition, the present study observes the rotating effect in the two-pass smooth duct to obtain the baseline data in comparison with the ribbed duct for various rib arrangements.     

Computational analysis of heat transfer and pressure drop performance for internally finned tubes with three different longitudinal wavy fins

Turbulent pressure drop and heat transfer characteristics in tubes with three different kinds of internally longitudinal fin patterns (interrupted wavy, sinusoidal wavy and plain) are numerically investigated for Re = 904–4,520. The channel velocity, temperature, and turbulence fields are obtained to discern the mechanisms of heat transfer enhancement. Numerical results indicate that the steady and spatially periodic growth and disruption of cross-sectional vortices occur near the tube/fin walls along the streamwise locations. The thermal boundary layers near the tube/fin surfaces are thereby periodically interrupted, with heat transfer near the recirculation zones being enhanced. The overall heat transfer coefficients in wavy channels are higher than those in a plain fin channel, while with larger pressure drop penalties. At the same waviness, the interrupted wavy fin tube could enhance heat transfer by 72–90%, with more than 2–4 times of pressure drop penalty. Among the fins studied, the sinusoidal wavy fin has the best comprehensive performance.     

Heat transfer and flow temperature measurements in a rotating triangular channel with various rib arrangements

In the present study, the regionally-averaged heat transfer coefficients and flow temperature distributions were measured in an equilateral triangular channel with three different rib arrangements (α = 45, 90 and 135°). To measure regionally-averaged heat transfer coefficients in the channel, two rows of copper blocks and a single heater were installed on two ribbed walls. The fluid temperature distributions were obtained using a thermocouple-array. The rotation number ranged from 0.0 to 0.1 with a fixed Reynolds number of 10,000. For the 90° ribs, the heat transfer coefficients on the pressure side surface were increased significantly with rotation, while the suction side surface had lower heat transfer coefficients than the stationary channel. For the angled ribs, rib-induced secondary flow dominated the heat transfer characteristics and high heat transfer rates were observed on the regions near the inner wall for the 45° angled ribs and near the leading edge for the 135° angled ribs.     

DNS of heat transfer in turbulent and transitional channel flow obstructed by rectangular prisms

Direct numerical simulation (DNS) of heat transfer in a channel flow obstructed by rectangular prisms has been performed for Re_τ = 80–20, where Re_τ is based on the friction velocity, the channel half width and the kinematic viscosity. The molecular Prandtl number is set to be 0.71. The flow remains unsteady down to Re_τ = 40 owing to the disturbance induced by the prism. For Re_τ = 30 and 20, the flow results in a steady laminar flow. In the vicinity of the prism, the three-dimensional complex vortices are generated and heat transfer is enhanced. The Reynolds number effect on the time-averaged vortex structure and the local Nusselt number are investigated. The mechanism of the heat transfer enhancement is discussed. In addition, the mean flow parameters such as the friction factor and the Nusselt number are examined in comparison with existing DNS and experimental data.     

Experimental investigation of the coherent structures in a spirally corrugated pipe

Previous numerical and theoretical results (Chen et al., 2019; Liu et al., 2018; Zhao et al., 2019) based on the optimization theory of convective heat transfer reveal that the optimized flow structures in a straight circular pipe enhancing convective heat transfer are multiple longitudinal vortices. This conclusion encourages us to find out whether such flow structures really exist in some enhanced heat transfer pipes by means of advanced experimental techniques. Therefore, a typical enhanced heat transfer pipe was selected, namely a spirally corrugated pipe, and stereoscopic particle image velocimetry (SPIV) was employed to measure its internal instantaneous flow field. Moreover, the proper orthogonal decomposition (POD) method was used to extract the large-scale coherent structures from the measured instantaneous velocity fields. Besides the spirally corrugated pipe, the fully developed turbulent flow in a straight pipe was also analyzed as benchmark of the enhanced heat transfer pipes. The results reveal that longitudinal whirling flow with multi-vortices is formed in both the fully developed turbulent flow field of the straight pipe and the spirally corrugated one. It is thus easy to explain the heat transfer enhancement mechanism of the above flow structures from the perspective of momentum transfer. The flow structures of the fully developed turbulent flow in a straight pipe are quite similar to the optimal flow pattern from the optimization theory. More specifically, multiple longitudinal vortices are spontaneously generated due to turbulence without external heat transfer enhancement techniques. Furthermore, the flow structures similar to multiple longitudinal vortices also exist in the spirally corrugated pipe, although these flow structures deviate from symmetric multiple vortices. Moreover, the flow structures in the spirally corrugated pipe are much more energetic than those in the fully developed turbulent flow in a straight pipe. This is probably the reason why a spirally corrugated pipe can enhance heat transfer compared with a straight circular pipe.     

Heat transfer in film-cooled turbulent boundary layers at different blowing ratios as affected by longitudinal vortices

Effects of embedded longitudinal vortices on heat transfer in film-cooled turbulent boundary layers at different blowing ratios are discussed. These results were obtained in boundary layers at free-stream velocities of 10 and 15 m/s. Film coolant was injected from a single row of holes at blowing ratios of 0.47–1.26. A single longitudinal vortex was induced upstream of the film-cooling holes using a half-delta wing attached to the wind tunnel floor. Heat transfer measurements were made downstream of injection using a constant heat flux surface with 126 thermocouples for surface temperature measurements. For all blowing ratios examined, the embedded vortices cause significant alterations to wall heat transfer and to film cooling distributions. Measurrments of mean temperatures and mean velocities in spanwise planes show that high wall heat transfer regions are associated with regions of high near-wall longitudinal velocity where very little film coolant is present. In addition to high heat transfer regions associated with the vortex downwash, there are also secondary heat transfer peaks. These secondary peaks develop due to shear layer mixing and interaction between the vortex and cooling jets and become higher in magnitude and more persistent with downstream distance as the blowing ratio increases from 0.47 to 1.26.     

Numerical simulation of fluid flow and heat transfer characteristics in channel with V corrugated plates

A detailed numerical study is carried out to investigate fluid flow and heat transfer characteristics in a channel with heated V corrugated upper and lower plates. The parameters studied include the Reynolds number (Re = 2,000–5,500), angles of V corrugated plates (θ = 20°, 40°, 60°), and constant heat fluxs (q″ = 580, 830, 1,090 W/m²). Numerical results have been validated using the experimented data reported by Naphon, and a good agreement has been found. The angles of V corrugated plates (θ) and the Reynolds number are demonstrated to significantly affect the fluid flow and the heat transfer rate. Increasing the angles of V corrugated plates can make the heat transfer performance become better. The increasing Reynolds number leads to a more complex fluid flow and heat transfer rate. The numerical calculations with a non-equilibrium wall function have a better accuracy than with a standard wall function for solving high Reynolds numbers or complex flow problems.     

A rotating elliptic airfoil in fluid at rest and in a parallel freestream

This paper reports results of DPIV measurements on a two-dimensional elliptic airfoil rotating about its own axis of symmetry in a fluid at rest and in a parallel freestream. In the former case, we examined three rotating speeds (Re _c,Ω = 400, 1,000 and 2,000), and in the later case, four rotating speeds (Ro _c,Ω = 2.4, 1.2, 0.6 and 0.4), together with two freestream velocities (Re _c,u  = 200 and 1,000) and two starting configurations of the airfoil (i.e., chord parallel to (α ₀ = 0°) or normal (α ₀ = 90°) to the freestream). Results show that a rotating airfoil in a stationary fluid produces two distinct types of vortex structures depending on the Reynolds number. The first type occurs at the lowest Reynolds number (Re _c,Ω = 400), where vortices shed from the two edges or tips of the airfoil dissipated quickly, resulting in the airfoil rotating in a layer of diffused vorticity. The second type occurs at higher Reynolds numbers (i.e., Re _c,Ω = 1,000 and 2,000), where the corresponding vortices rotated together with the airfoil. Due to the vortex suction effect, the torque characteristics are likely to be heavily damped for the first type because of the rapidly subsiding vortex shedding, and more oscillatory for the second type due to persistent presence of tip vortices. In a parallel freestream, increasing the tip-speed ratio (V/U) of the airfoil (i.e., decreasing the Rossby number, Ro _c,Ω) transformed the flow topology from periodic vortex shedding at Ro _c,Ω = 2.4 to the generation of a “hovering vortex” at Ro _c,Ω = 0.6 and 0.4. The presence of the hovering vortex, which has not been reported in literature before, is likely to enhance the lift characteristics of the airfoil. Freestream Reynolds number is found to have minimal effect on the vortex formation and shedding process, although it enhances shear layer instability and produces more small-scale flow structures that affect the dynamics of the hovering vortex. Likewise, initial starting configuration of the airfoil, while affecting the flow transient during the initial phase of rotation, has insignificant effect on the overall flow topology. Unfortunately, technical constraint of our apparatus prevented us from carrying out complimentary force measurements; nevertheless, the results presented herein, which are more extensive than those computed by Lugt and Ohring (1977), will provide useful benchmark data, from which more advanced numerical calculations can be carried out to ascertain the corresponding force characteristics, particularly for those conditions with the presence of hovering vortex.     

Synchronized analysis of an unsteady laminar flow downstream of a circular cylinder centred between two parallel walls using PIV and mass transfer probes

Experiments are carried out in the wake of a cylinder of d _c  = 10 mm diameter placed symmetrically between two parallel walls with a blockage ratio r = 1/3 and a Reynolds number varying between 75 ≤ Re ≤ 277. Particle image velocimetry is exerted to obtain the instantaneous velocity components in the cylinder wake. A snapshot proper orthogonal decomposition (POD) is also applied to these PIV results in order to extract the dominant modes through the implementation of an inhomogeneous filtering of these different snapshots, apart from an interpolation to estimate the wall shear rate at the lower wall downstream the cylinder. Mass transfer circular probes are placed at the lower wall downstream this obstacle so as to further determine the time evolution of the wall shear rate, by bringing the inverse method to bear on the convective-diffusion equation. Comparisons between the two synchronized techniques demonstrate that electrochemical method can give more accurate information about the coherent structures present in the flow and about the interaction of the von Kármán vortices with the walls of the tunnel as well. The comparison between the two measurement techniques in the flow regions concerns the spatiotemporal evolutions of the wall shear rate obtained from PIV measurements and the wall shear rate using mass transfer probes. Discrepancy between the PIV measurements and the electrochemical ones near the wall, where the secondary vortices P ₁′ are generated at wall, are caused by a PIV bias and a limitations of the singular mass transfer probes.     

Heat/Mass transfer in a two-pass rotating rectangular duct with and without 70°-angled ribs

The present study investigates convective heat/mass transfer and flow characteristics inside a cooling passage of rotating gas-turbine blades. The rotating duct with and without rib turbulators are used. The ribs of 70° attack angle are attached on leading and trailing surfaces in a staggered arrangement. A naphthalene sublimation technique is employed to determine detailed local heat transfer coefficients using the heat and mass transfer analogy. Additional numerical calculations are conducted to analyze the flow patterns in the cooling. The local heat/mass transfer and the flow pattern in the passage are changed significantly according to rib configurations, duct turning geometry and duct rotation speed. The results show that the duct rotation generates the heat transfer discrepancy between the leading and trailing walls due to the secondary flows induced by the Coriolis force. The heat/mass transfer on the ribbed duct shows 80% higher than the smooth duct because the ribs attached on the walls disturb the mainflow resulting in recirculation and secondary flows near the rib with the secondary flow generated by rotation. The overall heat transfer pattern on the leading and trailing walls for the first and second passes depend on the rotating speed and the turning geometry, but the local heat transfer trend is affected mainly by the rib arrangeements.     

Heat transfer and friction factor of turbulent flow through a horizontal semi-circular duct

Forced convection heat transfer and pressure drop characteristics of air flow inside a horizontal semi-circular duct are investigated experimentally. The experiments are carried out on a semi-circular duct of 23 mm inner radius, 2 mm thickness, and 2,000 mm length within a range of Reynolds number (8,242 ≤ Re ≤ 57,794)., under uniform wall heat flux conditions. The friction factor is determined by measuring the axial static pressure at different selected axial stations along the semi-circular duct. The variations of surface and mean air temperatures, local heat transfer coefficient, local Nusselt number, and the friction factor with the axial dimensionless distance are presented. It is observed that, for a given value of Reynolds number, each of the local heat transfer coefficient and the friction factor has a relatively high value near the entrance of the semi-circular duct then it decreases with increasing the dimensionless axial distance until it approaches a nearly constant value at the fully developed region. Also, it is found that, with increasing the Reynolds number the average heat transfer coefficient is increased and the friction factor is decreased. Moreover, empirical correlations for the heat transfer coefficient and friction factor as a function of the Reynolds number are obtained.     

Computation of heat transfer enhancement in a plate‐fin heat exchanger with triangular inserts and delta wing vortex generator

Longitudinal vortices disrupt the growth of the thermal boundary layer, thereby the vortex generators producing the longitudinal vortices are well known for the enhancement of heat transfer in compact heat exchangers. The present investigation determines the heat transfer characteristics with secondary flow analysis in plate fin triangular ducts with delta wing vortex generators. This geometrical configuration is investigated for various angles of attack of the wing i.e. 15°, 20°, 26° and 37° and Reynolds numbers 100 and 200. The constant wall temperature boundary condition is used. The solution of the complete Navier Stokes equation and the energy equation is carried out using the staggered grid arrangement. The performance of the combination of triangular secondary fins and delta wing with stamping on slant surfaces has also been studied. Copyright © 2009 John Wiley & Sons, Ltd.     

Local heat transfer around a wall-mounted cube at 45° to flow in a turbulent boundary layer

The flow and local heat transfer around a wall-mounted cube oriented 45° to the flow is investigated experimentally in the range of Reynolds number 4.2 × 10³–3.3 × 10⁴ based on the cube height. The distribution of local heat transfer on the cube and its base wall are examined, and it is clarified that the heat transfer distribution under the angled condition differs markedly to that for cube oriented perpendicular to the flow, particularly on the top face of the cube. The surface pressure distribution is also investigated, revealing a well-formed pair of leading-edge vortices extending from the front corner of the top face downstream along both front edges for Re>(1−2)×10⁴. Regions of high heat transfer and low pressure are formed along the flow reattachment and separation lines caused by these vortices. In particular, near the front corner of the top face, pressure suction and heat transfer enhancement are pronounced. The average heat transfer on the top face is enhanced at Re>(1−2)×10⁴ over that of a cube aligned perpendicular to the flow.