Measurements and high-speed visualizations of flow boiling of a dielectric fluid in a silicon microchannel heat sink

Experiments were conducted to investigate flow boiling heat transfer to a dielectric fluid in a silicon chip-integrated microchannel heat sink. Twenty-four microchannels, each 389 μm × 389 μm in cross-section, were fabricated into the 12.7 mm × 12.7 mm silicon substrate. High-speed visualizations (at 12,500 frames per second) were performed simultaneously with heat transfer and pressure drop measurements to investigate the physics of flow boiling in parallel microchannel arrays. At low heat fluxes, bubbly flow is dominant, with the bubbles coalescing to form vapor slugs as the heat flux is increased. At high heat fluxes, the flow regimes in the downstream portion of the microchannels are characteristic of alternating wispy-annular flow and churn flow, while flow reversal is observed in the upstream region near the microchannel inlet. Local heat transfer measurements, obtained at three flow rates ranging from 35 to 60 ml/min, show that at lower heat fluxes, the heat transfer coefficient increases with increasing heat flux. The heat transfer coefficient in fully developed boiling is seen to be independent of flow rate in this range. At higher heat fluxes (exceeding 542, 673, 730 kW/m², respectively, for flow rates of 35, 47 and 60 ml/min), this trend is reversed, and the heat transfer coefficient decreases with further increases in heat flux due to partial dryout in some of the microchannels. Heat fluxes at which fully developed boiling is achieved depend on the flow rate. The pressure drop in fully developed boiling increases with increasing heat flux and is independent of flow rate for the test conditions considered in this work.     

Freon R141b flow boiling in silicon microchannel heat sinks: experimental investigation

This paper presents experimental investigations on Freon R141b flow boiling in rectangular microchannel heat sinks. The main aim is to provide an appropriate working fluid for microchannel flow boiling to meet the cooling demand of high power electronic devices. The microchannel heat sink used in this work contains 50 parallel channels, with a 60 × 200 (W × H) μm cross-section. The flow boiling heat transfer experiments are performed with R141b over mass velocities ranging from 400 to 980 kg/(m² s) and heat flux from 40 to 700 kW/m², and the outlet pressure satisfying the atmospheric condition. The fluid flow-rate, fluid inlet/outlet temperature, wall temperature, and pressure drop are measured. The results indicate that the mean heat transfer coefficient of R141b flow boiling in present microchannel heat sinks depends heavily on mass velocity and heat flux and can be predicted by Kandlikar’s correlation (Heat Transf Eng 25(3):86–93, 2004). The two-phase pressure drop keeps increasing as mass velocity and exit vapor quality rise.     

Experimental investigation of flow boiling characteristics in a cross-linked microchannel heat sink

This paper experimentally investigates flow boiling characteristics in a cross-linked microchannel heat sink at low mass fluxes and high heat fluxes. The heat sink consists of 45 straight microchannels each with a hydraulic diameter of 248 μm and heated length of 16 mm. Three cross-links, of width 500 μm, are introduced in the present microchannel heat sink to achieve better temperature uniformity and to avoid flow mal-distribution. Flow visualization, flow instability, two-phase pressure drop, and two-phase heat transfer measurements are conducted using the dielectric coolant FC-72 over a range of heat flux from 7.2 to 104.2 kW/m², mass flux from 99 to 290 kg/m² s, and exit quality from 0.01 to 0.71. Thermochromic liquid crystals are used in the present study as full-field surface temperature sensors to map the temperature distribution on the heat sink surface. Flow visualization studies indicate that the observed flow regime is primarily slug. Visual observations of flow patterns in the cross-links demonstrate that bubbles nucleate and grow rapidly on the surface of the cross-links and in the tangential direction at the microchannels’ entrance due to the effect of circulations generated in those regions. The two-phase pressure drop strongly increases with the exit quality, at x_e,o < 0.3, and the two-phase frictional pressure drop increases by a factor of 1.6–2 compared to the straight microchannel heat sink. The flow boiling heat transfer coefficient increases with increasing exit quality at a constant mass flux, which is caused by the dominance of the nucleation boiling mechanism in the cross-link region.     

