Prediction of Dryout in Flat Heat Pipes at High Heat Fluxes From Multiple Discrete Sources

Author(s):  
Unnikrishnan Vadakkan ◽  
Suresh V. Garimella ◽  
Jayathi Y. Murthy

A three-dimensional model has been developed to analyze the transient and steady-state performance of flat heat pipes with discrete heat sources. Three-dimensional flow and energy equations are solved in the vapor and liquid regions, along with conduction in the wall. Saturated flow models are used for heat transfer and fluid flow through the wick. In the wick region, the analysis uses an equilibrium model for heat transfer and a Brinkman-Forchheimer extended Darcy model for fluid flow. Averaged properties weighted with the porosity are used for the wick analysis. The state equation is used in the vapor core to relate density change to the operating pressure. The density change due to pressurization of the vapor core is accounted for in the continuity equation. Vapor flow, temperature and hydrodynamic pressure fields are computed at each time step from coupled continuity/momentum and energy equations in the wick and vapor regions. The mass flow rate at the interface is obtained from the application of kinetic theory. Predictions are made for the magnitude of heat flux at which dryout would occur in a flat heat pipe. The input heat flux and the spacing between the discrete heat sources are studied as parameters. The location in the heat pipe at which dryout is initiated is found to be different from that of the maximum temperature. The location where the maximum capillary pressure head is realized also changes during the transient. Axial conduction through the wall and wick are seen to play a significant role in determining the axial temperature variation.

2004 ◽  
Vol 126 (3) ◽  
pp. 347-354 ◽  
Author(s):  
Unnikrishnan Vadakkan ◽  
Suresh V. Garimella ◽  
Jayathi Y. Murthy

A three-dimensional model has been developed to analyze the transient and steady-state performance of flat heat pipes subjected to heating with multiple discrete heat sources. Three-dimensional flow and energy equations are solved in the vapor and liquid regions, along with conduction in the wall. Saturated flow models are used for heat transfer and fluid flow through the wick. In the wick region, the analysis uses an equilibrium model for heat transfer and a Brinkman-Forchheimer extended Darcy model for fluid flow. Averaged properties weighted with the porosity are used for the wick analysis. The state equation is used in the vapor core to relate density change to the operating pressure. The density change due to pressurization of the vapor core is accounted for in the continuity equation. Vapor flow, temperature and hydrodynamic pressure fields are computed at each time step from coupled continuity/momentum and energy equations in the wick and vapor regions. The mass flow rate at the interface is obtained from the application of kinetic theory. Predictions are made for the magnitude of heat flux at which dryout would occur in a flat heat pipe. The input heat flux and the spacing between the discrete heat sources are studied as parameters. The location in the heat pipe at which dryout is initiated is found to be different from that of the maximum temperature. The location where the maximum capillary pressure head is realized also changes during the transient. Axial conduction through the wall and wick are seen to play a significant role in determining the axial temperature variation.


2007 ◽  
Vol 589 ◽  
pp. 1-31 ◽  
Author(s):  
JIN ZHANG ◽  
STEPHEN J. WATSON ◽  
HARRIS WONG

