Macro-to-microchannel transition in two-phase flow: Part 1 - Two-phase flow patterns and film thickness measurements

The classification of macroscale, mesoscale and microscale channels with respect to two-phase processes is still an open question. The main objective of this study focuses on investigating the macro-to-microscale transition during flow boiling in small scale channels of three different sizes with three different refrigerants over a range of saturation conditions to investigate the effects of channel confinement on two-phase flow patterns and liquid film stratification in a single circular horizontal channel (Part 2 covers the flow boiling heat transfer and critical heat flux). This paper presents the experimental two-phase flow pattern transition data together with a top/bottom liquid film thickness comparison for refrigerants R134a, R236fa and R245fa during flow boiling in small channels of 1.03, 2.20 and 3.04 mm diameter. Based on this work, an improved flow pattern map has been proposed by determining the flow patterns transitions existing under different conditions including the transition to macroscale slug/plug flow at a confinement number of Co ≈ 0.3-0.4. From the top/bottom liquid film thickness comparison results, it was observed that the gravity forces are fully suppressed and overcome by the surface tension and shear forces when the confinement number approaches 1, Co ≈ 1. Thus, as a new approximate rule, the lower threshold of macroscale flow is Co = 0.3-0.4 while the upper threshold of symmetric microscale flow is Co ≈ 1 with a transition (or mesoscale) region in-between.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

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.     

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.     

Subcooled flow boiling of R-134a in vertical channels of small diameter

Subcooled flow boiling heat transfer for refrigerant R-134a in vertical cylindrical tubes with 0.83, 1.22 and 1.70 mm internal diameter was experimentally investigated. The effects of the heat flux, q″ = 1–26 kW/m², mass flux, G = 300–700 kg/m² s, inlet subcooling, ΔT_sub,i = 5–15 °C, system pressure, P = 7.70–10.17 bar, and channel diameter, D, on the subcooled boiling heat transfer were explored in detail. The results are presented in the form of boiling curves and heat transfer coefficients. The boiling curves evidenced the existence of hysteresis when increasing the heat flux until the onset of nucleate boiling, ONB. The wall superheat at ONB was found to be essentially higher than that predicted with correlations for larger tubes. An increase of the mass flux leads, for early subcooled boiling, to an increase in the heat transfer coefficient. However, for fully developed subcooled boiling, increases of the mass flux only result in a slight improvement of the heat transfer. Higher inlet subcooling, higher system pressure and smaller channel diameter lead to better boiling heat transfer. Experimental heat transfer coefficients are compared to predictions from classical correlations available in the literature. None of them predicts the experimental data for all tested conditions.     

Flow boiling heat transfer of R134a in the multiport minichannel heat exchangers

The flow boiling heat transfer characteristics of R134a in the multiport minichannel heat exchangers are presented. The heat exchanger was designed as the counter flow tube-in-tube heat exchanger with refrigerant flowing in the inner tube and hot water in the gap between the outer and inner tubes. Two inner tubes were made from extruded multiport aluminium with the internal hydraulic diameter of 1.1 mm for 14 numbers of channels and 1.2 mm for eight numbers of channels. The outer surface areas of two inner test sections are 5979 mm² and 6171 m², while the inner surface areas are 13,545 mm² and 8856 mm² for 14 and eight numbers of channels, respectively. The outer tube of heat exchanger was made from circular acrylic tube with an internal hydraulic diameter of 25.4 mm. The experiments were performed at the heat fluxes between 15 and 65 kW/m², mass flux of refrigerant between 300 and 800 kg/m² s and saturation pressure ranging from 4 to 6 bar. For instance the boiling curve, average heat transfer coefficients are discussed. The comparison results of two test sections with different the number of channels are investigated. The results are also compared with nine existing correlations. The new correlation for predicting the heat transfer coefficient was also proposed.     

Saturated critical heat flux in a multi-microchannel heat sink fed by a split flow system

An extensive experimental campaign has been carried out for the measurement of saturated critical heat flux in a multi-microchannel copper heat sink. The heat sink was formed by 29 parallel channels that were 199 μm wide and 756 μm deep. In order to increase the critical heat flux and reduce the two-phase pressure drop, a split flow system was implemented with one central inlet at the middle of the channels and two outlets at either end. The base critical heat flux was measured using three HFC Refrigerants (R134a, R236fa and R245fa) for mass fluxes ranging from 250 to 1500 kg/m² s, inlet subcoolings from ?25 to ?5 K and saturation temperatures from 20 to 50 °C. The parametric effects of mass velocity, saturation temperature and inlet subcooling were investigated. The analysis showed that significantly higher CHF was obtainable with the split flow system (one inlet–two outlets) compared to the single inlet–single outlet system, providing also a much lower pressure drop. Notably several existing predictive methods matched the experimental data quite well and quantitatively predicted the benefit of higher CHF of the split flow.     

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.     

