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.     

Periodic boiling in parallel micro-channels at low vapor quality

Based on experimental investigations the present study evaluates instability and heat transfer phenomenon under condition of periodic flow boiling of water and ethanol in parallel triangular micro-channels. Tests were performed in the range of hydraulic diameter 100–220 μm, mass flux 32–200 kg/m² s, heat flux 120–270 kW/m², vapor quality x = 0.01–0.08. The period between successive events depends on the boiling number and decreases with an increase in the boiling number. The initial film thickness decreases with increasing heat flux. When the liquid film reached the minimum initial film thickness CHF regime occurred. Temporal variations of pressure drop, fluid and heater temperatures were periodic. Oscillation frequency is the same for the pressure drop, for the fluid temperature at the outlet manifold, and for the mean and maximum heater temperature fluctuations. All these fluctuations are in phase. The CHF phenomenon is different from that observed in a single channel of conventional size. A key difference between micro-channel heat sink and single conventional channel is amplification of parallel-channel instability prior to CHF. The dimensionless experimental values of the heat transfer coefficient are presented as the Nusselt number dependence on the Eotvos number and the boiling number.     

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.     

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].     

R410A condensation inside a commercial brazed plate heat exchanger

This paper presents the heat transfer coefficients and the pressure drop measured during HFC-410A condensation inside a commercial brazed plate heat exchanger: the effects of saturation temperature, refrigerant mass flux and vapour super-heating are investigated. The heat transfer coefficients show weak sensitivity to saturation temperature and great sensitivity to refrigerant mass flux and vapour super-heating. At low refrigerant mass flux (<20 kg/m² s) the saturated vapour condensation heat transfer coefficients are not dependent on mass flux and are well predicted by Nusselt [W. Nusselt, Die oberflachenkondensation des wasserdampfes, Energy 60 (1916) 541–546, 569–575] analysis for vertical surface: the condensation process is gravity controlled. For higher refrigerant mass flux (>20 kg/m²s) the saturated vapour condensation heat transfer coefficients depend on mass flux and are well predicted by Akers et al. [W.W. Akers, H.A. Deans, O.K. Crosser, Condensing heat transfer within horizontal tubes, Chem. Eng. Prog. Symp. Series 55 (1959) 171–176] equation: forced convection condensation occurs. In the forced convection condensation region the heat transfer coefficients show a 30% increase for a doubling of the refrigerant mass flux. The condensation heat transfer coefficients of super-heated vapour are 8–10% higher than those of saturated vapour and are well predicted by Webb [R.L. Webb, Convective condensation of superheated vapor, ASME J. Heat Transfer 120 (1998) 418–421] model. A simple linear equation based on the kinetic energy per unit volume of the refrigerant flow is proposed for the computation of the frictional pressure drop.     

Two-phase flow operational maps for multi-microchannel evaporators

The current paper presents new operational maps for several different multi-microchannel evaporators, with and without any inlet restrictions (micro-orifices), for the two-phase flow of refrigerants R245fa, R236fa, and R1234ze(E). The test fluids flowed in 67 parallel channels, each having a cross-sectional area of 100 × 100 μm². In order to emulate the power dissipated by active components in a 3D CMOS CPU chip, two aluminium microheaters were sputtered onto the back-side of the test section providing a 0.5 cm² each. Without any inlet restrictions in the micro-evaporator, significant parallel channel flow instabilities, vapor back flow, and flow maldistribution led to high-amplitude and high-frequency temperature and pressure oscillations. Such undesired phenomena were then prevented by placing restrictions at the inlet of each channel. High-speed flow visualization distinguished eight different operating regimes of the two-phase flow depending on the tested operating conditions. Therefore, the preferred operating regimes can be easily traced. In particular, flashing two-phase flow without back flow appeared to be the best operating regime without any flow and temperature instabilities.     

Periodic flow boiling in a non-uniformly heated microchannel heat sink

Hydrodynamic and thermal characteristics of flow boiling in a non-uniformly heated microchannel were studied. Experiments were performed with a single microchannel and a series of microheaters to study the microscale boiling of water under axially non-uniform heat input conditions. A simultaneous real time visualization of the flow pattern was performed with the measurement of experimental parameters. Tests were performed over a mass flux of 309.8 kg/m² s, and heat flux of 200–600 kW/m². Test results showed different fluctuations of heated wall temperature, pressure drop, and mass flux with variations of the heat input along the flow direction. The unique periodic flow boiling in a single microchannel was observed at all heat flux conditions except for the increasing heat input distribution case which is the nearly uniform effective heat input distribution condition. The instability is correlated with flow pattern transition. For the nearly uniform effective heating condition, no fluctuation of the wall temperature, pressure drop, or mass flux was observed. We can relieve the instability by increasing total heat input along the flow direction and predict the instability using the transition criteria and flow pattern map.     

