Flow boiling heat transfer in two-phase micro-channel heat sinks--I. Experimental investigation and assessment of correlation methods

This paper is the first of a two-part study concerning measurement and prediction of saturated flow boiling heat transfer in a water-cooled micro-channel heat sink. In this paper, new experimental results are discussed which provide new physical insight into the unique nature of flow boiling in narrow rectangular micro-channels. The micro-channel heat sink contained 21 parallel channels having a m cross-section. Tests were performed with deionized water over a mass velocity range of 135-402 kg/m² s, inlet temperatures of 30 and 60 °C, and an outlet pressure of 1.17 bar. Results indicate an abrupt transition to annular flow near the point of zero thermodynamic equilibrium quality, and reveal the dominant heat transfer mechanism is forced convective boiling corresponding to annular flow. Contrary to macro-channel trends, the heat transfer coefficient is shown to decrease with increasing thermodynamic equilibrium quality. This unique trend is attributed to appreciable droplet entrainment at the onset of annular flow regime development, and the increase in mass flow rate of the annular film by droplet deposition downstream. Eleven previous empirical correlations are assessed and deemed unable to predict the correct trend of heat transfer coefficient with quality because of the unique nature of flow boiling in micro-channels, and the operating conditions of water-cooled micro-channel heat sinks falling outside the recommended application range for most correlations. Part II of this study will introduce a new annular flow model as an alternative approach to heat transfer coefficient prediction for micro-channels.     

Single-phase and two-phase heat transfer characteristics of low temperature hybrid micro-channel/micro-jet impingement cooling module

This study examines the single-phase and two-phase cooling performance of a hybrid micro-channel/micro-jet impingement cooling scheme using HFE 7100 as working fluid. This scheme consists of supplying coolant from a series of jets that deposit liquid into the micro-channels. A single-phase numerical scheme that utilizes the k–ε turbulent model and a method for determining the extent of the laminarized wall layer shows very good predictions of measured wall temperatures. It is shown jet velocity has a profound influence on single-phase cooling performance. High jet velocities enable jet fluid to penetrate the axial micro-channel flow and produce a strong impingement effect at the wall. On the other hand, the influence of jets at low jet velocities is greatly compromised compared to the micro-channel flow. During nucleate boiling, vapor layer development along the micro-channel in the hybrid module is fundamentally different from that encountered in conventional micro-channels. Here, subcooled jet fluid produces repeated regions of bubble growth followed by bubble collapse, rather than the continuous growth common to conventional micro-channel flow. By reducing void fraction along the micro-channel, the hybrid scheme contributes greater wall temperature uniformity. Increasing subcooling and/or flow rate delay the onset of boiling to higher heat fluxes and higher wall temperatures, and also increase critical heat flux considerably. A nucleate boiling heat transfer coefficient correlation is developed that fits the present data with a mean absolute error of 6.10%.     

Consolidated method to predicting pressure drop and heat transfer coefficient for both subcooled and saturated flow boiling in micro-channel heat sinks

Published studies concerning transport phenomena in micro-channel heat sinks can be divided into those concerning saturated boiling versus those focused on subcooled boiling, with the vast majority related to the former. What has been lacking is a single generalized method to tackle both boiling regimes. The primary objective of the present paper is to construct a consolidated method to predicting transport behavior of micro-channel heat sinks incurring all possible heat transfer regimes. First, a new correlation is developed for subcooled flow boiling pressure drop that accounts for inlet subcooling, micro-channel aspect ratio, and length-to-diameter ratio. This correlation shows excellent predictive capability against subcooled HFE 7100 pressure drop data corresponding to four different micro-channel geometries. Next, a consolidated method is developed for pressure drop that is capable of tackling inlet single-phase liquid, subcooled boiling, saturated boiling, and single-phase vapor regimes as well as inlet contraction and outlet expansion. A similar consolidated method is developed to predict the heat transfer coefficient that is capable of tackling all possible combinations of heat transfer regimes. The new consolidated method is shown to be highly effective at reproducing both data and trends for HFE 7100, water and R134a.     

Flow condensation in parallel micro-channels – Part 2: Heat transfer results and correlation technique

