Role of Melt Convection on the Thermal Performance of Heat Sinks with Phase Change Material

The role of melt convection oil the performance of beat sinks with Phase Change Material (PCM) is presented in this paper. The beat sink consists of aluminum plate fins embedded in PCM and heat flux is supplied from the bottom. The design of such a heat sink requires optimization with respect to its geometrical parameters. The objective of the optimization is to maximize the heat sink operation time for the prescribed heat flux and the critical chip temperature. The parameters considered for optimization are fin number and fill thickness. The height and base plate thickness of heat sink are kept constant in the present analysis. An enthalpy based CFD model is developed, which is capable Of Simulating phase change and associated melt convection. The CFD model is Coupled with Genetic Algorithm (GA) for carrying out the optimization. Two cases are considered, one without melt convection (conduction regime) and the other with convection. It is found that the geometrical optimizations of heat sinks are different for the two cases, indicating the importance of inch convection in the design of heat sinks with PCMs.

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. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of Melt Convection on The Optimum Thermal Design Of Heat Sinks With Phase Change Material

In this paper, the role of melt convection on the performance of heat sinks with phase change material (PCM) is investigated numerically. The heat sink consists of aluminum plate fins embedded in PCM, and is subjected to heat flux supplied from the bottom. A single-domain enthalpy-based CFD model is developed, which is capable of simulating the phase change process and the associated melt convection. The CFD model is coupled with a genetic algorithm for carrying out the optimization. Two cases are considered, namely, one without melt convection (i.e., conduction heat transfer analysis), and the other with convection. It is found that the geometrical optimizations of heat sinks are different for the two cases, indicating the importance of melt convection in the design of heat sinks with PCMs. In the case of conduction analysis, the optimum width of half fin (i.e., sum of half pitch and half fin thickness) is a constant, which is in good agreement with results reported in the literature. On the other hand, if melt convection is considered, the optimum half fin width depends on the effective thermal diffusivity due to conduction and convection. With melt convection, the optimized design results in a significant improvement of operational time.

Performance Analysis of Heat Sinks With Phase-Change Materials Subjected to Transient and Cyclic Heating

Phase-change cooling technique is a suitable method for thermal management of electronic equipment subjected to transient or cyclic heat loads. The thermal performance of a phase-change based heat sink under cyclic heat load depends on several design parameters, namely, applied heat flux, cooling heat transfer coefficient, thermophysical properties of phase-change materials (PCMs), and physical dimensions of phase-change storage system during melting and freezing processes. A one-dimensional conduction heat transfer model is formulated to evaluate the effectiveness of preliminary design of practical PCM-based energy storage units. In this model, the phase-change process of the PCM is divided into melting and solidification subprocesses, for which separate equations are written. The equations are solved sequentially and an explicit closed-form solution is obtained. The efficacy of analytical model is estimated by comparing with a finite-volume-based numerical solution for both transient and cyclic heat loads.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Are the available boiling heat transfer coefficients suitable for silicon microchannel heat sinks?

The typical MEMS fabrication of micro evaporators ensures the perfect smooth wall surface that is lack of nucleation sites, significantly decreasing the heat transfer coefficients compared with miniature evaporators fabricated using copper or stainless steel. In the present paper, we performed the boiling heat transfer experiment in silicon triangular microchannel heat sink over a wide parameter range for 102 runs. Acetone was used as the working fluid. The measured boiling heat transfer coefficients versus the local vapor mass qualities are compared with the classical Chen’s correlation and other correlations for macro and miniature capillary tubes. It is found that most of these correlations significantly over-predict the measured heat transfer coefficients. New correlations are given. There are many reasons for such deviations. The major reason is coming from the perfect smooth silicon surface that lowers the heat transfer performances. New theory is recommended for the silicon microchannel heat sink that should be different from metallic capillary tubes.

Numerically Investigating the Effects of Cross-Links in Scaled Microchannel Heat Sinks

Thermal management as a method of heightening performance in miniaturized electronic devices using microchannel heat sinks has recently become of interest to researchers and the industry. One of the current challenges is to design heat sinks with uniform flow distribution. A number of experimental studies have been conducted to seek appropriate designs for microchannel heat sinks. However, pursuing this goal experimentally can be an expensive endeavor. The present work investigates the effect of cross-links on adiabatic two-phase flow in an array of parallel channels. It is carried out using the three dimensional mixture model from the computational fluid dynamics software, FLUENT 6.3. A straight channel and two cross-linked channel models were simulated. The cross-links were located at 1/3 and 2/3 of the channel length, and their widths were one and two times larger than the channel width. All test models had 45 parallel rectangular channels, with a hydraulic diameter of 1.59 mm. The results showed that the trend of flow distribution agrees with experimental results. A new design, with cross-links incorporated, was proposed and the results showed a significant improvement of up to 55% on flow distribution compared with the standard straight channel configuration without a penalty in the pressure drop. Further discussion about the effect of cross-links on flow distribution, flow structure, and pressure drop was also documented.