Experimental study on condensation heat transfer in vertical minichannels for new refrigerant R1234ze(E) versus R134a and R236fa

Experimental condensation heat transfer data for the new refrigerant R1234ze(E), trans-1,3,3,3-tetrafluoropropene, are presented and compared with refrigerants R134a and R236fa for a vertically aligned, aluminum multi-port tube. Local condensation heat transfer measurements with such a multi-microchannel test section are very challenging due to the large uncertainties related to the heat flux estimation. Presently, a new experimental test facility was designed with a test section to directly measure the wall temperature along a vertically aligned aluminum multi-port tube with rectangular channels of 1.45 mm hydraulic diameter. Then, a new data reduction process was developed to compute the local condensation heat transfer coefficients accounting for the non-uniform distribution of the local heat flux along the channels. The condensation heat transfer coefficients showed the expected decrease as the vapor quality decreased (1.0-0.0) during the condensation process, as the mass velocity decreased (260-50 kg m⁻² s⁻¹) and as the saturation temperature increased (25-70 °C). However, the heat transfer coefficients were not affected by the condensing heat flux (1-62 kW m⁻²) or by the entrance conditions within the tested range. It was found that the heat transfer performance of R1234ze(E) was about 15-25% lower than for R134a but relatively similar to R236fa. The experimental data were then compared with leading prediction methods from the literature for horizontal channels. In general, the agreement was poor, over-predicting the high Nusselt number data and under-predicting the low Nusselt number data, but capturing the mid-range quite well. A modified correlation was developed and yielded a good agreement with the current database for all three fluids over a wide range of operating conditions.     

Natural convection boiling of water and surfactants in narrow horizontal annular channels

Natural convection boiling of water and surfactants at atmospheric pressure in narrow horizontal annular channels was studied experimentally in the range of Bond numbers Bo = 0.185–1.52. The flow pattern was visualized by high-speed video recording to identify the different regimes of boiling of water and surfactants. The channel length was 24 mm and 36 mm, the gap size was 0.45, 1.2, 2.2, and 3.7 mm. The heat flux was in the range of 20–500 kW/m², the concentration of surfactant solutions was varied from 10 to 600 ppm. For water boiling at Bond numbers Bo < 1 the CHF in restricted space is lower than that in unconfined space. This effect increases with increasing the channel length. For water at Bond number Bo = 1.52, boiling can almost be considered as unconfined. Additive of surfactant led to enhancement of heat transfer compared to water boiling in the same gap size, however, this effect decreased with decreasing gap size. For the same gap size, CHF in surfactant solutions was significantly lower than that in water. Hysteresis was observed for boiling in degraded surfactant solutions.     

Nucleate boiling and critical heat flux of HFE-7100 in horizontal narrow spaces

The pool boiling heat transfer and critical heat flux CHF of saturated HFE-7100 at atmospheric pressure on a confined smooth copper surface were experimentally studied. The horizontal upward boiling surface was confined by a face-to-face parallel unheated surface. We analysed the effects obtained by changing the diameter of the unheated surface and the gap between the boiling surface and the adiabatic surface. The gap values investigated were s = 0.5, 1.0, 2.0, 3.5 mm. To confine the circular boiling surface (d = 30 mm), two different Plexiglas discs were used: one with a diameter D = 30 mm, equal to that of the copper boiling surface, and the other with a diameter D = 60 mm, equal to that of the overall test section support. For each configuration, boiling curves were obtained up to the thermal crisis. For both configurations, it was observed that, at low wall superheat, the effect of confinement was not significant if Bo > 1, while for Bo ? 1 the heat transfer coefficient increased as the channel width s decreased. By contrast, at high wall superheat, a drastic reduction in both heat transfer and CHF was seen when the channel width s decreased; this reduction was less pronounced when the smaller confinement disc (D = 30 mm) was used. CHF data were also compared with the values predicted by literature correlations.     

Investigation of saturated critical heat flux in a single, uniformly heated microchannel

