Bubble dynamics at boiling incipience in subcooled upward flow boiling

Bubble dynamics in water subcooled flow boiling was investigated through visualization using a high-speed camera. The test section was a vertical rectangular channel, and a copper surface of low contact angle was used as a heated surface. Main experimental parameters were the pressure, mass flux and liquid subcooling. Although all the experiments were conducted under low void fraction conditions close to the onset of nucleate boiling, no bubbles stayed at the nucleation sites at which they were formed. Depending on the experimental conditions, the following two types of bubble behavior were observed after nucleation: (1) lift-off from the heated surface followed by collapsing rapidly in subcooled bulk liquid due to condensation, and (2) sliding along the vertical heated surface for a long distance. Since the bubble lift-off was observed only when the wall superheat was high, the boundary between the lift-off and the sliding could be determined in terms of the Jakob number. Based on the present experimental results, discussion was made for the possible mechanisms governing the bubble dynamics.     

On the rise paths of single vapor bubbles after the departure from nucleation sites in subcooled upflow boiling

A photographic study was carried out for the subcooled flow boiling of water to elucidate the rise characteristics of single vapor bubbles after the departure from nucleation sites. The test section was a transparent glass tube of 20 mm in inside diameter and the flow direction was vertical upward; liquid subcooling was parametrically changed within 0–16 K keeping system pressure and liquid velocity at 120 kPa and 1 m/s, respectively. The bubble rise paths were analyzed from the video images that were obtained at the heat flux slightly higher than the minimum heat flux for the onset of nucleate boiling. In the present experiments, all the bubbles departed from their nucleation sites immediately after the inception. In low subcooling experiments, bubbles slid upward and consequently were not detached from the vertical heated wall; the bubble size was increased monotonously with time in this case. In moderate and high subcooling experiments, bubbles were detached from the wall after sliding for several millimeters and migrated towards the subcooled bulk liquid. The bubbles then reversed the direction of lateral migration and were reattached to the wall at moderate subcooling while they collapsed due to the condensation at high subcooling. It was hence considered that the mechanisms of the heat transfer from heated wall and the axial growth of vapor volume were influenced by the difference in bubble rise path. It was observed after the inception that bubbles were varied from flattened to more rounded shape. This observation suggested that the bubble detachment is mainly caused by the change in bubble shape due to the surface tension force.     

A BUBBLE MECHANISTIC MODEL FOR SUBCOOLED BOILING FLOW PREDICTIONS

Population balance equations combined with a three-dimensional two-fluid model are employed to predict subcooled boiling flow at low pressure in a vertical annular channel. The MUSIG (Multiple-Size-Group) model implemented in CFX4.4 is extended to account for the wall nucleation and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements is made for the void fraction, bubble Sauter diameter, interfacial area concentration, bubble population density, and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Good agreement is achieved with the local radial void fraction, bubble Sauter diameter, interfacial area concentration, bubble population density, and liquid velocity profiles against measurements. However, further improvement is needed for the accurate prediction of the vapor velocity using the present bubble mechanistic model. A proposal to include an algebraic slip model to account for bubble separation in the MUSIG boiling model is presented.     

Fundamental consideration of wall heat partition of vertical subcooled boiling flows

Improved wall heat flux partitioning accounting sliding bubbles and a mechanistic model that incorporates the fundamental consideration of bubble frequency during low-pressure subcooled flow boiling is presented. A model considering the forces acting on departing bubbles at the heated surface is employed. Coupled with a three-dimensional two-fluid and population balance equations based on the modified MUSIG (MUltiple-SIze-Group) model, the behavior of an upward forced convective subcooled boiling flows in a vertical annular channel is simulated. Comparison of model predictions against local and axial measurements (heat fluxes ranged from 152.9 to 705.0 kW/m²) is made for the void fraction, Sauter mean bubble diameter and interfacial area concentration covering a range of different mass and heat fluxes and inlet subcoolings. Good agreement is achieved between the predicted and measured profiles. Reasonable agreement with recent experimental measurements is also attained for the predicted growth and waiting times of bubble frequency at particular local wall superheat and subcooling temperatures.     

Micro‐bubble emission boiling from horizontal and vertical surfaces to subcooled parallel flow water

Heat removal of more than 10 MW/m² in heat flux has been required in high‐heat‐generation equipment in nuclear fusion reactors. In some conditions of water subcooling and velocity, there appears an extraordinary high heat flux boiling in the transition boiling region. This boiling regime is called micro‐bubble emission boiling (MEB) because many micro‐bubbles are spouted from the heat transfer surface accompanying a huge sound. The study intent is to obtain heat transfer performance of MEB in horizontal and vertical heated surfaces to parallel flow of subcooled water, comparing with CHF of this system. Three types of MEB with different heat transfer performance and bubble behavior are observed according to the flow velocity and liquid subcooling. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(2): 130–140, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10077     

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.     

