Effects of Inlet Temperature Uniformity and Nonuniformity on the Tip Leakage Flow and Rotor Blade Tip and Casing Heat Transfer Characteristics

High thermal load appears at the blade tip and casing of a gas turbine engine. It becomes a significant design challenge to protect the turbine materials from this severe situation. As a result of geometric complexity and experimental limitations, computational fluid dynamics tools have been used to predict blade tip leakage flow aerodynamics and heat transfer at typical engine operating conditions. In this paper, the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer has been studied numerically. Uniform low (444 K) and high (800 K) inlet temperatures and nonuniform (parabolic) temperature profiles have been considered at a fixed rotor rotation speed (9500 rpm). The results showed that the change of flow properties at a higher inlet temperature yields significant variations in the leakage flow aerodynamics and heat transfer relative to the lower inlet temperature condition. Aerodynamic behavior of the tip leakage flow varies significantly with the distortion of turbine inlet temperature. For more realistic inlet condition, the velocity range is insignificant at all the time instants. At a high inlet temperature, reverse secondary flow is strongly opposed by the tip leakage flow and the heat transfer fluctuations are reduced greatly.

Tip leakage flow and heat transfer on turbine blade tip and casing, Part 1: Effect of tip clearance height and rotation speed

Steady simulations were performed to investigate tip leakage flow and heat transfer characteristics on the rotor blade tip and casing in a single-stage gas turbine engine. A typical high-pressure gas turbine stage was modeled with a pressure ratio of 3.2. The predicted isentropic Mach number and adiabatic wall temperature on the casing showed good agreement with available experimental data under similar operating condition. The present numerical study focuses extensively on the effects of tip clearance heights and rotor rotational speeds on the blade tip and casing heat transfer characteristics. It was observed that the tip leakage flow structure is highly dependent on the height of the tip gap and the speed of the rotor. In all cases, the tip leakage flow was seen to separate and recirculate just around the corner of the pressure side of the blade tip. This region of re-circulating flow enlarges with increasing clearance heights. The separated leakage flow reattaches afterwards on the tip surface. Leakage flow reattachment was shown to enhance surface heat transfer at the tip. The interaction between tip leakage flow and secondary flows that is induced by the relative casing motion is found to significantly influence the blade tip and casing heat transfer distribution. A region of critical heat transfer exists on the casing near the blade tip leading edge and along the pressure-side edge for all the clearance heights that were investigated. At high rotation speed, the region of critical heat transfer tends to move towards the trailing edge due to the change in inflow angle.

Tip Leakage Flow and Heat Transfer on Turbine Stage Tip and Casing: Effect of Unsteady Stator-Rotor Interactions

Unsteady simulations were performed to investigate time dependent behaviors of the leakage flow structures and heat transfer on the rotor blade tip and casing in a single stage gas turbine engine. This paper mainly illustrates the unsteady nature of the leakage flow and heat transfer, particularly, that caused by the stator–rotor interactions. In order to obtain time-accurate results, the effects of varying the number of time steps, sub iterations, and the number of vane passing periods was firstly examined. The effect of tip clearance height and rotor speeds was also examined. The results showed periodic patterns of the tip leakage flow and heat transfer rate distribution for each vane passing. The relative position of the vane and vane trailing edge shock with respect to time alters the flow conditions in the rotor domain, and results in significant variations in the tip leakage flow structures and heat transfer rate distributions. It is observed that the trailing edge shock phenomenon results in a critical heat transfer region on the blade tip and casing. Consequently, the turbine blade tip and casing are subjected to large fluctuations of Nusselt number (about Nu = 2000 to 6000 and about Nu = 1000 to 10000, respectively) at a high frequency (coinciding with the rotor speed).

Effect of turbine inlet temperature on blade tip leakage flow and heat transfer

One of the most critical gas turbine engine components, rotor blade tip and casing, are exposed to high thermal load. It becomes a significant design challenge to protect the turbine materials from this severe situation. As a result of geometric complexity and experimental limitations, Computational Fluid Dynamics (CFD) tools have been used to predict blade tip leakage flow aerodynamics and heat transfer at typical engine operating conditions. In this paper, the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer has been studied numerically. Uniform low (LTIT: 444 K) and high (HTIT: 800 K) turbine inlet temperature have been considered. The results showed the higher turbine inlet temperature yields the higher velocity and temperature variations in the leakage flow aerodynamics and heat transfer. For a given turbine geometry and on-design operating conditions, the turbine power output can be increased by 1.48 times, when the turbine inlet temperature increases 1.80 times. Whereas the averaged heat fluxes on the casing and the blade tip become 2.71 and 2.82 times larger, respectively. Therefore, about 2.8 times larger cooling capacity is required to keep the same turbine material temperature. Furthermore, the maximum heat flux on the blade tip of high turbine inlet temperature case reaches up to 3.348 times larger than that of LTIT case. The effect of the interaction of stator and rotor on heat transfer features is also explored using unsteady simulations.

Experimental investigation and numerical simulation of heat transfer in quenching fluidised beds

