Bubble growth and horizontal coalescence in saturated pool boiling on a titanium foil,investigated by high-speed IR thermography

Growth of an isolated bubble and horizontal coalescence events between bubbles of dissimilar size were examined during pool nucleate boiling of water on a horizontal, electrically-heated titanium foil 25 μm thick. Wall temperature measurements on the back of the foil by high-speed IR camera, synchronized with high-speed video camera recordings of the bubble motion, improved the temporal and spatial resolution of previous observations by high-speed liquid crystal thermography to 1 ms and 40 μm, respectively, leading to better detailed maps of the transient distributions of wall heat flux. The observations revealed complex behaviour that disagreed with some other observations and current modelling assumptions for the mechanisms of heat transfer over the wall contact areas of bubbles and interactions between bubbles. Heat transfer occurred from the entire contact area and was not confined to a narrow peripheral triple-contact zone. There was evidence of an asymmetrical interaction between bubbles before coalescence. It was hypothesised that a fast-growing bubble pushed superheated liquid under a slow-growing bubble. Contact of this liquid with regions of the wall that had been pre-cooled during bubble growth caused local reductions in the wall heat flux. During coalescence, movement of liquid under both bubbles caused further changes in the wall heat flux that also depended on pre-cooling. Contraction of the contact area caused a peripheral reduction in the heat flux and there was no evidence of a large increase in heat flux during detachment. Boiling on very thin foils imposes special conditions. Sensitivity to the thermal history of the wall must be taken into account when applying the observations and hypotheses to other conditions.     

On the Relation between Nucleation Site Density and Critical Heat Flux of Pool Boiling

ABSTRACT

It is traditionally accepted that the critical heat flux (CHF) decreases with increasing nucleation site density (NSD). However, such a CHF-NSD relation was no longer observed in the BETA-B experiment performed on nano-film heaters; instead the increase of NSD resulted in a gain in CHF. To address this seeming contradiction in the relation between critical heat flux and nucleation site density, the present work employed probabilistic analysis to reveal the different tendencies. A concept of effective NSD was proposed, which concerns the active nucleation sites appear within a bubble lifetime, and the resulting bubbles have the chance of direct interaction. We assumed that the boiling crisis on a heater surface is mainly induced by two mechanisms: dry spot expanding in isolated bubble regime for low-NSD surface, coalescence of dry spots under multiple bubbles in fully developed nucleate boiling regime for high-NSD surface, or a combination of the two in the transition regime for medium-NSD surface. Accordingly, we estimated the critical heat flux of each boiling regime at which the boiling crisis occurs. The result indicated that there is a threshold of nucleation site density below which the increase of NSD is contributing to CHF enhancement, while the trend is inverted beyond the threshold.     

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.     

Dynamic characteristics of transient boiling on a square platinum microheater under millisecond pulsed heating

A series of experimental investigations of boiling incipience and bubble dynamics of water under pulsed heating conditions for various pulse durations ranging from 1 ms to 100 ms were conducted. Using a very smooth square platinum microheater, 100 μm on a side, and a high-speed digital camera, the boiling incipience was observed and investigated as a function of the bulk temperature of the microheater, pulse power level, and pulse duration. Given a specific pulse duration, for low pulse power levels, there would be no bubble nucleation or bubble mergence, for moderate pulse power levels, individual bubbles generated on the heater merged to form a single large bubble, while for high pulse power levels, the rapid growth of the individual bubbles and subsequent bubble interaction, resulted in a reduction in bubble coalescence into a single larger bubble, referred to as bubble splash. The transient heat flux range at which bubble coalescence occurs was identified experimentally, along with the temporal variations of bubble size, bubble interface velocity and interface acceleration.     

Critical Heat Flux of Boiling Heat Transfer in a Moderate Narrow Space

