LES Investigation of Flow and Heat Transfer in a Channel With Dimples and Protrusions

Author(s):  
Mohammad A. Elyyan ◽  
Danesh K. Tafti

LES 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, with form losses increasing from 45% to 80% with an increase in the 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.

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.


Author(s):  
X. Zhang ◽  
D. K. Tafti

High-resolution time-dependent numerical simulations are used to investigate the effect of thickness ratio on fluid flow and heat transfer performance in multi-louvered fins. Results for three fin thickness ratios, two louver angles, and a fin pitch to louver pitch ratio of one are reported for Reynolds number ranging from 50 to 1200. Thickness ratio is found to have a significant effect on flow efficiency, especially in geometries with small louver angles. For small louver angles, increasing thickness to louver pitch ratio from 0.05 to 0.15, decreases the flow efficiency by as much as 35–40%. As expected, increasing thickness ratio increases total pressure drop, most of which results from an increase in form drag. Heat transfer coefficient, on the other hand, is not influenced strongly by the thickness ratio. The increase in flow acceleration and local Reynolds number with increase in thickness ratio, on one hand, is offset by low flow efficiencies and recirculation zones on the other. As a consequence, some heat transfer degradation is found at low Reynolds numbers, however the degradation diminishes as the Reynolds number increases beyond 300. In general, larger thickness ratios lead to a lower ratio of j/f.


Author(s):  
Cody Dowd ◽  
Danesh Tafti

The focus of this research is to predict the flow and heat transfer in a rotating two-pass duct geometry with staggered ribs using Large-Eddy Simulations (LES). The geometry consists of a U-Bend with 17 ribs in each pass. The ribs are staggered with an e/Dh = 0.1 and P/e = 10. LES is performed at a Reynolds number of 100,000, a rotation number of 0.2 and buoyancy parameters (Bo) of 0.5 and 1.0. The effects of Coriolis forces and centrifugal buoyancy are isolated and studied individually. In all cases it is found that increasing Bo from 0.5 to 1.0 at Ro = 0.2 has little impact on heat transfer. It is found that in the first pass, the heat transfer is quite receptive to Coriolis forces which augment and attenuate heat transfer at the trailing and leading walls, respectively. Centrifugal buoyancy, on the other hand has a bigger effect in augmenting heat transfer at the trailing wall than in attenuating heat transfer at the leading wall. On contrary, it aids heat transfer in the second half of the first pass at the leading wall by energizing the flow near the wall. The heat transfer in the second pass is dominated by the highly turbulent flow exiting the bend. Coriolis forces have no impact on the augmentation of heat transfer on the leading wall till the second half of the passage whereas it attenuates heat transfer at the trailing wall as soon as the flow exits the bend. Contrary to phenomenological arguments, inclusion of centrifugal buoyancy augments heat transfer over Coriolis forces alone on both the leading and trailing walls of the second pass.


2009 ◽  
Vol 131 (11) ◽  
Author(s):  
Y.-C. Shih ◽  
J. M. Khodadadi ◽  
K.-H. Weng ◽  
A. Ahmed

The periodic state of laminar flow and heat transfer due to an insulated or isothermal rotating cylinder object in a square cavity is investigated computationally. A finite-volume-based computational methodology utilizing primitive variables is used. Various rotating objects (circle, square, and equilateral triangle) with different sizes are placed in the middle of a square cavity. A combination of a fixed computational grid and a sliding mesh was utilized for the square and triangle shapes. For the insulated and isothermal objects, the cavity is maintained as differentially heated and isothermal enclosures, respectively. Natural convection heat transfer is neglected. For a given shape of the object and a constant angular velocity, a range of rotating Reynolds numbers are covered for a Pr=5 fluid. The Reynolds numbers were selected so that the flow fields are not generally affected by the Taylor instabilities (Ta<1750). The periodic flow field, the interaction of the rotating objects with the recirculating vortices at the four corners, and the periodic channeling effect of the traversing vertices are clearly elucidated. The simulations of the dynamic flow fields were confirmed against experimental data obtained by particle image velocimetry. The corresponding thermal fields in relation to the evolving flow patterns and the skewness of the temperature contours in comparison to the conduction-only case were discussed. The skewness is observed to become more marked as the Reynolds number is lowered. Transient variations of the average Nusselt numbers of the respective systems show that for high Re numbers, a quasiperiodic behavior due to the onset of the Taylor instabilities is dominant, whereas for low Re numbers, periodicity of the system is clearly observed. Time-integrated average Nusselt numbers of the insulated and isothermal object systems were correlated with the rotational Reynolds number and shape of the object. For high Re numbers, the performance of the system is independent of the shape of the object. On the other hand, with lowering of the hydraulic diameter (i.e., bigger objects), the triangle and the circle exhibit the highest and lowest heat transfers, respectively. High intensity of the periodic channeling and not its frequency is identified as the cause of the observed enhancement.