Surrogate-based analysis and optimization for the design of heat sinks with jet impingement

Heat sinks are widely used for cooling electronic devices and systems. Their thermal performance is usually determined by the material, shape, and size of the heat sink. With the assistance of computational fluid dynamics (CFD) and surrogate-based optimization, heat sinks can be designed and optimized to achieve a high level of performance. In this paper, the design and optimization of a plate-fin-type heat sink cooled by impingement jet is presented. The flow and thermal fields are simulated using the CFD simulation; the thermal resistance of the heat sink is then estimated. A Kriging surrogate model is developed to approximate the objective function (thermal resistance) as a function of design variables. Surrogate-based optimization is implemented by adaptively adding infill points based on an integrated strategy of the minimum value, the maximum mean square error approach, and the expected improvement approaches. The results show the influence of design variables on the thermal resistance and give the optimal heat sink with lowest thermal resistance for given jet impingement conditions. 

Surrogate-based analysis and optimization for the design of heat sinks with jet impingement

Heat sinks are widely used for cooling electronic devices and systems. Their thermal performance is usually determined by the material, shape, and size of the heat sink. With the assistance of computational fluid dynamics (CFD) and surrogate-based optimization, heat sinks can be designed and optimized to achieve a high level of performance. In this paper, the design and optimization of a plate-fin-type heat sink cooled by impingement jet is presented. The flow and thermal fields are simulated using the CFD simulation; the thermal resistance of the heat sink is then estimated. A Kriging surrogate model is developed to approximate the objective function (thermal resistance) as a function of design variables. Surrogate-based optimization is implemented by adaptively adding infill points based on an integrated strategy of the minimum value, the maximum mean square error approach, and the expected improvement approaches. The results show the influence of design variables on the thermal resistance and give the optimal heat sink with lowest thermal resistance for given jet impingement conditions.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

Numerical Study Of Heat Transfer From Pin Fin Heat Sink Using Steady And Pulsated Impinging Jets

The work reported in this thesis is an attempt to enhance heat transfer in electronic devices with the use of impinging air jets on pin-finned heat sinks. The cooling per-formance of electronic devices has attracted increased attention owing to the demand of compact size, higher power densities and demands on system performance and re-liability. Although the technology of cooling has greatly advanced, the main cause of malfunction of the electronic devices remains overheating. The problem arises due to restriction of space and also due to high heat dissipation rates, which have increased from a fraction of a W/cm2to 100s of W /cm2. Although several researchers have at-tempted to address this at the design stage, unfortunately the speed of invention of cooling mechanism has not kept pace with the ever-increasing requirement of heat re- moval from electronic chips. As a result, efficient cooling of electronic chip remains a challenge in thermal engineering. Heat transfer can be enhanced by several ways like air cooling, liquid cooling, phase change cooling etc. However, in certain applications due to limitations on cost and weight, eg. air borne application, air cooling is imperative. The heat transfer can be increased by two ways. First, increasing the heat transfer coefficient (forced convec- tion), and second, increasing the surface area of heat transfer (finned heat sinks). From previous literature it was established that for a given volumetric air flow rate, jet im-pingement is the best option for enhancing heat transfer coefficient and for a given volume of heat sink material pin-finned heat sinks are the best option because of their high surface area to volume ratio. There are certain applications where very high jet velocities cannot be used because of limitations of noise and presence of delicate components. This process can further be improved by pulsating the jet. A steady jet often stabilizes the boundary layer on the surface to be cooled. Enhancement in the convective heat transfer can be achieved if the boundary layer is broken. Disruptions in the boundary layer can be caused by pulsating the impinging jet, i.e., making the jet unsteady. Besides, the pulsations lead to chaotic mixing, i.e., the fluid particles no more follow well defined streamlines but move unpredictably through the stagnation region. Thus the flow mimics turbulence at low Reynolds number. The pulsation should be done in such a way that the boundary layer can be disturbed periodically and yet adequate coolant is made available. So, that there is not much variation in temperature during one pulse cycle. From previous literature it was found that square waveform is most effective in enhancing heat transfer. In the present study the combined effect of pin-finned heat sink and impinging slot jet, both steady and unsteady, has been investigated for both laminar and turbulent flows. The effect of fin height and height of impingement has been studied. The jets have been pulsated in square waveform to study the effect of frequency and duty cycle. This thesis attempts to increase our understanding of the slot jet impingement on pin-finned heat sinks through numerical investigations. A systematic study is carried out using the finite-volume code FLUENT (Version 6.2) to solve the thermal and flow fields. The standard k-ε model for turbulence equations and two layer zonal model in wall function are used in the problem Pressure-velocity coupling is handled using the SIMPLE algorithm with a staggered grid. The parameters that affect the heat transfer coefficient are: height of the fins, total height of impingement, jet exit Reynolds number, frequency of the jet and duty cycle (percentage time the jet is flowing during one complete cycle of the pulse). From the studies carried out it was found that: a) beyond a certain height of the fin the rate of enhancement of heat transfer becomes very low with further increase in height, b) the heat transfer enhancement is much more sensitive to any changes at low Reynolds number than compared to high Reynolds number, c) for a given total height of impingement the use of fins and pulsated jet, increases the effective heat transfer coefficient by almost 200% for the same average Reynolds number, d) for all the cases it was observed that the optimum frequency of impingement is around 50 − 100 Hz and optimum duty cycle around 25-33.33%, e) in the case of turbulent jets the enhancement in heat transfer due to pulsations is very less compared to the enhancement in case of laminar jets.