Low-GWP refrigerants flow boiling heat transfer in a 5 PPI copper foam

This paper reports an experimental investigation of the heat transfer performance of the new low-GWP refrigerants, R1234yf and R1234ze(E), during flow boiling heat transfer inside a horizontal high porosity copper foam with 5 Pores Per Inch (PPI). Metal foams are a class of cellular structured materials consisting of a stochastic distribution of interconnected pores; these materials have been proposed as effective solutions for heat transfer enhancement during both single and two-phase heat transfer. R1234yf and R1234ze(E) refrigerants are appealing alternatives of the more traditional R134a by virtue of their negligible values of GWP and normal boiling temperatures close to that of R134a, which make them suitable solution in several different applications, such as: refrigeration and air conditioning and electronic thermal management. This work compares the two-phase heat transfer behaviour of these new HFO refrigerants, studying the boiling process inside a porous medium and permitting to understand their effective heat transfer capabilities. The experimental measurements were carried out by imposing three different heat fluxes: 50, 75, and 100 kW m⁻², at a constant saturation temperature of 30 °C; the refrigerant mass velocity was varied between 50 and 200 kg m⁻² s⁻¹, whilst the mean vapour quality varied from 0.2 to 0.95. The two-phase heat transfer and pressure drop performance of the two new HFO refrigerants is compared against that of the more traditional R134a.     

Pressure drop and heat transfer characteristics of boiling water in sub-hundred micron channel

The current work focuses on the pressure drop, heat transfer and stability in two phase flow in microchannels with hydraulic diameter of less than one hundred microns. Experiments were conducted in smooth microchannels of hydraulic diameter of 45, 65 μm, and a rough microchannel of hydraulic diameter of 70 μm, with deionised water as the working fluid. The local saturation pressure and temperature vary substantially over the length of the channel. In order to correctly predict the local saturation temperature and subsequently the heat transfer characteristics, numerical techniques have been used in conjunction with the conventional two phase pressure drop models. The Lockhart–Martinelli (liquid–laminar, vapour–laminar) model is found to predict the two phase pressure drop data within 20%. The instability in two phase flow is quantified; it is found that microchannels of smaller hydraulic diameter have lesser instabilities as compared to their larger counterparts. The experiments also suggest that surface characteristics strongly affect flow stability in the two phase flow regime. The effect of hydraulic diameter and surface characteristics on the flow characteristics and stability in two phase flow is seldom reported, and is of considerable practical relevance.     

Experimental evaluation of cooling performance of circular heat sinks for heat dissipation from electronic chips using nanofluid

In this work, an experimental investigation on cooling performance of using nanofluid to replace the pure water as the coolant in a minichannel heat sink is conducted. The heat sink comprises of four circular channels with hydraulic diameter of 6 mm. Thermal and hydraulic performances of the nanofluid cooled minichannel heat sink are evaluated from the results obtained for the Nusselt number, friction factor, thermal resistance and pumping power, with the volume flow rate ranging from 0.3 to 1.5 L/min. The experimental results show that the nanofluid cooled heat sink outperforms the water-cooled one, having significantly higher average heat transfer coefficient. Despite the marked increase in dynamic viscosity due to dispersing the nanoparticles in water, the friction factor for the nanofluid-cooled heat sink is found slightly increased only.     

Thermal performance of a carbon fiber composite material heat sink in an FC-72 thermosyphon

This study analyzes the use of a carbon fiber epoxy heat sink for evaporator surface enhancement in a FC-72 thermosyphon. The pin-fin heat sink features 945 small-cross-section (1.27 mm by 0.965 mm) fins fabricated with an integral base plate. These fins have a high thermal conductivity (500 W/m K) along the length of the fin. The influence of heat load, thermosyphon fill volume, and condenser operating temperature on the overall thermal performance is examined. The results of this experiment provide significant insight into the possible implementation and potential benefits of carbon-fiber heat sink technology in two-phase flow leading to significant improvements in thermal management strategies for advanced electronics.     

Flow boiling heat transfer study of R-134a/R-290/R-600a mixture in 9.52 and 12.7 mm smooth horizontal tubes: Experimental investigation