The influence of bubble size on void fraction distribution in subcooled flow boiling at low pressure

The effect of bubble size on void fraction distribution in subcooled flow boiling in a vertical annular channel at low pressure is studied numerically. It is found that a simple linear formula used by Anglart and Nylund [1] and adopted in Reference [2] for calculation of bubble size as a function of local subcooling lacks a physical and experimental basis limiting the general application of the model for predicting subcooled flow boiling. A bubble size correlation proposed by Zeitoun and Shoukri [3] has been employed in this study. The predictions of void fraction profiles and the bubble size distributions, after incorporating the above bubble size correction, show very good agreement with the experimental data.     

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.     

Liquid–vapor structure near heating surface at high heat flux in subcooled pool boiling

Liquid–vapor behavior close to a heating surface was measured using two conductance probes with tip diameters smaller than 5 μm. Measurements were carried out for water boiling on an upward-facing copper surface under subcooling from 0 to 30 K. The probe signals and the void fraction distributions showed that there is little difference in the liquid–vapor structure beneath large vapor masses in saturated and subcooled boiling, that a macrolayer remains on the heating surface, and that in subcooled boiling it does not dry out even at heat fluxes far higher than CHF for saturated boiling. The thickness of the macrolayer forming beneath large vapor masses was determined from the location where the probe signals corresponding to the large vapor masses disappear. It was found that the thicknesses of the macrolayer formed in subcooled boiling are comparable to or thicker than those near the CHF in saturated boiling, and it is considered that this is most likely to be one of the causes why the CHF increases with the increasing subcooling.     

Subcooled flow boiling heat transfer and associated bubble characteristics of R-134a in a narrow annular duct

Experiments are conducted here to investigate how the channel size affects the subcooled flow boiling heat transfer and associated bubble characteristics of refrigerant R-134a in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. From the measured boiling curves, the temperature undershoot at ONB is found to be relatively significant for the subcooled flow boiling of R-134a in the duct. The R-134a subcooled flow boiling heat transfer coefficient increases with a reduction in the gap size, but decreases with an increase in the inlet liquid subcooling. Besides, raising the imposed heat flux can cause a substantial increase in the subcooled boiling heat transfer coefficient. However, the effects of the refrigerant mass flux and saturated temperature on the boiling heat transfer coefficient are small in the narrow duct. Visualization of the subcooled flow boiling processes reveals that the bubbles are suppressed to become smaller and less dense by raising the refrigerant mass flux and inlet subcooling. Moreover, raising the imposed heat flux significantly increases the bubble population, coalescence and departure frequency. The increase in the bubble departure frequency by reducing the duct size is due to the rising wall shear stress of the liquid flow, and at a high imposed heat flux many bubbles generated from the cavities on the heating surface tend to merge together to form big bubbles. Correlation for the present subcooled flow boiling heat transfer data of R-134a in the narrow annular duct is proposed. Additionally, the present data for some quantitative bubble characteristics such as the mean bubble departure diameter and frequency and the active nucleation site density are also correlated.     

Temperature measurements near a heating surface at high heat fluxes in subcooled pool boiling

In previous papers (Int J Heat Mass Transfer, 2008;50:3481–3489, 2009;52: 814–821), the authors conducted measurements of liquid–vapor structures in the vicinity of a heating surface for subcooled pool boiling on an upward‐facing copper surface by using a conducting probe method. We reported that the macrolayer dryout model is the most appropriate model of the CHF and that the reason why the CHF increases with increasing subcooling is most likely that a thick macrolayer is able to form beneath large vapor masses and the lowest heat flux of the vapor mass region shifts towards the higher heat flux. To develop a mechanistic model of the CHF for subcooled boiling, therefore, it is necessary to elucidate the effects of local subcooling on boiling behaviors in the vicinity of a heating surface. This paper measured local temperatures close to a heating surface using a micro‐thermocouple at high heat fluxes for water boiling on an upward‐facing surface in the 0 to 40 K range of subcooling. A value for the effective subcooling, defined as the local subcooling during the period while vapor masses are being formed was estimated from the detected bottom peaks of the temperature fluctuations. It was established that the effective subcooling adjacent to the surface remains at considerably lower values than the bulk liquid subcooling. This suggests that, from nucleation to coalescence, the subcooling of a bulk liquid has a smaller effect on the behavior of primary bubbles than the extent of the subcooling would appear to suggest. An empirical correlation of the effective subcooling is proposed to provide a step towards quantitative modeling of the CHF for subcooled boiling. © 2009 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.20277     

