Local heat transfer distribution and effect of instabilities during flow boiling in a silicon microchannel heat sink

Flow boiling of the perfluorinated dielectric fluid FC-77 in a silicon microchannel heat sink is investigated. The heat sink contains 60 parallel microchannels each of 100 μm width and 389 μm depth. Twenty-five evenly distributed temperature sensors in the substrate yield local heat transfer coefficients. The pressure drop across the channels is also measured. Experiments are conducted at five flow rates through the heat sink in the range of 20–80 ml/min with the inlet subcooling held at 26 K in all the tests. At each flow rate, the uniform heat input to the substrate is increased in steps so that the fluid experiences flow regimes from single-phase liquid flow to the occurrence of critical heat flux (CHF). In the upstream region of the channels, the flow develops from single-phase liquid flow at low heat fluxes to pulsating two-phase flow at high heat fluxes during flow instability that commences at a threshold heat flux in the range of 30.5–62.3 W/cm² depending on the flow rate. In the downstream region, progressive flow patterns from bubbly flow, slug flow, elongated bubbles or annular flow, alternating wispy-annular and churn flow, and wall dryout at highest heat fluxes are observed. As a result, the heat transfer coefficients in the downstream region experience substantial variations over the entire heat flux range, based on which five distinct boiling regimes are identified. In contrast, the heat transfer coefficient midway along the channels remains relatively constant over the heat flux range tested. Due to changes in flow patterns during flow instability, the heat transfer is enhanced both in the downstream region (prior to extended wall dryout) and in the upstream region. A previous study by the authors found no effect of instabilities during flow boiling in a heat sink with larger microchannels (each 300 μm wide and 389 μm deep); it appears therefore that the effect of instabilities on heat transfer is amplified in smaller-sized channels. While CHF increases with increasing flow rate, the pressure drop across the channels has only a minimal dependence on flow rate once boiling is initiated in the microchannels, and varies almost linearly with increasing heat flux.     

New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes

A new flow boiling heat transfer model and a new flow pattern map based on the flow boiling heat transfer mechanisms for horizontal tubes have been developed specifically for CO₂. Firstly, a nucleate boiling heat transfer correlation incorporating the effects of reduced pressure and heat flux at low vapor qualities has been proposed for CO₂. Secondly, a nucleate boiling heat transfer suppression factor correlation incorporating liquid film thickness and tube diameters has been proposed based on the flow boiling heat transfer mechanisms so as to capture the trends in the flow boiling heat transfer data. In addition, a dryout inception correlation has been developed. Accordingly, the heat transfer correlation in the dryout region has been modified. In the new flow pattern map, an intermittent flow to annular flow transition criterion and an annular flow to dryout region transition criterion have been proposed based on the changes in the flow boiling heat transfer trends. The flow boiling heat transfer model predicts 75.5% of all the CO₂ database within ±30%. The flow boiling heat transfer model and the flow pattern map are applicable to a wide range of conditions: tube diameters (equivalent diameters for non-circular channels) from 0.8 to 10 mm, mass velocities from 170 to 570 kg/m² s, heat fluxes from 5 to 32 kW/m² and saturation temperatures from −28 to 25 °C (reduced pressures from 0.21 to 0.87).     

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.     

High heat flux flow boiling in silicon multi-microchannels – Part I: Heat transfer characteristics of refrigerant R236fa

This article is the first in a three part study on flow boiling of refrigerants R236fa and R245fa in a silicon multi-microchannel heat sink. The heat sink was composed of 67 parallel channels, which are 223 μm wide, 680 μm high and 20 mm long with 80 μm thick fins separating the channels. The base heat flux was varied from 3.6 to 221 W/cm², the mass velocity from 281 to 1501 kg/m² s and the exit vapour quality from 2% to 75%. The working pressure and saturation temperature were set nominally at 273 kPa and 25 °C, respectively. The present database includes 1217 local heat transfer coefficient measurements, for which three different heat transfer trends were identified, but in most cases the heat transfer coefficient increased with heat flux and was almost independent of vapour quality and mass velocity. Importantly, it was found for apparently the first time that the heat transfer coefficient as a function of vapour quality reaches a maximum at very high heat fluxes and then decreases with further increase of heat flux.     