A detailed experimental investigation is carried out to study the flow boiling heat transfer behavior of R-134a/R-290/R-600a (91%/4.068%/4.932% by mass) refrigerant mixture in smooth horizontal tubes of diameter 9.52 and 12.7 mm. The heat transfer coefficients of the mixture are experimentally measured under varied heat flux conditions for stratified flow patterns using a coaxial counter-current heat exchanger test section. The tests are conducted for refrigerant inlet temperatures between ?9 and 5 °C and mass flow rates ranging from 3 to 5 g s^?1. Kattan–Thome–Favrat maps are used to confirm the flow patterns for the tested conditions. The magnitude of the heat transfer coefficient with respect to flow patterns and different mechanisms of boiling are discussed. The heat transfer coefficient of the refrigerant mixture is also compared with that of R-134a for selected working conditions. The significance of nucleate boiling in the overall heat transfer process under these testing conditions is highlighted.     

Flow pattern characterization for R-245fa in minichannels: Optical measurement technique and experimental results

An optical measurement method using image processing for two-phase flow pattern characterization in minichannel is developed. The bubble frequency, the percentage of small bubbles as well as their velocity are measured. A high-speed high-definition video camera is used to measure these parameters and to identify the flow regimes and their transitions. The tests are performed in a 3.0 mm glass channel using saturated R-245fa at 60 °C (4.6 bar). The mass velocity is ranging from 100 to 1500 kg/m² s, the heat flux is varying from 10 to 90 kW/m² and the inlet vapor quality from 0 to 1. Four flow patterns (bubbly flow, bubbly–slug flow, slug flow and annular flow) are recognized. The comparison between the present experimental intermittent/annular transition lines and five transition lines from macroscale and microscale flow pattern maps available in the literature is presented. Finally, the influence of the flow pattern on the heat transfer coefficient is highlighted.     

Flow boiling of halogenated refrigerants at high saturation temperature in a horizontal smooth tube

Several correlations are available in the open literature for computing the heat transfer coefficient during flow boiling inside plain channels. With respect to halogenated refrigerants, these correlations are usually compared to data taken in a limited range of evaporation temperature and reduced pressure. More recently, the adoption of new refrigerants, such as high pressure HFCs and carbon dioxide, requires to largely extend the pressure range of application of such correlations. Besides, the design of evaporators for some heat pumping applications, where temperatures are set at higher values as compared to usual evaporating temperatures in air-conditioning equipment, requires proper validation of the computing methods.The present paper aims at comparing four well-known predicting models to a new database collected during flow boiling of HCFC (R22) and HFC refrigerants (R134a, R125 and R410A) in a horizontal 8 mm internal diameter tube. This database is characterized by saturation temperature ranging between 25 °C and 45 °C, reduced pressure spanning between 0.19 and 0.53. Mass velocity ranges between 200 and 600 kg m^?2 s^?1 and heat flux between 9 and 53 kW m^?2.Evaporating heat transfer coefficients of halogenated refrigerants at such high temperatures have not been reported in the open literature so far. The discussion of the results will enlighten some similarities with experimental trends presented in the literature for evaporating carbon dioxide.Two models tested here show good prediction capabilities of the present experimental data, but not for all the data sets in the same way. For the purpose of practical use, a simple modification of the correlation by Gungor and Winterton [1] is proposed, showing that this is able to catch the experimental trends of the present database with good agreement.     

Experimental study of the blockage boundary for dense-phase pneumatic conveying of powders through a horizontal slit

The estimation of the blockage boundary for pneumatic conveying through a slit is of significant importance. In this paper, we investigate the characteristics for blockage of powder (48 μm average diameter) through a horizontal slit (1.6 m × 0.05 m × 0.002 m). The results show that the required critical solid mass flow rate increases as the superficial air velocity increases superficial air velocity. The solid loading ratio and superficial air velocity displayed a decreasing power law relationship. This finding agrees with existing theory and experimental results. However, a minimum inlet solid loading ratio exists. When the air velocity is greater than the corresponding air velocity of the minimum solid loading ratio, the solid loading ratio exhibits an increasing trend in power law. We also found that when the inlet conveying pressure increased, the critical solid mass flow rate required for blockage, the inlet solid loading ratio, and the minimum inlet solid loading ratio increased.     

Impact of aspect ratio on flow boiling of water in rectangular microchannels