This second part of a two-part study concerns heat transfer characteristics for FC-72 condensing along parallel, square micro-channels with a hydraulic diameter of 1 mm, which were formed in the top surface of a solid copper plate. Heat from the condensing flow was rejected to a counter flow of water through channels brazed to the underside of the copper plate. The FC-72 condensation heat transfer coefficient was highest near the channel inlet, where the annual liquid film is thinnest. The heat transfer coefficient decreased along the micro-channel because of the film thickening and eventual collapse of the annular regime. Notable heat transfer enhancement was observed for annular flow regions of the micro-channel associated with interfacial waves. Comparing the present data to predictions of previous annular condensation heat transfer correlations shows correlations intended for macro-channels generally provide better predictions than correlations intended specifically for mini/micro-channels. A new condensation heat transfer coefficient correlation is proposed for annular condensation heat transfer in mini/micro-channels. The new correlation shows excellent predictive capability based on both the present FC-72 data and a large database for mini/micro-channel flows amassed from eight previous sources.     

Critical heat flux for subcooled flow boiling in micro-channel heat sinks

Critical heat flux (CHF) was measured and examined with high-speed video for subcooled flow boiling in micro-channel heat sinks using HFE 7100 as working fluid. High subcooling was achieved by pre-cooling the working fluid using a secondary low-temperature refrigeration system. The high subcooling greatly reduced both bubble departure diameter and void fraction, and precluded flow pattern transitions beyond the bubbly regime. CHF was triggered by vapor blanket formation along the micro-channel walls despite the presence of abundant core liquid, which is consistent with the mechanism of Departure from Nucleate Boiling (DNB). CHF increased with increasing mass velocity and/or subcooling and decreasing hydraulic diameter for a given total mass flow rate. A pre-mature type of CHF was caused by vapor backflow into the heat sink’s inlet plenum at low mass velocities and small inlet subcoolings, and was associated with significant fluctuations in inlet and outlet pressure, as well as wall temperature. A systematic technique is developed to modify existing CHF correlations to more accurately account for features unique to micro-channel heat sinks, including rectangular cross-section, three-sided heating, and flow interaction between micro-channels. This technique is shown to be successful at correlating micro-channel heat sink data corresponding to different hydraulic diameters, mass velocities and inlet temperatures.     

垂直窄缝通道内环状流沸腾传热特性

本研究基于液膜和蒸汽的质量、动量和能量方程,建立了均匀热流垂直窄缝通道内环状流沸腾传热模型,通过相关文献估算环状流起始点处液膜厚度,利用有限差分法对环状流模型方程组进行数值求解,得到沿流道环状流区域的液膜厚度,并进一步预测了局部沸腾传热系数,结果表明:环状流区域的局部沸腾传热系数随质量流量和干度的增加而增加,与Kenning关联式对比,模型预测沸腾传热系数较关联式计算值偏低。将不同工况下的226组两相环状流实验数据与模型预测结果进行对比,平均绝对误差为18.2%。     

微通道和小通道内流动沸腾换热的研究

介绍了小通道、微通道内流动沸腾换热机理的最新进展，其中包括国内外学者对各种形状的单通道小通道内流动沸腾现象的解释及其对换热系数影响的各种不同看法，以及微通道内流动沸腾换热现象研究的概况；还介绍了小通道和微通道内流型转变的研究现状，指出了各研究者研究结果的差异，并提出要进一步探索小通道微通道内流动沸腾换热现象，必须借助更先进的观测手段等。     

Flow Pattern and Heat Transfer Behavior of Boiling Two—Phase flow in Inclined Pipes

Movable Electrical Conducting Probe (MECP), a kind of simple and reliable measuring transducer, used for predicting full-flow-path flow pattern in a boiling vapor/liquid two-phase flow is introduced in this paper. When the test pipe is set at different inclination angles, several kinds of flow patterns, such as bubble, slug, churn, intermittent, and annular flows, may be observed in accordance with the locations of MECP. By means of flow pattern analysis, flow field numerical calculations have been carried out, and heat transfer coefficient correlations along full-flow-path derived. The results show that heat transfer performance of boiling two-phase flow could be significantly augmented as expected in some flow pattern zones.The results of the investigation, measuring techniques and conclusions contained in this paper would be a useful reference in foundational research for prediction of flow pattern and heat transfer behavior in boiling two-phase flow, as well as for turbine vane liquid-cooling design.     

Saturated flow boiling heat transfer and associated bubble characteristics of R-134a in a narrow annular duct

