Contribution of thin-film evaporation during flow boiling inside microchannels

Flow boiling through microchannels is characterized by nucleation and growth of vapor bubbles that fill the entire channel cross-sectional area. As the bubbles nucleate and grow inside the microchannel, a thin film of liquid or a microlayer gets trapped between the bubbles and the channel walls. The heat transfer mechanism present at the channel walls during flow boiling is studied numerically. It is then compared to the heat transfer mechanisms present during nucleate pool boiling and in a moving evaporating meniscus. Increasing contact angle improved wall heat transfer in case of nucleate boiling and moving evaporating meniscus but not in the case of flow boiling inside a microchannel. It is shown that the thermal and the flow fields present inside the microchannel around a bubble are fundamentally different as compared to nucleate pool boiling or in a moving evaporating meniscus. It is explained why thin-film evaporation is the dominant heat transfer mechanism and is responsible for creating an apparent nucleate boiling effect inside a microchannel.     

Similarities and Differences Between Flow Boiling in Microchannels and Pool Boiling

Recent literature indicates that under certain conditions the heat transfer coefficient during flow boiling in microchannels is quite similar to that under pool boiling conditions. This is rather unexpected, as microchannels are believed to provide significant heat transfer enhancement under single-phase as well as flow boiling conditions. This article explores the underlying heat transfer mechanisms and illustrates the similarities and differences between the two processes. Formation of elongated bubbles and their passage over the microchannel walls have similarities to the bubble ebullition cycle in pool boiling. During the passage of elongated bubbles, the longer duration between two successive liquid slugs leads to wall dryout and a critical heat flux that may be lower than that under pool boiling conditions. A clear understanding of these phenomena will help in overcoming these limiting factors and in developing strategies for enhancing heat transfer during flow boiling in microchannels.     

Effects of heater size and working fluids on nucleate boiling heat transfer

The present study is an experimental investigation of nucleate boiling heat transfer mechanism in pool boiling from wire heaters immersed in saturated FC-72 coolant and water. The vapor volume flow rate departing from a wire during nucleate boiling was determined by measuring the volume of bubbles from the wire utilizing the consecutive-photo method. The effects of the wire size on heat transfer mechanism during a nucleate boiling were investigated, varying 25 μm, 75 μm, and 390 μm, by measuring vapor volume flow rate and the frequency of bubbles departing from a wire immersed in saturated FC-72. One wire diameter of 390 μm was selected and tested in saturated water to investigate the fluid effect on the nucleate boiling heat transfer mechanism. Results of the study showed that an increase in nucleate boiling heat transfer coefficients with reductions in wire diameter was related to the decreased latent heat contribution. The latent heat contribution of boiling heat transfer for the water test was found to be higher than that of FC-72. The frequency of departing bubbles was correlated as a function of bubble diameters.     

Effect of aspect ratio on saturated boiling flow in microchannels with nonuniform heat flux

Microchannel two‐phase flow is an effective cooling method used in microelectronics, in which the heat flux density is unevenly distributed usually. The paper is focused on numerical study the effect of aspect ratio on the flow boiling of microchannels with nonuniform heat flux. The heat source is a three‐dimensional (3D) integrated circuit. 3D microchannel model and volume of fluid method are coupled in numerical simulation. The results show that the aspect ratio has no relationship with the two‐phase pressure drop of the microchannel. It has a certain influence on the distribution of bubble shape. In terms of the heat transfer coefficient, the aspect ratio has a certain influence on a section of the inlet. Due to the nucleate boiling, the convective heat transfer in the remaining areas is the dominant factor and the average heat transfer coefficient is mainly determined by the heat flux at the bottom of the channel.     

On the transition from partial to fully developed subcooled flow boiling

