Experimental study of pool temperature effects on nucleate pool boiling

Nucleate pool boiling experiments with constant wall temperature were performed using pure R113 for subcooled, saturated, and superheated pool conditions. A microscale heater array and Wheatstone bridge circuits were used to maintain the constant wall temperature and to measure the instantaneous heat flow rate accurately with high temporal and spatial resolutions. Images of bubble growth were taken at 5000 frames per second using a high-speed CCD camera synchronized with the heat flow rate measurements. The bubble geometry was obtained from the captured bubble images. The effect of the pool conditions on the bubble growth behavior was analyzed using dimensionless parameters for the initial and thermal growth regions. The effect of the pool conditions on the heat flow rate behavior was also examined. The bubble growth behaviors during subcooled, saturated, and superheated pool boiling were analyzed using a modified Jakob number that we newly defined. Dimensionless time and bubble radius parameters with the modified Jakob number characterized the bubble growth behavior well. These phenomena require further analysis for various pool temperature conditions, and this study will provide good experimental data with precise constant wall temperature boundary condition for such works.     

Single bubble growth in saturated pool boiling on a constant wall temperature surface

Nucleate pool boiling experiments with constant wall temperatures were performed using R11 and R113 for saturated pool boiling conditions. A microscale heater array and Wheatstone bridge circuits were used to maintain a constant wall temperature condition and to obtain measurements with high temporal and spatial resolution. Accurate heat flow rate data were obtained from microscale heater array by controlling the surface conditions at a high temporal resolution. Images of the bubble growth were captured using a high-speed CCD camera synchronized with the heat flow rate measurements. The geometry of the bubble was obtained from the images. In the asymptotic growth region, the bubble showed a growth rate that was proportional to t^1/5, which was slower than the growth rate proposed in previous analytical analyses. The bubble growth behavior was analyzed using a new dimensionless parameter to permit comparisons with previous results at the same scale. The comparisons showed good agreement in the asymptotic growth region. A non-dimensional correlation for the bubble radius that can predict the bubble growth and the heat flow rate simultaneously, was suggested. The required heat flow rate for the volume change of the observed bubble was estimated to be larger than the instantaneous heat flow rate measured from the wall. Heat, other than the instantaneous heat supplied from the wall, is estimated to be transferred through the interface between bubble and liquid, even with saturated pool conditions. This phenomenon under a saturated pool condition needs to be analyzed and the data from this study can supply the good experimental data with the precise boundary condition (constant wall temperature).     

Studies on heat transfer and flow characteristics in subcooled flow boiling—Part 2. Flow characteristics

Distributions of fluid temperature and its fluctuation are measured across a R-113 subcooled boiling flow channel with heat fluxes up to the CHF. A microthermocouple probe associated with an electric compensation circuit for the time constant is used for this purpose. Applying statistical treatments to the recorded temperature fluctuation, the heat transfer process in the flow and the characteristics of the bubbles flowing close to the heated surface are investigated. For high heat fluxes nearby the CHF, some bubbles adjacent to the heated surface show a clear trend to coalesce to large volume bubbles with relatively long passing periods, suggesting a mechanism of departure from nucleate boiling by periodical wall temperature rise due to momentary liquid film dryout underneath the large bubbles.     

Experimental investigation of micro-scale temperature transients in sub-cooled flow boiling on a horizontal heater

