Saturated 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 saturated 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. The measured heat transfer data indicate that the saturated flow boiling heat transfer coefficient increases with a decrease in the gap of the duct. Besides, raising the imposed heat flux 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 milder. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases slightly at increasing R-134a mass flux. Moreover, the bubble departure frequency increases at reducing duct size mainly due to the rising shear stress of the liquid flow, and at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Correlation for the present saturated flow boiling heat transfer data of R-134a in the narrow annular duct is proposed. Additionally, data for some quantitative bubble characteristics such as the mean bubble departure diameter and frequency and the active nucleation site density are also correlated.     

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

Experiments are conducted here to investigate how the channel size affects the R-407C saturated flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. The measured data indicate that the saturated flow boiling heat transfer coefficient increases with a decrease in the duct gap. Besides, raising the imposed heat flux 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 milder. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases noticeably at increasing R-407C mass flux. Moreover, the bubble departure frequency increases at reducing duct size and at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Meanwhile comparison of the present heat transfer data for R-407C with R-134a in the same duct and with some existing correlations is conducted. Furthermore, correlation for the present R-407C saturated flow boiling heat transfer data 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.     

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.     

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.     

Modeling of Particulate Fouling on Heat Exchanger Surfaces: Influence of Bubbles on Iron Oxide Deposition

Studies of iron oxide deposition on Alloy-800 heat exchanger tubes have been part of a continuing research program at the University of New Brunswick (UNB); the present work formulates mechanisms for the effect of bubbles on deposition in water under boiling conditions. To supplement results from earlier deposition experiments in a fouling loop at UNB, measurements of bubble frequency and departure diameter as a function of heat flux were performed. High-speed movies of bubbling air/water systems indicated that a pumping action moved particles from adjacent areas at the surface to bubble nucleation sites. To explain the observations, the model considers deposition and concomitant removal. Deposition includes microlayer evaporation and filtration through the porous deposit. The deposit is sparse in the first stage, when the dominant process is microlayer evaporation including particle trapping and pumping, creating spots of deposit. Filtration becomes more important as the deposit thickens to a stage when microlayer evaporation becomes negligible. Chimney effects then control. Turbulence due to detaching and collapsing bubbles affects removal. In subcooled boiling, collapsing bubbles generate enough turbulence to maintain much of the deposit labile, while in bulk boiling bubble detachment from the nucleation site is dominant and a smaller portion of the deposit is labile and subject to removal. Model predictions are presented and shown to agree quite well with experimental data.     

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.     

Unified theoretical prediction of fully developed nucleate boiling and critical heat flux based on a dynamic microlayer model

A new dynamic microlayer model has been proposed to predict theoretically the heat flux in fully developed nucleate boiling regions including critical heat flux (CHF). In this model, the heat transfer with boiling is mainly attributed to the evaporation of the microlayers which are periodically formed while the individual bubbles are forming. Since the initial microlayer thickness becomes thinner with the increase of wall superheat, both the local evaporation and the partial dryout speed of the microlayer increase. As a result, the time-averaged heat flux during the period of individual bubble has a maximum point, the CHF, at the predicted continuous boiling curve.     

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.     

The effect of fouling on nucleate pool boiling of small wires

Pool boiling experiments with small diameter wires were conducted in earth gravity and microgravity conditions. Bare wire and fouled wire with a scale deposition of calcium carbonate was used. The wettability on the scale wire was higher than that on the bare wire. Though more vigorous bubbling was observed on the scale wire when compared to that on the bare wire at the same heat flux, the boiling curve for the scale wire was approximately the same as that on the bare wire. However, the critical heat flux (CHF) on the scale wire was higher than that of the bare wire. On the scale wire, the departure diameter of bubbles was relatively smaller than that on the bare wire. The smaller diameter of bubbles detaching from the scale wire is considered to be due to the high wettability and high nucleation site density. As the result, the coalescence of bubbles near the wire was prevented, and the CHF was delayed and increased on the scale wire when compared to that on the bare wire. © 2004 Wiley Periodicals, Inc. Heat Trans Asian Res, 33(5): 316–329, 2004; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.20016     

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.     

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.     

Numerical prediction of critical heat flux in pool boiling with the two-fluid model

