Effect of surfactant addition on boiling heat transfer in a liquid film flowing in a diverging open channel

An experimental study was conducted to investigate how the addition of small amounts of a surfactant influences the heat transfer characteristics in a thin boiling liquid film flowing in a diverging open channel. Heat transfer experiments were conducted with fluid inlet temperatures from 40 °C to 92 °C. The flow field on the plate included thin film supercritical flow upstream of a hydraulic jump and thick film subcritical flow downstream of a hydraulic jump. Nusselt numbers for the non-boiling heat transfer without surfactant addition scaled linearly with the film Reynolds number. The boiling heat transfer produced higher Nusselt numbers with a weaker dependence on the Reynolds number. Experimental results showed that a boiling surfactant solution created a thick foam layer with high heat transfer rates and Nusselt numbers that are very weakly dependent on the inlet flow rate or the inlet Reynolds number.     

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.     

Contribution of thin-film evaporation during flow boiling inside microchannels

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

Analysis of pool boiling heat transfer: effect of bubbles sliding on the heating surface

Pool boiling on surfaces where sliding bubble mechanism plays an important role has been studied. The heat transfer phenomenon for such cases has been analysed. The model considers different mechanisms such as latent heat transfer due to microlayer evaporation, transient conduction due to thermal boundary layer reformation, natural convection and heat transfer due to the sliding bubbles. Both microlayer evaporation and transient conduction take place during the sliding of bubbles, which occurs in geometries such as inclined surfaces and horizontal tubes. The model has been validated against experimental results from literature for water, refrigerant R134a and propane. The model was found to agree well for these fluids over a wide range of pressures. The model shows the importance of the contributions of the different mechanisms for different fluids, wall superheats and pressures.     

The effect of liquid film evaporation on flow boiling heat transfer in a micro tube

Flow boiling in micro channels is attracting large attention since it leads to large heat transfer area per unit volume. Generated vapor bubbles in micro channels are elongated due to the restriction of channel wall, and thus slug flow becomes one of the main flow regimes. In slug flow, sequential bubbles are confined by the liquid slugs, and thin liquid film is formed between tube wall and bubble. Liquid film evaporation is one of the main heat transfer mechanisms in micro channels and liquid film thickness is a very important parameter which determines heat transfer coefficient. In the present study, liquid film thickness is measured by laser focus displacement meter under flow boiling condition and compared with the correlation proposed for an adiabatic flow. The relationship between liquid film thickness and heat transfer coefficient is also investigated. Initial liquid film thickness under flow boiling condition can be predicted well by the correlation proposed under adiabatic condition. Under flow boiling condition, liquid film surface fluctuates due to high vapor velocity and shows periodic pattern against time. Frequency of periodic pattern increases with heat flux. At low quality, heat transfer coefficients calculated from measured liquid film thickness show good accordance with heat transfer coefficients obtained directly from wall temperature measurements.     

Two-phase cooling characteristics of a saturated free falling circular jet of HFE7100 on a heated disk: Effect of jet length

Direct jet impingement boiling heat transfer operating at low flow rates is of great interest for the localized moderate heat fluxes from the targets with delicate mechanical structure, where the aggressive techniques such as high-speed jets are not suitable. Boiling heat transfer from an upward facing disk targeted by a falling jet was studied experimentally at different volumetric flow rates and various jet lengths. The working fluid was chosen to be the dielectric liquid HFE7100 and the heated spot was an 8-mm diameter disk. Using previous CHF correlations in their original form, valid at very low volumetric flow rates, results in large disagreements since it was found that variation in the jet length changes the boiling characteristics. It is demonstrated that although the circular hydraulic jump formation within the heater diameter may suppress the heat transfer under certain conditions, moving the jet closer to the target may significantly improve the boiling curves at the critical heat flux (CHF) regime. At low flow rates, the CHF increases as the jet length decreases while for moderate and high flow rates the boiling curves show approximately a universal behavior for different jet lengths. For such low flow rates, the effect of jet length on boiling curves was shown to be related to the variation of the cross section of the falling jet and the formation of hydraulic jump at radial distances smaller than the heater diameter. The current CHF results for different jet lengths are correlated by including the effect of jet length in the previous correlation proposed by Sharan and Lienhard.     

