Heat transfer correlations for PCM-based heat sinks with plate fins

Characterization of melting process in a Phase Change Material (PCM)-based heat sink with plate fin type thermal conductivity enhancers (TCEs) is numerically studied in this paper. Detailed parametric investigations are performed to find the effect of aspect ratio of enclosure and the applied heat flux on the thermal performance of the heat sinks. Various non-dimensional numbers, such as Nusselt number (Nu), Rayleigh number (Ra), Stefan number (Ste) and Fourier number (Fo) based on a characteristic length scale, are identified as important parameters. The half fin thickness and the fin height are varied to obtain a wide range of aspect ratios of an enclosure. It is found that a single correlation of Nu with Ra is not applicable for all aspect ratios of enclosure with melt convection taken into account. To find appropriate length scales, enclosures with different aspect ratios are divided into three categories, viz. (a) shallow enclosure, (b) rectangular enclosure and (c) tall enclosure. Accordingly, an appropriate characteristic length scale is identified for each type of enclosure and correlation of Nu with Ra based on that characteristic length scale is developed.     

Natural convection heat transfer in horizontal and vertical closed narrow enclosures with heated rectangular finned base plate

Natural convection heat transfer and fluid flow characteristics in horizontal and vertical narrow enclosures with heated rectangular finned base plate have been experimentally investigated at a wide range of Rayleigh number (Ra) for different fin spacings and fin lengths. Quantitative comparisons of finned surface effectiveness (ε) and heat transfer rate between horizontal and vertical enclosures have been reported. In comparison with enclosure of a bare base plate, insertion of heat conducting fins always enhances heat transfer rate. Optimization of fin-array geometry has been addressed. The results gave an optimum fin spacing at which Nusselt number (Nu_H) and finned surface effectiveness (ε) are maximum. It has been found that: (1) increasing fin length increases Nu_H and ε; (2) increasing Ra increasesNu_H for any fin-array geometries and (3) for any fin-array geometry and at Ra > 10,000, increasing Ra decreases ε while for fin-array geometries of large fin spacing and at Ra < 10,000, increasing Ra increases ε. Useful design guidelines have been suggested. Correlations of Nu_H have been developed for horizontal and vertical enclosures. Correlations predictions have been compared with previous data and good agreement was found.     

Natural convection heat transfer of molten salt nanofluids around vertical array of heated horizontal cylinders

Abstract

To analyze the natural convection heat transfer of molten salt nanofluids around vertical array of horizontal cylinders, a numerical simulation is performed with cylinder numbers in the range of N?=?2–8 and pitches in the range of S/D?=?5–10. The results show that the heat transfer of nanofluids around each cylinder is affected by the position in the tube row and its distance from the adjacent cylinder. The average Nusselt number (Nu_a) of the natural convection heat transfer over the whole array of cylinders is determined by cylinder spacing S/D, cylinder number N, and Ra. When S/D?=?5, Nu_a decreases as the cylinder number increases. When S/D?=?10, Nua increases as the cylinder number increases. Compared with the molten salts, the natural convection heat transfer of nanofluids is enhanced. This study provides a theoretical basis for the design of a single-tank energy storage system.     

Heat generation effects in natural convection inside a porous annulus

The interplay between internal heat generation and externally driven natural convection inside a porous medium annulus is studied in detail using numerical methods. The axisymmetric domain is bounded with adiabatic top and bottom walls and differentially heated side walls sustaining steady natural convection of a fluid with Prandtl number, Pr = 5, through a porous matrix of volumetric porosity, ? = 0.4. The generalized momentum equation with Brinkman–Darcy–Forchheimer terms and the local thermal non-equilibrium based two-energy equation model are solved to determine the flow and the temperature distribution. Beyond a critical heat generation value defined using an internal Rayleigh number, Ra_I,cr^?, the convection transits from unicellular to bicellular mode, as the annulus T_max becomes higher than the fixed hot-wall temperature. The Ra_I,cr^? increases proportionately when the permeability based external Rayleigh number Ra_E^? and the solid–fluid thermal conductivity ratio γ are independently increased. A correlation is proposed to predict the overall annulus Nu as a function of Ra_E^?, Ra_I^?, Da and γ. It predicts the results within ± 20% accuracy.     

Visualization of heat flow due to natural convection within triangular cavities using Bejan’s heatline concept