Surface temperature fluctuations that occur locally underneath departing bubbles in pool boiling are shown to result in local heat transfer coefficients ranging from 1 to 10 kW/cm². These estimates were reported in the literature involved both numerical and experimental approaches. Significantly higher heat fluxes are associated with flow boiling than pool boiling under similar conditions of wall superheat and liquid subcooling (e.g. at boiling inception and at critical heat flux). These enhancements are primarily caused by the convective transport, acceleration/distortion of the bubble departure process as well as the resultant potential enhancement of the local surface temperature fluctuations.In this study we measure the surface temperature fluctuations using temperature micro/nano-sensors fabricated on a silicon wafer during flow boiling on the silicon wafer which is heated from below. The silicon wafer is clamped on a constant heat flux type calorimeter consisting of a vertical copper cylinder with embedded cartridge heaters and K-type thermocouples. Micro/nano-thermocouples (thin film thermocouples or “TFT”) are fabricated on the surface of the silicon wafer. High speed data acquisition apparatus is used to record temperature data from the TFT at 1 kHz. A fluorinert was used as the test fluid (PF-5060, manufacturer: 3M Co.). The calorimeter and surface temperature measurement apparatus is housed in a test section with glass walls for visual observation. The liquid is pumped from a constant temperature bath to maintain a fixed subcooling during the experiments under steady state conditions. The transient temperature data from the FFT array during flow boiling on the silicon wafer is analyzed using fast Fourier transform (FFT). The FFT data is analyzed as a function of the wall heat flux and wall superheat. The number of temperature peaks in the FFT data is observed to increase with increase in wall heat flux and the peaks are found to cover a wider spectrum with peaks at higher frequencies with enhancement of heat flux. The surface temperature fluctuations, especially at small length and time scales, are perturbed potentially by the coupled hydrodynamic and thermal transport processes, resulting in enhanced local and global heat flux values. Boiling incipience condition and the flow boiling data are compared with correlations reported in the literature.     

Experimental study on bubble departure characteristics in subcooled nucleate pool boiling

The boiling models use departure diameter and frequency in closure relations for the calculation of nucleate boiling heat flux. These parameters are normally derived from empirical correlations which depend heavily on experiments. While these parameters are studied mostly for saturated conditions, there is not sufficient data for the values of departure diameter and frequency in subcooled boiling. In this work, the bubble departure characteristics, i.e. the departure diameters and frequency have been measured using high speed visualization experiments with subcooled demineralized water at atmospheric pressure for nucleate pool boiling conditions. The water pool dimensions were 300 mm × 135 mm × 250 mm with four different heating elements to carry out the parametric studies of bubble departure behavior. The considered parameters were heater surface roughness, heater geometry and heater inclination along with the experimental conditions like degree of subcooling (ΔT_sub = 5−20 K), superheat (ΔT_sat = 1−10 K) and the heat flux. The departure diameters and frequencies were directly measured from the images captured. It was intended to generate the subcooled nucleate pool boiling data under a wide range of conditions which are not present in the literature. The departure diameter was found to increase with the wall superheat, heater size and the inclination angle while the liquid subcooling and surface roughness produced a damping effect on the diameter. The departure frequency was found to increase with the wall superheat and the inclination angle, but decreases with an increase in the heater size. The frequency increases with the degree of subcooling except very close to the saturation, and is unaffected by the surface roughness beyond a certain superheat value.     

Heat transfer at the phase interface of bubbles collapsing in subcooled liquids and during subcooled boiling

Heat transfer phenomena with bubble condensation in subcooled liquids and during subcooled boiling are discussed, which were studied by using holographic interferometry. This inertialess and non-intrusive measuring technique allows interesting insights into fluid dynamics and thermodynamics of phase change processes. The resistance for the heat transport through the phase boundaries is mainly on the liquid side and therefore the heat transfer coefficient can expressed by simple correlations following the procedure in forced convective flow. For the sake of simple boundary conditions first bubble condensation in subcooled liquids with constant and homogenous temperature are studied and than the more complex situation with large temperature gradients in the wall near boundary during subcooled boiling are discussed.Dedicated to Prof. Dieter Mewes on the occasion of his 65th birthday.     

Enhanced heat transfer in confined pool boiling

We report the results of an experimental investigation of the heat transfer during nucleate boiling on a spatially confined boiling surface. The heat flux as a function of the boiling surface temperature was measured in pool boiling pots with diameters ranging from 15 mm down to 4.5 mm. It was found that a reduction of the pool diameter leads to an enhancement of the nucleate boiling heat flux for most of the boiling curve. Our experimental results indicate that this enhancement is not affected by the depth of the boiling pot, the material of the bounding wall, or the diameter of the inlet water supply. High-speed camera imaging shows that the heat transfer enhancement for the spatially confined pool boiling occurs in conjunction with a stable circulating flow, which is in contrast to the chaotic and mainly upward motion for boiling in larger pool diameters. An explanation for the enhancement of the heat transfer and the associated change in flow pattern is found in the singularisation of the nucleate boiling process.     