The subject of the present study is to relate the boiling heat transfer process with experimentally observed bubble behaviour during subcooled flow boiling of water in a vertical heated annulus. It presents an attempt to explain the transition from partial to fully developed flow boiling with regard to bubble growth rates and to the time that individual bubbles spend attached to the heater surface.Within the partial nucleate boiling region bubbles barely change in size and shape while sliding a long distance on the heater surface. Such behaviour indicates an important contribution of the microlayer evaporation mechanism in the overall heat transfer rate. With increasing heat flux, or reducing flow rate at constant heat flux, bubble growth rates increase significantly. Bubbles grow while sliding, detach from the heater, and subsequently collapse in the bulk fluid within a distance of 1-2 diameters parallel to the heater surface. This confirms that bubble agitation becomes a leading heat transfer mode with increasing heat flux. There is however, a sharp transition between the two observed bubble behaviours that can be taken as the transition from partial to fully developed boiling. Hence, this information is used to develop a new model for the transition from partial to fully developed subcooled flow boiling.     

Subcooled Flow Boiling in a Minichannel

It has been considered that dry-out occurs easily in boiling heat transfer for a small channel, a mini- or microchannel, because the channel was easily filled with coalescing vapor bubbles. In the present study, the experiments of subcooled flow boiling of water were performed under atmospheric conditions for a horizontal rectangular channel for which the size is 1 mm height and 1 mm width, with a flat heating surface of 10 mm length and 1 mm width placed on the bottom of the channel. The heating surface has a top of copper heating block and is heated by ceramic heaters. In the high heat flux region of nucleate boiling, about 70–80% of the heating surface was covered with a large coalescing bubble and the boiling reached critical heat flux as observed by high-speed video. In the beginning of transition boiling, coalescing bubbles were collapsed to many fine bubbles and microbubble emission boiling was observed at liquid subcooling higher than 30 K. The maximum heat flux obtained was 8 MW/m² (800 W/cm²) at liquid subcooling of higher than 40 K and a liquid velocity of 0.5 m/s. However, the surface temperature was very much higher than that of a centimeter-scale channel. The high-speed video photographs indicated that microbubble emission boiling occurs in the deep transition boiling region.     

Flow and Heat Transfer Characteristics of Liquid Nitrogen in Mini-/Microchannels

Flow and heat transfer of liquid nitrogen in mini-/microchannels have many particular characteristics and are very important for many cooling applications. In this study, the investigation of flow and heat transfer characteristics of liquid nitrogen in mini-/microchannels is presented by summarizing the experimental studies carried out in the author's group. In addition, some recent results about flow and heat transfer of liquid nitrogen in microchannel heat sink are also presented. It is found that small viscosity of liquid nitrogen enables the single-phase liquid flow in mini-/microchannels to be turbulent state, which proves that the classical theory for pressure drop is still valid if the surface roughness of the passage is properly taken into consideration. Experiments of flow boiling of liquid nitrogen are conducted under both adiabatic and diabatic conditions. It is shown that confinement number can be applicable in classifying the heat transfer characteristics of liquid nitrogen in macro- and microchannels. Flow visualization in microchannels at low temperatures poses big challenges on experimental aspects, which have been subtly overcome and clear images have been obtained. The flow patterns and flow regimes of two-phase flow of liquid nitrogen exhibit different features from the room-temperature fluids. Furthermore, three-dimensional (3D) flow visualization by only one high-speed camera is conducted to obtain more detailed information of flow patterns. Finally, the experiments of flow boiling of liquid nitrogen in microchannel heat sink are also presented and discussed.     

Numerical investigation of subcooled flow boiling in segmented finned microchannels

Bubble growth behavior and heat transfer characteristics during subcooled flow boiling in segmented finned microchannels have been numerically investigated. Simulations have been performed for a single row of segmented finned microchannel and predicted results are compared with experimental investigations. Onset of nucleation, formation of bubbles, their growth and movements have been investigated for different values of applied heat flux. Mechanism of bubble expansion without clogging resulting in enhanced heat transfer in segmented finned microchannels has been explained. Temperature and pressure fluctuations during subcooled flow boiling condition have been investigated. It is observed that at high heat flux, thin liquid film trapped between the bubble and channel wall is evaporated leading to localized heating effect. Predicted flow patterns are similar to experimental results. However, simulations over predict the bubble growth rate and heat transfer coefficient.     

