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

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

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

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

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

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

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

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

Flow pattern transition instability during flow boiling in a single microchannel

We studied the unique characteristics of flow boiling in a single microchannel, including the periodic pressure drop, mass flow rate, and temperature fluctuations, in terms of a long time period. Experiments were conducted using a single horizontal microchannel and deionized water to study boiling instabilities at very small mass and heat flow rate conditions. A Polydimethylsiloxane (PDMS) rectangular single microchannel had a hydraulic diameter of 103.5 μm and a length of 40 mm. A series of piecewise serpentine platinum microheaters were fabricated on the inner bottom wall of the rectangular microchannel to supply thermal energy to the test fluid. Real-time flow visualizations of the flow pattern inside the microchannel were performed simultaneously with measurements of the experimental parameters. Tests were performed for mass fluxes of 170 and 360 kg/m² s and heat fluxes of 200–530 kW/m². The test results showed that the heated wall temperature, pressure drop, and mass flux all fluctuated with a long period and large amplitude. These periodic fluctuations exactly matched the transition of two alternating flow patterns inside the microchannel: a bubbly/slug flow and an elongated slug/semi-annular flow. Therefore, the flow pattern transition instability in the single microchannel caused a cyclic behavior of the wall temperature, pressure drop, and mass flux, and this behavior had a very long period (100–200 s) and large amplitude.     

Effects of surfactant additive on flow boiling over a microheater under pulse heating

An experimental investigation has been carried out to study effects of surfactant additive on microscale boiling under pulse heating over a Pt microheater (140 × 100 μm²) fabricated in a trapezoidal microchannel (600 μm in width and 150 μm in depth). Experiments are carried out for six different surfactant concentrations of Triton X-100 ranging from 47 ppm to 2103 ppm, for mass flux in the range from 45 kg/m² s to 225 kg/m² s, pulse width in the range from 50 μs to 2 ms, and heat flux in the range from 3 MW/m² to 65 MW/m². As in existing work on pool boiling under steady heating, it is found that nucleate boiling becomes more vigorous and heat transfer is enhanced greatly with the addition of surfactant with maximum boiling heat transfer occurs at the critical micelle concentration (cmc). Furthermore, these maximum values of boiling heat transfer coefficient increase with decreasing pulse width. When concentration is below cmc, the heat flux needed for nucleation increases with increasing concentration and the nucleation temperature is reduced. When concentration is higher than cmc, the boiling heat transfer coefficient decreases and nucleation temperature is higher than that of pure water.     

Hydrodynamics and heat transfer during flow boiling instabilities in a single microchannel

Boiling in microchannels is widely considered as one of the front runners in process intensification heat removal. Flow boiling heat transfer in microchannel geometry and the associated flow instabilities are not well understood, further research is necessary into the flow instabilities adverse effect on heat transfer.Boiling is induced in microchannel geometry (hydraulic diameter 727 μm) to investigate several flow instabilities. A transparent, metallic, conductive deposit has been developed on the exterior of rectangular microchannels, allowing simultaneous heating and visualisation.Presented in this paper is data for a particular case with a uniform heat flux of 4.26 kW/m² applied to the microchannel and inlet liquid mass flowrate, held constant at 1.13 × 10^?5 kg/s. In conjunction with obtaining high-speed images, a sensitive infrared camera is used to record the temperature profiles on the exterior wall of the microchannel, and a data acquisition system is used to record the pressure fluctuations over time. Various phenomena are apparent during the flow instabilities; these can be characterised into timescales occurring at 100’s seconds, 10’s seconds, several seconds and finally milliseconds. Correlation of pressure oscillations with temperature fluctuations as a function of the heat flux applied to the microchannel is possible.From analysis of our results, images and video sequences with the corresponding physical data obtained, it is possible to follow simultaneously particular flow, pressure and temperature conditions leading to nucleate boiling, flow instabilities and transition regimes during flow boiling in a microchannel. The investigation allowed us to quantify and characterise the timescales of various observed instabilities during flow boiling in a microchannel. High speed imaging revealed some of the controlling physical mechanisms responsible for the observed instabilities.     