Natural convection and flow circulation within a cavity has received significant attention in recent times. The wide range of applicability of flow inside a cavity (food processing industries, molten metal industries, etc.) requires thorough understanding for cost efficient processes. This paper is based on comprehensive analysis of heat flow pattern using Bejan’s heatline concept. The key parameters for our study are the Prandtl number, Rayleigh number and Nusselt number. The values of Prandtl number (0.015, 0.026, 0.7 and 1000) have been chosen based on wide range of applicability. The Rayleigh number has been varied from 10² to 10⁵. Interesting results were obtained. For low Rayleigh number, it is found that the heatlines are smooth and perfectly normal to the isotherms indicating the dominance of conduction. But as Ra increases, flow slowly becomes convection dominant. It is also observed that multiple secondary circulations are formed for fluids with low Pr whereas these features are absent in higher Pr fluids. Multiple circulation cells for smaller Pr also correspond multiple cells of heatlines which illustrate less thermal transport from hot wall. On the other hand, the dense heatlines at bottom wall display enhanced heat transport for larger Pr. Further, local heat transfer (Nu_l, Nu_t) are explained based on heatlines. The comprehensive analysis is concluded with the average Nusselt number plots. A correlation for average heat transfer rate and Ra has been developed and the range of Rayleigh number is also found, to depict the conduction dominant heat transfer.     

Passive Laminar Heat Transfer Across Porous Cavities Using the Thermal Non-Equilibrium Model

This work presents a study on laminar free convection within a square cavity filled with a fluid saturated porous medium. Macroscopic flow equations are obtained by volume-averaging local instantaneous continuity and momentum equations. The so-called “two-energy equation model” is used, in which distinct macroscopic equations are applied to the working fluid and the solid material. Transport equations are discretized using the control-volume method and the system of algebraic equations is relaxed via the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm. The effect of Ra_m on Nu_w correctly predicted the enhancement of passive heat transfer across the cavity for increasing Ra_m. Increasing k_s/k_f enhances the conduction transport through the solid material and, consequently, dampens the overall Nusselt number, defined here as the ratio between conduction and convection mechanisms over conduction transport only. Further, results indicate that by increasing the void space within the porous material the overall Nusselt number is reduced rather than increased. Individual contributions to the average Nusselt number indicate that, although convection is enhanced with increasing porosity, the reduction of conduction heat transfer through the solid material is the controlling mechanics for Nu_w as porosity increases. The results herein might contribute to design and optimization of passive heat transfer systems.     

Hydromagnetic effect on mixed convection in a lid-driven cavity with sinusoidal corrugated bottom surface

The present numerical simulation is conducted to analyze the mixed convection flow and heat transfer in a lid-driven cavity with sinusoidal wavy bottom surface in presence of transverse magnetic field. The enclosure is saturated with electrically conducting fluid. The cavity vertical walls are insulated while the wavy bottom surface is maintained at a uniform temperature higher than the top lid. In addition, the transport equations are solved by using the finite element formulation based on the Galerkin method of weighted residuals. The implications of Reynolds number (Re), Hartmann number (Ha) and number of undulations (λ) on the flow structure and heat transfer characteristics are investigated in detail while, Prandtl number (Pr) and Rayleigh number (Ra) are considered fixed. The trend of the local heat transfer is found to follow a wavy pattern. The results of this investigation illustrate that the average Nusselt number (Nu) at the heated surface increases with an increase of the number of waves as well as the Reynolds number, while decreases with increasing Hartmann number.     

Laminar heat transfer in a porous channel simulated with a two-energy equation model

Laminar heat transfer in a porous channel is numerically simulated with a two-energy equation model for conduction and convection. Macroscopic equations for continuity, momentum and energy transport for the fluid and solid phases are presented. The numerical methodology employed is based on the control volume approach with a boundary-fitted non-orthogonal coordinate system. Fully developed forced convection in a porous channel bounded by parallel plates is considered. Solutions for Nusselt numbers along the channel are presented for laminar flows. Results simulate the effects Reynolds number Re, porosity, particle size and solid-to-fluid thermal conductivity ratio on Nusselt sumber, Nu, which is defined for both the solid and fluid phases. High Re, low porosities, low particle diameters and low thermal conductivity ratios promote thermal equilibrium between phases leading to higher values of Nu.     

Natural convection slip flow in a vertical microchannel heated at uniform heat flux