Simulation on nucleate boiling in micro-channel

Using the VOF multiphase flow model, numerical simulations are conducted to investigate the nucleate boiling of water in micro-channels. The Marangoni heat transfer through the bubble surface is analyzed, and is compared with the incipient heat flux at the onset of nucleate boiling in micro-channels. The bubble growth in the channel is divided into two stages. At the initial stage, bubble growth is controlled by surface tension, while at the second stage the incipient heat transfer dominated the boiling process. In the results, the full process of bubble generating, growing, departing, combining, and shrinking in the channel is displayed. The simulated results with similar condition are agreed well with some experimental results in references. The method and discussion in the paper are helpful to the investigation of the mechanism of micro-scale two-phase flow and heat transfer.     

Numerical Investigation of Saturated Flow Boiling in Microchannels by the Lattice Boltzmann Method

Bubble formation in saturated flow boiling in 2D microchannels, generated from a microheater under constant wall heat flux or constant wall temperature conditions, is studied numerically based on a newly developed lattice Boltzmann model for liquid-vapor phase change. Simulations are carried out to study effects of inlet velocity, contact angle, and heater size on saturated flow boiling of water under constant wall heat flux conditions. Important information, such as effects of static contact angle on nucleation time and nucleation temperature, which was unable to be obtained by other numerical simulation methods, is obtained. Furthermore, effects of inlet velocity, contact angle, and superheat on nucleate boiling heat transfer in steady flow boiling of water under constant wall temperature conditions are also presented. It is found that the nucleate boiling heat transfer at the microheater is higher if the heater surface is more hydrophilic, because the superheated vapor at the hydrophilic wall has a thinner thermal boundary layer and a larger thermal conductivity.     

Bubble Dynamics,Flow, and Heat Transfer during Flow Boiling in Parallel Microchannels

Significant efforts have recently been made to investigate flow boiling in microchannels, which is considered an effective cooling method for high-power microelectronic devices. However, a fundamental understanding of the bubble motion and flow reversal observed during flow boiling in parallel microchannels is lacking in the literature. In this study, complete numerical simulations are performed to further clarify the boiling process by using the level-set method for tracking the liquid–vapor interface which is modified to treat an immersed solid surface. The effects of contact angle, wall superheat, and the number of channels on the bubble growth, reverse flow, and heat transfer are analyzed.     

Subcooled flow boiling heat transfer of R-134a and the associated bubble characteristics in a vertical plate heat exchanger

Subcooled flow boiling heat transfer characteristics of refrigerant R-134a in a vertical plate heat exchanger (PHE) are investigated experimentally in this study. Besides, the associated bubble characteristics are also inspected by visualizing the boiling flow in the vertical PHE. In the experiment two vertical counterflow channels are formed in the exchanger by three plates of commercial geometry with a corrugated sinusoidal shape of a chevron angle of 60°. Upflow boiling of subcooled refrigerant R-134a in one channel receives heat from the downflow of hot water in the other channel. The effects of the boiling heat flux, refrigerant mass flux, system pressure and inlet subcooling of R-134a on the subcooled boiling heat transfer are explored in detail. The results are presented in terms of the boiling curves and heat transfer coefficients. The measured data showed that the slopes of the boiling curves change significantly during the onset of nucleate boiling (ONB) especially at low mass flux and high saturation temperature. Besides, the boiling hysteresis is significant at a low refrigerant mass flux. The subcooled boiling heat transfer coefficient is affected noticeably by the mass flux of the refrigerant. However, increases in the inlet subcooling and saturation temperature only show slight improvement on the boiling heat transfer coefficient.The photos from the flow visualization reveal that at higher imposed heat flux the plate surface is covered with more bubbles and the bubble generation frequency is substantially higher, and the bubbles tend to coalesce to form big bubbles. But these big bubbles are prone to breaking up into small bubbles as they move over the corrugated plate, producing strong agitating flow motion and hence enhancing the boiling heat transfer. We also note that the bubbles nucleated from the plate are suppressed to a larger degree for higher inlet subcooling and mass flux. Finally, empirical correlations are proposed to correlate the present data for the heat transfer coefficient and the bubble departure diameter in terms of boiling, Froude, Reynolds and Jakob numbers.     