Experimental investigation and flow visualization to determine the optimum dimension range of microgap heat sinks

The rapid increase of heat flux in high performance electronic devices has necessitated the development of high capacity thermal management techniques that can support extremely high heat transfer rates. Flow boiling in microgap is very promising for this purpose due to its high heat transfer rate and ease of fabrication. However, the effects of microgap dimension on heat transfer and pressure drop characteristics along with flow visualization have not been investigated extensively. This paper focuses on flow boiling experiments of deionized water in silicon microgap heat sink for ten different microgap dimensions from a range of 80 μm–1000 μm to determine the most effective and efficient range of microgap dimensions based on heat transfer and pressure drop performance. High speed flow visualization is conducted simultaneously along with experiments to illustrate the bubble characteristics in the boiling flow in microgap. The results of this study show that confinement in flow boiling occurs for microgap sizes 500 μm and below and confined slug/annular flow is the main dominant regime whereas physical confinement does not occur for microgap sizes 700 μm and above and bubbly flow is the dominant flow regime. The microgap is ineffective below 100 μm as partial dryout strikes very early and the wall temperature is much higher for a fixed heat flux as microgap size increases above 500 μm. In addition, results show that pressure drop and pressure fluctuation decrease with the increases of gap size whereas wall temperature and wall temperature fluctuation increase with the increases of gap size. A strong dependence of heat transfer coefficient on microgap sizes is observed for microgap sizes 500 μm and below. However, the heat transfer coefficient is independent of microgap size for microgap sizes 700 μm and above.     

Subcooled flow boiling heat transfer and associated bubble characteristics of FC-72 on a heated micro-pin-finned silicon chip

Experiments are conducted here to investigate subcooled flow boiling heat transfer and associated bubble characteristics of FC-72 on a heated micro-pin-finned silicon chip flush-mounted on the bottom of a horizontal rectangular channel. In the experiments the mass flux is varied from 287 to 431 kg/m² s, coolant inlet subcooling from 2.3 to 4.3 °C, and imposed heat flux from 1 to 10 W/cm². Besides, the silicon chips contain three different geometries of micro-structures, namely, the smooth, pin-finned 200 and pin-finned 100 surfaces. The pin-finned 200 and 100 surfaces, respectively, contain micro-pin-fins of size 200 μm × 200 μm × 70 μm (width × length × height) and 100 μm × 100 μm × 70 μm. The measured data show that the subcooled flow boiling heat transfer coefficient is reduced at increasing inlet liquid subcooling but is little affected by the coolant mass flux. Besides, adding the micro-pin-fin structures to the chip surface can effectively raise the single-phase convection and flow boiling heat transfer coefficients. Moreover, the mean bubble departure diameter and active nucleation site density are reduced for rises in the FC-72 mass flux and inlet liquid subcooling. Increasing coolant mass flux or reducing inlet liquid subcooling results in a higher mean bubble departure frequency. Furthermore, larger bubble departure diameter, higher bubble departure frequency, and higher active nucleation site density are observed as the imposed heat flux is increased. Finally, empirical correlations for the present data for the heat transfer and bubble characteristics in the FC-72 subcooled flow boiling are proposed.     

Dryout characteristics during flow boiling of R134a in vertical circular minichannels

In this paper, the experimental results of dryout during flow boiling in minichannels are reported and analysed. Experiments were carried out in vertical circular minichannels with internal diameters of 1.22 mm and 1.70 mm and a fixed heated length of 220 mm. R134a was used as working fluid. Mass flux was varied from 50 kg/m² s to 600 kg/m² s and experiments were performed at two different system pressures corresponding to saturation temperatures of 27 °C and 32 °C. Experimental results show that the dryout heat flux increases with mass flux and decreases with tube diameter while system pressure has no clear effect for the range of experimental conditions covered. Finally, the prediction capabilities of the well known critical heat flux (CHF) correlations are also tested.     

