Large Eddy Simulation of Flow and Heat Transfer Mechanism in Matrix Cooling Channel

Author(s):  
Yigang Luan ◽  
Lianfeng Yang ◽  
Bo Wan ◽  
Tao Sun

Gas turbine engines have been widely used in modern industry especially in the aviation, marine and energy fields. The efficiency of gas turbines directly affects the economy and emissions. It’s acknowledged that the higher turbine inlet temperatures contribute to the overall gas turbine engine efficiency. Since the components are subject to the heat load, the internal cooling technology of turbine blades is of vital importance to ensure the safe and normal operation. This paper is focused on exploring the flow and heat transfer mechanism in matrix cooling channels. In order to analyze the internal flow field characteristics of this cooling configuration at a Reynolds number of 30000 accurately, large eddy simulation method is carried out. Methods of vortex identification and field synergy are employed to study its flow field. Cross-sectional views of velocity in three subchannels at different positions have been presented. The results show that the airflow is strongly disturbed by the bending part. It’s concluded that due to the bending structure, the airflow becomes complex and disordered. When the airflow goes from the inlet to the turning, some small-sized and discontinuous vortices are formed. Behind the bending structure, the size of the vortices becomes big and the vortices fill the subchannels. Because of the structure of latticework, the airflow is affected by each other. Airflow in one subchannel can exert a shear force on another airflow in the opposite subchannel. It’s the force whose direction is the same as the vortex that enhances the longitudinal vortices. And the longitudinal vortices contribute to the energy exchange of the internal airflow and the heat transfer between airflow and walls. Besides, a comparison of the CFD results and the experimental data is made to prove that the numerical simulation methods are reasonable and acceptable.

2021 ◽  
Vol 11 (15) ◽  
pp. 7167
Author(s):  
Liang Xu ◽  
Xu Zhao ◽  
Lei Xi ◽  
Yonghao Ma ◽  
Jianmin Gao ◽  
...  

Swirling impinging jet (SIJ) is considered as an effective means to achieve uniform cooling at high heat transfer rates, and the complex flow structure and its mechanism of enhancing heat transfer have attracted much attention in recent years. The large eddy simulation (LES) technique is employed to analyze the flow fields of swirling and non-swirling impinging jet emanating from a hole with four spiral and straight grooves, respectively, at a relatively high Reynolds number (Re) of 16,000 and a small jet spacing of H/D = 2 on a concave surface with uniform heat flux. Firstly, this work analyzes two different sub-grid stress models, and LES with the wall-adapting local eddy-viscosity model (WALEM) is established for accurately predicting flow and heat transfer performance of SIJ on a flat surface. The complex flow field structures, spectral characteristics, time-averaged flow characteristics and heat transfer on the target surface for the swirling and non-swirling impinging jets are compared in detail using the established method. The results show that small-scale recirculation vortices near the wall change the nearby flow into an unstable microwave state, resulting in small-scale fluctuation of the local Nusselt number (Nu) of the wall. There is a stable recirculation vortex at the stagnation point of the target surface, and the axial and radial fluctuating speeds are consistent with the fluctuating wall temperature. With the increase in the radial radius away from the stagnation point, the main frequency of the fluctuation of wall temperature coincides with the main frequency of the fluctuation of radial fluctuating velocity at x/D = 0.5. Compared with 0° straight hole, 45° spiral hole has a larger fluctuating speed because of speed deflection, resulting in a larger turbulence intensity and a stronger air transport capacity. The heat transfer intensity of the 45° spiral hole on the target surface is slightly improved within 5–10%.


Author(s):  
Sunil Patil ◽  
Danesh Tafti

Large eddy simulations of swirling flow and the associated convective heat transfer in a gas turbine can combustor under cold flow conditions for Reynolds numbers of 50,000 and 80,000 with a characteristic Swirl number of 0.7 are carried out. A precursor Reynolds averaged Navier-Stokes (RANS) simulation is used to provide the inlet boundary conditions to the large-eddy simulation (LES) computational domain, which includes only the can combustor. A stochastic procedure based on the classical view of turbulence as a superposition of the coherent structures is used to simulate the turbulence at the inlet plane of the computational domain using the mean flow velocity and Reynolds stress data from the precursor RANS simulation. To further reduce the overall computational resource requirement and the total computational time, the near wall region is modeled using a zonal two layer model (WMLES). A novel formulation in the generalized co-ordinate system is used for the solution of effective tangential velocity and temperature in the inner layer virtual mesh. The WMLES predictions are compared with the experimental data of Patil et al. (2011, “Experimental and Numerical Investigation of Convective Heat Transfer in Gas Turbine Can Combustor,” ASME J. Turbomach., 133(1), p. 011028) for the local heat transfer distribution on the combustor liner wall obtained using robust infrared thermography technique. The heat transfer coefficient distribution on the liner wall predicted from the WMLES is in good agreement with experimental values. The location and the magnitude of the peak heat transfer are predicted in very close agreement with the experiments.


Author(s):  
Sourabh Kumar ◽  
Ryoichi S. Amano

