Effects of 3-D Local Unsteadiness on Heat Transfer Prediction in Turbulated Passages

Author(s):  
Pamela A. McDowell ◽  
William D. York ◽  
D. Keith Walters ◽  
James H. Leylek

A newly developed unsteady turbulence model was used to predict heat transfer in a turbulated passage typical of turbine airfoil cooling applications. Comparison of fullyconverged computational solutions to experimental measurements reveal that accurate prediction of heat transfer coefficient requires the effects of local small-scale unsteadiness to be captured. Validation was accomplished through comparison of the time- and area-averaged Nusselt number on the passage wall between adjacent ribs with experimental data from the open literature. The straight channel had a square cross-sectional area with multiple rows of staggered and rounded-edge ribs on opposite walls that were orthogonal to the flow. Simulations were run for Reynolds numbers of 5500, 16500, and 25000. Computational solutions were obtained on a multi-block, multi-topology, unstructured, and adaptive grid, using a pressure-correction based, fully-implicit Navier-Stokes solver. The computational results include two-dimensional (2-D) and three-dimensional (3-D) steady and unsteady simulations with viscous sublayers resolved (y+ ≤ 1) on all the walls in every case. Turbulence closure was obtained using a new turbulence model developed in-house for the unsteady simulations, and a realizable k-ε turbulence model was used for the steady simulations. The results obtained from the unsteady simulations show greatly improved agreement with the experimental data, especially at realistically high Reynolds numbers. The key 3-D physics mechanisms responsible for the successful outcome include: (1) shear layer roll-up over the turbulators; (2) recirculation zones both upstream and downstream of the rib faces; and (3) reattachment regions between each rib pair. Results from the unsteady case are superior to those of the steady because they capture the aforementioned mechanisms, and therefore more accurately predict the heat transfer.

Author(s):  
Alexander Kayne ◽  
Ramesh Agarwal

In recent years Computational Fluid Dynamics (CFD) simulations are increasingly used to model the air circulation and temperature environment inside the rooms of residential and office buildings to gain insight into the relative energy consumptions of various HVAC systems for cooling/heating for climate control and thermal comfort. This requires accurate simulation of turbulent flow and heat transfer for various types of ventilation systems using the Reynolds-Averaged Navier-Stokes (RANS) equations of fluid dynamics. Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS) of Navier-Stokes equations is computationally intensive and expensive for simulations of this kind. As a result, vast majority of CFD simulations employ RANS equations in conjunction with a turbulence model. In order to assess the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for accurate simulations, it is critical to validate the calculations against the experimental data. For this purpose, we use three well known benchmark validation cases, one for natural convection in 2D closed vertical cavity, second for forced convection in a 2D rectangular cavity and the third for mixed convection in a 2D square cavity. The simulations are performed on a number of meshes of different density using a number of turbulence models. It is found that k-epsilon two-equation turbulence model with a second-order algorithm on a reasonable mesh gives the best results. This information is then used to determine the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for flows in 3D enclosures with different ventilation systems. In particular two cases are considered for which the experimental data is available. These cases are (1) air flow and heat transfer in a naturally ventilated room and (2) airflow and temperature distribution in an atrium. Good agreement with the experimental data and computations of other investigators is obtained.


2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.


Author(s):  
Abdulrahman Alenezi ◽  
Abdulrahman Almutairi ◽  
Hamad Alhajeri ◽  
Abdulaziz Gamil ◽  
Faisal Alshammari

Abstract A detailed heat transfer numerical study of a three-dimensional impinging jet on a roughened isothermal surface is presented and is investigated from flow physics vantage point under the influence of different parameters. The effects of the Reynolds number, roughness location, and roughness dimension on the flow physics and heat transfer parameters are studied. Additionally, the relations between average heat transfer coefficient (AHTC) and flow physics including pressure, wall shear and flow vortices with thermodynamic nonequilibrium are offered. This paper studies the effect of varying both location and dimension of the roughness element which took the shape of square cross-sectional continuous ribs to deliver a favorable trade-off between total pressure loss and heat transfer rate. The roughness element was tested for three different radial locations (R/D) = 1, 1.5, and 2 and at each location its height (i.e., width) (e) was changed from 0.25 to 1 mm in incremental steps of 0.25. The study used a jet angle (α) of 90 deg, jet-to-target distance (H/D = 6), and Re ranges from 10,000 to 50,000, where H is the vertical distance between the target plate and jet exit. The results show that the AHTC can be significantly affected by changing the geometry and dimensions of the roughness element. This variation can be either an augmentation of, or decrease in, the (HTC) when compared with the baseline case. An enhancement of 12.9% in the AHTC was achieved by using optimal location and dimensions of the roughness element at specific Reynolds number. However, a diminution between 10% and 30% in (AHTC) was attained by the use of rib height e = 1 mm at Re = 50k. The variation of both rib location and height showed better contribution in increasing heat transfer for low-range Reynolds numbers.