Boiling two-phase flow and efficiency of co- and counter-current microchannel heat exchangers with gas heating

To learn how to utilize the exhaust heat from a high-temperature gas product of a methanol reformer, the present study experimentally investigates the boiling two-phase flow in co- and counter-current microchannel heat exchangers (MCHE) with gas heating. Boiling two-phase flow patterns, two-phase flow instability, and efficiency are explored. The working fluid on the hot and cold sides are helium and liquid methanol, respectively. The silicon-based MCHE, which has dimensions of 20 mm × 20 mm, is designed with 18 parallel microchannels on both sides and is prepared using microfabrication processes. Four types of two-phase flow patterns – bubbly-elongated slug flow, annular flow, annular flow with liquid film breakup, and dryout are identified in both types of MCHE that are studied. A flow pattern map is then constructed on the plane of the methanol mass flux versus heat flux for both types of MCHE. In the counter-current MCHE, the efficiency increases significantly with an increase in the mass flux in both the single- and two-phase flow regions, while the effect of mass flux is insignificant in the co-current MCHE. In the two-phase flow region, the efficiency of both types of MCHEs gradually increases with an increase in the hot-side thermal power until the CHF is approached. The highest efficiency obtained in the present study is about 0.85 and 0.90 for the co- and counter-current MCHEs, respectively.     

High heat flux flow boiling in silicon multi-microchannels – Part III: Saturated critical heat flux of R236fa and two-phase pressure drops

New experimental critical heat flux results for saturated boiling conditions have been obtained for R236fa flowing in a silicon multi-microchannel heat sink composed of 67 parallel channels, 223 μm wide, 680 μm high and with 80 μm thick fins separating the channels. The microchannel length was 20 mm. The footprint critical heat fluxes measured varied from 112 to 250 W/cm² and the wall critical heat fluxes from 21.9 to 52.2 W/cm² for mass velocities from 276 to 992 kg/m²s. When increasing the mass velocity, the wall critical heat flux was observed to increase. The inlet saturation temperatures (20.31 ? T_sat,in ? 34.27 °C) and the inlet subcoolings (0.4 ? Δ T_sub ? 15.3 K) were found to have a negligible influence on the saturated CHF. The best methods for predicting the data were those of Wojtan et al. [L. Wojtan, R. Revellin, J. R. Thome, Investigation of critical heat flux in single, uniformly heated microchannels, Exp. Therm. Fluid Sci. 30 (2006) 765–774] and Revellin and Thome [R. Revellin, J. R. Thome, A theoretical model for the prediction of the critical heat flux in heated microchannels, Int. J. Heat Mass Transfer 50 (in press)]. They both predict the experimental CHF results with a mean absolute error of around 9%. Using the critical vapour quality, an annular-to-dryout transition is also proposed as a limit in a diabatic microscale flow pattern map. Pressure drop measurements were measured and analysed, showing that the homogeneous model could correctly predict the observed trends.     

The importance of turbulence during condensation in a horizontal circular minichannel

Three-dimensional simulations of condensation of refrigerant R134a in a horizontal minichannel are presented. Mass fluxes ranging from 50 kg m^?2 s^?1 up to 1000 kg m^?2 s^?1 are considered in a circular minichannel of 1 mm diameter, and uniform wall and vapour–liquid interface temperatures are imposed as boundary conditions. The Volume of Fluid (VOF) method is used to track the vapour–liquid interface; the effects of interfacial shear stress, gravity and surface tension are taken into account. The influence of turbulence in the condensate film is analysed and compared against the assumption of laminar condensate flow by employing different computational approaches and validating the results against experimental data. Under the assumption of laminar condensate flow, experimental heat transfer coefficient values at low mass fluxes can be predicted, but the computed heat transfer coefficient is found to be almost independent of mass flux and vapour quality. Only when turbulence in the condensate film is taken into account does the numerical model capture the influence of mass flux that is observed in the experimental measurements.     