INTanDUCTI0NBoilingheattransferandcriticalheatflux(CHF)inaconfinednarrowspacehavebeenstudiedexperi-melltallybyanumberofinvestigatorsinthepastfewdecades.However,thereisnoanypopularlyacceptedmodelintheheattransferinnarrowspaceboiling,althoughsomepopularknowledgeabouttheboilingheattransferinthenarrowspacehavebeenacceptedbymanyresearchers.Theknowledgecanbecon-cludedasthatthenucleateboilingheattransferisenhancedatlowheatfluxregionanddeterioratedathighheatfiuxregi0nespeciallyatCHF.Theenhanceme…     

Time and space resolved wall temperature and heat flux measurements during nucleate boiling with constant heat flux boundary conditions

The lack of time and space resolved measurements under nucleating bubbles has complicated efforts to fully explain pool-boiling phenomena. In this work, time and space resolved temperature and heat flux distributions under nucleating bubbles on a constant heat flux surface were obtained using a 10 × 10 microheater array with 100 μm resolution along with high-speed images. A numerical simulation was used to compute the substrate conduction, which was then subtracted from the heater power to obtain the wall-to-liquid heat transfer. The data indicated that most of the energy required for bubble growth came from the superheated layer around the bubble. Microlayer evaporation and contact line heat transfer accounted for not more than 23% of the total heat transferred from the surface. The dominant heat transfer mechanism was transient conduction into the liquid during bubble departure. Bubble coalescence was not observed to transfer a significant amount of heat.     

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.     

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.     

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.     

微纳结构铜表面低液位池沸腾实验研究

采用两步电镀法,在改变电流密度的情况下制备出具有不同微纳结构和润湿特性的A、B两个表面,并应用于低液位饱和池沸腾的实验研究中.通过与铜表面对比,发现两个表面在低热流密度情况下,传热系数要高于铜表面,但液位降低时传热系数提升幅度较小,原因在于铜表面沸腾气泡较大,液位降低气泡脱离有很大影响,而表面A、B沸腾气泡较小,液位降...     

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

An experiment is conducted here to investigate how the channel size affects the subcooled flow boiling heat transfer and the associated bubble characteristics of refrigerant R-407C in a horizontal narrow annular duct with the gap of the duct fixed at 1.0 and 2.0 mm. The measured boiling curves indicate that the temperature overshoot at ONB is relatively significant for the subcooled flow boiling of R-407C in the duct. Besides, the subcooled flow boiling heat transfer coefficient increases with a reduction in the duct gap, but decreases with an increase in the inlet liquid subcooling. Moreover, raising the heat flux imposed on the duct can cause a significant increase in the boiling heat transfer coefficients. However, the effects of the refrigerant mass flux and saturated temperature on the boiling heat transfer coefficient are slighter. Visualization of the subcooled flow boiling processes in the duct reveals that the bubbles are suppressed to become smaller and less dense by raising the refrigerant mass flux and inlet subcooling. Raising the imposed heat flux, however, produces positive effects on the bubble population, coalescence and departure frequency. Meanwhile, the present heat transfer data for R-407C are compared with the R-134a data measured in the same duct and with some existing correlations. We also propose empirical correlations for the present data for the R-407C subcooled flow boiling heat transfer and some quantitative bubble characteristics such as the mean bubble departure diameter and frequency and the active nucleation site density.     

Fluid flow and heat transfer characteristics of low temperature two-phase micro-channel heat sinks – Part 2. Subcooled boiling pressure drop and heat transfer

This second part of a two-part study explores the performance of a new cooling scheme in which the primary working fluid flowing through a micro-channel heat sink is indirectly cooled by a refrigeration cooling system. The objective of this part of study is to explore the pressure drop and heat transfer characteristics of the heat sink. During single-phase cooling, pressure drop decreased with increasing heat flux because of decreased liquid viscosity. However, pressure drop began increasing with increasing heat flux following bubble departure. These opposite trends produced a minimum in the variation of pressure drop with heat flux. Increasing liquid subcooling decreased two-phase pressure drop because of decreased void fraction caused by strong condensation at bubble interfaces as well as decreased likelihood of bubble coalescence. It is shown macro-channel subcooled boiling pressure drop and heat transfer correlations are unsuitable for micro-channel flows. However, two new modified correlations produced good predictions of the present heat transfer data.     