Experiments are conducted here to investigate how the channel size affects the saturated flow boiling heat transfer and associated bubble characteristics of refrigerant R-134a in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. The measured heat transfer data indicate that the saturated flow boiling heat transfer coefficient increases with a decrease in the gap of the duct. Besides, raising the imposed heat flux can cause a significant increase in the boiling heat transfer coefficients. However, the effects of the refrigerant mass flux and saturated temperature on the boiling heat transfer coefficient are milder. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases slightly at increasing R-134a mass flux. Moreover, the bubble departure frequency increases at reducing duct size mainly due to the rising shear stress of the liquid flow, and at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Correlation for the present saturated flow boiling heat transfer data of R-134a in the narrow annular duct is proposed. Additionally, data for some quantitative bubble characteristics such as the mean bubble departure diameter and frequency and the active nucleation site density are also correlated.     

Assessment of the effectiveness of nanofluids for single-phase and two-phase heat transfer in micro-channels

Experiments were performed to explore the micro-channel cooling benefits of water-based nanofluids containing small concentrations of Al₂O₃. The high thermal conductivity of nanoparticles is shown to enhance the single-phase heat transfer coefficient, especially for laminar flow. Higher heat transfer coefficients were achieved mostly in the entrance region of micro-channels. However, the enhancement was weaker in the fully developed region, proving that nanoparticles have an appreciable effect on thermal boundary layer development. Higher concentrations also produced greater sensitivity to heat flux. Despite this enhancement, the overall cooling effectiveness of nanoparticles was quite miniscule because of the large axial temperature rise associated with the decreased specific heat for the nanofluid compared to the base fluid. For two-phase cooling, nanoparticles caused catastrophic failure by depositing into large clusters near the channel exit due to localized evaporation once boiling commenced. These and other practical disadvantages bring into question the overall merit of using nanofluids in micro-channel heat sinks.     

Effects of heat flux,vapour quality,channel hydraulic diameter on flow boiling heat transfer in variable aspect ratio micro-channels using transparent heating

Experiments on flow boiling heat transfer in high aspect ratio micro-channels with FC-72 were carried out. Three channels with different hydraulic diameters (571, 762 and 1454 μm) and aspect ratios (20, 20 and 10) were selected. The tested mass fluxes were 11.2, 22.4 and 44.8 kg m^?2 s^?1 and heat fluxes ranging from 0–18.6 kW m^?2. In the present study, boiling curves with obvious temperature overshoots are presented. Average heat transfer coefficient and local heat transfer coefficient along stream-wise direction are measured as a function of heat flux and vapour quality respectively. Slug-annular flow and annular flow are the main flow regimes. Convective boiling is found to be the dominant heat transfer mechanism. Local heat transfer coefficient increases with decreasing hydraulic diameter. Moreover, the effect of hydraulic diameter is more significant when mass flux is higher. The unique channel geometry is considered as the decisive reason of the flow regimes as well as heat transfer mechanisms.     

The effect of heating direction on flow boiling heat transfer of R134a in micro-channels

This paper presents effects of heating directions on heat transfer performance of R134a flow boiling in micro-channel heat sink.The heat sink has 30 parallel rectangular channels with cross-sectional dimensions of 500μm width 500μm depth and 30mm length.The experimental operation condition ranges of the heat flux and the mass flux were 13.48 to 82.25 W/cm2 and 373.3 to 1244.4 kg/m2s respectively.The vapor quality ranged from 0.07 to 0.93.The heat transfer coefficients of top heating and bottom heating both were up to 25 kW/m2 K.Two dominate transfer mechanisms of nucleate boiling and convection boiling were observed according to boiling curves.The experimental results indicated that the heat transfer coefficient of bottom heating was 13.9％ higher than top heating in low heat flux,while in high heat flux,the heat transfer coefficient of bottom heating was 9.9％.higher than the top heating,because bubbles were harder to divorce the heating wall.And a modified correlation was provided to predict heat transfer of top heating.     

Saturated flow boiling heat transfer of R-410A and associated bubble characteristics in a narrow annular duct

Experiments are conducted here to investigate how the channel size affects the R-410A saturated flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. The measured data indicate that the saturated flow boiling heat transfer coefficient increases with increasing refrigerant mass flux and saturated temperature and with a decrease in the gap size. Besides, raising the imposed heat flux can cause a significant increase in the boiling heat transfer coefficient. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases slightly with increasing refrigerant mass flux and saturated temperature. Moreover, the bubble departure frequency increases at reducing duct size and increasing mass flux. And at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Meanwhile, comparisons of the present heat transfer data for R-410A with R-407C and R-134a in the same duct and with some existing correlations are conducted. Furthermore, an empirical correlation for the present R-410A saturated flow boiling heat transfer data is proposed.     

Saturated flow boiling heat transfer of R-407C and associated bubble characteristics in a narrow annular duct