Gas turbine plays a significant role throughout the industrial world. Aircraft propulsion, land-based power generation, and marine propulsion are most notable sectors where gas turbines are extensively used. The power output in these applications can be increased by raising the temperature of the gas entering the turbines. Turbine blades and vanes constrain the temperature of hot gases. For internal cooling design, techniques for heat extraction from the surfaces exposed to hot stream are based on increasing heat transfer areas and the promotion of turbulence of the cooling flow. Heat transfer is enhanced for example due to ribs, bends, rotation and buoyancy effects; all characterizes flow within the channels. Computational Fluid Dynamics (CFD) simulations are carried out using turbulence models like Large Eddy Simulation (LES) and Reynolds stress model (RSM). These CFD simulations were based on advanced computing technology to improve the accuracy of three-dimensional metal temperature prediction that can be applied routinely in the design stage of turbine cooled vanes and blades. The present work is done to study the effect of secondary flow due to the presence of ribs on heat transfer. In this paper, it is obtained by casting repeated continuous V and broken V-shaped ribs on one side of the two passes square channel into the core of blade. Two different combinations of 60° V and Broken 60° V-ribs in the channel are considered. This work is an attempt to collect information about Nusselt number inside the ribbed duct. Large Eddy Simulation (LES) is carried out on the Inlet V and Inverted V outlet continuous and Broken Inlet V and Inverted V-rib arrangements to analyze the flow pattern inside the channel. Hybrid LES/Reynolds Averaged Navier-Strokes (RANS) modeling is used to modify Reynolds stresses using Algebraic Stress Model (ASM), and a CFD strategy is proposed to predict heat transfer across the cooling channel.


Author(s):  
Samer Abdel-Wahab ◽  
Danesh K. Tafti

This paper presents results from large eddy simulation (LES) of fully developed flow in a 90° ribbed duct with rib pitch-to-height ratio P/e = 10 and a rib height to hydraulic diameter ratio e/Dh = 0.1. Three rotation numbers Ro = 0.18, 0.35 and 0.67 are studied at a nominal Reynolds number based on bulk velocity of 20,000. Mean flow and turbulent quantities, together with heat transfer and friction augmentation data are presented. Turbulence and heat transfer are augmented on the trailing surface and reduced at the leading surface. The heat transfer augmentation ratio on the trailing surface asymptotes to a value of 3.7 ± 5% and does not show any further increasing trend as the rotation number increases beyond 0.2. On the other hand, augmentation ratios on the leading surface keep decreasing with an increase in rotation number with values ranging from 1.7 at Ro = 0.18 to 1.2 at Ro = 0.67. Secondary flow cells augment the heat transfer coefficient on the smooth walls by 20% to 30% over a stationary duct. An increase in rotation number from 0.35 to 0.67 decreases the frictional losses from an augmentation ratio of 9.6 to 8.75 and is a consequence of decrease in form drag and wall shear. Overall augmentation compared with a non-rotating duct ranges from +15% to +20% for heat transfer, and +10% to +15% for friction over the range of rotation numbers studied. Comparison of heat transfer augmentation with previous experimental results in the literature shows very good agreement.


2008 ◽  
Vol 130 (4) ◽  
Author(s):  
Mohammad A. Elyyan ◽  
Danesh K. Tafti

Large eddy simulation calculations are conducted for flow in a channel with dimples and protrusions on opposite walls with both surfaces heated at three Reynolds numbers, ReH=220, 940, and 9300, ranging from laminar, weakly turbulent, to fully turbulent, respectively. Turbulence generated by the separated shear layer in the dimple and along the downstream rim of the dimple is primarily responsible for heat transfer augmentation on the dimple surface. On the other hand, augmentation on the protrusion surface is mostly driven by flow impingement and flow acceleration between protrusions, while the turbulence generated in the wake has a secondary effect. Heat transfer augmentation ratios of 0.99 at ReH=220,2.9 at ReH=940, and 2.5 at ReH=9300 are obtained. Both skin friction and form losses contribute to pressure drop in the channel. Form losses increase from 45% to 80% with increasing Reynolds number. Friction coefficient augmentation ratios of 1.67, 4.82, and 6.37 are obtained at ReH=220, 940, and 9300, respectively. Based on the geometry studied, it is found that dimples and protrusions may not be viable heat transfer augmentation surfaces when the flow is steady and laminar.


2012 ◽  
Vol 152-154 ◽  
pp. 1319-1324
Author(s):  
Tao Lu ◽  
Xing Guo Zhu ◽  
Ping Wang ◽  
Wei Yyu Zhu

In the present paper, large-eddy simulation (LES) based on commercial computational fluid dynamics (CFD) software FLUENT for prediction of flow and heat transfer in a mixing T-junction was completed. Mean and root mean square (RMS) temperature and velocity were defined to describe the distributions and fluctuations of temperature and velocity. Numerical results indicate that profiles between symmetrical planes are almost same and the root mean square temperature and velocity close to the center of the main duct in the downstream are larger than those near the main duct wall. The prediction of the fluctuations of temperature and velocity is significant to understand the knowledge of the cause of thermal fatigue in a mixing T-junction.


Sign in / Sign up

Export Citation Format

Share Document