Nucleate and film boiling on a microheater under pulse heating in a microchannel

Using MEMS technology, a Pt microheater (60 × 100 µm²) fabricated on a glass wafer is placed in a silicon-based microchannel of trapezoidal cross section. With the aid of a high-speed CCD and based on Pt's linear temperature-resistance characteristic, flow boiling phenomena and temperature response on the surface of the microheater in the microchannel under pulse heating are observed and recorded. At a given mass flux, nucleate boiling and film boiling begin to appear on the microheater with increasing heat flux. A flow boiling map, showing the effects of heat and mass flux on nucleate and film boiling regimes on the microheater at a pulse heating width of 2 ms, is presented. It is found that nucleate boiling is changed to film boiling as the heat flux supplied to the microheater is increased. Furthermore, increasing mass flux increases the heat flux required for the incipience of nucleate boiling and film boiling on the microheater in the microchannel.     

Heat transfer characteristics during subcooled flow boiling of a kerosene kind hydrocarbon fuel in a 1 mm diameter channel

The subcooled flow boiling heat transfer characteristics of a kerosene kind hydrocarbon fuel were investigated in an electrically heated horizontal tube with an inner diameter of 1.0 mm, in the range of heat flux: 20–1500 kW/m², fluid temperature: 25–400 °C, mass flux: 1260–2160 kg/m² s, and pressure: 0.25–2.5 MPa. It was proposed that nucleate boiling heat transfer mechanism is dominant, as the heat transfer performance is dependent on heat flux imposed on the channel, rather than the fuel flow rate. It was found that the wall temperatures along the test section kept constant during the fully developed subcooled boiling (FDSB) of the non-azeotropic hydrocarbon fuel. After the onset of nucleate boiling, the temperature differences between inner wall and bulk fluid begin to decrease with the increase of heat flux. Experimental results show that the complicated boiling heat transfer behavior of hydrocarbon fuel is profoundly affected by the pressure and heat flux, especially by fuel subcooling. A correlation of heat transfer coefficients varying with heat fluxes and fuel subcooling was curve fitted. Excellent agreement is obtained between the predicted values and the experimental data.     

Experimental investigation and modelling of heat transfer during convective boiling in a minichannel

Flow boiling heat transfer characteristics of water are experimentally studied in a circular minichannel with an inner diameter of 1500 μm. The fluid flows upwards and the test section, made of the nickel alloy Inconel 600, is directly electrically heated. Thus, the evaporation takes place under the defined boundary condition of constant heat flux. Mass fluxes between 50 and 100 kg/(m² s) and heat fluxes from 10 to 115 kW/m² at an inlet pressure of 3 bar are examined.Infrared thermography is applied to measure the outer wall temperatures of the minichannel. This experimental method permits the identification of different boiling regions, boiling mechanisms and the determination of local heat transfer coefficients. Measurements are carried out in single-phase flow, subcooled and saturated boiling regions. The experimental heat transfer coefficients in the region of saturated boiling are compared with correlations available in literature and with a physically founded model developed for convective boiling.     

On the eruptive boiling in silicon-based microchannels

This study investigates experimentally eruptive boiling in a silicon-based rectangular microchannel with a hydraulic diameter of 33.7 μm, a width of 99.8 μm and a depth-to-width ratio of 0.203. The microchannel is made of SOI wafer and prepared using bulk micro-machining and anodic bonding. The surface roughness for both the bottom and the side walls was measured using an atomic force microscope. The evolution of the eruptive boiling of water in the smooth microchannel was clearly observed using an ultra high-speed video camera (up to 50,000 frames/s) at mass fluxes of 417 and 625 kg/m² s and a heat flux from 14.9 to 372 kW/m². It is confirmed that eruptive boiling is a form of rapid bubble nucleation after which the bubble merges with a slug bubble downstream in a short distance or evolve to a slug bubble. The bubble frequency in all of the cases studied is provided. Eruptive boiling may be predicted classically with nano-sized cavities that are consistent with the measured surface roughness.     

Refrigerant flow boiling heat transfer in parallel microchannels as a function of local vapor quality