Bubble Dynamics and Heat Transfer During Pool Boiling on Wettability Patterned Surfaces

Surfaces with spatial wettability patterns have been proven to enhance heat transfer coefficient and critical heat flux in pool boiling. To understand the physical mechanism behind this phenomenon and obtain the correlation among some critical parameters (bubble departure frequency, bubble size, nucleation site density, surface tension), pool boiling experiments were conducted. A Pyrex glass with a layer of indium-tin-oxide was used as the substrate. Hydrophobic patterns will serve as nucleation sites. Experiments were conducted in deionized water under atmospheric pressure at a relatively low heat flux. The processes of nucleation, growth, and departure of individual bubbles were visualized by using a high speed camera through the bottom of the heater surface. It has been found that the patterned surface performed the best in heat transfer for subcooled pool boiling when compared with hydrophilic and hydrophobic surfaces. The nucleation site density of the biphilic surface was much higher, when compared with that of the homogeneous surface. The individual bubbles always nucleate on the edge of the hydrophobic and hydrophilic area, and then move onto the hydrophobic pattern. Most of the individual bubbles detach from the wettability patterned surface in the diameter range from 300 µm to 450 µm (around 77.3%). The bubble departure periods scatter in the range from 80 ms to 1500 ms.     

Heat transfer during pool boiling based on evaporation from micro and macrolayer

An analytical model of heat transfer based on evaporation from the micro and macrolayers to the vapor bubble during pool boiling is developed. Evaporation of microlayer and macrolayer during the growth of individual bubbles is taken care of by using temporal and spatial variation of temperature in the liquid layer. Change of bubble shape during the entire cycle of bubble growth and departure is meticulously considered to find out the rate of heat transfer from the solid surface to the boiling liquid. Continuous boiling curve is developed by considering the bubble dynamics and decreasing thickness of liquid layer along with the increase of dry spot radius. Transient variation of macrolayer and microlayer thickness is predicted along with their effect on CHF. Present model exhibits a good agreement with reported experimental data as well as theories.     

Numerical study of patterns and influencing factors on flow boiling in vertical tubes by thermal LBM simulation

Heat transfer and flow pattern of flow boiling in vertical tube are investigated numerically based on the phase-change Lattice Boltzmann method (LBM) which includes an improved pseudo-potential LB model and a thermal LB model. Two-dimensional numerical simulations are carried out under constant heat flux conditions for the first time. The processes of growth, slippage, detachment and coalescence of the bubbles are captured to verify the correctness of the model. The effects of gravity, contact angle and wall superheat on bubble departure diameter and nucleation waiting time are illustrated. The multiple flow patterns, single-phase flow, bubble flow, slug flow and DNB have been illustrated with the behaviors of bubble nucleation, growth, departure, and coalescence. Some basic features of flow boiling have been clearly observed in the simulation. The influences of several factors such as heat flux, Reynolds number, the width of flow channel, and the width of nucleation point on flow boiling especially on the point of DNB are investigated. The numerical results show that the DNB could be avoided by reducing the heating density, increasing the Reynolds number, increasing the width of the tube and reducing the heating concentration.     

An experimental investigation of bubble nucleation of a refrigerant in pressurized boiling flows

An experimental investigation is performed to determine the effect of system pressure and heat flux on flow boiling and associated bubble characteristics of a refrigerant in a narrow vertical duct. A high-pressure flow boiling test loop was built and TLC (thermo-chromic liquid crystal) was applied to the back of the heater foil for high resolution and accurate measurement of heater surface temperature. Refrigerant R-134a is used as the test fluid at different pressures ranging from 690 to 827 kPa and different heat fluxes to quantify their influence in bubble characteristics such as bubble nucleation, growth, departure, and coalescence. Two synchronized high resolution and high-speed cameras are used to simultaneously capture TLC images as well as bubbling activities at high frame rates. By varying flow rate and system pressure, TLC and bubble images were captured and analyzed. Results show that the bubble generation frequency and size increase with heat flux. An increase in pressure from 690 to 827 kPa increased the bubble frequency and size by about 32 Hz and 20 μm, respectively. Bubble coalescence was also observed after departure from the nucleation site.     