Convection and evaporation of microlayer underneath bubbles under microgravity

The multidimensional heat transfer and fluid flow in the microlayer region below a vapor bubble formed during boiling in microgravity are investigated by numerically solving the Navier–Stokes equations with the energy equation. The flow is driven by Marangoni flow due to the surface tension gradient along the bubble surface that results from the temperature gradient. The model also includes condensation and evaporation at the bubble surface. The flow field and heat transfer are calculated for microlayer thicknesses from 0.01 mm to 10 mm to investigate the effect of microlayer thickness. The results show that the velocities are small and have only a small effect on the temperature distribution as compared to the solution for pure conduction in the liquid. Natural convection is shown to have a negligible effect on heat transfer. For less than ideal evaporative heat transfer at the bubble interface, Marangoni convection caused the heat transfer to increase several percent. The flow in the microlayer is shown to agree with the lubrication analogy only for thin, relatively flat interfaces. © 2000 Scripta Technica, Heat Trans Asian Res, 30(1): 1–10, 2001     

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.     

The influence of bubbly flow on boiling from a tube in a bundle

The forms of bubbly flow occurring within a tube bundle are discussed and the boiling process in the bundle is notionally divided into mechanisms due to liquid forced convection, sliding bubbles and nucleation. A novel experimental analysis of heat transfer from a tube in a bundle indicates the predominance of the sliding bubble part. There is a virtual absence of nucleation in a bundle except at the lowest tubes indicating that, once enough bubbles have been produced, the other mechanisms are sufficient to transfer the heat from the tubes.     

Numerical investigation of subcooled flow boiling in segmented finned microchannels

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

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.     

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.     

Nanofluid stabilizes and enhances convective boiling heat transfer in a single microchannel

The flow boiling heat transfer in a single microchannel was investigated with pure water and nanofluid as the working fluids. The microchannel had a size of 7500 × 100 × 250 μm, which was formed by two pyrex glasses and a silicon wafer. A platinum film with a length of 3500 μm and a width of 80 μm was deposited at the bottom channel surface, acting as the heater and temperature sensor. The nanofluid had a low weight concentration of 0.2%, consisting of de-ionized water and 40 nm Al₂O₃ nanoparticles. The nanoparticle deposition phenomenon was not observed. The boiling flow displays chaotic behavior due to the random bubble coalescence and breakup in the milliseconds timescale at moderate heat fluxes for pure water. The flow instability with large oscillation amplitudes and long cycle periods was observed with further increases in heat fluxes. The flow patterns are switched between the elongated bubbles and isolated miniature bubbles in the timescale of 100 s. It is found that nanofluid significantly mitigate the flow instability without nanoparticle deposition effect. The boiling flow is always stable or quasi-stable with significantly reduced pressure drop and enhanced heat transfer. Miniature bubbles are the major flow pattern in the microchannel. Elongated bubbles temporarily appear in the milliseconds timescale but isolated miniature bubbles will occupy the channel shortly. The decreased surface tension force acting on the bubble accounts for the smaller bubble size before the bubble departure. The inhibition of the dry patch development by the structural disjoining pressure, and the enlarged percentage of liquid film evaporation heat transfer region with nanoparticles, may account for the heat transfer enhancement compared to pure water.     

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.     

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.     

Flow Boiling in an Inclined Channel With Downward-Facing Heated Upper Wall