Author(s):  
Robert Pearce ◽  
Peter Ireland ◽  
Ed Dane ◽  
Janendra Telisinghe

Leading edge impingement systems are increasingly being used for high pressure turbine blades in gas turbine engines, in regions where very high heat loads are encountered. The flow structure in such systems can be very complex and high resolution experimental data is required for engine-realistic systems to enable code validation and optimal design. This paper presents spatially resolved heat transfer distributions for an engine-realistic impingement system for multiple different hole geometries, with jet Reynolds numbers in the range of 13,000–22,000. Following this, Reynolds-averaged Navier-Stokes computational fluid dynamics simulations are compared to the experimental data. The experimental results show variation in heat transfer distributions for different geometries, however average levels are primarily dependent on jet Reynolds number. The computational simulations match the shape of the distributions well however with a consistent over-prediction of around 10% in heat transfer levels.


Author(s):  
A. A. Boretti

The paper presents a computer code for steady and unsteady, three-dimensional, compressible, turbulent, viscous flow simulations. The mathematical model is based on the Favre-averaged Navier-Stokes conservation equations, closed by a statistical model of turbulence. Turbulence effects are represented by using a low Reynolds number K-ω model. The numerical model makes use of a finite difference approximation in generalized coordinates for space discretization. The solution of time-dependent, three-dimensional, non-homogeneous, partial differential equations is obtained by solving, in a prescribed, symmetric pattern, three time-dependent, one-dimensional, homogeneous partial differential equations, representing convection and diffusion along each generalized coordinate direction, and one ordinary differential equation, representing generation and destruction. An explicit, multi-step, dissipative, Runge-Kutta scheme is finally adopted for time discretization. The code is applied to simulate the flow through a linear cascade of turbine rotor blades, where detailed experimental data are available. Blade aerodynamic and heat transfer are computed, at variable Reynolds and Mach numbers and turbulence levels, and compared with experimental data. While the aerodynamic prediction is relatively unaffected by the properties of both mathematical and numerical models, the heat transfer prediction proves to be extremely sensitive to models details. Low Reynolds number K-ω turbulence models theoretically reproduce laminar, turbulent and transitional boundary layers. However, their practical use in a Navier-Stokes code does not allow to entirely capture the effects of turbulence intensity and Mach and Reynolds numbers on blade heat transfer.


Mathematics ◽  
2022 ◽  
Vol 10 (2) ◽  
pp. 272
Author(s):  
Chenyu Wu ◽  
Haoran Li ◽  
Yufei Zhang ◽  
Haixin Chen

The accuracy of an airfoil stall prediction heavily depends on the computation of the separated shear layer. Capturing the strong non-equilibrium turbulence in the shear layer is crucial for the accuracy of a stall prediction. In this paper, different Reynolds-averaged Navier–Stokes turbulence models are adopted and compared for airfoil stall prediction. The results show that the separated shear layer fixed k−v2¯−ω (abbreviated as SPF k−v2¯−ω) turbulence model captures the non-equilibrium turbulence in the separated shear layer well and gives satisfactory predictions of both thin-airfoil stall and trailing-edge stall. At small Reynolds numbers (Re~105), the relative error between the predicted CL,max of NACA64A010 by the SPF k−v2¯−ω model and the experimental data is less than 3.5%. At high Reynolds numbers (Re~106), the CL,max of NACA64A010 and NACA64A006 predicted by the SPF k−v2¯−ω model also has an error of less than 5.5% relative to the experimental data. The stall of the NACA0012 airfoil, which features trailing-edge stall, is also computed by the SPF k−v2¯−ω model. The SPF k−v2¯−ω model is also applied to a NACA0012 airfoil, which features trailing-edge stall and an error of CL relative to the experiment at CL>1.0 is smaller than 3.5%. The SPF k−v2¯−ω model shows higher accuracy than other turbulence models.