The critical role of channel cross-sectional area in microchannel flow boiling heat transfer

Experiments are conducted with a perfluorinated dielectric fluid, Fluorinert FC-77, to identify the critical geometric parameters that affect flow boiling heat transfer and flow patterns in microchannels. In recent work by the authors (Harirchian and Garimella, 2009), seven different silicon test pieces containing parallel microchannels of widths ranging from 100 to 5850 μm, all with a depth of 400 μm were tested and it was shown that for a fixed channel depth, the heat transfer coefficient was independent of channel width for microchannels of widths 400 μm and larger, with the flow regimes in these microchannels being similar; nucleate boiling was also found to be dominant over a wide range of heat fluxes. In the present study, experiments are performed with five additional microchannel test pieces with channel depths of 100 and 250 μm and widths ranging from 100 to 1000 μm. Flow visualizations are performed using a high-speed digital video camera to determine the flow regimes, with simultaneous local measurements of the heat transfer coefficient and pressure drop. The aim of the present study is to investigate as independent parameters the channel width and depth as well as the aspect ratio and cross-sectional area on boiling heat transfer in microchannels, based on an expanded database of experimental results. The flow visualizations and heat transfer results show that the channel cross-sectional area is the important governing parameter determining boiling mechanisms and heat transfer in microchannels. For channels with cross-sectional area exceeding a specific value, nucleate boiling is the dominant mechanism and the boiling heat transfer coefficient is independent of channel dimensions; below this threshold value of cross-sectional area, vapor confinement is observed in all channels at all heat fluxes, and the heat transfer rate increases as the microchannel cross-sectional area decreases before premature dryout occurs due to channel confinement.     

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

Visualisation of flow boiling heat transfer in a microtube

This paper presents the results of the flow boiling patterns and heat transfer coefficients of FC-72 in a small tube. The internal diameter of the tube is 0.48 mm, with a heated length of 73 mm. The mass flow rate varies from 50 to 3,000 kg/m²-s. The microtube is made of Pyrex in order to obtain the visualisation of the flow pattern along the heated channel. Different types of flow pattern have been observed: bubbly flow, deformed bubbly flow, bubbly/slug flow, slug flow, slug/annular flow, and annular flow. The experiments show the presence of flow instabilities in a large portion of the tests at low mass flow rates and low subcooling. Flow patterns in presence of flow instabilities are mainly characterized by bubbly/slug flow and slug/annular flow. Heat transfer rates have been studied in all flow pattern conditions. The two groups of results, with flow instabilities and without flow instabilities, show similar heat transfer behaviour. The heat transfer characteristics of the pipes have been studied in comparison with mass flux and vapour quality.     

An experimental investigation of flow boiling in an asymmetrically heated rectangular microchannel

By using unique experimental techniques and carefully constructed experimental apparatus, the characteristics of flow boiling of water in microscale were investigated using a single horizontal rectangular microchannel. A polydimethylsiloxane rectangular microchannel (D_h = 103.5 and 133 μm) was fabricated by using the replica molding technique, a kind of soft lithography. A piecewise serpentine platinum microheater array on a Pyrex substrate was fabricated with the surface micromachining MEMS technique. Real time flow visualization of the phase change phenomena inside the microchannel was performed using a high speed CCD camera with microscope. The experimental local boiling heat transfer coefficients were studied, and single bubble inception, growth, and departure, as well as elongated bubble behavior were analyzed to elucidate the microscale heat transfer mechanisms. Tests were performed for mass fluxes of 77.5, 154.9, and 309.8 kg/m² s and heat fluxes of 180–500 kW/m². The effects of mass flux, heat flux, and vapor qualities on flow boiling heat transfer in a microchannel were studied.     

Effects of heat flux,mass flux,vapor quality,and saturation temperature on flow boiling heat transfer in microchannels