Unsteady Natural Convective Boiling in a Narrow Vertical Channel

The purpose of this work is to describe the two-phase flow structure and heat transfer of unsteady natural convective boiling in a narrow vertical channel. The experiments are performed with saturated n-pentane at a pressure of 1 bar. An unheated plate is placed parallel to the heating surface and lateral sides are closed. The distance between the heated surface and the confinement plate is 800 μm. Void fraction measurements are performed using capacitive sensors. The void fraction increases with heat flux and reaches a maximum of 0.80 in the mid-height of the channel when the heat flux is equal to 90% of the critical heat flux. Flow observations using a high-speed video camera show an unsteady thermo-hydraulic behavior. The frequency of the cycles increases with the wall temperature during nucleate and transient boiling. Local velocities of the bubble meniscus developing within the confined space are determined during the boiling cycles. The time-averaged liquid flow rate increases significantly with heat flux and reaches a maximum for heat flux close to the critical heat flux.     

Experimental investigation of bubble behaviors in subcooled flow boiling

The experimental investigation on vapor bubble growth is performed for analyzing subcooled boiling in a vertical annular channel with inner heating surface and upward water flow under atmospheric pressure. Bulk liquid mass flux ranges from 79 kg/m2s to 316 kg/m2s, and subcooling is from 40 K to 60 K. The bubble behaviors from inception to collapse are captured by High-speed photography. The performance of bubble growth recorded by the high-speed photography is given in this paper. The bubble behaviors, effect of the bubble slippage on the heat transfer, and various forces acting on the bubble are discussed.     

Axial developments of interfacial area and void concentration profiles in subcooled boiling flow of water

Axial developments of the local void fraction, interfacial area concentration and bubble Sauter mean diameter were measured in subcooled boiling flow of water in a vertical internally heated annulus using the double-sensor conductivity probe technique. Measurements were performed under varying conditions of heat flux, inlet liquid velocity and inlet liquid temperature. A total of 10 data sets were acquired. Based on these measurements with the previous data obtained in the present test loop, the influence of flow condition on the profiles of local two-phase flow parameters was discussed. The measured average void fraction and interfacial area concentration were compared with the predictions by existing correlations for drift-flux parameters and interfacial area concentration. Also, the recently proposed bubble layer thickness model in subcooled boiling was evaluated for the measurement data.     

Numerical analysis of the dynamics of moving vapor bubbles

Bubbles have been observed rapidly sweeping along very fine heated wires during subcooled nucleate boiling with jet flows emanating from the tops of the vapor bubbles. This paper analyzes the physical mechanisms driving the bubble and the jet flows from the tops of these moving bubbles. The flows are analyzed by numerically solving the governing equations for the velocity and temperature distributions around the bubble and the heated wire as the bubble moves along the wire. The bubble motion is due to the non-uniform temperature distribution in the liquid and in the wire caused by the bubble as it moves along the wire. The flow is driven by the horizontal Marangoni flow induced by the temperature difference across the bubble which thrusts the bubble forward. Comparisons with experimental observations suggest that the condensation heat transfer at the bubble interface is restricted by non-condensable gases that increases the surface temperature gradient and the resulting Marangoni flow.     

Vapor bubble dynamics and heat transfer characteristics over the spherical surface during saturated pool boiling

The present study deals with numerical investigations of the boiling phenomena over a spherical surface at different degrees of superheat (ΔT), varying from 10 to 500 K. Various phenomena like vapor sliding, bubble formation, pinch-off, induced vorticity have been illustrated for a deep understanding of the boiling process over a spherical surface. The effect of the degree of superheat on the bubble pinch-off time and volume is also investigated. Further, reported the spatial observation of vapor sliding and retention over the surface with time scale, overall and average characteristics. The fast Fourier transform of a spaced average void fraction of liquid and Nusselt number showed the dominance of film boiling with respect to the degree of superheat. As the degree of superheat increases, the vapor generation rate also increases, which produces a more vapor–liquid interface. Further, with an increase in the degree of superheat, the vapor generation progression shifted from linear to nonlinear patterns. A sphere with ΔT = 500 K generated 32.59 times more vapor than a sphere with ΔT = 10 K. It is found that the vapor generation is dependent on the degree of superheating and exposed time for heating. Thus, a correlation and artificial neural network model have been developed to predict vapor generation during boiling over the spherical surface as a function of time and degree of superheat.     