Author(s):  
Vijay K. Garg ◽  
Raymond E. Gaugler

In order to study the effect of film cooling on the flow and heat transfer characteristics of actual turbine blades, a three-dimensional Navier-Stokes code has been developed. An existing code (Chima and Yokota, 1990) has been modified for the purpose. The code is an explicit finite difference code with an algebraic turbulence model. The thin-layer Navier-Stokes equations are solved using a general body-fitted coordinate system. The effects of film cooling have been incorporated into the code in the form of appropriate boundary conditions at the hole locations on the blade surface. Each hole exit is represented by several control volumes, thus providing an ability to study the effect of hole shape on the film-cooling characteristics. Comparison with experimental data is fair. Further validation of the code is required, however, and in this respect, there is an urgent need for detailed experimental data on actual turbine blades.


1996 ◽  
Vol 118 (2) ◽  
pp. 307-313 ◽  
Author(s):  
A. A. Ameri ◽  
A. Arnone

The effect of transition modeling on the heat transfer predictions from rotating turbine blades was investigated. Three-dimensional computations using a Reynolds-averaged Navier–Stokes code were performed. The code utilized the Baldwin–Lomax algebraic turbulence model, which was supplemented with a simple algebraic model for transition. The heat transfer results obtained on the blade surface and the hub endwall were compared with experimental data for two Reynolds numbers and their corresponding rotational speeds. The prediction of heat transfer on the blade surfaces was found to improve with the inclusion of the transition length model and wake-induced transition effects over the simple abrupt transition model.


Author(s):  
David L. Rigby

Numerical results are presented for flow in a rotating internal passage with a 180 degree turn and ribbed walls. Reynolds numbers ranging from 5200 to 7900, and Rotation numbers of 0.0 and 0.24 were considered. The straight sections of the channel have a square cross section, with square ribs spaced one hydraulic diameter (D) apart on two opposite sides. The ribs have a height of 0.1D and are not staggered from one side to the other. The full three dimensional Reynolds Averaged Navier-Stokes equations are solved combined with the Wilcox k-ω turbulence model. By solving an additional equation for mass transfer, it is possible to isolate the effect of buoyancy in the presence of rotation. That is, heat transfer induced buoyancy effects can be eliminated as in naphthalene sublimation experiments. Heat transfer, mass transfer and flow field results are presented with favorable agreement with available experimental data. It is shown that numerically predicting the reattachment between ribs is essential to achieving an accurate prediction of heat/mass transfer. For the low Reynolds numbers considered, the standard turbulence model did not produce reattachment between ribs. By modifying the wall boundary condition on ω, the turbulent specific dissipation rate, much better agreement with the flow structure and heat/mass transfer was achieved. It is beyond the scope of the present work to make a general recommendation on the ω wall boundary condition. However, the present results suggest that the ω boundary condition should take into account the proximity to abrupt changes in geometry.


Author(s):  
R. J. Boyle ◽  
Forrest E. Ames ◽  
P. W. Giel

An approach to predicting the effects of freestream turbulence on turbine vane and blade heat transfer is described. Four models for predicting the effects of freestream turbulence were incorporated into a Navier-Stokes CFD analysis. Predictions were compared with experimental data in order to identify an appropriate model for use across a wide range of flow conditions. The analyses were compared with data from five vane geometries and from four rotor geometries. Each of these nine geometries had data for different Reynolds numbers. Comparisons were made for twenty four cases. Steady state calculations were done because all experimental data were obtained in steady state tests. High turbulence levels often result in suction surface transition upstream of the throat, while at low to moderate Reynolds numbers the pressure surface remains laminar. A two-dimensional analysis was used because the flow is predominantly two-dimensional in the regions where freestream turbulence significantly augments surface heat transfer. Because the evaluation of models for predicting turbulence effects can be affected by other factors, the paper discusses modeling for transition, relaminarization, and near wall damping. Quantitative comparisons are given between the predictions and data.


Sign in / Sign up

Export Citation Format

Share Document