The Influence of System Pressure on Bubble Coalescence in Nucleate Boiling

Boiling is one of the most effective heat transfer mechanisms. In spite of a long time of research, the physical fundamentals are still not sufficiently understood. Pursuing the objective to predict heat transfer based on physical and geometrical properties, experimental and numerical investigations are conducted at the institute of the authors. The focus of the presented research is the coalescence of two single bubbles under varying pressure conditions. In the experiment a thin stainless-steel foil is used as a Joule heater. The experiments were performed in a pressure range of 300–1000 mbar using FC72 as working fluid. Two types of heaters with a distance between two artificial nucleation sites of 300 μm (type 3) and 500 μm (type 5) were used. The experimental results indicate a strong dependence of the occurrence of bubble coalescence on pressure. For the type 5 heater, a Gaussian distribution for the coalescence frequency when plotted over pressure is observed. Experimental results with the type 3 heater show a similar distribution of the frequency with a shifted maximum. Further, it is shown that during bubble coalescence a small droplet can remain inside the bubble and enhance the heat transfer, which is attributed to an additional thin film region. The formation of this remaining droplet is sensitive to system pressure. Numerical investigations of bubble coalescence were conducted with the computational fluid dynamics (CFD) software OpenFOAM. In OpenFOAM, dynamic mesh handling allows high spatial resolution at the phase boundary, which is captured with the volume-of fluid method. Evaporation and a subgrid microscale model were implemented in the flow solver to account for evaporation at the phase boundary and the three-phase contact line. The results show a strong dependence of bubble dynamics and coalescence on contact angle and bubble growth rate. Although it was possible to observe the creation of the residual droplet, more effort needs to be put into finding appropriate initial conditions.     

Experimental study and heat transfer analysis on the boiling of saturated liquid nitrogen under transient pulsed laser irradiation

The boiling behavior of the liquid nitrogen (LN2) under the transient high heat flux urgently needs to be researched systematically. In this paper, the high power short pulse duration laser was used to heat the saturated LN2 rapidly, and the high-speed photography aided by the spark light system was employed to take series of photos which displayed the process of LN2's boiling behavior under such conditions. Also, a special temperature measuring system was applied to record the temperature variation of the heating surface. The experiments indicated that an explosive boiling happened within LN2 by the laser heating, and a conventional boiling followed up after the newly-defined changeover time. By analyzing the temperature variation of the heating surface, it is found that the latent heat released by the crack of the bubbles in the bubble cluster induced by the explosive boiling is an important factor that greatly influences the boiling heat transfer mechanism.     

Bubble characteristics in time periodic evaporation flow of R-134a in a narrow annular pipe due to heat flux oscillation

In this work, bubble characteristics of periodic evaporation flow with refrigerant R-134a in a horizontal narrow annular pipe were examined experimentally in details. Attention is focused on the time periodic evaporation flow characteristics affected by the mean levels, amplitudes, and periods of the heat flux oscillation. The photos of the R-134a time periodic evaporating flow taken from the duct side are presented to show the change of the dominant two-phase flow pattern in the duct with the experimental parameters. The results show that at the low vapor quality, the bubbles get smaller with time and become less crowded in the duct in the first half of the cycle in which the R-134a heat flux decreases. The changes of the bubble characteristics with the instantaneous heat flux become more pronounced for an increase in the amplitude of the heat flux oscillation. At the very high mean vapor quality the bubble nucleation can be barely seen in the entire periodic cycle since the liquid film covering the heating surface is very thin. In addition, the duct flow is dominated by the annular two-phase flow at all time.     

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.