A series of tests have been performed to determine the saturated critical heat flux (CHF) in 0.5 and 0.8 mm internal diameter microchannel tubes as a function of refrigerant mass velocity, heated length, saturation temperature and inlet liquid subcooling. The tested refrigerants were R-134a and R-245fa and the heated length of microchannel was varied between 20 and 70 mm. The results show a strong dependence of CHF on mass velocity, heated length and microchannel diameter but no influence of liquid subcooling (2–15 °C) was observed. The experimental results have been compared to the well-known CHF single-channel correlation of Y. Katto and H. Ohno [An improved version of the generalized correlation of critical heat flux for the forced convective boiling in uniformly heated vertical tubes, Int. J. Heat and Mass Transfer 27 (9) (1984) 1641–1648] and the multichannel correlation of W. Qu and I. Mudawar [Measurement and correlation of critical heat flux in two-phase microchannel heat sinks, Int. J. Heat and Mass Transfer 47 (2004) 2045–2059]. The comparison shows that the correlation of Katto–Ohno predicts microchannel data with a mean absolute error of 32.8% with only 41.2% of the data falling within a ±15% error band. The correlation of Qu and Mudawar shows the same trends as the CHF data but significantly overpredicts them. Based on the present experimental data, a new microscale version of the Katto–Ohno correlation for the prediction of CHF during saturated boiling in microchannels has been proposed.     

Flow boiling enhancement of FC-72 from microporous surfaces in minichannels

The microporous coatings can remarkably enhance the liquid boiling heat transfer. Therefore, they are promising to be introduced into minichannels in the design of the cooling system of high-power microchips. However, the flow boiling heat transfer characteristics from microporous surfaces in the minichannels have not been extensively studied, and the pertinent knowledge is rather fragmentary. The present research is an experimental investigation on flow boiling of a dielectric fluid FC-72 from microporous coating surfaces in horizontal, rectangular minichannels of 0.49, 0.93 and 1.26 mm hydraulic diameter. Effects of coating structural parameters, such as the particle diameter and coating thickness, were investigated to identify the optimum microporous coating for heat transfer enhancement. All microporous surfaces in this paper were found to significantly enhance FC-72 flow boiling heat transfer in minichannels. With the optimum coating, the heat transfer coefficients could be 7-10 times those of the uncoated surface, and the boiling wall temperature was reduced by about 10 K. The flow boiling phenomena in the present minichannels were distinctly different from those in conventional-sized channels, due to the wall confinement effect on vapor bubbles. The confinement effect was evaluated by taking the contributions of the liquid mass flux and channel size into consideration. It was found that the very strong confinement effect was unfavorable with respect to flow boiling enhancement of the microporous coatings in the minichannels.     

Visualization study of the effects of nanoparticles surface deposition on convective flow boiling CHF from a short heated wall

Enhancements of nucleate boiling critical heat flux (CHF) using nanofluids in a pool boiling are well known. Considering importance of flow boiling heat transfer in various practical applications, an experimental study on CHF enhancements of nanofluids under convective flow conditions was performed. Changing flow velocity from 0 m/s to 4 m/s, the water boiling on nanoparticles-coated heater was conducted and CHF increased at a given velocity. To understand clearly the mechanism of flow boiling CHF enhancement in nanofluid, the visualization of the nucleate boiling and CHF phenomenon was conducted using the high-speed video camera. It was found that the boiling heat transfer on the nanoparticles-coated heater was lower than that on bare heater, which induced the different flow regime at same heat flux. The different wetting zone on bare and nanoparticles-coated heaters was observed by visualization study. Based the wetting zone fraction, there was brief that the nucleate boiling fraction on heater would be related with the surface wettability. A new concept of flow boiling model was proposed based on the wetting zone fraction. Finally, the effect of nanoparticles deposition layer on the heater was interpreted with the physical mechanisms to increase CHF.     

Experimental study of critical heat flux enhancement during forced convective flow boiling of nanofluid on a short heated surface

Enhancements of nucleate boiling critical heat flux (CHF) using nanofluids in a pool boiling are well-known. Considering importance of flow boiling heat transfer in various practical applications, an experimental study on CHF enhancements of nanofluids under convective flow conditions was performed. A rectangular flow channel with 10-mm width and 5-mm height was used. A 10 mm-diameter disk-type copper surface, heated by conduction heat transfer, was placed at the bottom surface of the flow channel as a test heater. Aqueous nanofluids with alumina nanoparticles at the concentration of 0.01% by volume were investigated. The experimental results showed that the nanofluid flow boiling CHF was distinctly enhanced under the forced convective flow conditions compared to that in pure water. Subsequent to the boiling experiments, the heater surfaces were examined with scanning electron microscope and by measuring contact angle. The surface characterization results suggested that the flow boiling CHF enhancement in nanofluids is mostly caused by the nanoparticles deposition of the heater surface during vigorous boiling of nanofluids and the subsequent wettability enhancements.     

Effects of channel dimension,heat flux,and mass flux on flow boiling regimes in microchannels