Numerical solutions for steady state developing natural convection flow in air, in vertical parallel-plate microchannels are accomplished. An asymmetric heating is considered and the walls are assumed to be at uniform heat flux. A first-order model is used for slip and jump boundary conditions and an analytical solution for the fully developed flow is also given. Results are performed for air, for the heat flux ratio in the 0.0–1.0 range, for Rayleigh, Ra, and Knudsen, Kn, numbers from 10^?1 to 8 × 10³ and from 0.0 to 0.10, respectively. The maximum mass flow rate is always obtained for the highest considered Kn value, whereas the average Nusselt number, Nu, increases for lower Ra (<10) and decreases for Ra > 100. Wall temperature profiles have the lowest values for highest considered Kn value at lower Ra, whereas for the developing flow, they present opposite trends. For developing flow, velocity profiles for asymmetric and symmetric heating are completely different. In developing flow velocity profiles along the wall present the highest increases for asymmetric heating and the highest considered Kn value.     

Analysis of microchannel heat sinks for electronics cooling

This paper presents an analytical and numerical study on the heat transfer characteristics of forced convection across a microchannel heat sink. Two analytical approaches are used: the porous medium model and the fin approach. In the porous medium approach, the modified Darcy equation for the fluid and the two-equation model for heat transfer between the solid and fluid phases are employed. Firstly, the effects of channel aspect ratio (α_s) and effective thermal conductivity ratio (k?) on the overall Nusselt number of the heat sink are studied in detail. The predictions from the two approaches both show that the overall Nusselt number (Nu) increases as α_s is increased and decreases with increasing k?. However, the results also reveal that there exists significant difference between the two approaches for both the temperature distributions and overall Nusselt numbers, and the discrepancy becomes larger as either α_s or k? is increased. It is suggested that this discrepancy can be attributed to the indispensable assumption of uniform fluid temperature in the direction normal to the coolant flow invoked in the fin approach. The effect of porosity (ε) on the thermal performance of the microchannel is subsequently examined. It is found that whereas the porous medium model predicts the existence of an optimal porosity for the microchannel heat sink, the fin approach predicts that the heat transfer capability of the heat sink increases monotonically with the porosity. The effect of turbulent heat transfer within the microchannel is next studied, and it is found that turbulent heat transfer results in a decreased optimal porosity in comparison with that for the laminar flow. A new concept of microchannel cooling in combination with microheat pipes is proposed, and the enhancement in heat transfer due to the heat pipes is estimated. Finally, two-dimensional numerical calculations are conducted for both constant heat flux and constant wall temperature conditions to check the accuracy of analytical solutions and to examine the effect of different boundary conditions on the overall heat transfer.     

Heat transfer characteristics of a slot jet impinging on a semi-circular convex surface

Surface heat transfer characteristics of a heated slot jet impinging on a semi-circular convex surface have been investigated by using the transient heating liquid crystal technique. Free jet velocity, turbulence and temperature characteristics have been determined by using a combination of an X-wire and a cold wire anemometry. The parametric effects of jet Reynolds number (Re_W) ranging from 5600 to 13,200 and the dimensionless slot nozzle-to-impingement surface distance (Y/W) ranges from 2 to 10 on the local circumferential heat transfer have been studied. Local circumferential Nusselt number (Nu_S) decreases with increasing the dimensionless circumferential distance (S/W) from its maximum value at the stagnation point up to S/W=3.1. The transition in the wall jet from laminar to turbulent flow was completed by about 3.3?S/W?4.2 which coincided with a secondary peak in heat transfer. Correlations of local and average Nusselt numbers with Re_W, Y/W and S/W have been established for the stagnation point and the circumferential distribution. The rate of decay of average circumferential Nusselt numbers around the semi-circular convex surface is much faster than that which occurs laterally along the flat surface. As Y/W increases, the effect of surface curvature becomes apparent and the difference between the flat surface correlation and the convex surface becomes more pronounced.     

Laminar and turbulent free convection in a composite enclosure

Turbulent natural convection in a two-dimensional horizontal composite square cavity, isothermally heated at the left side and cooled from the opposing surface, is numerically analyzed using the finite volume method. The composite square cavity is formed by three distinct regions, namely, clear, porous and solid region. The development of a numerical tool able to treat all these regions as one computational domain is of advantage for engineering design and analysis of passive thermal control systems. Governing equations are written in terms of primitive variables and are recast into a general form. It was found that the fluid begins to permeate the porous medium for values of Ra greater than 10⁶. Nusselt number values show that for the range of Ra analyzed there is no significant variation between the laminar and turbulent model solution. When comparing the effects of Ra, k_s/k_f and Da on Nu, results indicate that the solid phase properties have a greater influence in enhancing the overall heat transferred trough the cavity.     