Three-dimensional numerical simulations of the atmospheric saturated pool boiling are performed with the aim of predicting the critical heat flux. The two-phase mixture in pool boiling is described with the transient two-fluid model. The transient heat conduction in the horizontal heated wall is also solved. Dynamics of vapor generation on the heated wall is modeled through the density of nucleation sites and the bubble residence time on the wall. The heater’s surface is divided into zones, which number per unit area equals the density of nucleation sites, while the location of nucleation site within each zone is determined by a random function. The results show a replenishment of the heater’s surface with water and surface wetting for lower heat fluxes, while heater’s surface dry-out is predicted at critical heat flux values. Also, it is shown that the decrease of nucleation site density leads to the reduction of critical heat flux values. Obtained results of critical heat flux are in good agreement with available measured data. The presented approach is original regarding both the application of the two-fluid two-phase model for the prediction of boiling crisis in pool boiling and the defined boundary conditions at the heated wall surface.     

Time periodic flow boiling heat transfer of R-134a and associated bubble characteristics in a narrow annular duct due to flow rate oscillation

An experiment is conducted here to investigate the effects of the imposed time periodic refrigerant flow rate oscillation in the form of nearly a triangular wave on refrigeriant R-134a flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct with the duct gap fixed at 2.0 mm. The results indicate that when the imposed heat flux is close to that for the onset of stable flow boiling, intermittent flow boiling appears in which nucleate boiling on the heated surface does not exist in an entire periodic cycle. At somewhat higher heat flux persistent boiling prevails. Besides, the refrigerant flow rate oscillation only slightly affects the time-average boiling curves and heat transfer coefficients. Moreover, the heated wall temperature, bubble departure diameter and frequency, and active nucleation site density are found to oscillate periodically in time as well and at the same frequency as the imposed mass flux oscillation. Furthermore, in the persistent boiling the resulting heated wall temperature oscillation is stronger for a longer period and a larger amplitude of the mass flux oscillation. And for a larger amplitude of the mass flux oscillation, stronger temporal oscillations in the bubble characteristics are noted. The effects of the mass flux oscillation on the size of the departing bubble and active nucleation site density dominate over the bubble departure frequency, causing the heated wall temperature to decrease and heat transfer coefficient to increase at reducing mass flux in the flow boiling, opposing to that in the single-phase flow. But they are only mildly affected by the period of the mass flux oscillation. However, a short time lag in the wall temperature oscillation is also noted. Finally, a flow regime map is provided to delineate the boundaries separating different boiling regimes for the R-134a flow boiling in the annular duct.     

Saturated flow boiling heat transfer of R-410A and associated bubble characteristics in a narrow annular duct

Experiments are conducted here to investigate how the channel size affects the R-410A saturated flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm in this study. The measured data indicate that the saturated flow boiling heat transfer coefficient increases with increasing refrigerant mass flux and saturated temperature and with a decrease in the gap size. Besides, raising the imposed heat flux can cause a significant increase in the boiling heat transfer coefficient. The results from the flow visualization show that the mean diameter of the bubbles departing from the heating surface decreases slightly with increasing refrigerant mass flux and saturated temperature. Moreover, the bubble departure frequency increases at reducing duct size and increasing mass flux. And at a high imposed heat flux many bubbles generated from the cavities in the heating surface tend to merge together to form big bubbles. Meanwhile, comparisons of the present heat transfer data for R-410A with R-407C and R-134a in the same duct and with some existing correlations are conducted. Furthermore, an empirical correlation for the present R-410A saturated flow boiling heat transfer data is proposed.     

A Visualization Study of Bubble Behavior in Saturated Flow Boiling through a Vertical Mini-Tube

Forced convection saturated R-134a boiling experiments were conducted in a vertical mini-quartz tube coated with a transparent heater; the inner diameter of the tube was 1.33 mm and the heated length 235.5 mm. The dynamics of the saturated boiling process, bubble characteristics, and behavior were studied using a high-speed CCD camera at different mass fluxes in up-flow at 6.425 bar. The heat fluxes were 5 and 20 kW/m². The flow visualization results show that the bubble departure frequency generally increases as the heat flux increases. In addition, bubble lift-off diameter, growth rate, and velocity after bubble lift-off were determined by analyzing the images.     

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.     