Enhanced flow boiling heat transfer of FC-72 on micro-pin-finned surfaces

For the purpose of cooling electronic components with high heat flux efficiently, some experiments were conducted to study the flow boiling heat transfer performance of FC-72 on silicon chips. Micro-pin-fins were fabricated on the chip surface using a dry etching technique to enhance boiling heat transfer. Three different fluid velocities (0.5, 1 and 2 m/s) and three different liquid subcoolings (15, 25 and 35 K) were performed, respectively. A smooth chip (chip S) and four micro-pin-finned chips with the same fin thickness of 30 μm and different fin heights of 60 μm (chip PF30–60) and 120 μm (chip PF30–120), respectively, were tested. All the micro-pin-finned surfaces show a considerable heat transfer enhancement compared to the smooth one, and the critical heat flux increases in the order of chip S, PF30–60 and PF30–120. For a lower ratio of fin height to fin pitch and/or higher fluid velocity, the fluid velocity has a positive effect on the nucleate boiling curves for the micro-pin-finned surfaces. At the velocities lower than 1 m/s, the micro-pin-finned surfaces show a sharp increase in heat flux with increasing wall superheat, and the wall temperature at the critical heat flux (CHF) is less than the upper limit, 85 °C, for the reliable operation of LSI chips. The CHF values for all surfaces increase with fluid velocity and subcooling. The maximum CHF can reach nearly 150 W/cm² for chip PF30–120 at the fluid velocity of 2 m/s and the liquid subcooling of 35 K.     

Wall heat flux partitioning during subcooled forced flow film boiling of water on a vertical surface

Subcooled flow film boiling experiments were conducted on a vertical flat plate, 30.5 cm in height, and 3.175 cm wide with forced convective upflow of subcooled water at atmospheric pressure. Data have been obtained for mass fluxes ranging from 0 to 700 kg/m²s, inlet subcoolings ranging from 0 to 25 °C and wall superheats ranging from 200 to 400 °C. Correlations for wall heat transfer coefficient and wall heat flux partitioning were developed as part of this work. These correlations derive their support from simultaneous measurements of the wall heat flux, fluid temperature profiles, liquid side heat flux and interfacial wave behavior during steady state flow film boiling. A new correlation for the film collapse temperature was also deduced by considering the limiting case of heat flux to the subcooled liquid being equal to the wall heat flux. The premise of this deduction is that film collapse under subcooled conditions occurs when there is no net vapor generation. These correlations have also been compared with the data and correlations available in the literature.     

Correlation of critical heat flux and two-phase friction factor for subcooled convective boiling in structured surface microchannels

Critical heat flux (CHF) and pressure drop of subcooled flow boiling are measured for a microchannel heat sink containing 75 parallel 100 μm × 200 μm structured surface channels. The heated surface is made of a Cu metal sheet with/without 2 μm thickness diamond film. Tests and measurements are conducted with de-ionized water, de-ionized water +1 vol.% MCNT additive solution, and FC-72 fluids over a mass velocity range of 820–1600 kg/m² s, with inlet temperatures of 15(8.6)°C, 25(13.6)°C, 44(24.6)°C, and 64(36.6)°C for DI water (FC-72), and heat fluxes up to 600 W/cm². The CHF of subcooled flow boiling of the test fluids in the microchannels is measured parametrically. The two-phase pressure drop is also measured. Both CHF and the two-phase friction factor correlation for one-side heating with two other side-structured surface microchannels are proposed and developed in terms of the relevant parameters.     

Unstable and stable flow boiling in parallel microchannels and in a single microchannel