Lattice Boltzmann Simulation of Natural Convection in an Inclined Square Cavity with Spatial Temperature Variation

In this paper, natural convection heat transfer in an inclined square cavity filled with pure air (Pr = 0.71) was numerically analyzed with the lattice Boltzmann method. The heat source element is symmetrically embedded over the center of the bottom wall, and its temperature varies sinusoidally along the length. The top and the rest part of the bottom wall are adiabatic while the sidewalls are fixed at a low temperature. The influences of heat source length, inclination angle, and Rayleigh number (Ra) on flow and heat transfer were investigated. The Nusselt number (Nu) distributions on the heat source surface, the streamlines, and the isotherms were presented. The results show that the inclination angle and heat source length have a significant impact on the flow and temperature fields and the heat transfer performance at high Rayleigh numbers. In addition, the average Nu firstly increases with γ and reaches a local maximum at around γ = 45°, then decreases with increasing γ and reaches minimum at γ = 180° in the cavity with ? = 0.4. Similar behaviors are observed for ? = 0.2 at Ra = 10⁴. Moreover, nonuniform heating produces a significant different type of average Nu and two local minimum average Nu values are observed at around γ = 45° and γ = 180° for Ra = 10⁵ in the cavity with ? = 0.2.     

Aspect ratio effect on natural convection in a square enclosure with a sinusoidal active thermal wall using a thermal non-equilibrium model

ABSTRACT

A numerical investigation of the aspect ratio effect on natural convection in a square enclosure is carried out by adopting the local thermal non-equilibrium model. The top and bottom walls of the enclosure are adiabatic, the left vertical wall is partially heated and cooled by the sinusoidal thermal boundary condition, and the right vertical wall is maintained at uniform thermal boundary condition. The results show the value of periodicity parameter increasing. The streamlines vary in different patterns, rotating clockwise and counterclockwise simultaneously when N > 1, and the number of clockwise and counterclockwise rotating cells increases with the increase of N and equals the value of N. The sinusoidal local Nusselt number profiles are observed and the wave amplitude of local Nusselt number decreases with the increase of aspect ratio, and the absolute values of average Nusselt number at left wall of porous cavity reach maximum when Ar = 1. The absolute value of solid-to-fluid temperature differences decreases as the inter-phase heat transfer coefficient (H) increases and it increases as the value of aspect ratio increases. The total heat transfer of porous cavity can be enhanced by increasing the aspect ratio and the thermal conductivity ratio.     

Developing laminar mixed convection with heat and mass transfer in horizontal and vertical tubes

The effects of the solutal and thermal Grashof numbers on the flow, temperature and concentration fields in tubes with uniform heat flux and concentration at the fluid-solid interface have been investigated numerically using a three-dimensional axially parabolic model. Results show a complex development of the flow field which is strongly influenced by the values of the two Grashof numbers and by the tube inclination. For vertical tubes the flow field is also influenced by the relative direction of the flow and the buoyancy forces. In general, very close to the tube inlet forced convection boundary layer development dominates. Further downstream, the effects of solutal buoyancy predominate while those of thermal origin determine the flow field far downstream and, in particular, the fully developed conditions. The axial evolution of the wall shear stress τ_z, the Nusselt number Nu_z and the Sherwood number Sh_z in both horizontal and vertical tubes are presented for different combinations of the two Grashof numbers. For horizontal tubes and vertical tubes with upward flow these three variables are greater than the corresponding ones for forced convection. The opposite is true for downward flow in vertical tubes.     

Natural convection from narrow horizontal plates at moderate Rayleigh numbers

The present work deals with the natural convection flow and heat transfer from a horizontal plate cooled from above. Experiments are carried out for rectangular plates having aspect ratios between φ=0.036 and 0.43 and Rayleigh numbers in the range 290?Ra_w?3.3×10⁵. These values of Ra_w and φ have been selected below those commonly considered in previous research in view of a future application to the design of printed circuit boards. The plates are made of two different metals, copper and steel. The choice of a metal is relevant to the present problem because the plates are heated by means of an electric current. Important variations of the surface temperature are observed along the transverse direction for the steel plates. The surface of the copper plates is almost isothermal because of the high thermal conductivity of the metal.Calculations for a semi-infinite plate are carried out to predict the transverse profiles of the surface temperature and heat flux and to visualize the structure of the flow. Three-dimensional calculations are also used at a qualitative level to observe the changes in the flow structure due to the finite length of the plate. Present results are compared with both previous experimental work and analyses that are based on boundary layer theory. It is shown that analyses for an infinite boundary layer are not completely applicable to the present problem because of its different physics. The most relevant feature of the natural convection flow, which is not predicted by boundary layer analyses, is a thermal plume rising near the center of the plate.Present heat transfer results differ from previous experimental work because of the lower Rayleigh numbers and aspect ratios investigated here. The Nusselt number is found to depend on Ra_wⁿ, with the exponent n=0.17 being lower than most of the values reported in the literature. This comparatively low value is related to the transverse conduction of heat through the air, which becomes increasingly significant as Ra_w approaches zero. It is shown that such a low-Ra_w effect can be accounted for in a physically consistent manner by adding a constant term to the heat transfer correlation. On the other hand, it is found that the Nusselt number does not significantly depend on the aspect ratio in the range of φ investigated contrary to what has been previously reported for wider plates.     