In this paper we focus on the impact of varying the aspect ratio of rectangular microchannels, on the overall pressure drop involving water boiling. An integrated system comprising micro-heaters, sensors and microchannels has been realized on (1 1 0) silicon wafers, following CMOS compatible process steps. Rectangular microchannels were fabricated with varying aspect ratios (width [W] to depth [H]) but constant hydraulic diameter of 142 ± 2 μm and length of 20 mm. The invariant nature of the hydraulic diameter is confirmed through two independent means: physical measurements using profilometer and by measuring the pressure drop in single-phase fluid flow. The experimental results show that the pressure drop for two-phase flow in rectangular microchannels experiences minima at an aspect ratio of about 1.6. The minimum is possibly due to opposing trends of frictional and acceleration pressure drops, with respect to aspect ratio. In a certain heat flux and mass flux range, it is observed that the two-phase pressure drop is lower than the corresponding single-phase value. This is the first study to investigate the effect of aspect ratio in two-phase flow in microchannels, to the best of our knowledge. The results are in qualitative agreement with annular flow model predictions. These results improve the possibility of designing effective heat-sinks based on two-phase fluid flow in microchannels.     

Two-phase flow of refrigerants in 85 μm-wide multi-microchannels: Part I – Pressure drop

This article is the first part of a study on flow boiling of R236fa and R245fa. This part presents pressure drop measurements obtained on a silicon multi-microchannel evaporator with 85 μm wide and 560 μm high channels separated by 46 μm wide fins. The 135 microchannels were 12.7 mm long. Dielectric refrigerants R236fa and R245fa were used as the evaporating test fluids. The inlet saturation temperature was maintained at 30.5 °C while the mass fluxes were varied from 499 to 1100 kg/m² s and the base heat flux was tested from 130 to 1400 kW/m². A new experimental technique was developed to measure the outlet pressure losses, which represented up to 30% of the total pressure drop and thus cannot be neglected. The microchannel pressure drop measurements were very well predicted by the method of Cioncolini et al. (2009).     

Experimental research on heat transfer of natural convection in vertical rectangular channels with large aspect ratio

This work presents the experimental research on the steady laminar natural convection heat transfer of air in three vertical thin rectangular channels with different gap clearance. The much higher ratio of width to gap clearance (60–24) and the ratio of length to gap clearance (800–320) make the rectangular channels similar with the coolant flow passage in plate type fuel reactors. The vertical rectangular channels were composed of two stainless steal plates and were heated by electrical heating rods. The wall temperatures were detected with the K-type thermocouples which were inserted into the blind holes drilled in the steal plates. Also the air temperatures at the inlet and outlet of the channel were detected. The wall heat fluxes added to the air flow were calculated by the Fourier heat conduction law. The heat transfer characteristics were analyzed, and the average Nusselt numbers in all the three channels could be well correlated with the Rayleigh number or the modified Rayleigh number in a uniform correlation. Furthermore, the maximum wall temperatures were investigated, which is a key parameter for the fuel’s integrity during some accidents. It was found that even the wall heat flux was up to 1500 W/m², the maximum wall temperature was lower than 350 °C. All this work is valuable for the plate type reactor’s design and safety analysis.     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

Wavelet decomposition method decoupled boiling/evaporation oscillation mechanisms over two to three timescales: A study for a microchannel with pin fin structure

Boiling/evaporation heat transfer in a microchannel with pin fin structure was performed with water as the working fluid. Simultaneous measurements of various parameters were performed. The chip wall temperatures were measured by a high spatial-time resolution IR image system, having a sensitivity of 0.02 °C. The flow pattern variations synchronously changed wall temperatures due to ultra-small Bi number. The wavelet decomposition method successfully identified the noise signal and decoupled various temperature oscillations with different amplitudes and frequencies. Three types of temperature oscillations were identified according to heat flux q and mass flux G. The first type of oscillation occurred at q/G < 0.62 kJ/kg. The approximation coefficient of wavelet decomposition decided the dominant cycle period which was ∼3 times of the fluid residence time in the microchannel, behaving the density wave oscillation characteristic. The detail coefficients of wavelet decomposition decided the dominant cycle period, which matched the flow pattern transition determined value well, representing the flow pattern transition induced oscillation. For the second type of oscillation, the wavelet decomposition decoupled the three oscillation mechanisms. The pressure drop oscillation caused the temperature oscillation amplitudes of 5–10 °C and cycle periods of 10–15 s. The density wave oscillation and flow pattern transition induced oscillation are embedded with both the pressure rise and decrease stages of the pressure drop oscillation. The third type of oscillation happened at q/G > 1.13 kJ/kg, having the density wave oscillation coupled with the varied liquid film evaporation induced oscillation. The liquid island, retention bubble induced nucleation sites and cone-shape two-phase developing region are unique features of microchannel boiling with pin fin structure. This study illustrated that pressure drop oscillation and density wave oscillation, usually happened in large size channels, also take place in microchannels. The flow pattern transition and varied liquid film evaporation induced oscillations are specific to microchannel boiling/evaporation flow.