1997 ◽  
Vol 119 (1) ◽  
pp. 61-67 ◽  
Author(s):  
J.-X. Chen ◽  
X. Gan ◽  
J. M. Owen

A superposed radial outflow of air is used to cool two disks that are rotating at equal and opposite speeds at rotational Reynolds numbers up to 1.2 × 106. One disk, which is heated up to 100°C, is instrumented with thermocouples and fluxmeters; the other disk, which is unheated, is made from transparent polycarbonate to allow the measurement of velocity using an LDA system. Measured Nusselt numbers and velocities are compared with computations made using an axisymmetric elliptic solver with a low-Reynolds-number k–ε turbulence model. Over the range of flow rates and rotational speeds tested, agreement between the computations and measurements is mainly good. As suggested by the Reynolds analogy, the Nusselt numbers for contrarotating disks increase strongly with rotational speed and weakly with flow rate; they are lower than the values obtained under equivalent conditions in a rotor–stator system.


Author(s):  
Y.-C. Shih ◽  
J. M. Khodadadi ◽  
K.-H. Weng ◽  
H. F. Oztop

Computational analysis of transient phenomenon followed by the periodic state of laminar flow and heat transfer due to an insulated rotating object in a square cavity is investigated. A finite-volume-based computational methodology utilizing primitive variables is used. Various rotating objects (circle, square and equilateral triangle) with different sizes are placed in the middle of the cavity. A combination of a fixed computational grid with a sliding mesh was utilized for the square and triangle shapes. The cavity is maintained as a differentially-heated enclosure and the motionless insulated object is set in rotation at time t = 0. Natural convection heat transfer is neglected. For a given shape of the object and a constant angular velocity, a range of rotating Reynolds numbers are covered for a Pr = 5 fluid. The Reynolds numbers were selected so that the flow fields are not generally affected by the Taylor instabilities (Ta &lt; 1750). The evolving flow field and the interaction of the rotating objects with the recirculating vortices at the four corners are elucidated. The corresponding thermal fields in relation to the evolving flow patterns and the skewness of the temperature contours in comparison to conduction-only case were discussed. The skewness is observed to become more marked as the Reynolds number is lowered. At the same time, similarity of the thermal fields for various shapes for the same Reynolds number varifies the appropriate selection of the hydraulic diameter. Transient variations of the average Nusselt numbers on the two walls show that for high Re numbers, a quasi-periodic behavior due to the onset of the Taylor instabilities is dominant, whereas for low Re numbers, periodicity of the system is clearly observed. Time-integrated average Nusselt number of the cavity is correlated to the rotational Reynolds number and shape of the object. The triangle object clearly gives rise to high heat transfer followed by the square and circle objects.


2000 ◽  
Vol 123 (2) ◽  
pp. 347-358 ◽  
Author(s):  
P. Bagchi ◽  
M. Y. Ha ◽  
S. Balachandar