Wall boiling and bubble population balance equations combined with a two-fluid model are employed to predict boiling two-phase flow in an inclined channel with a downward-facing heated upper wall. In order to observe the boiling behavior on the inclined, downward-facing heated wall, a visualization experiment was carried out with a 100 mm × 100 mm of the cross section, 1.2-m-long rectangular channel, inclined by 10° from the horizontal plane. The size of the heated wall was 50 mm by 750 mm and the heat flux was provided by Joule heating using DC electrical current. The temperatures of the heater surface were measured and used in calculating heat transfer coefficients. The wall superheat for 100 kW/m² heat flux and 200 kg/m²s mass flux ranged between 9.3°C and 15.1°C. High-speed video images showed that bubbles were sliding, continuing to grow, and combining with small bubbles growing at their nucleation sites in the downstream. Then large bubbles coalesced together when the bubbles grew too large to have a space between them. Finally, an elongated slug bubble formed and it continued to slide along the heated wall. For these circumstances of wall boiling and two-phase flow in the inclined channel, the existing wall boiling model encompassing bubble growth and sliding was improved by considering the influence of large bubbles near the heated wall and liquid film evaporation under the large slug bubbles. With this improved model, the predicted wall superheat agreed well with the experimental data, while the RPI model largely overpredicted the wall superheat.     

Bubble growth and transient heat transfer at the onset of boiling

An experimental study was conducted to investigate transient local heat transfer around a bubble at onset of boiling on a thin glass heating plate immersed in saturated n-hexane at low pressure. Eight rapid response Cu-Ni thermocouples consisting of a vacuum deposited thin film were used to measure the temperature change of the heating surface. Simultaneous high-speed video photographs were also obtained. The surface temperatures near a nucleation site decreased rapidly owing to the evaporation of a thin layer (microlayer) of liquid formed beneath the bubble in the early period and the rate of bubble growth increased with increasing incipient boiling superheat (ΔT_IB). The thickness of the microlayer decreased markedly with increasing ΔT_IB. © 1998 Scripta Technica, Heat Trans Jpn Res, 26(7): 484–492, 1997     

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.     

Microlayer thickness under a vapor bubble on a cylindrical probe

The thickness of a liquid microlayer underneath a vapor bubble on a heated, cylindrical probe was determined by simultaneously solving the fourth‐order differential equation for the microlayer thickness that incorporates the momentum and energy equations in the microlayer in conjunction with the pressure distribution in the microlayer and the evaporative heat flux at the interface. The analysis also considers the temperature gradient along the probe due to heat transfer in the probe. The results show that the microlayer on a cylindrical surface is very thin and short except for very low probe surface temperatures, superheated less than 1 K. The microlayer size and the evaporative heat flux both decrease rapidly as the surface temperature increases. The results show that most of the evaporation occurs along the curved portion of the interface. © 2000 Scripta Technica, Heat Trans Asian Res, 29(3): 193–203, 2000     

Boiling of binary mixtures—study of individual bubbles

Experiments are reported in which individual bubbles of vapour are grown at a plane wall in initially stagnant isothermal liquid. Such tests had already been done and satisfactorily analysed and understood for pure liquids, but here the liquids are binary mixtures of hexane and octane of varying composition.The chief result is that ciné observations of the general behaviour of such bubbles (rate of growth, change of shape, time and size at departure) presents no new problems; it is identical to behaviour in pure liquids, provided one change of parameter is made. That change is to substitute the correct temperature for the evaporating binary interface in place of the saturation temperature for the pure liquid. That interface temperature was determined by reference to recent analytic solutions for evaporation of semi-infinite binary liquid by diffusion of heat and mass.Detailed rapid thermometric observations were also made, of the temperatures at the wall below the bubble and inside the bubble. They showed some aspects which differ markedly from pure liquids (though not such as to affect general behaviour). Those aspects were explained by detailed analysis of diffusion of heat and mass in a thin layer of liquid (microlayer) left beneath the growing bubble.Implications for boiling heat transfer are also considered and it is shown that this analysis of individual bubbles does not, of itself, explain the known reduction in heat transfer coefficient between pure and binary liquids. Several existing methods of explaining that reduction are examined and shown to have a common basis close to that determined by our analysis of single bubbles. But on extending from that basis in a simple way, to predict boiling heat transfer rates, they fail, except where an additional pressure-dependent multiplier is introduced, on purely empirical grounds.