Coalescence of bubbles in nucleate boiling on microheaters

In this paper, an experiment was performed which is based on a heating surface consisting of microheaters where the temperature of each heater can be individually controlled by an electronic feedback loop. The power consumed by the heaters throughout the cycle of individual bubble growth, coalescence, detachment and departure was measured at high frequencies, thus the heat flux and its variation were obtained. By a careful timing and control of two individual microheaters, we were able to produce two individual bubbles side-by-side. The coalescence would takes place when they grow to a certain size that allows them to touch each other. We have recorded two major heat flux spikes for a typical cycle of boiling with coalescence. The first one corresponds to the nucleation of bubbles; the second one is for the coalescence of the two bubbles. We found that the heat flux variation is closely related to the bubble dynamics and bubble-bubble interaction. By comparing with the single bubble results without coalescence, we also found that the heat transfer is highly enhanced due to the coalescence.     

Semi-analytical model for boiling from enhanced structures

Understanding of the fundamental mechanism of boiling from enhanced structures is currently incomplete. The main focus of the present study was to develop a semi-analytical model to predict the bubble departure diameter, frequency, and nucleation site density for a boiling enhancement structure. Existing models in literature were used as a framework to develop this. Salient feature of the model is the improvements on sub-models for bubble departure diameter, evaporation within the channels and convective heat transfer from the external surfaces of the enhanced structure. The model was used to calculate the total heat dissipated from the structures. Comparisons to experiments revealed that the bubble departure is predicted within ±10%, frequency within ±30% for all data points except a couple, nucleation site density within ±40% and the heat flux within ±50%. Sample calculations have been included to show its use in optimizing the geometrical parameters for maximizing heat transfer.     

Study of bubble growth in water pool boiling through synchronized,infrared thermometry and high-speed video

High-speed video and infrared thermometry were used to obtain time- and space-resolved information on bubble nucleation and heat transfer in pool boiling of water. The bubble departure diameter and frequency, growth and wait times, and nucleation site density were directly measured for a thin, electrically-heated, indium–tin-oxide surface, laid on a sapphire substrate. These data are very valuable for validation of two-phase flow and heat transfer models, including computational fluid dynamics with interface tracking methods. Here, detailed experimental bubble-growth data from individual nucleation sites were used to evaluate simple, commonly-used, but poorly-validated, bubble-growth and nucleate-boiling heat-transfer models. The agreement between the data and the models was found to be reasonably good. Also, the heat flux partitioning model, to which our data on nucleation site density, bubble departure diameter and frequency were directly fed, suggests that transient conduction following bubble departure is the dominant contribution to nucleate-boiling heat transfer.     

Bubble growth,departure and the following flow pattern evolution during flow boiling in a mini-tube

In the present study, the bubble growth, departure and the following flow pattern evolution during flow boiling in the mini-tube were visualized and quantitatively investigated, along with the simultaneous measurement of the local heat transfer coefficient around a specified nucleation site. Liquid nitrogen was employed as the working fluid and the test section was a segment of vertically upward quartz glass tube with the inner diameter range of 1.3–1.5 mm, which was coated by a layer of transparent ITO film as the heater on the outer surface. The growth rates of bubbles had similar and constant growth rate in two periods of time, i.e., before and after the bubbles departing from the nucleation site, which indicated the bubble growth was primarily governed by the inertial force. The bubble departure diameter and bubble period were investigated and the corresponding correlation was obtained based on the experimental data, which showed that the tube size of the mini-tube had no notable effect on the bubble departure and the trend of the bubble departure was similar to that in macro-tubes. Whereas the following flow pattern evolution was apparently confined due to the size effect, which presented desirable heat transfer performance in mini-tubes. The heat transfer coefficients for different flow patterns along the mini-tube were obtained in terms of bubbly, slug, annular flow and the flow regimes of flow reversal and post dryout. It was found that the dominant heat transfer mechanism was the liquid film evaporation which offered desirable heat transfer capability. The heat transfer performance would be deteriorated in the post dryout regime, while flow reversal could somewhat enhance the heat transfer upstream of the nucleation site. Boiling curves around the specified nucleation site were recorded and analyzed based on the recorded flow patterns.