A simultaneous visualization and measurement study have been carried out to investigate flow boiling instabilities of water in microchannels at various heat fluxes and mass fluxes. Two separate flow boiling experiments were conducted in eight parallel silicon microchannels (with flow interaction from neighboring channels at headers) and in a single microchannel (without flow interaction), respectively. These microchannels, at a length of 30 mm, had an identical trapezoidal cross-section with a hydraulic diameter of 186 μm. At a given heat flux and inlet water temperature, it was found that stable and unstable flow boiling regimes existed, depending on the mass flux. A flow boiling map, in terms of heat flux vs mass flux, showing stable flow boiling regime and unstable flow boiling regime is presented for parallel microchannels as well as for a single microchannel, respectively, at an inlet water temperature of 35 °C. In the stable flow boiling regime, isolated bubbles were generated and were pushed away by the incoming subcooled liquid. Two unstable flow boiling regimes, with long-period oscillation (more than 1 s) and short-period oscillation (less than 0.1 s) in temperature and pressure, were identified. The former was due to the expansion of vapor bubble from downstream while the latter was owing to the flow pattern transition from annular to mist flow. A comparison of results of flow boiling in parallel microchannels and in a single microchannel shows that flow interaction effects from neighboring channels at the headers are significant.     

Analysis of flow patterns emerging during evaporation in parallel microchannels

The evaporation processes of 2-propanol and water in cyclo olefin polymer (COP) and silicon microchannels of square cross-section are studied with a high-speed camera. The COP channels with a cross-section of 50 μm × 50 μm are rather smooth, whereas the 30 μm × 30 μm silicon channels have comparatively rough surfaces. For the COP channels, two different evaporation modes are identified, both with oscillating liquid–vapor menisci. One of these modes is characterized by an extremely rapid evaporation and a corresponding discontinuous shift of the meniscus. In the silicon channels four different evaporation modes are observed. Oscillatory motion of the liquid fronts also dominates here, and depending on the total mass flow and the wall temperature the oscillations in different channels are synchronized or desynchronized. Besides the flow patterns also the velocity trajectories of the evaporating liquid fronts are analyzed in detail and show a rather good reproducibility over different channels and different cycles. Compared to most other studies reported in this field, bubble nucleation is found to be of secondary importance for the evaporation processes.     

Infrared micro-particle image velocimetry measurements and predictions of flow distribution in a microchannel heat sink

The flow distribution in a silicon microchannel heat sink was studied using infrared micro-particle image velocimetry (IR μPIV). The microchannel test piece consisted of seventy-six 110 μm wide × 371 μm deep channels etched into a silicon substrate. Inlet and outlet manifolds, also etched into the substrate, were fed by 1.4 mm inner-diameter tubing ports. An image-processing algorithm was developed that significantly improves the quality of IR μPIV recordings in low signal-to-noise ratio environments. A general expression for the PIV measurement depth is presented, which is valid for PIV images that have undergone a threshold image-processing operation. Experiments were performed at two different flow rates: 10 ml/min (Re = 10.2) and 100 ml/min (Re = 102). Little flow maldistribution was observed at the lower flow rate. However, significant flow maldistribution was observed at Re = 102, with the channels near the centerline having an approximately 30% greater mass flux than the channels near the lateral edges of the heat sink. Numerical simulations carried out for flow in the microchannel heat sink agreed very well with the experimental measurements, validating the use of a computational approach for studying the effect of manifold design on flow distribution in microchannel heat sinks.     

Saturated flow boiling heat transfer and pressure drop in silicon microchannel arrays

Flow boiling in arrays of parallel microchannels is investigated using a silicon test piece with imbedded discrete heat sources and integrated local temperature sensors. The microchannels considered range in width from 102 μm to 997 μm, with the channel depth being nominally 400 μm in each case. Each test piece has a footprint of 1.27 cm by 1.27 cm with parallel microchannels diced into one surface. Twenty five microsensors integrated into the microchannel heat sinks allow for accurate local temperature measurements over the entire test piece. The experiments are conducted with deionized water which enters the channels in a purely liquid state. Results are presented in terms of temperatures and pressure drop as a function of imposed heat flux. The experimental results allow a critical assessment of the applicability of existing models and correlations in predicting the heat transfer rates and pressure drops in microchannel arrays, and lead to the development of models for predicting the two-phase pressure drop and saturated boiling heat transfer coefficient.     