Hydrodynamics and heat transfer prior to onset of nucleate boiling in a rectangular microchannel heat sink

The conjugate heat transfer of flow boiling in a rectangular microchannel heat sink (MCHS) was modelled numerically to investigate the hydrodynamics and thermal responses of flow prior to the onset of nucleate boiling (ONB). Local hydrodynamics and thermal conditions leading to ONB are analysed numerically for different heat flux. The flow patterns of different modes of microconvection and mixed convective flows including the circulating flow, wavy flow and seeping flow were demonstrated and discussed. The numerical study proposes the mechanism leading to the first bubble nucleation which cover the initiation of fluid instability until the ONB. This work provides better understanding of the superheat induced flow instability and the progressive fluid convection under transient heating.     

Fundamental Issues Related to Flow Boiling and Two-Phase Flow Patterns in Microchannels – Experimental Challenges and Opportunities

ABSTRACT

Flow boiling heat transfer in microchannels is used today in many diverse applications. The previous studies addressing the effect of channel size, heat flux, vapor quality, and mass flux on heat transfer during flow boiling are reviewed in the present paper. The relationship between flow characteristics and flow boiling heat transfer was studied experimentally for refrigerant R-C318 at moderate reduced pressures where the contribution of nucleate boiling is decisive. Flow boiling mechanisms were identified using an annular microchannel with transparent outer wall for successive visualization of boiling. The considerable suppression of nucleate boiling heat transfer was observed at transition to annular flow and explained by formation of a liquid flow with thin film and dry spots. A general equation for prediction of two-phase flow boiling heat transfer inside the circular, annular, and rectangular microchannels is proposed and verified using the experimental data. This equation accounts for the nucleate boiling suppression, forced convection, and thin film evaporative heat transfer in the form that allows to distinguish more clearly the contribution of each mechanism of heat transfer under the conditions, when it is predominant. A new approach for prediction of transition to the annular flow is proposed and verified, using the experimental data.     

A Study of Bubble Behavior and Boiling Heat Transfer Enhancement under Electric Field

The effect of a D.C. electric field on nucleate boiling heat transfer for refrigerants, R11, R113, and FC72, was investigated experimentally in a single-tube shell/tube heat exchanger by using the temperature control method of wall superheat. Also the behavior of bubble under nonuniform electric field produced by wire electrodes was studied by numerical calculation. For R11, the electrohydrodynamic (EHD) enhancement for boiling heat transfer was observed for all ranges of wall superheat tested. However, the enhancement in boiling heat transfer disappeared if the wall superheat exceeded 13°C for R113, and no electric field effect on the boiling heat transfer was observed for FC72. An application of approximately 5 kV was enough to eliminate the boiling hysteresis for R11 and R113. Numerical study of the electric field in a single medium has hinted that the bubbles are forced away from the heating surface and toward the electrostatic stagnation point by the dielectrophoretic force. Such modified bubble motion turns out to promote the boiling heat transfer if one uses proper electrode configuration.     

Numerical investigation of the bubble growth in horizontal rectangular microchannels

To explore the mechanism of flow boiling in microchannels, the processes of a single-vapor bubble evaporating and two lateral bubbles merging in a 2D microchannel are investigated. The temperature recovery model based on volume of fluid method is adopted to perform the flow boiling phenomena. The effects of wall superheat, Reynolds number, contact angle, surface tension, and two-bubble merger on heat transfer are discussed. Wall superheat dominates the bubble growth and is roughly proportional to wall heat flux. The update of velocity and temperature fields’ distribution in the channel increases with increasing inflow Reynolds number, which improves the wall heat flux markedly. Besides, the area of thin liquid film between the wall and the bubble is enlarged by reducing the contact angle, thus, expanding the wall heat flux several times compared with the single-phase cross section. However, variation of surface tension (0.0589, 0.1?N/m) is found to be insignificant.     