Flow boiling heat transfer with the refrigerants R-134a and R-245fa in copper microchannel cold plate evaporators is investigated. Arrays of microchannels of hydraulic diameter 1.09 and 0.54 mm are considered. The aspect ratio of the rectangular cross section of the channels in both test sections is 2.5. The heat transfer coefficient is measured as a function of local thermodynamic vapor quality in the range −0.2 to 0.9, at saturation temperatures ranging from 8 to 30 °C, mass flux from 20 to 350 kg m⁻² s⁻¹, and heat flux from 0 to 22 W cm⁻². The heat transfer coefficient is found to vary significantly with heat flux and vapor quality, but only slightly with saturation pressure and mass flux for the range of values investigated. It was found that nucleate boiling dominates the heat transfer. In addition to discussing measurement results, several flow boiling heat transfer correlations are also assessed for applicability to the present experiments.     

Flow boiling heat transfer and pressure drop of R-134a in a mini tube: an experimental investigation

This paper presents the results of an experimental study carried out with R-134a during flow boiling in a horizontal tube of 2.6 mm ID. The experimental tests included (i) heat fluxes in the range from 10 to 100 kW/m², (ii) the refrigerant mass velocities set to the discrete values in the range of 240-930 kg/(m² s) and (iii) saturation temperature of 12 and 22 °C. The study analyzed the heat transfer, through the local heat transfer coefficient along of flow, and pressure drop, under the variation of these different parameters. It was possible to observe the significant influence of heat flux in the heat transfer coefficient and mass velocity in the pressure drop, besides the effects of saturation temperature. In the low quality region, it was possible to observe a significant influence of heat flux on the heat transfer coefficient. In the high vapor quality region, for high mass velocities, this influence tended to vanish, and the coefficient decreased. The influence of mass velocity in the heat transfer coefficient was detected in most tests for a threshold value of vapor quality, which was higher as the heat flux increased. For higher heat flux the heat transfer coefficient was nearly independent of mass velocity. The frictional pressure drop increased with the increase in vapor quality and mass velocity. Predictive models for heat transfer coefficient in mini channels were evaluated and the calculated coefficient agreed well with measured data within a range 35% for saturation temperature of 22 °C. These results extend the ranges of heat fluxes and mass velocities beyond values available in literature, and add a substantial contribution to the comprehension of boiling heat transfer phenomena inside mini channels.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

Experimental study on the heat transfer coefficient of water flow boiling in mini/microchannels

The heat transfer coefficients of the evaporative water flow in mini/microchannels are studied experimentally to explore the novel heat dissipation for high power electronics. Two sets of parallel channels which are 61 channels with hydraulic diameter of 0.293 mm and 20 channels with hydraulic diameter of 1.2 mm are investigated respectively. The inlet and outlet temperatures of fluids, and the temperatures beneath the channels are measured to calculate the heat dissipation of the evaporative water in channels. The experiments are carried out with the mass flow rates range from 11.09 kg/(m² s) to 44.36 kg/(m² s) for minichannels and 49.59 kg/(m² s) to 198.37 kg/(m² s) for microchannels. The effective heat flux range from 5 W/cm² to 50 W/cm², and the resulted outlet vapor qualities range from 0 to 0.8. The relations of the heat transfer coefficient with heat flux and vapor quality are analyzed according to the results. The experimental heat transfer coefficients are compared with the prediction of latest developed correlations. A new correlation takes the effect of Bond number is proposed, and be verified that it is effective to predict the heat transfer coefficient of both minichannels and microchannels in a large range of vapor qualities.     

Two-phase flow of R-134a refrigerant during flow boiling through a horizontal circular mini-channel

This research focuses on heat transfer to R-134a during flow boiling in a 1.75 mm internal diameter tube. Flow visualisation and heat transfer experiments are conducted to obtain heat transfer coefficients for different flow patterns. The measured data in each flow regime are compared with predictions from a three-zone flow boiling model. The calculations are in fair agreement with the experimental results which correspond in particular to slug flow, throat-annular flow and churn flow regimes under conditions of low heat flux.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

Two-phase flow of refrigerants in 85 μm-wide multi-microchannels: Part II – Heat transfer with 35 local heaters