Flow boiling of refrigerant HFC-134a in a multi-microchannel copper cold plate evaporator is investigated. The heat transfer coefficient is measured locally for the entire range of vapor qualities starting from subcooled liquid to superheated vapor. The test piece contains 17 parallel, rectangular microchannels (0.762 mm wide) of hydraulic diameter 1.09 mm and aspect ratio 2.5. The design of the test facility is validated by a robust energy balance as well as a comparison of single-phase heat transfer coefficients with results from the literature. Results are presented for four different mass fluxes of 20.3, 40.5, 60.8, and 81.0 kg m^?2 s^?1, which correspond to refrigerant mass flow rates of 0.5–2.0 g s^?1, and at three different pressures 400, 550 and 750 kPa corresponding to saturation temperatures of 8.9, 18.7, and 29 °C. The wall heat flux varies from 0 to 20 W/cm² in the experiments. The heat transfer coefficient is found to vary significantly with refrigerant inlet quality and mass flow rate, but only slightly with saturation pressure for the range of values investigated. The peak heat transfer coefficient is observed for a vapor quality of approximately 20%.     

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

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

A four-zone model for saturated flow boiling in a microchannel of rectangular cross-section

A four-zone flow boiling model is presented to describe saturated flow boiling heat transfer mechanisms in a microchannel of rectangular cross-section. The boiling process in the microchannel is assumed to be a cyclic passage of four zones: (i) liquid-slug zone, (ii) elongated bubble zone, (iii) partially-dryout zone, and (iv) fully-dryout zone. The existence of the partially-dryout zone in this model is proposed to take into consideration of corner effects on boiling heat transfer in the microchannel. To verify this new model, an experimental study was carried out to investigate flow boiling heat transfer of water in a microchannel having a rectangular cross-section with a hydraulic diameter of 137 μm (202 μm in width and 104 μm in depth) with a length of 30 mm under three-side heating condition. The data for bubble nucleation frequency was correlated in terms of the Boiling number, which was used to determine the heat transfer coefficient. It is found that the present four-zone flow boiling model successfully predicts trends of boiling heat transfer data in a microchannel with a rectangular cross-section, having a sharp peak at low vapor quality depending on the mass flow rate. The predictions of flow boiling heat transfer coefficient in the microchannel are found in good agreement with experimental data with a MAE of 13.9%.     

Correlation of critical heat flux for flow boiling of water in mini-channels

In view of practical significance of a correlation of critical heat flux (CHF) in the aspects of engineering design and prediction, this study is aiming at evaluation of existing CHF correlations for flow boiling of water with available databases taken from small-diameter tubes, and then development of a new, simple CHF correlation. Available CHF databases in the literature for flow boiling of water in small-diameter tubes (0.33 < D_h < 6.22 mm) are collected, covering wide parametric ranges. Three correlations by Bowring, Katto and Shah are evaluated with the CHF data for saturated flow boiling, and three correlations by Inasaka–Nariai, Celata et al. and Hall–Mudawar evaluated with the CHF data for subcooled flow boiling. The Hall–Mudawar correlation and the Shah correlation seem to be the most reliable tools for CHF prediction in the subcooled and saturated flow boiling regions, respectively. In order to avoid the defect of predictive discontinuities often encountered when applying previous correlations, a simple, nondimensional, inlet conditions dependent CHF correlation for saturated flow boiling has been formulated. Its functional form is determined by the application of the artificial neural network and parametric trend analyses to the collected database. Superiority of this correlation has been verified by the database. The new correlation has a mean deviation of 16.8% for this collected databank, smallest among all tested correlations. Compared to many inordinately complex correlations, this new correlation consists only of a single equation.     

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

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

Enhanced boiling of HFE-7100 dielectric liquid on porous graphite