Critical heat flux for subcooled flow boiling in micro-channel heat sinks

Critical heat flux (CHF) was measured and examined with high-speed video for subcooled flow boiling in micro-channel heat sinks using HFE 7100 as working fluid. High subcooling was achieved by pre-cooling the working fluid using a secondary low-temperature refrigeration system. The high subcooling greatly reduced both bubble departure diameter and void fraction, and precluded flow pattern transitions beyond the bubbly regime. CHF was triggered by vapor blanket formation along the micro-channel walls despite the presence of abundant core liquid, which is consistent with the mechanism of Departure from Nucleate Boiling (DNB). CHF increased with increasing mass velocity and/or subcooling and decreasing hydraulic diameter for a given total mass flow rate. A pre-mature type of CHF was caused by vapor backflow into the heat sink’s inlet plenum at low mass velocities and small inlet subcoolings, and was associated with significant fluctuations in inlet and outlet pressure, as well as wall temperature. A systematic technique is developed to modify existing CHF correlations to more accurately account for features unique to micro-channel heat sinks, including rectangular cross-section, three-sided heating, and flow interaction between micro-channels. This technique is shown to be successful at correlating micro-channel heat sink data corresponding to different hydraulic diameters, mass velocities and inlet temperatures.     

Effects of jet pattern on two-phase performance of hybrid micro-channel/micro-circular-jet-impingement thermal management scheme

This paper explores the two-phase cooling performance of a hybrid cooling scheme in which a linear array of micro-jets deposits liquid gradually along each channel of a micro-channel heat sink. The study also examines the benefits of utilizing differently sized jets along the micro-channel. Three micro-jet patterns, decreasing-jet-size (relative to center of channel), equal-jet-size and increasing-jet-size, were tested using HFE 7100 as working fluid. It is shown feeding subcooled coolant into the micro-channel in a gradual manner greatly reduces vapor growth along the micro-channel. Void fraction increased between jets but decreased sharply beneath each jet, creating a repeated pattern of growth followed by coalesce, and netting only a mild overall increase in void fraction along the flow direction with predominantly liquid flow at outlet. Unlike most flow boiling situations, where pressure drop increases with increasing heat flux, pressure drop in the hybrid configurations actually decreased and reached a minimum just before CHF. This behavior is closely related to the low void fraction and predominantly liquid flow. Pressure drop in the two-phase region is highest for the equal-jet-size pattern, followed by the decreasing-jet-size and increasing-jet-size patterns, respectively. Low void fraction increased the effectiveness of the hybrid cooling schemes in utilizing bulk liquid subcooling and therefore helped achieve high CHF values. The decreasing-jet-size pattern, which had the highest outlet subcooling, achieved the highest CHF. A single correlation was constructed for the three jet patterns, which relates the two-phase heat transfer coefficient to heat flux and wall superheat.     

Pulsed photothermal heating of the media during bubble generation around gold nanoparticles

Pulsed laser induced heating of the transparent media with light-absorbing nanoparticles (NP) in water and in living cells was experimentally studied with photothermal microscopy and thermal lens methods. At low laser-induced temperatures of the NPs and without vapor bubble generation the media was heated by NPs due to thermal diffusion. At higher laser-induced temperatures of the NPs that have caused the generation of vapor bubbles around NPs no heating of the media outside the bubbles and after their collapse was observed. This effect of “cold” generation of the bubbles is associated with NP properties and was not observed during bubble generation in homogeneously heated light-absorbing media under equal conditions.     

A visualized study of micro-bubble emission boiling

Experiments were designed, and experimental equipment was built, to study the characteristics of micro-bubble emission boiling (MEB) in water contacting a copper heating surface 10 mm in diameter. The behavior of bubbles on the heating surface was captured by a high-speed video camera. The results of these experiments indicated that after subcooling exceeded 25 K, MEB occurred and was accompanied by the emission of numerous extremely small bubbles. During the initial stage of MEB, two different bubble behaviors were observed: a film of vapor on the heating surface expanded and shrank periodically, emitting micro-bubbles, and the film of vapor expanded unevenly before condensing or collapsing into many micro-bubbles. During fully developed MEB, the film of vapor exhibited irregular changes at its surface and partially collapsed in several milliseconds. Nearly simultaneously, a new vapor film layer formed on the heating surface. MEB never occurred during water subcooling when the heating surface was embedded 0.5 mm within a ceramic thimble.