Saturated flow boiling heat transfer of refrigerant R-410A in a horizontal annular finned duct

An experiment is conducted here to investigate the saturated flow boiling heat transfer characteristics of ozone friendly refrigerant R-410A in a horizontal annular finned duct. Meanwhile the associated bubble characteristics in the duct are also inspected from the flow visualization. The experimental data are presented in terms of saturated flow boiling curves, boiling heat transfer coefficients and flow photos. In addition, empirical correlation equations for the saturated flow boiling heat transfer coefficient and mean bubble departure diameter are proposed. The saturated flow boiling curves show that boiling hysteresis is insignificant in the flow and the wall superheat needed for the onset of nucleate boiling is slightly affected by the refrigerant mass flux. Besides, the boiling curves are mainly affected by the imposed heat flux and refrigerant mass flux. Moreover, the measured saturated flow boiling heat transfer coefficient increases with the imposed heat flux and refrigerant mass flux. Furthermore, at a higher refrigerant mass flux the departing bubbles are smaller.     

Correlation of the Flow Pattern and Flow Boiling Heat Transfer in Microchannels

Flow boiling in microchannels is characterized by the considerable influence of capillary forces and constraint effects on the flow pattern and heat transfer. In this article we utilize the features of gas–liquid flow patterns in rectangular microchannels under adiabatic conditions to explain the regularities of refrigerants flow boiling heat transfer. The flow-pattern maps for the upward and horizontal nitrogen–water flow in a microchannel with the size of 1500 × 720 μm were determined via dual-laser flow scanning and compared with corrected Mishima and Ishii prediction. Flow boiling heat transfer was studied for vertical and horizontal microchannel heat sink with similar channels using refrigerants R-21 and R-134a. The data on local heat transfer coefficients were obtained in the range of mass flux from 33 to 190 kg/m²-s, pressure from 1.5 to 11 bar, and heat flux from 10 to 160 kW/m². The nucleate and convective flow boiling modes were observed for both refrigerants. It was found that heat transfer deterioration occurred for annular flow when the film thickness became small to suppress nucleate boiling. The mechanism of heat transfer deterioration was discussed and a model of heat transfer deterioration was applied to predict the experimental data.     

Falling film hydrodynamics of black liquor under evaporative conditions

Bubbles were observed in a thin, evaporating, falling film of black liquor (a fluid mixture generated during the pulp production) on the exterior wall of a research evaporator. Because the presence of bubbles could not be explained by nucleate boiling, a combination of turbulent vapour entrainment and effects due to surface-active compounds – surfactants – is proposed. Black liquor contains numerous surfactants, which are likely to enhance bubble formation and stabilization in the fluid and on the film interface. One observed important effect of bubble formation was fluid loss due to bubble-bursting aerosolization (sputtering). Also, bubbles and bubble processes probably alter the film velocity-profile and heat transfer resistance, thereby affecting heat transfer across the film and hence evaporator efficiency.     

Bubble sweeping and jet flows during nucleate boiling of subcooled liquids

An experimental investigation was conducted to investigate nucleate boiling on a very fine heating wire. By using zoom routine and CCD camera system, the dynamical process of nucleate boiling was visually observed. Sweeping bubbles and several modes of jet flows were described and discussed. For some cases, big bubbles, small bubbles, sweeping bubbles and jet flows coexisted in boiling system, and greatly enhanced heat transfer. These phenomena are quite different from usual observation of nucleate boiling. In this paper, the process of bubble sweeping phenomenon is described in detail and the effect induced by sweeping bubbles is argued. And also, several jet flows are illustrated and discussed, as well as the interaction between bubble sweeping and jet flows.