Studies on heat transfer and flow characteristics in subcooled flow boiling—Part 1. Boiling characteristics

Measurement of wall temperature profile and photographic observation are performed for R-113 subcooled boiling flow in a channel with heat fluxes up to the CHF. The incipient boiling superheats measured are little affected by mass velocity and liquid subcooling. Hysteresis in boiling observed by increasing and decreasing heat flux seems to be ascribed to variation in size of active nucleation cavities on the wall. Increasing heat flux up to the CHF, the bubble density on the heated surface increases and remarkably large coalescent bubbles appear periodically near the heating section outlet.     

Experimental and numerical study on nucleate bubble deformation in subcooled flow boiling

Experimental and numerical study was conducted to investigate the bubble behaviors in subcooled flow nucleate boiling. The bubble behaviors in subcooled flow boiling in an upward annular channel were investigated in the range of subcooling degree 5–30 K by visualization with high spatial and temporal resolutions using a high speed video camera and Cassegrain tele-microscope. Obvious deformation on the upstream side surface of the bubble during its growth process was frequently observed. This deformation phenomenon was caused by the condensation occurring at the upstream side bottom of the bubble, which results from the Marangoni flow along the bubble surface from the bubble bottom to the top. Since the Marangoni flow cannot be directly observed by the current experiments because it occurs in a very thin interface along the bubble surface, the numerical simulations of bubble growth and departure behaviors in subcooled flow boiling were carried out. As a result, it was confirmed that the bubble deformation was caused by the Marangoni flow along the bubble surface. Moreover, the phenomenon of wave propagation on the bubble surface during the condensation process was observed, and it can enhance the heat transfer between the bubble and the surrounding subcooled liquid.     

Electric field effects during nucleate boiling from an artificial nucleation site

An experimental study of saturated pool boiling from a single artificial nucleation site on a polished copper surface has been performed. Isolated bubbles grow and depart from the artificial cavity and the bubble dynamics are recorded with a high speed camera. Experimental results are obtained for bubble growth, departure and vertical rise both with and without the application of an electric field between an upper electrode and the boiling surface. As detailed in a previous paper from the same research group the high spatial and temporal resolution of the video sequences facilitated the development of a baseline experimental bubble growth law which predicts the bubble volumetric growth characteristics for a range of surface superheats at atmospheric pressure. The presence of an electric field has been found to positively augment the convective heat transfer over that of buoyant natural convection. Further to this, for high electric field strengths, the bubble shape, volumetric growth characteristics and bubble rise are different from that of the baseline cases. These results provide compelling evidence that electric fields can be implemented to alter the bubble dynamics and subsequent heat transfer rates during boiling of dielectric liquids.     

A study of nucleate bubble growth on microstructured surface through high speed and infrared visualization

We studied bubble growth on a microstructured surface during nucleate boiling using optical high-speed and infrared (IR) cameras. The effects of structured surfaces on bubble growth and dynamics were examined and their role analyzed with the use of simple models. A smooth, bare surface was prepared, and four microstructured test sample surfaces were fabricated with microscale gaps ranging from 5 to 80 µm. The optical high-speed camera was used to observe the bubble growth profile with high temporal resolution; the IR camera was focused on the underside of the sample for direct visualization of the boiling process. Overall, the microstructured surfaces produced more bubbles, a lower frequency and nucleation site density than the bare surface for the low heat flux range (100–300 kW/m²), corresponding to the isolated bubble growth regime. The liberated bubble size was dependent on the size of the microstructure gap. Analysis of the high-speed images revealed that the liquid between the microstructures did not evaporate during bubble growth; however, during the initial growth stage, there was a brief period in which the liquid at the nucleation site evaporated. The large surface area and relatively high number of nucleation points contributed to enhanced bubble growth on the structured surfaces.     

Analytical calculation of DNB-superheating by a postulated thermo-mechanical effect of nucleate boiling