Fluidisation characteristics at different surfaces of a work-piece of complex geometry are conducted in a fluidised bed at various conditions including fluidising number, bed temperature and fluidising medium. The quenching of the work-piece is performed experimentally. In particular, the major frequency and energy of the pressure fluctuations are measured as a function of either fluidising velocity or heat transfer position and the results are used to develop a mathematic model. A computational model is developed to simulate gas dynamics and heat transfer between the fluidised bed and the work-piece surface, as well as simulating the temperature within the work-piece. The predicted cooling curves are in good agreement with the experimental results. Based on the simulation results, the flow characteristics of the gas and the temperature of the dense gas-solid phase near the work-piece surface are analysed to understand the heat transfer mechanism in the fluidised bed.

Alternating heat transfer between objects of large geometrical difference

The heat transfer on the surface of an object in a gas fluidised bed is sequentially and alternately induced by particle-packet and gas bubble. This phenomenon is studied with computational simulation. The particle-packet and bubble are modelled by a double particle-layers and porous medium model and a hemispherical model, respectively. The heat transfer to and within the object is simulated concurrently. Different grid schemes are applied and different grid sizes are used in meshing the particle-packet and the object as there is a very large difference in their geometrical sizes. Based on theoretical analysis, an approximate method is developed to calculate the heat flux at the surface of the object. The simulation is implemented in a CFD package and the results are compared with experiments.

Numerical analysis of heat transfer and the optimisation of regenerators

To optimize a regenerator’s structure and its operation parameters and, consequently, to increase the efficiency of heat recovery and to save energy, a computational approach is used to study the unsteady three-dimensional flow and heat transfer. The simulation is performed in two steps. In the first step, the gas flow and heat transfer in a typical sphere-bed unit is simulated to deduce a dimensionless equation of heat transfer between gas and sphere. In the second step, a model is developed to simplify the prototype and to simulate the gas flow and heat transfer in the whole regenerator. The heat exchange process in regenerators and the effects of the regenerator’s structure and operation parameters, such as gas mass flux, reversal time, regenerator height, sphere diameter, and thermophysical properties of the spheres, are studied with the model to determine efficiency of heat recovery.

Influence of Goertler vortices spanwise wavelength on heat transfer rates

The boundary layer over concave surfaces can be unstable due to centrifugal forces, giving rise to Goertler vortices. These vortices create two regions in the spanwise direction—the upwash and downwash regions. The downwash region is responsible for compressing the boundary layer toward the wall, increasing the heat transfer rate. The upwash region does the opposite. In the nonlinear development of the Goertler vortices, it can be observed that the upwash region becomes narrow and the spanwise–average heat transfer rate is higher than that for a Blasius boundary layer. This paper analyzes the influence of the spanwise wavelength of the Goertler the heat transfer. The equation is written in vorticity-velocity formulation. The time integration is done via a classical fourth-order Runge-Kutta method. The spatial derivatives are calculated using high-order compact finite difference and spectral methods. Three different wavelengths are analyzed. The results show that steady Goertler flow can increase the heat transfer rates to values close to the values of turbulence, without the existence of a secondary instability. The geometry (and computation domain) are presented

Implementation of the conjugate heat transfer code in KIVA-4

KIVA is an open Computational Fluid Dynamics (CFD) source code that is capable to compute the transient two and three-dimensional chemically reactive fluid flows with spray. The latest version in the family of KIVA codes is the KIVA-4 which is capable of handling the unstructured mesh. This project focuses on the implementation of the Conjugate Heat Transfer code (CHT) in KIVA-4. The previous version of KIVA code with conjugate heat transfer code has been developed at Michigan Technological University by Egel Urip and is be used in this project. During the first phase of the project, the difference in the code structure between the previous version of KIVA and the KIVA-4 has been studied, which is the most challenging part of the project. The second phase involves the reverse engineering where the CHT code in previous version is extracted and implemented in KIVA-4 according to the new code structure. The validation of the implemented code is performed using a 4-valve Pentroof engine case. A solid cylinder wall has been developed using GRIDGEN which surrounds 3/4th of the engine cylinder and heat transfer to the solid wall during one engine cycle (0-720 Crank Angle Degree) is compared with that of the reference result. The reference results are nothing but the same engine case run in the previous version with the original code developed by Egel. The results of current code are very much comparable to that of the reference results which verifies that successful implementation of the CHT code in KIVA-4.

A compact difference scheme for numerical solutions of second order dual-phase-lagging models of microscale heat transfer

Dual-phase-lagging (DPL) models constitute a family of non-Fourier models of heat conduction that allow for the presence of time lags in the heat flux and the temperature gradient. These lags may need to be considered when modeling microscale heat transfer, and thus DPL models have found application in the last years in a wide range of theoretical and technical heat transfer problems. Consequently, analytical solutions and methods for computing numerical approximations have been proposed for particular DPL models in different settings. In this work, a compact difference scheme for second order DPL models is developed, providing higher order precision than a previously proposed method. The scheme is shown to be unconditionally stable and convergent, and its accuracy is illustrated with numerical examples.

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Heat transfer by impinging jets on a moving strip

This paper deals with heat transfer on a moving plate by mean of an impinging jet. Three different turbulence models are used and it turns out that Lam-Bremhorst model is in good agreement with measurements when Re is lower that 5000. In case of moving strip (ratio m=V strip/V jet lower than 1/3), there is almost no effect of m on Nusselt distribution in the stagnation region.