Micro heat pipes have been used to cool micro electronic devices, but their heat transfer coefficients are low compared with those of conventional heat pipes. In this work, a dual-wet pipe is proposed as a model to study heat transfer in micro heat pipes. The dual-wet pipe has a long and narrow cavity of rectangular cross-section. The bottom-half of the horizontal pipe is made of a wetting material, and the top-half of a non-wetting material. A wetting liquid fills the bottom half of the cavity, while its vapour fills the rest. This configuration ensures that the liquid–vapour interface is pinned at the contact line. As one end of the pipe is heated, the liquid evaporates and increases the vapour pressure. The higher pressure drives the vapour to the cold end where the vapour condenses and releases the latent heat. The condensate moves along the bottom half of the pipe back to the hot end to complete the cycle. We solve the steady-flow problem assuming a small imposed temperature difference between the two ends of the pipe. This leads to skew-symmetric fluid flow and temperature distribution along the pipe so that we only need to focus on the evaporative half of the pipe. Since the pipe is slender, the axial flow gradients are much smaller than the cross-stream gradients. Thus, we can treat the evaporative flow in a cross-sectional plane as two-dimensional. This evaporative motion is governed by two dimensionless parameters: an evaporation number E defined as the ratio of the evaporative heat flux at the interface to the conductive heat flux in the liquid, and a Marangoni number M. The motion is solved in the limit E→∞ and M→∞. It is found that evaporation occurs mainly near the contact line in a small region of size E−1W, where W is the half-width of the pipe. The non-dimensional evaporation rate Q* ~ E−1 ln E as determined by matched asymptotic expansions. We use this result to derive analytical solutions for the temperature distribution Tp and vapour and liquid flows along the pipe. The solutions depend on three dimensionless parameters: the heat-pipe number H, which is the ratio of heat transfer by vapour flow to that by conduction in the pipe wall and liquid, the ratio R of viscous resistance of vapour flow to interfacial evaporation resistance, and the aspect ratio S. If HR≫1, a thermal boundary layer appears near the pipe end, the width of which scales as (HR)−1/2L, where L is the half-length of the pipe. A similar boundary layer exists at the cold end. Outside the boundary layers, Tp varies linearly with a gradual slope. Thus, these regions correspond to the evaporative, adiabatic and condensing regions commonly observed in conventional heat pipes. This is the first time that the distinct regions have been captured by a single solution, without prior assumptions of their existence. If HR ~ 1 or less, then Tp is linear almost everywhere. This is the case found in most micro-heat-pipe experiments. Our analysis of the dual-wet pipe provides an explanation for the comparatively low effective thermal conductivity in micro heat pipes, and points to ways of improving their heat transfer capabilities.


2021 ◽  
Vol 321 ◽  
pp. 04010
Author(s):  
Ramazan Aykut Sezmen ◽  
Barbaros Çetin ◽  
Zafer Dursunkaya

Heat pipes are phase change heat transfer devices used in wide range of heat transport applications due to their high thermal transport capacities with low temperature differences. Heat pipes are especially preferred for electronic cooling applications and aerospace avionics to satisfy high heat transfer rate requirements. In this study, heat transfer and phase change mechanisms of working fluid are investigated and modeled using a 3-D thermal resistance network for multichannel flat grooved heat pipes. First, heat transfer and fluid flow are modeled in half of a single grooved structure due to symmetry, and is subjected to uniform heat flux. Radius of meniscus curvature and temperature distribution along the groove are calculated. Results are compared with experiments in the literature and show good agreement. The validated heat transfer and fluid flow models are extended to a multichannel model to observe performance of grooved heat pipes with localized heat sources, not covering the entire width, a vital feature for realistic simulation of operational devices. Predictions of the temperature distribution along the multichannel of the heat pipe are provided and the effect of the distribution of heat sources on the heat pipe is discussed.


2021 ◽  
Vol 2119 (1) ◽  
pp. 012088
Author(s):  
A. A. Litvintceva ◽  
N. I. Volkov ◽  
N. I. Vorogushina ◽  
V. A. Moskovskikh ◽  
V. V. Cheverda

Abstract Heat pipes are a good solution for temperature stabilization, for example, of microelectronics, because these kinds of systems are without any moving parts. Experimental research of the effect of operating parameters on the heat transfer in a cylindrical heat pipe has been conducted. The effect of the working fluid properties and the porous layer thickness on the heat flux and temperature difference in the heat pipe has been investigated. The temperature field of the heat pipe has been investigated using the IR-camera and K-type thermocouples. The data obtained by IR-camera and K-type thermocouples have been compared. It is demonstrated the power transferred from the evaporator to the condenser is a linear function of the temperature difference between them.


1990 ◽  
Vol 112 (3) ◽  
pp. 602-607 ◽  
Author(s):  
B. R. Babin ◽  
G. P. Peterson

A computer model was developed to aid in the design of a flexible bellows heat pipe for cooling small discrete heat sources or arrays of small heat sources. This model was used to evaluate the operational characteristics and performance limitations of these heat pipes and to formulate and optimize a series of conceptual designs. Three flexible bellows heat pipes approximately 40 mm in length and 6 mm in diameter were constructed and tested using three different wick configurations. The test pipes were found to be boiling limited over most of the operating temperature range tested. Heat fluxes in excess of 200 W/cm2 were obtained and thermal resistance values of less than 0.7 °C/W were measured. Although the computer model slightly underestimated the experimentally determined transport limit for one of the wicking configurations, the remaining transport predictions were consistently within 8 percent of the experimental values.