Thermal Management of Electronics Using PCM-Based Heat Sink Subjected to Cyclic Heat Load

Designing a heat sink based on a phase change material (PCM) under cyclic loading is a critical issue. For cyclic operation, it is required that the fraction of the PCM melting during the heating cycle should completely resolidify during the cooling period, so that that thermal storage unit can be operated for an unlimited number of cycles. Accordingly, studies are carried out to find the parameters influencing the behavior of a PCM under cyclic loading. A number of parameters are identified in the process, the most important ones being the duty cycle and heat transfer coefficient (h) for cooling. The required h or the required cooling period for complete resolidification for infinite cyclic operation of a conventional PCM-based heat sink is found to be very high and unrealistic with air cooling from the surface. To overcome this problem, the conventional design is modified where h and the area exposed to heat transfer can be independently controlled. With this arrangement, the enhanced area provided for cooling keeps h within realistic limits. Analytical investigation is carried out to evaluate the thermal performance of this modified PCM-based heat sink in comparison to those with conventional designs. Experiments are also performed on both the conventional and the modified PCM-based heat sinks to validate the new findings.

Role of melt convection on optimization of PCM-based heat sink under cyclic heat load

In this article, we study the thermal performance of phase-change material (PCM)-based heat sinks under cyclic heat load and subjected to melt convection. Plate fin type heat sinks made of aluminum and filled with PCM are considered in this study. The heat sink is heated from the bottom. For a prescribed value of heat flux, design of such a heat sink can be optimized with respect to its geometry, with the objective of minimizing the temperature rise during heating and ensuring complete solidification of PCM at the end of the cooling period for a given cycle. For given length and base plate thickness of a heat sink, a genetic algorithm (GA)-based optimization is carried out with respect to geometrical variables such as fin thickness, fin height, and the number of fins. The thermal performance of the heat sink for a given set of parameters is evaluated using an enthalpy-based heat transfer model, which provides the necessary data for the optimization algorithm. The effect of melt convection is studied by taking two cases, one without melt convection (conduction regime) and the other with convection. The results show that melt convection alters the results of geometrical optimization.

Numerical simulation of heat transfer in rectangular microchannel

Numerical simulation of heat transfer in a high aspect ratio rectangular microchannel with heat sinks has been conducted, similar to an experimental study. Three channel heights measuring 0.3 mm, 0.6mmand 1mmare considered and the Reynolds number varies from 300 to 2360, based on the hydraulic diameter. Simulation starts with the validation study on the Nusselt number and the Poiseuille number variations along the channel streamwise direction. It is found that the predicted Nusselt number has shown very good agreement with the theoretical estimation, but some discrepancies are noted in the Poiseuille number comparison. This observation however is in consistent with conclusions made by other researchers for the same flow problem. Simulation continues on the evaluation of heat transfer characteristics, namely the friction factor and the thermal resistance. It is found that noticeable scaling effect happens at small channel height of 0.3 mm and the predicted friction factor agrees fairly well with an experimental based correlation. Present simulation further reveals that the thermal resistance is low at small channel height, indicating that the heat transfer performance can be enhanced with the decrease of the channel height.