Experiments are conducted with a perfluorinated dielectric fluid, Fluorinert FC-77, to investigate the effects of channel size and mass flux (225–1420 kg/m²s) on microchannel flow boiling regimes by means of high-speed photography. Seven different silicon test pieces with parallel microchannels of widths ranging from 100 to 5850 μm, all with a depth of 400 μm, are considered. Flow visualizations are performed with a high-speed digital video camera while local measurements of the heat transfer coefficient are simultaneously obtained. The visualizations and the heat transfer data show that flow regimes in the microchannels of width 400 μm and larger are similar, with nucleate boiling being dominant in these channels over a wide range of heat flux. In contrast, flow regimes in the smaller microchannels are different and bubble nucleation at the walls is suppressed at a relatively low heat flux for these sizes. Two types of flow regime maps are developed and the effects of channel width on the flow regime transitions are discussed.     

Status of prediction methods for critical heat fluxes in mini and microchannels

Saturated critical heat flux (CHF) is an important issue during flow boiling in mini and microchannels. To determine the best prediction method available in the literature, 2996 data points from 19 different laboratories have been collected since 1958. The database includes nine different fluids (R-134a, R-245fa, R-236fa, R-123, R-32, R-113, nitrogen, CO₂ and water) for a wide range of experimental conditions. This database has been compared to 6 different correlations and 1 theoretically based model. For predicting the non-aqueous fluids, the theoretical model by Revellin and Thome [Revellin, R., Thome, J.R., 2008. A theoretical model for the prediction of the critical heat flux in heated microchannels. Int. J. Heat Mass Transfer 51, 1216–1225] is the best method. It predicts 86% of the CHF data for non-aqueous fluids within a 30% error band. The data for water are best predicted by the correlation by Zhang et al. [Zhang, W., Hibiki, T., Mishima, K., Mi, Y., 2006. Correlation of critical heat flux for flow boiling of water in minichannels. Int. J. Heat Mass Transfer 49, 1058–1072]. This method predicts 83% of the CHF data for water within a 30% error band. Some suggestions have also been proposed in this paper for the future studies.     

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.     

Influence of the aspect ratio on boiling flows in rectangular mini-channels

This article presents experiments conducted with two single rectangular mini-channels of same hydraulic diameter (1.4 mm) and different aspect ratios for conditions of horizontal boiling flow. The Forane® 365 HX used was subcooled (ΔT_sub = 15 °C) for all the boiling curves presented in the paper. Local heat transfer coefficients were measured for heat flux ranging from 25 to 62 kW m⁻² and mass flux from 200 kg m⁻² s⁻¹ to 400 kg m⁻² s⁻¹. The boiling flows were observed with two different cameras (depending on the flow velocity) through a visualization window. The flow patterns in the two channels were compared for similar conditions. The results show that the boiling heat transfer coefficient and the pressure drop values are different for the two single mini-channels. For low heat flux condition, the channel with lowest aspect ratio (H/W = 0.143) has a higher heat transfer coefficient. On the other hand, for high heat flux condition, the opposite situation occurs, namely the heat transfer coefficient becomes higher for the channel with highest aspect ratio (H/W = 0.43). This is probably due to the earlier onset of dryout in the channel with lowest aspect ratio. For the two cases of heating, the pressure drop for the two-phase flow remains lower for the channel with lowest aspect ratio. These results show that the aspect ratio plays a substantial role for boiling flows in rectangular channels. As for single-phase flows, the heat transfer characteristics are significantly influenced (even though the hydraulic diameter remains the same) by this parameter.     

Critical heat flux in forced convective subcooled boiling with a plane wall jet (effect of subcooling on CHF)

Critical heat flux (CHF) during forced convective subcooled boiling with a plane jet has been yet made insufficient except for saturation condition when comparing CHF with impinging jet system including multiple jets. The present experiment has measured the CHF with plane jet on a rectangular heated surface of 40, 60, and 80 mm in length and 10 and 20 mm in width. Subcooled liquid being supplied through the plane jet with a different thickness of 1 and 2 mm, covers the heated surface where rigorous boiling takes place. The experiment varies a jet velocity from 3 to 15 m/s, a subcooling from 0 to 60 K, and system pressure at 0.1 MPa for water and at 1.5 to 3.0 MPa for R22. It is found that the existing correlation for saturation can be applied to the CHF of water, but hardly to the CHF of R22 in spite of saturation condition. After the effects of jet velocity and subcooling on the CHF can be elucidated, a new correlation including the effect of subcooling is proposed to predict most of the CHF data within an accuracy of ±20 percent. This correlation for saturated condition is found to interestingly agree with that theoretically derived by applying the Katto and Haramura criterion to this system. Received on 8 January 1997