2013 ◽  
Vol 597 ◽  
pp. 3-8
Author(s):  
Lahoucine Belarche ◽  
Btissam Abourida ◽  
Slawomir Smolen ◽  
Touria Mediouni

Natural convection in inclined cubic cavity, discretely heated, is studied numerically using a three-dimensional finite volume formulation. Two heating square portions are placed on the vertical wall of the enclosure, while the rest of the considered wall is adiabatic. These sections, similar to the integrated electronic components, generate a heat flux q". The opposite vertical wall is maintained at a cold uniform temperature Tc and the other walls are adiabatic. The fluid flow and heat transfer in the cavity are studied for different sets of the governing parameters, namely the Rayleigh number Ra (103 ≤ Ra ≤ 107), the cavity inclination γ (- 45° ≤ γ ≤ 45°) and the position of the heating sections λ (0.3 ≤ λ ≤ 0.7). The dimensions of the heater sections, ε = D / H and the longitudinal aspect ratio of the cavity Ax = H / L are respectively fixed to 0.35 and 1.


Author(s):  
Hongmin Li ◽  
Minel J. Braun ◽  
G.-X. Wang ◽  
Edward A. Evans

Hydrothermal growth is the industry method of preference to obtain high quality single crystals. Due to the high pressure and high temperature growth conditions, growth process is carried out in closed containers. During a growth run, the only flow and heat transfer that control crystal growers have is the outside heating. An inverse algorithm, used to obtain the heating distribution for an autoclave with a two-dimensional flow, is further developed and used to determine the heating distribution for an industry autoclave with three-dimensional flows. A cross-section area average temperature distribution is set as a target. With the three steps, including CFD simulation of the fluid flow, heat conduction in the metal wall, and heat conduction in the insulation layer, the heater heat flux distribution is determined. The distributions appear close to linear from the median height to the top/bottom with small magnitude deviation in the circumferential direction. Linearly distributed heaters, based on the determined heat flux distribution, are then used and heat transfer and fluid flow is numerically simulated with a conjugate model. The achieved temperature agrees well with the targeted one. The distribution and heating rates of linearly distributed heaters can be applied to industry autoclaves.


1991 ◽  
Vol 113 (3) ◽  
pp. 728-734 ◽  
Author(s):  
A. Faghri ◽  
M. Buchko

A numerical analysis and experimental verification of the effects of heat load distribution on the vapor temperature, wall temperature, and the heat transport capacity for heat pipes with multiple heat sources is presented. A numerical solution of the elliptic conjugate mass, momentum and energy equations in conjunction with the thermodynamic equilibrium relations and appropriate boundary conditions for the vapor region, wick structure, and the heat pipe wall are given. The experimental testing of a copper-water heat pipe with multiple heat sources was also made showing excellent agreement with the numerical results. An optimization of the heat distribution for such heat pipes was performed and it was concluded that by redistribution of the heat load, the heat capacity can be increased.


Author(s):  
Chemseddine Maatki ◽  
Kaouther Ghachem ◽  
Mohammed A. Almeshaal ◽  
Nidhal Ben Khedher ◽  
Lioua Kolsi

Abstract The Three-dimensional natural convection with isothermal discrete heat sources in a cubical cavity has been carefully studied using the 3D vector potential-vorticity formulation. Based on the finite volume method, the governing equations are solved with a home-made computational code (written in Fortran). Assuming that all cavity vertical walls are adiabatic, the upper wall of the cavity is kept at a cold temperature. However, in the bottom face, heat sources are placed under different configurations. The size of the discrete sources, their positions, and their numbers are varied for different Rayleigh numbers. The Prandtl number is fixed at 0.71. Three-dimensional distribution of the temperature iso-surfaces, the heat transfer rate, and entropy generation are evaluated. It is found that heat transfer and entropy generation are strongly affected by the arrangement of the discrete heated sources. In conclusion, the heat transfer rate is maximized, and the entropy generation is minimized for the inline arrangement of more than two heaters compared to the diagonal one.


Sign in / Sign up

Export Citation Format

Share Document