Experiments are conducted here to investigate how the channel size affects the R-407C saturated flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. The measured data indicate that the saturated flow boiling heat transfer coefficient increases with a decrease in the duct gap. Besides, raising the imposed heat flux can cause a significant increase in the boiling heat transfer coefficients. However, the effects of the refrigerant mass flux and saturated temperature on the boiling heat transfer coefficient are milder. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases noticeably at increasing R-407C mass flux. Moreover, the bubble departure frequency increases at reducing duct size and at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Meanwhile comparison of the present heat transfer data for R-407C with R-134a in the same duct and with some existing correlations is conducted. Furthermore, correlation for the present R-407C saturated flow boiling heat transfer data is proposed. Additionally, the present data for some quantitative bubble characteristics such as the mean bubble departure diameter and frequency and the active nucleation site density are also correlated.     

Two-phase flow in high-heat-flux micro-channel heat sink for refrigeration cooling applications: Part II—heat transfer characteristics

This paper is the second of a two-part study concerning two-phase flow and heat transfer characteristics of R134a in a micro-channel heat sink incorporated as an evaporator in a refrigeration cycle. Boiling heat transfer coefficients were measured by controlling heat flux (q″ = 15.9 − 93.8 W/cm²) and vapor quality (x_e = 0.26 − 0.87) over a broad range of mass velocity. While prior studies point to either nucleate boiling or annular film evaporation (convective flow boiling) as dominant heat transfer mechanisms in small channels, the present study shows heat transfer is associated with different mechanisms for low, medium and high qualities. Nucleate boiling occurs only at low qualities (x_e < 0.05) corresponding to very low heat fluxes, and high fluxes produce medium quality (0.05 < x_e < 0.55) or high quality (x_e > 0.55) flows dominated by annular film evaporation. Because of the large differences in heat transfer mechanism between the three quality regions, better predictions are possible by dividing the quality range into smaller ranges corresponding to these flow transitions. A new heat transfer coefficient correlation is recommended which shows excellent predictions for both R134a and water.     

单面加热竖直窄通道内环状流沸腾传热特性

基于汽芯的动量方程和液膜的质量和动量方程,建立了单面均匀热流竖直窄通道内环状流沸腾传热模型,利用数值法对方程组进行求解,得出了环状流区域的液膜厚度,并进一步预测了环状流两相沸腾传热系数。研究表明:模型预测的两相沸腾传热系数比Mahmound关联式计算值偏小;将不同工况下的291组环状流两相沸腾传热系数实验值与模型预测值进行对比,平均绝对误差为12.7%。     

Prediction and measurement of incipient boiling heat flux in micro-channel heat sinks

Experiments were performed to measure the incipient boiling heat flux in a heat sink containing 21 rectangular (231 μm wide and 713 μm deep) micro-channels. Tests were performed using deionized water with inlet liquid velocities of 0.13-1.44 m/s, inlet temperatures of 30, 60 and 90 °C, and an outlet pressure of 1.2 bar. Using a microscope, boiling incipience was identified when the first bubbles were detected growing at, and departing from the micro-channel wall near the outlet. A comprehensive model was developed to predict the incipient boiling heat flux, accounting for the complexities of bubble formation along the flat and corner regions of a rectangular flow channel, as well as the likelihood of bubbles growing sufficiently large to engulf the entire flow area of a micro-channel. The model is based on a bubble departure criterion, which combines both mechanical considerations (force balance on a bubble) and thermal considerations (superheating entire bubble interface). The model shows good agreement with the experimental results.     

Saturated flow boiling heat transfer of refrigerant R-410A in a horizontal annular finned duct

An experiment is conducted here to investigate the saturated flow boiling heat transfer characteristics of ozone friendly refrigerant R-410A in a horizontal annular finned duct. Meanwhile the associated bubble characteristics in the duct are also inspected from the flow visualization. The experimental data are presented in terms of saturated flow boiling curves, boiling heat transfer coefficients and flow photos. In addition, empirical correlation equations for the saturated flow boiling heat transfer coefficient and mean bubble departure diameter are proposed. The saturated flow boiling curves show that boiling hysteresis is insignificant in the flow and the wall superheat needed for the onset of nucleate boiling is slightly affected by the refrigerant mass flux. Besides, the boiling curves are mainly affected by the imposed heat flux and refrigerant mass flux. Moreover, the measured saturated flow boiling heat transfer coefficient increases with the imposed heat flux and refrigerant mass flux. Furthermore, at a higher refrigerant mass flux the departing bubbles are smaller.