Enhanced boiling of HFE-7100 dielectric liquid on porous graphite measuring 10 mm × 10 mm is investigated, and results are compared with those for smooth copper (Cu) of the same dimensions. Although liquid is out-gassed for hours before performing the pool boiling experiments, air entrapped in re-entrant type cavities, ranging in size from tens to hundreds of microns, not only enhanced the nucleate boiling heat transfer and the critical heat flux (CHF), but also, the mixing by the released tiny air bubbles from the porous graphite prior to boiling incipience enhanced the natural convection heat transfer by ∼19%. No temperature excursion is associated with the nucleate boiling on porous graphite, which ensues at very low surface superheat of 0.5–0.8 K. Conversely, the temperature overshoot at incipient boiling on Cu is as much as 39.2, 36.6, 34.1 and 32.8 K in 0 (saturation), 10, 20 and 30 K subcooled boiling, respectively. Nucleate boiling ensues on Cu at a surface superheat of 11.9, 10.9, 9.5 and 7.5 K in 0 (saturation), 10, 20 and 30 K subcooled boiling, respectively. The saturation nucleate boiling heat flux on porous graphite is 1700% higher than that on Cu at a surface superheat of ∼10 K and decreases exponentially with increased superheat to ∼60% higher near CHF. The CHF values of HFE-7100 on porous graphite of 31.8, 45.1, 55.9 and 66.4 W/cm² in 0 (saturation), 10, 20 and 30 K subcooled boiling, are 60% higher and the corresponding superheats are 25% lower than those on Cu. In addition, the rate of increase in CHF with increased liquid subcooling is 50% higher than that on Cu.     

Flow boiling CHF enhancement with surfactant solutions under atmospheric pressure

Surfactant effect on CHF (critical heat flux) was determined during water flow boiling at atmospheric pressure in closed loop filled with solution of tri-sodium phosphate (TSP, Na₃PO₄ · 12H₂O). TSP was added to the containment sump water to adjust pH level during accident in nuclear power plants. CHF was measured for four different water surfactant solutions in vertical tubes, at different mass fluxes (100–500 kg/m² s) and two inlet subcooling temperatures (50 °C and 75 °C). Surfactant solutions (0.05–0.2%) at low mass flux (～100 kg/m² s) showed the best CHF enhancement. CHF was decreased at high mass flux (500 kg/m² s) compared to the reference plain water data. Maximum increase in CHF was about 48% as compared to the reference data. Surfactant caused a decrease in contact angle associated with an increase of CHF from surfactant addition.     

Effects of pulse width and mass flux on microscale flow boiling under pulse heating

A Pt microheater (140 × 100 μm²) is fabricated on a glass wafer and enclosed in a silicon microchannel of trapezoidal cross section by MEMS technology. With the aid of a high-speed CCD and data acquisition system, subcooled flow boiling phenomena and temperature response on the surface of the microheater under pulse heating are observed and recorded. Experiments are conducted for six pulse widths (50 μs, 100 μs, 200 μs, 600 μs, 1 ms, and 2 ms) under different mass and heat fluxes. With increasing heat flux at a fixed pulse width and different mass fluxes, four flow regimes including single phase, nucleate boiling, film boiling and dry out are identified. Since flow boiling regimes are relatively independent of mass flux, correlation equations based on experimental data for the transitional heat flux of different flow boiling regimes are obtained in terms of pulse width only. It is also found that pulse width and mass flux have little influence on boiling inception time, and the classical analytical solution for the nucleation inception time in terms of heat flux is verified experimentally.     

Study of the critical heat flux condition with water and R-123 during flow boiling in microtubes. Part I: Experimental results and discussion of parametric effects

Extensive experimentation was performed to obtain flow boiling critical heat flux data in single stainless steel microtubes with diameters from 0.286 to 0.700 mm over a wide range of mass fluxes, inlet subcoolings, and exit pressures for two different working fluids (water and R-123). The effect of different operating parameters – mass flux, inlet subcooling, exit quality, heated length and diameter – were assessed in detail (Part I of the paper). The conventional DNB-type behavior is observed in the high subcooled region, and the typical dryout type behavior is seen in the high-quality saturated region when the flow is completely annular. The flow in transitional flow patterns (churn–annular or slug–annular) causes a peculiar increase of CHF with exit quality. Also, the increased void fraction near the saturated region in subcooled boiling results in increased subcooled CHF values. Part II of the paper deals with comparison of data with existing correlations and development of a new correlation to predict the CHF condition in the subcooled liquid region.