This article is the second part of a study on flow boiling of R236fa and R245fa. This part presents the heat transfer coefficients obtained in a 12.7 mm silicon evaporator composed of 135 microchannels with 85 μm wide and 560 μm high channels separated by 46 μm wide fins. There were 35 local heaters and temperature measurements arranged in a 5 × 7 array. The heat transfer results were uniform in the lateral direction to the flow (attributable to the inlet restriction) and a function of the heat flux, vapor quality and mass flux. The steady-state standard deviation of the local base temperature was less than 0.2 °C, inferring that the boiling process was very stable. For wall heat fluxes over 45 kW/m², the heat transfer coefficient curves were V-shaped, decreasing for intermittent flow regimes and increasing for annular flow. The three-zone model of Thome et al. (2004) was the best heat transfer prediction method when setting the dryout thickness equal to the channel roughness.     

Boiling heat transfer,pressure drop,and flow pattern in a horizontal square minichannel

This study experimentally investigated the flow boiling heat transfer, pressure drop, and flow pattern in a horizontal square minichannel with a hydraulic diameter of 2.0 mm, and the effects of mass flux, vapor quality, heat flux, and refrigerant properties on the flow boiling characteristics were clarified. The heat transfer coefficient and pressure drop of R32 and R1234yf were measured in a mass flux range of 50–400 kgm⁻²s⁻¹ at a saturation temperature of 15 °C. The flow pattern of the square minichannel outlet was observed and was classified as plug, wavy, churn, and annular flows. The heat transfer coefficients in the square minichannel were larger than those in the circular minichannel with a similar hydraulic diameter at low mass flux conditions. The heat transfer coefficients of R32 indicated higher values compared with those of R1234yf at same mass flux and qualities. An empirical heat transfer model taking into account the forced convection, nucleate boiling, and thin liquid film evaporation was developed for horizontal square and circular minichannels. The frictional pressure drop of R32 was 1.5–2 times higher than that of R1234yf at same mass flux and vapor quality condition, and the effect of channel shape on the frictional pressure drop was small unlike the boiling heat transfer.     

Adiabatic and diabatic two-phase venting flow in a microchannel

The growth and advection of the vapor phase in two-phase microchannel heat exchangers increase the system pressure and cause flow instabilities. One solution is to locally vent the vapor formed by capping the microchannels with a porous, hydrophobic membrane. In this paper we visualize this venting process in a single 124 μm by 98 μm copper microchannel with a 65 μm thick, 220 nm pore diameter hydrophobic Teflon membrane wall to determine the impact of varying flow conditions on the flow structures and venting process during adiabatic and diabatic operation. We characterize liquid velocities of 0.14, 0.36 and 0.65 m/s with superficial air velocities varying from 0.3 to 8 m/s. Wavy-stratified and stratified flow dominated low liquid velocities while annular type flows dominated at the higher velocities. Gas/vapor venting can be improved by increasing the venting area, increasing the trans-membrane pressure or using thinner, high permeability membranes. Diabatic experiments with mass flux velocities of 140 and 340 kg/s/m² and exit qualities up to 20% found that stratified type flows dominate at lower mass fluxes while churn-annular flow became more prevalent at the higher mass-flux and quality. The diabatic flow regimes are believed to significantly influence the pressure-drop and heat transfer coefficient in vapor venting heat exchangers.     

Flow boiling heat transfer characteristics of R123 and R134a in a micro-channel

Flow boiling heat transfer in a single circular micro-channel of 0.19 mm ID has been experimentally investigated with R123 and R134a for various experimental conditions: heat fluxes (10, 15, 20 kW/m²), mass velocities (314, 392, 470 kg/m² s), vapor qualities (0.2–0.85) and different saturation pressures (158, 208 kPa for R123; 900, 1100 kPa for R134a). The heat transfer trends between R123 and R134a are clearly distinguished. Whether nucleate boiling is suppressed at low vapor quality or not determines the heat transfer trend and mechanism in the flow boiling of micro-channels. High convective heat transfer coefficients in the two-phase flow of micro-channels enables nucleate boiling to be suppressed even at low vapor quality, depending on the wall superheat requirement for nucleate boiling. In the case of early suppression of nucleate boiling, specifically R123, heat transfer is dominated by evaporation of thin liquid films around elongated bubbles. In the contrary case, namely R134a, nucleate boiling is dominant heat transfer mechanism until its suppression at high vapor quality and then two-phase forced convection heat transfer becomes dominant. It is similar to the heat transfer characteristic of macro-channels except the enhancement of nucleate boiling and the short forced convection region.     