Observations of stationary subcooled nucleate boiling with forced convection in a glass-pipe at a pressure of 0.1 MPa lead to a new interpretation of the wall. This interpretation applied to the rules of linear non-equilibrium thermodynamics results in framework equations with non-linear heat transfer as indicated by several empirical correlations. Momentum balances at the singular surfaces of a bubble have been combined with a sonic limit for the mass transfer through the interfaces serving as a maximum condition with respect to heat transfer (DNB). The calculation predicts the corresponding wall temperature from properties of state without using empirical coefficients and therefore must be valid independent of coolants' geometries and surface conditions. Documented measurements at pool boiling, forced and free convection with cryogenic liquids, water and liquid metals emphasize this if data are properly selected corresponding to the precondition that void fraction remains low.     

Convective boiling heat transfer in a concentric annular gap

Heat transfer and pressure drop characteristics were measured in both subcooled and saturated concentric narrow gap forced convection boiling. Four boiling regimes were observed: isolated bubbles, coalesced bubbles, coalesced bubbles with small dry patches, and liquid deficient flow. Data were obtained to qualitatively identify the effects of gap size, pressure, flow rate, inlet quality and wall superheat on the boiling regimes and the transitions between the various regimes. Some significant differences from unconfined forced convection boiling were found.     

Study of transient heat transfer and synchronized flow visualizations during sub-cooled flow boiling in a small aspect ratio microchannel

The challenges that microchannel flow boiling technology faces are the lack of understanding of underlying mechanisms of heat transfer during various flow boiling regimes and a dearth of analytical models that can predict heat transfer. This paper aims to understand flow boiling heat transfer mechanisms by analyzing results obtained by synchronously captured high-speed flow visualizations with local, transient temperature data. Using Inverse Heat Conduction Problem (IHCP) solution methodology, the transient wetted surface heat flux and temperature as well as heat transfer coefficient are calculated. These are then correlated with the visual data. Experiments are performed on a single microchannel embedded with fast response temperature sensors located (630 µm) below the wetted surface. The height, width and length of the microchannel are 0.42 mm, 2.54 mm and 25.4 mm respectively. De-ionized, de-gassed water is used as the working fluid. Two heat fluxes are tested at each of the mass fluxes of 182 kg/(m²s) and 380 kg/(m²s). Because of vapor confinement, slug flow is observed for the tested conditions. The present study provides detailed insights into the effect of various events such as passage of vapor slug, 3-phase contact line, partial-dry-out and liquid slug on transient heat transfer coefficient. Transient heat transfer coefficient peaks when thin film evaporation mechanism is prevalent. The peak value is influenced by the distance of bubble incipience as well as downstream events obstructing the flow. Heat transfer coefficient during the passage of liquid slug and 3-phase contact line were relatively lower for the tested experimental conditions.     

Boiling regimes in uncoated polydimethylsiloxane microchannels with a fine wire heater

Polydimethylsiloxane (PDMS) is a type of gas permeable media widely used in microfluidic applications. In this work multiphase patterns and boiling curves in PDMS square microchannels were experimentally investigated. Very fine platinum wires with diameter of 50 μm were embedded through the microchannels and serving as heater. The multiphase patterns were visualized by means of high speed CCD camera with microscope. Curves of temperature versus heat flux on the wire heaters were plotted. Based on the evolution of multiphase patterns, five boiling regimes were classified, that is, single phase, bubble formation, slug formation, slug dominated and dry out. Interestingly, the bubbles were generated from the channel walls rather than the heater surface, and so-called “droplets-in-bubble” phenomenon drew attention in which bunches of microdroplets kept forming, growing, and disappearing within the big bubbles. The boiling curves were plotted and compared to boiling in open space and in glass tubes. The heat transfer in the PDMS microchannels got deteriorated when the bubbles formed.     

Single bubble growth in saturated pool boiling of binary mixtures