Direct numerical solution for flow and heat transfer past a sphere in a uniform flow is obtained using an accurate and efficient Fourier-Chebyshev spectral collocation method for Reynolds numbers up to 500. We investigate the flow and temperature fields over a range of Reynolds numbers, showing steady and axisymmetric flow when the Reynolds number is less than 210, steady and nonaxisymmetric flow without vortex shedding when the Reynolds number is between 210 and 270, and unsteady three-dimensional flow with vortex shedding when the Reynolds number is above 270. Results from three-dimensional simulation are compared with the corresponding axisymmetric simulations for Re>210 in order to see the effect of unsteadiness and three-dimensionality on heat transfer past a sphere. The local Nusselt number distribution obtained from the 3D simulation shows big differences in the wake region compared with axisymmetric one, when there exists strong vortex shedding in the wake. But the differences in surface-average Nusselt number between axisymmetric and three-dimensional simulations are small owing to the smaller surface area associated with the base region. The shedding process is observed to be dominantly one-sided and as a result axisymmetry of the surface heat transfer is broken even after a time-average. The one-sided shedding also results in a time-averaged mean lift force on the sphere.


Author(s):  
M. Elyyan ◽  
A. Rozati ◽  
D. K. Tafti

Flow field and heat transfer for parallel fins with dimples and protrusions are predicted with large-eddy simulations at a nominal Reynolds number based on fin pitch of 15,000. Dimple and protrusion depth and imprint diameter to channel height ratio are 0.4 and 2.0, respectively. The results show that on the dimple side, the flow and heat transfer is dominated by unsteady vorticity generated and ejected out by the separated shear layer in the dimple. The high turbulent energy which results from the unsteady dynamics is mostly responsible for heat transfer augmentation on the dimple side. A maximum augmentation of about 4 occurs in the reattachment zone of the dimple and immediately downstream of it. On the protrusion side, however, the augmentation in heat transfer is dominated by flow impingement at the front of the protrusion, which results in a maximum augmentation of 5.2. The overall heat transfer and friction coefficient augmentations of 2.34 and 6.35 are calculated for this configuration. Pressure drag from the dimple cavity and protrusion contribute 82% of the total pressure drop.


1989 ◽  
Vol 111 (3) ◽  
pp. 301-309 ◽  
Author(s):  
J. Moore ◽  
J. G. Moore ◽  
G. S. Henry ◽  
U. Chaudhry

The effects of Reynolds number on flow through a tip gap are investigated by performing laminar flow calculations for an idealized two-dimensional tip gap geometry. The results of the calculations aid in understanding and reconciliation of low Much number turbine tip gap measurements, which range in tip gap Reynolds number from 100 to 10,000. For the higher Reynolds numbers, both the calculations and the measurements show a large separation off the sharp edge of the blade tip corner. For a high Reynolds number, fully turbulent flow calculations were also made. These also show a large separation and the results are compared with heat transfer measurements. At high Mach numbers, there are complex shock structures in the tip gap. These are modeled experimentally using a water table.


2002 ◽  
Vol 124 (6) ◽  
pp. 1056-1063 ◽  
Author(s):  
Xundan Shi ◽  
J. M. Khodadadi

A finite-volume-based computational study of steady laminar flow and heat transfer (neglecting natural convection) within a lid-driven square cavity due to a single thin fin is presented. The lid moves from left to right and a fixed thin fin can be positioned perpendicular to any of the three stationary walls. Three fins with lengths equal to 5, 10, and 15 percent of the side, positioned at 15 locations were examined for Re=500, 1000, 2000, and Pr=1 (total of 135 cases). Placing a fin on the right wall brings about multi-cell recirculating vortices compared to the case without a fin that exhibits a primary vortex and two small corner cells. A fin slows the flow near the anchoring wall and reduces the temperature gradients, thus degrading heat transfer capacity. A fin positioned near the top right corner of the cavity can reduce heat transfer most effectively in cases with all three different Reynolds numbers and lengths. Regardless of the Reynolds number, placing a fin on the right wall—compared to putting a fin on the left and bottom walls—can always enhance heat transfer on the left wall and at the same time, reduce heat transfer on the bottom, right and top walls. A long fin has the most marked effect on the system’s heat transfer capabilities. Mean Nusselt number was successfully correlated to the Reynolds number, length of the fin and its position.


Sign in / Sign up

Export Citation Format

Share Document