Two-phase pressure drop and flow visualization of FC-72 in a silicon microchannel heat sink

The rapid development of two-phase microfluidic devices has triggered the demand for a detailed understanding of the flow characteristics inside microchannel heat sinks to advance the cooling process of micro-electronics. The present study focuses on the experimental investigation of pressure drop characteristics and flow visualization of a two-phase flow in a silicon microchannel heat sink. The microchannel heat sink consists of a rectangular silicon chip in which 45 rectangular microchannels were chemically etched with a depth of 276 μm, width of 225 μm, and a length of 16 mm. Experiments are carried out for mass fluxes ranging from 341 to 531 kg/m² s and heat fluxes from 60.4 to 130.6 kW/m² using FC-72 as the working fluid. Bubble growth and flow regimes are observed using high speed visualization. Three major flow regimes are identified: bubbly, slug, and annular. The frictional two-phase pressure drop increases with exit quality for a constant mass flux. An assessment of various pressure drop correlations reported in the literature is conducted for validation. A new general correlation is developed to predict the two-phase pressure drop in microchannel heat sinks for five different refrigerants. The experimental pressure drops for laminar-liquid laminar-vapor and laminar-liquid turbulent-vapor flow conditions are predicted by the new correlation with mean absolute errors of 10.4% and 14.5%, respectively.     

Stratified flow boiling heat transfer study of a HFC/HC refrigerant mixture in smooth horizontal tubes

The flow boiling heat transfer coefficients of R-134a/R-290/R-600a (91%:4.068%:4.932% by mass) refrigerant mixture are experimentally arrived in two tubes of diameter 9.52 and 12.7 mm. The tests are conducted to target the varied heat flux condition and stratified flow pattern found in evaporators of refrigerators and deep freezers. The varied heat flux condition is imposed on the refrigerant using a coaxial counter-current heat exchanger test section. The experiments are performed for mass flow rates of the refrigerant mixture between 3 and 5 g s⁻¹ and entry temperature between −8.59 and 5.33°C which are bubble temperatures corresponding to a pressure of 3.2 and 5 bar. The influences of heat flux, mass flow rate, pressure, flow pattern, tube diameter on the heat transfer coefficient are discussed. The profound effects of nucleate boiling prevailing even at higher vapor qualities in evaporators are highlighted. The heat transfer coefficient of the refrigerant mixture is also compared with that of R-134a.     

Flow boiling heat transfer in a vertical spirally internally ribbed tube

 Experiments of flow boiling heat transfer and two-phase flow frictional pressure drop in a spirally internally ribbed tube (φ22×5.5 mm) and a smooth tube (φ19×2 mm) were conducted, respectively, under the condition of 6×10⁵ Pa (absolute atmosphere pressure). The available heated length of the test sections was 2500 mm. The mass fluxes were selected, respectively, at 410, 610 and 810 kg/m² s. The maximum heat flux was controlled according to exit quality, which was no more than 0.3 in each test run. The experimental results in the spirally internally ribbed tube were compared with that in the smooth tube. It shows that flow boiling heat transfer coefficients in the spirally internally ribbed tube are 1.4–2 times that in the smooth tube, and the flow boiling heat transfer under the condition of smaller temperature differences can be achieved in the spirally internally ribbed tube. Also, the two-phase flow frictional pressure drop in the spirally internally ribbed tube increases a factor of 1.6–2 as compared with that in the smooth tube. The effects of mass flux and pressure on the flow boiling heat transfer were presented. The effect of diameters on flow boiling heat transfer in smooth tubes was analyzed. Based on the fits of the experimental data, correlations of flow boiling heat transfer coefficient and two-phase flow frictional factor were proposed, respectively. The mechanisms of enhanced flow boiling heat transfer in the spirally internally ribbed tube were analyzed. Received on 1 December 1999