Nucleate pool boiling experiments for binary mixtures, which are consisted of R11 and R113, were performed with constant wall temperature condition. Results for binary mixtures were also compared with pure fluids. A microscale heater array and Wheatstone bridge circuits were used to maintain the constant temperature of the heating surface and to obtain heat flow rate measurements with high temporal and spatial resolutions. Bubble growth images were captured using a high-speed CCD camera synchronized with the heat flow rate measurements.The departure time for binary mixtures was longer than that for pure fluids, and binary mixtures had a higher onset of nucleate boiling (ONB) temperature than pure fluids. In the asymptotic growth region, the bubble growth rate was proportional to a value between t^1/6 and t^1/4. The bubble growth behavior was analyzed to permit comparisons with binary mixtures and pure fluids at the same scale using dimensionless parameters. There was no discernible difference in the bubble growth behavior between binary mixtures and pure fluids for a given ONB temperature. And the departure radius and time were well predicted within a ±30% error.The minimum heat transfer coefficient of binary mixtures occurred near the maximum |y−x| value, and the average required heat flux during bubble growth did not depend on the mass fraction of R11 as more volatile component in binary mixtures. Finally, the results showed that for binary mixtures, a higher ONB temperature had the greatest effect on reducing the heat transfer coefficient.     

Review of nucleate pool boiling bubble heat transfer mechanisms

Enhanced convection, transient conduction, microlayer evaporation, and contact line heat transfer have all been proposed as mechanisms by which bubbles transfer energy during boiling. Models based on these mechanisms contain fitting parameters that are used to fit them to the data, resulting a proliferation of “validated” models. A review of the recent experimental, analytical, and numerical work into single bubble heat transfer is presented to determine the contribution of each of the above mechanisms to the overall heat transfer. Transient conduction and microconvection are found to be the dominant heat transfer mechanisms. Heat transfer through the microlayer and at the three-phase contact line do not contribute more than about 25% of the overall heat transfer.     

Bubble nucleation characteristics in pool boiling of a wetting liquid on smooth and rough surfaces

Quantitative measurements are obtained from high-speed visualizations of pool boiling at atmospheric pressure from smooth and roughened surfaces, using a perfluorinated hydrocarbon (FC-77) as the working fluid. The boiling surfaces are fabricated from aluminum and prepared by mechanical polishing in the case of the smooth surface, and by electrical discharge machining (EDM) in the case of the roughened surface. The roughness values (R_a) are 0.03 and 5.89 μm for the polished and roughened surfaces, respectively. The bubble diameter at departure, bubble departure frequency, active nucleation site density, and bubble terminal velocity are measured from the monochrome movies, which have been recorded at 8000 frames per second with a digital CCD camera and magnifying lens. Results are compared to predictions from existing models of bubble nucleation behavior in the literature. Wall superheat, heat flux, and heat transfer coefficient are also reported.     

A method of concurrent thermographic–photographic visualization of flow boiling in a minichannel

A method is developed to capture the distribution of surface temperature while simultaneously imaging the bubble motions in diabatic flow boiling in a horizontal minichannel. Liquid crystal thermography is used to obtain highly resolved surface temperature measurements on the uniformly heated upper surface of the channel. High-speed images of the flow field are acquired simultaneously and are overlaid with the thermal images. The local surface temperature and heat transfer coefficient can be analyzed with the knowledge of the nucleation site density and location, and bubble motion and size evolution. The horizontal channel is 1.2 mm high × 23 mm wide × 357 mm long, and the working fluids are Novec 649 and R-11. Optical access is through a machined glass plate which forms the bottom of the channel. The top surface is an electrically heated 76 μm-thick Hastelloy foil held in place by a water-cooled aluminum and glass frame. The heat loss resulting from this construction is computed using a conduction model in Fluent. The model is driven by temperature measurements on the foil, glass plate and aluminum frame. This model produces a corrected value for the local surface heat flux and enables the computation of the bulk fluid temperature and heat transfer coefficient along the channel. The streamwise evolution of the heat transfer coefficient for single-phase laminar flow is compared to theoretical values for a uniform-flux boundary condition. Examples of the use of the facility for visualizing subcooled two-phase flows are presented. These examples include measurements of the surface temperature distribution around active nucleation sites and the construction of boiling curves for locations along the test surface. Points on the curve can be associated with specific image sequences so that the role of mechanisms such as nucleation and the sliding of confined bubbles may be discerned.