Microchannel size effects on local flow boiling heat transfer to a dielectric fluid

Heat transfer with liquid–vapor phase change in microchannels can support very high heat fluxes for use in applications such as the thermal management of high-performance electronics. However, the effects of channel cross-sectional dimensions on the two-phase heat transfer coefficient and pressure drop have not been investigated extensively. In the present work, experiments are conducted to investigate the local flow boiling heat transfer of a dielectric fluid, Fluorinert FC-77, in microchannel heat sinks. Experiments are performed for mass fluxes ranging from 250 to 1600 kg/m² s. Seven different test pieces made from silicon and consisting of parallel microchannels with nominal widths ranging from 100 to 5850 μm, all with a nominal depth of 400 μm, are considered. An array of temperature sensors on the substrate allows for resolution of local temperatures and heat transfer coefficients. The results of this study show that for microchannels of width 400 μm and greater, the heat transfer coefficients corresponding to a fixed wall heat flux as well as the boiling curves are independent of channel size. Also, heat transfer coefficients and boiling curves are independent of mass flux in the nucleate boiling region for a fixed channel size, but are affected by mass flux as convective boiling dominates. A strong dependence of pressure drop on both channel size and mass flux is observed. The experimental results are compared to predictions from a number of existing correlations for both pool boiling and flow boiling heat transfer.     

Investigation of heat transfer and pressure drop of CO2 two-phase flow in a horizontal minichannel

An innovative cooling system based on evaporative CO₂ two-phase flow is under investigation for the tracker detectors upgrade at CERN (European Organization for Nuclear Research). The radiation hardness and the excellent thermodynamic properties emphasize carbon dioxide as a cooling agent in the foreseen minichannels. A circular stainless steel tube in horizontal orientation with an inner diameter of 1.42 mm and a length of 0.3 m has been used as a test section to perform the step-wise scanning of the vapor quality in the entire two-phase region. To characterize the heat transfer and the pressure drop depending on the vapor quality in the tube, measurements have been performed by varying the mass flux from 300 to 600 kg/m² s, the heat flux from 7.5 to 29.8 kW/m² and the saturation temperature from ?40 to 0 °C (reduced pressures from 0.136 to 0.472). Heat transfer coefficients between 4 kW/m² K and 28 kW/m² K and pressure gradients up to 75 kPa/m were registered. The measured data was analyzed corresponding to the dependencies on heat flux, mass flux and saturation temperature. A database has been established containing about 2000 measurement points. The experimental data was compared with common models recently developed by Cheng et al. [1], [2] to cross check their applicability. The overall trends and experimental data were reproduced as predicted by the models before the dryout onset, and deviations have been analyzed. A modified friction factor for the pressure drop model [1] in mist flow has been proposed based on the experimental data.     

Experimental investigation of single-phase turbulent flow of R-134a in a multiport microchannel heat sink

This experimental study aims to investigate the heat transfer characteristics of single-phase turbulent flow of R-134a refrigerant in a rectangular multi-micro channel heat sink having 27 channels where each channel has a hydraulic diameter of 421 μm. Experimental results were obtained for inlet temperatures ranging from 24 to 33 °C, mass fluxes from 1485 to 2784 kg m^− 2 s^− 1 and wall heat fluxes from 3 to 24 kW m^− 2. The results indicate that the heat transfer coefficients are found to be higher at lower inlet temperatures than those at higher ones. In addition, when equal amount of heat supplied to the heat sink, the heat transfer coefficients increase with increasing the mass flux of refrigerant. They were also compared with 12 well-known correlations and it was seen that 4 of 12 were in good agreement with each other with the average deviation < 10%. The findings demonstrate that well-known correlations in fundamental sources can be used to predict the heat transfer coefficient of R-134a during its single phase flow in a multiport microchannel heat sink under turbulent regime.