Carbon dioxide flow boiling in a single microchannel - Part I: Pressure drops

Carbon dioxide two-phase flow pressure drops have been investigated in a single horizontal stainless-steel micro-tube having a 0.529 mm inner diameter. Experiments were carried out in adiabatic conditions for four saturation temperatures of −10; −5; 0; 5 °C and mass fluxes ranging from 200 to 1400 kg/m² s, for inlet qualities up to unity. Measurements have been compared to the predictions of well-known methods. The Müller-Steinhagen and Heck correlation and the Friedel correlation gave the best fit as well as the homogeneous model when the liquid viscosity is used to represent the apparent two-phase viscosity. Further, an analysis based on the homogeneous model has not shown any clear appearance of the laminar or the transition regimes in any given range of Reynolds number. The apparent viscosity of the two-phase mixture was found larger than the liquid viscosity at low vapour qualities, namely at the lowest temperatures. Hence, a new expression to determine the equivalent viscosity was suggested as a function of the reduced pressure. Lastly, the Chisholm parameter from the Lockhart-Martinelli correlation was found lower than expected and also mainly dependent on the saturation temperature.     

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.     

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.     

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.     

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.     

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 and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included ground-based tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.     

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.     

Flow boiling behaviors in hydrophilic and hydrophobic microchannels

Surface wettability is a critical parameter in small scale phenomena, especially two-phase flow, since the surface force becomes dominant as size decreases. In present study, experiments of water flow boiling in hydrophilic and hydrophobic rectangular microchannels were conducted to investigate the wettability effect on flow boiling in rectangular microchannels. The rectangular microchannels were fabricated with a photosensitive glass to visualize flow pattern. The hydrophilic bare photosensitive glass microchannel was chemically treated to obtain a hydrophobic microchannel. And, visualization of flow patterns was carried out. And boiling heat transfer and two-phase pressure drop was analyzed with visualization results. The boiling heat transfer coefficient in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with nucleation site density and liquid film motion. And the pressure drop in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with unstable motions of bubble and liquid film. Finally, we find out the wettability is important parameter on the flow pattern, which were highly related with two-phase heat and mass transfer.     

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.     

A fractal analysis of subcooled flow boiling heat transfer

A fractal model for the subcooled flow boiling heat transfer is proposed in this paper. The analytical expressions for the subcooled flow boiling heat transfer are derived based on the fractal distribution of nucleation sites on boiling surfaces. The proposed fractal model for the subcooled flow boiling heat transfer is found to be a function of wall superheat, liquid subcooling, bulk velocity of fluid (or Reynolds number), fractal dimension, the minimum and maximum active cavity size, the contact angle and physical properties of fluid. No additional/new empirical constant is introduced, and the proposed model contains less empirical constants than the conventional models. The proposed model takes into account all the possible mechanisms for subcooled flow boiling heat transfer. The model predictions are compared with the existing experimental data, and fair agreement between the model predictions and experimental data is found for different bulk flow rates.     

Two-phase annular flow and evaporative heat transfer in a microchannel

A physical and mathematical model has been developed to predict the two-phase flow and heat transfer in a microchannel with evaporative heat transfer. Sample solutions to the model were obtained for both constant wall temperature and constant wall heat flux conditions. Results are provided for evaporation rate, liquid film thickness, liquid and vapor phase pressure and temperature distributions. In addition to the sample calculations that were used to illustrate the transport characteristics, computations based on the current model were performed to generate results for comparisons with the experimental results of Qu and Mudawar (2004) where two different mass flow rates of the working fluid were used in the experiment. The comparisons of total pressure drops with the experimental data of Qu and Mudawar (2004) cover the wall heat flux range of 142.71-240 W/cm² with a total channel mass flux of 400.1 kg/m² s and also the wall heat flu range of 99.54-204.39 W/cm² with total channel mass flux of 401.9 kg/m² s. The calculated results from the current model match closely with those of Qu and Mudawar (2004).     

Heat transfer in pool boiling of ammonia/water mixture

The nonazeotropic binary mixtures such as, methanol/water, ethanol/water and ammonia/water, have variable boiling and dew points, depending on the combination of substance and those mass fraction. It is expected to have a higher performance as a result of decreasing the thermodynamically irreversible loss, when there is a relevant mass fraction. Therefore, ammonia/water mixture is expected to use as working fluid in small temperature difference power generation cycles and absorption refrigeration cycles. However, few experiments were carried out for measuring heat transfer coefficient for ammonia/water mixture in the world. An experimental study has been carried out to measure boiling heat transfer coefficient of an ammonia/water mixture on a horizontal heated surface at low pressure of 0.2, 0.4 and 0.7 MPa and at low mass fraction of 0 < C < 0.27 and at high pressure 0.7, 1.0 and 1.5 MPa and at mass fraction of 0.5 < C < 1.0 and at heat flux under critical heat flux the heat transfer coefficient are compared with existing correlations prediction and a revised correlation can be proposed to predict them well.     

Numerical simulation of compressible gas flow and heat transfer in a microchannel surrounded by solid media

In this study, a numerical model is developed to investigate the coupled compressible gas flow and heat transfer in a microchannel surrounded by solid media. To accommodate the varying flow cross-section, the compressible gas flow model is established in a non-orthogonal curvilinear coordinate system. An iterative numerical procedure is employed to solve the coupled heat transfer and gas flow equations. The computer code for the compressible gas flow is first validated against two test problems, and then extended by including the heat conduction in the solid media. The effect of the inlet Mach number on the Nusselt number is examined. It is found that the pressure difference from the pyrolysis front to the heated surface is induced essentially by the gas addition from the channel wall, instead from the pyrolysis front. The necessity of accounting for the gas compressibility is clearly demonstrated when severe heating is applied. The pressure distribution obtained along the channel axial direction is useful for further structural analysis of composite materials.     

Enhanced heat transfer in confined pool boiling

We report the results of an experimental investigation of the heat transfer during nucleate boiling on a spatially confined boiling surface. The heat flux as a function of the boiling surface temperature was measured in pool boiling pots with diameters ranging from 15 mm down to 4.5 mm. It was found that a reduction of the pool diameter leads to an enhancement of the nucleate boiling heat flux for most of the boiling curve. Our experimental results indicate that this enhancement is not affected by the depth of the boiling pot, the material of the bounding wall, or the diameter of the inlet water supply. High-speed camera imaging shows that the heat transfer enhancement for the spatially confined pool boiling occurs in conjunction with a stable circulating flow, which is in contrast to the chaotic and mainly upward motion for boiling in larger pool diameters. An explanation for the enhancement of the heat transfer and the associated change in flow pattern is found in the singularisation of the nucleate boiling process.     

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.     

Flow pattern and boiling heat transfer of CO2 in horizontal small-bore tubes

Increasing attention has been focused on carbon dioxide (CO₂) heat pump system where the temperature level is rather low, while the operating pressure is rather high. In this system, the density difference between vapor and liquid becomes rather small, which significantly affects flow patterns. Low surface tension and latent heat also have significant influence on two-phase flow patterns and heat transfer. This paper describes experimental and numerical investigation on flow patterns and heat transfer characteristics of boiling flow CO₂ at high pressure in horizontal small-bore tubes ranging from 1.0 mm to 3.0 mm I.D. Even though the density difference is rather small at high pressure, phase stratification takes place, which leads to the intermittent dryout at the upper wall. So far developed discrete bubble model by the authors for vertical flows is modified so as to include horizontal flow mechanisms. The predicted flow patterns with this new model agree on the whole with the experimental observation.     

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.     

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

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.