Unsteady analysis of natural convection in a carbon nanotube-water filled cavity with an inclined heater

ABSTRACT

A finite element solution has been performed in this work to solve unsteady governing equations of natural convection in a carbon nanotube–water-filled cavity with inclined heater. The temperature of ceiling and left vertical walls is lower than that of the heater while the other walls are adiabatic. The main governing parameters are nanofluid volume fraction and Rayleigh number (Ra). It is found that the heat transfer rate shows different trends based on Rayleigh number and it increases with increase in nanoparticle volume fraction. It has been estimated that average Nusselt number (Nu) is dependent onRa through power regression models with strong positive linear correlation between ln (Nu) and ln (Ra). In particular, for the maximum time, when the solid volume fraction is varied from 0 to 0.1 the dependence between average Nu and linear Ra, on a logarithmic scale, is very high.     

Thermal description of pseudosteady-state natural convection inside a vertical cylinder

The heat transfer correlation, minimum temperature, and mean temperature are presented for pseudosteady-state natural convection heat transfer to a fluid inside a vertical cylinder. The SIMPLER numerical method was used for calculation in the range 0.25 < H/D <2, Ra < 10⁷, and Pr = 7. This range includes conduction to weak turbulence. The overall heat transfer for the convection-dominated range was found to be correlated by Nu = 0.519Ra^0.255 where the temperature difference for both the Nusselt and Rayleigh numbers was the center temperature minus the inside wall temperature. Correlations using other temperature differences are also presented and provide a method for prediction of the mean temperature, minimum temperature, or center temperature by knowing any one of them.     

Natural convection heat transfer from a horizontal cylinder embedded in a porous medium

This paper presents the results of an experimental investigation of heat transfer by natural convection from a horizontal cylinder embedded in porous media consisting of randomly packed glass spheres saturated by either water or silicone oil. It is shown that the overall range of the Rayleigh number, Ra, can be divided into two subregions, called ‘low’ and ‘high’, in each of which the Nusselt number, Nu, behaves differently. It is demonstrated that the low Ra region corresponds to Darey flow and the high to Forchheimer flow. Correlation equations for Nu for the Darcy regime are presented that account for viscous dissipation, and others for the Forchheimer regime that involve the first and second Forchheimer coefficients. The variation of properties with temperature and the wall effect on porosity (and consequently on heat transfer) are considered. The paper includes information concerning the resistance to flow in porous media that was obtained in conjunction with the heat transfer study.     

Experimental study of forced and free convective heat transfer in the thermal entry region of horizontal concentric annuli

Forced and free convective heat transfer for thermally developing and thermally fully developed laminar air flow inside horizontal concentric annuli in the thermal entrance length has been experimentally investigated. The experimental setup consists of a stainless steel annulus having a radius ratio of 2 and an inner tube with a heated length of 900 mm subjected to a constant wall heat flux boundary condition and an adiabatic outer annulus. The investigation covers Reynolds number range from 200 to 1000, the Grashof number was ranged from 6.2 × 10⁵ to 1.2 × 10⁷. The entrance sections used were long tube with length of 2520 mm (L/D_h = 63) and short tube with length of 504 mm (L/D_h = 12.6). The surface temperature distribution along the inner tube surface, and the local Nusselt number distribution versus dimensionless axial distance Z_t were presented and discussed. It is inferred that the free convection effects tended to decrease the heat transfer at low Re number while to increase the heat transfer for high Re number. This investigation reveals that the Nusselt number values were considerably greater than the corresponding values for fully developed combined convection over a significant portion of the annulus. The average heat transfer results were correlated in terms of the relevant dimensionless variables with an empirical correlation. The local Nusselt number results were compared with available literature and show similar trend and satisfactory agreement.