Bubble sweeping and interactions on wires during subcooled boiling

Multi-bubble dynamics on a heated wire are modeled with a phenomenological model of the temperature distributions around the bubbles and the forces acting on the bubbles. The model calculates the heat transfer from the wire to the bubble and the temperature distribution in the wire. The model balances the Marangoni force, the drag force, and the contact line force acting on the bubble to predict the bubble velocities. The predicted velocities and the predicted interactions agree well with experimentally observed bubble dynamics. The predictions show that when a moving bubble approaches a stationary bubble at moderate superheats, the reduced wire temperatures around the stationary bubble cause the moving bubble to slow and reverse direction before colliding with the stationary bubble. At higher superheats, the bubbles coalesce. The model also shows that when two bubbles approach each other from opposite directions, they will collide and coalesce at lower superheats than when only one bubble is moving because the temperature gradients in front of the moving bubbles are much steeper than in front of a stationary bubble; thus, moving bubbles do not slow much before coalescing.     

Microbubble return phenomena during subcooled boiling on small wires

An experimental investigation was conducted to explore the characteristics of subcooled boiling on microwires of 25 and 100 μm diameter. Microbubbles were observed to return to the wire surface after detachment, with two types of bubble return identified, i.e., isolated bubble return, and bubble return with liquid–vapor trailing jets. The former mode of bubble return occurred when isolated small bubbles (of less than 50 μm diameter) were generated from bubble collapse, while in the latter mode, a larger bubble (of up to 200 μm in diameter) at the end of a liquid–vapor jet issuing from the wire departed and then returned to the wire surface. The numerical simulations conducted show that the isolated bubble return is caused by large temperature gradients in the vicinity of the wire which lead to Marangoni flows and result in a strong thrust force driving the bubble back to the wire. Existence of large temperature gradients close to the microwire surface was demonstrated by experimental measurements, confirming numerical predictions. The numerical model accounts for the influence of noncondensable gas on the vapor saturation temperature as well as the interfacial condensation coefficient. The presence of noncondensable gas facilitates bubble return.     

Mixed Thermocapillary and Forced Convection Heat Transfer Around a Hemispherical Bubble in a Miniature Channel—A 3D Numerical Study

It has been established that for certain conditions, such as microgravity boiling, thermocapillary Marangoni flow has associated with it a significant enhancement of heat transfer. Typically, this phenomenon was investigated for the idealized case of an isolated and stationary bubble resting atop a heated solid that is immersed in a semi-infinite quiescent fluid or within a two-dimensional cavity. This article presents a three-dimensional numerical study that investigates the influence of thermal Marangoni convection on the fluid dynamics and heat transfer around a bubble during laminar flow of water in a minichannel. This mixed thermocapillary and forced convection problem is investigated for channel liquid inlet velocity of 0.01 m/s to 0.03 m/s and Marangoni numbers in the range of 10 to 300 under microgravity conditions. Three-dimensional effects become particularly important on the side and rear regions of the bubble. The thermocapillary forces accelerate the flow along almost the entire bubble interface. The hot core fluid from the heated bottom wall region is forced inward and propelled upward into the thermocapillary jet above the bubble. It can be quantified that the influence of thermocapillary flow on heat transfer enhancement shows an average increase by 40% at the downstream of the bubble and by 60% at the front and rear regions. This heat transfer enhancement depends mainly on the temperature differential as the driving potential for thermocapillary flow and bulk liquid velocity.     

Bubble sweeping and jet flows during nucleate boiling of subcooled liquids

An experimental investigation was conducted to investigate nucleate boiling on a very fine heating wire. By using zoom routine and CCD camera system, the dynamical process of nucleate boiling was visually observed. Sweeping bubbles and several modes of jet flows were described and discussed. For some cases, big bubbles, small bubbles, sweeping bubbles and jet flows coexisted in boiling system, and greatly enhanced heat transfer. These phenomena are quite different from usual observation of nucleate boiling. In this paper, the process of bubble sweeping phenomenon is described in detail and the effect induced by sweeping bubbles is argued. And also, several jet flows are illustrated and discussed, as well as the interaction between bubble sweeping and jet flows.     

Heat transfer for Marangoni‐driven boundary layer flow

Marangoni convection induced by variation of the surface tension with temperature along a surface influences crystal growth melts and other processes with liquid–vapor interfaces, such as boiling in both microgravity and normal gravity in some cases. This paper presents the Nusselt number for Marangoni flow over a flat surface calculated using a similarity solution for both the momentum equations and the energy equation assuming developing boundary layer flow along a surface. Solutions are presented for the surface velocity, the total flow rate, and the Nusselt number for various temperature profiles, Marangoni numbers, and Prandtl numbers. For large bubbles, the predicted boundary layer thickness would be less than the bubble diameter, so the curvature effects could be neglected and this analysis could be used as a first estimate of the effect of Marangoni flow around a vapor bubble. © 2002 Scripta Technica, Heat Trans Asian Res, 31(2): 105–116, 2002; DOI 10.1002/htj.10019     

MARANGONI EFFECT IN MICROBUBBLE SYSTEMS

This work explores the application of the Marangoni effect in micro systems involving small gas or vapor bubbles in a liquid environment subjected to a temperature gradient. The Marangoni effect characterizes the variation of surface tension along the bubble surface resulting from the temperature gradient around the bubble, thus driving the bubble toward the higher temperature region. This phenomenon is more pronounced as the bubble becomes smaller and the temperature gradient becomes steeper, both of which can be achieved in microbubble systems. Potential applications based on the Marangoni effect include linear bubble actuators, dynamic microvalves, and hot-spot locators. The optimum bubble size for these applications is expected to be of the order of 10 mu m. A smaller bubble may be difficult to introduce into the working system and maintain its size. Presented for illustration is a feasibility analysis for both a noncondensable gas bubble and a vapor bubble situated above a microheater. The analysis yields results for the temperature field around the bubble surface and the Marangoni driving pressure, which is a key element of the performance evaluation for all Marangoni-effect-based applications. The findings demonstrate clearly that the Marangoni effect on microbubbles is very significant and shows great promise for applications in microelectromechanical systems (MEMS).     

The Influence of the Magnitude of Gravitational Acceleration on Marangoni Convection About an Isolated Bubble under a Heated Wall

Thermocapillary or Marangoni convection is the liquid motion caused by surface tension variation in the presence of a temperature gradient along a gas–liquid or vapor–liquid interface. This work numerically investigates the effect of the magnitude of gravitational acceleration on the flow and temperature fields resulting from the presence of a hemispherical air bubble of constant radius of 1.0 mm, situated on a heated wall immersed in a liquid silicone oil layer of constant depth of 5.0 mm. The model is oriented such that the Marangoni and gravitational forces act to oppose one another. To elucidate the effect of gravity on Marangoni flow and heat transfer, the simulations were carried out for a silicone oil of Prandtl number 83, at a Marangoni number of 915. The gravity levels tested were 0g, 0.01g, 0.1g, 0.25g, 0.5g, 0.75g, and 1g, where g represents the earth gravitational acceleration of 9.81 m/s ² . The influence of the magnitude of gravitational acceleration on the velocity profile along the bubble interface and on the location of maximum velocity was analyzed. It was found that the gravity level affects the velocity profile by influencing the interfacial temperature gradient, but that the location of maximum velocity was almost independent of gravity level. The increase in heat flux on the wall to which the bubble is attached was calculated and it has been determined that local heat transfer enhancement of up to nearly 1.7 times that of the conduction only case can be achieved for the parameter range tested. Furthermore, local enhancement was observed to occur up to a distance of seven bubble radii for the zero-gravity case, but increased gravity levels cause a reduction in the effective radius of enhancement. The influence of the Marangoni flow on the heat transfer for the opposite cooled wall has also been analyzed.     

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     

Dynamics of discrete bubble in nucleate pool boiling on thin wires in microgravity

A space experiment on bubble behavior and heat transfer in subcooled pool boiling phenomenon has been performed utilizing the temperature-controlled pool boiling (TCPB) device both in normal gravity in the laboratory and in microgravity aboard the 22^nd Chinese recoverable satellite. The fluid is degassed R113 at 0.1 MPa and subcooled by 26°C nominally. A thin platinum wire of 60 μm in diameter and 30 mm in length is simultaneously used as heater and thermometer. Only the dynamics of the vapor bubbles, particularly the lateral motion and the departure of discrete vapor bubbles in nucleate pool boiling are reported and analyzed in the present paper. It’s found that these distinct behaviors can be explained by the Marangoni convection in the liquid surrounding vapor bubbles. The origin of the Marangoni effect is also discussed.     

竖直通道内相邻气泡对上升的直接数值模拟

采用Level Set方法和耦合表面张力模型的Navier-Stokes方程,结合ALE数值算法,直接模拟了竖直通道内两个相邻气泡的上升。重点研究不同空间布置的8mm气泡对后面的尾迹流及其相互作用。数值模拟准确再现气泡对的变形、吸引及排斥行为,气泡上升速度计算结果与经验式吻合。模拟结果表明,两个气泡后面的尾迹流及其相互作用决定了上升气泡对的行为。并排上升的气泡对,由于尾流区被一个射流流动分隔,气泡对没有聚并;然而当垂直上升气泡对中的尾随气泡有超过50%的投影面积进入到前头气泡的尾流区,聚并现象发生。     

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.     

Interface Shape and Marangoni Effect Around aBubble Within the Thermal Boundary Layer

INTRoDUCTIoNDuetohighheattransferperformancecharacter-izedbysmalltemPeraturedifferencesandhighheatfluxes,transportprocesseswithphasechange,espe-ciallyboilingandcondensationprocessesarewidelyemployedinnumerousenergyconversionandtrans-portsystems,heatingand/orcoolingdevices,andaerospaceaPplications.Priortotheutilizationofboil-ingprocessesinspaceapplications,suchasspacecraftthermalcontrol,additionalunderstandingofboilingheattransferbehaviorisneeded.Becauselargedmer-encesekistinthefiuiddensiti…     

Mixed Convective Heat Transfer Due to Forced and Thermocapillary Flow Around Bubbles in a Miniature Channel: A 2D Numerical Study

Marangoni thermocapillary convection and its contribution to heat transfer during boiling has been the subject of some debate in the literature. Currently, for certain conditions, such as microgravity boiling, it has been shown that Marangoni thermocapillary convection has a significant contribution to heat transfer. Typically, this phenomenon is investigated for the idealized case of an isolated and stationary bubble resting on a heated surface, which is immersed in a semi-infinite quiescent fluid or within a two-dimensional cavity. However, little information is available with regard to Marangoni heat transfer in miniature confined channels in the presence of a cross flow. As a result, this article presents a two-dimensional (2D) numerical study that investigates the influence of steady thermal Marangoni convection on the fluid dynamics and heat transfer around a bubble during laminar flow of water in a miniature channel with the view of developing a refined understanding of boiling heat transfer for such a configuration. This mixed convection problem is investigated under microgravity conditions for channel Reynolds numbers in the range of 0 to 500 at liquid inlet velocities between 0.01 m/s and 0.0 5m/s and Marangoni numbers in the range of 0 to 17,114. It is concluded that thermocapillary flow may have a significant impact on heat transfer enhancement. The simulations predict an average increase of 35% in heat flux at the downstream region of the bubble, while an average 60% increase is obtained at the front region of the bubble where mixed convective heat transfer takes place due to forced and thermocapillary flow.     

Interaction of two cavitation bubbles in a tube and its effects on heat transfer

When two cavitation bubbles exist in a confined space, the interaction between the bubbles significantly affects the characteristics of bubble dynamic behaviors. In this paper, a three-dimensional (3D) model is established to study the growth and collapse of two cavitation bubbles in a heated tube and its effects on heat transfer. The liquid and gas phases throughout the calculation domain are solved by a set of Navier-Stokes equations. It is assumed that the gas inside the bubble is compressible vapor, and the surrounding liquid is incompressible water. The mass transfer between two phases is ignored. The calculated bubble profiles were compared to the available experimental data, and a good agreement has been achieved. Then, the relationship among the bubble motion, flow field and pressure distributions was analyzed. On this basis, the effects of bubble interaction on the heat transfer between the wall surface and sounding liquid were discussed. It is found that heat transfer in the centre wall region is enhanced owing to the vortex flow and micro-jet induced by the bubble contraction and collapse. In contrast, the highest surface temperature appears in the surrounding region, which is mainly attributed to the thermal resistance induced by the bubble. The present study is helpful to understand the heat transfer phenomenon with cavitation in the liquid.     

Fundamental consideration of wall heat partition of vertical subcooled boiling flows

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

Bubble behavior during boiling of subcooled binary mixtures on very thin wires

Experimental investigations were conducted to explore the characteristics of subcooled nucleate boiling of binary mixtures composed of water and ethanol on very thin wires, and an emphasis was addressed to observe bubble movement and bubble jet flow. The experiments exhibited many different bubble jet flows, such as continuous, discontinuous, diverged, bifurcate jet flows, and small bubble spraying, corresponding to different vapor bubble sizes and bubble‐moving modes such as bubble initial momentum and bubble‐generating frequency. The bubble movement within the jet flow also was analyzed. © 2007 Wiley Periodicals, Inc. Heat Trans Asian Res, 36(2): 105–111, 2007; Published online in Wiley InterScience ( www.interscience. wiley.com ). DOI 10.1002/htj.20142     

The flow characteristics of an upward gas‐liquid two‐phase flow in a vertical tube with a wire coil: Part 2. Effect on void fraction and liquid film thickness

An experimental study has been carried out to clarify the characteristics of the void fraction and the liquid film thickness of the air‐water two‐phase flow in vertical tubes of 25‐mm inside diameter with wire coils of varying wire diameter and pitch. The flow pattern in the experiment on the average void fraction and the local void fraction distribution in cross section was a bubble flow, and the liquid film thickness was in the region of semiannular and annular flows. It is clarified from these experiments that the average void fraction in tubes with wire coils is lower than that in a smooth tube and decreases with the wire diameter owing to the centrifugal force of the swirl flow which concentrates bubbles at the center of the tube, that the local liquid film thickness becomes more uniform with a decrease in the pitch of the wire coil, and that the liquid film becomes thicker after the passage through the wire coil with an increase in the wire diameter. © 2002 Wiley Periodicals, Inc. Heat Trans Asian Res, 31(8): 652–664, 2002; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10067     

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

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

Premature transition to film boiling with stepwise heat generation 2nd report: Effect of wall material and surface condition

The effects of wall material and surface condition on the behavior of an initial boiling bubble of R113 subjected to transient heating were investigated using a heater with a large heat capacity. The behavior of the initial bubble is closely related to premature transition to film boiling of liquids with high wettability. An initial bubble, which is peculiarly shaped like a “straw hat” and leads to premature transition in saturated liquid nitrogen (as reported in a previous paper), also appears on a heated wall with large heat capacity and grows rapidly to cover the entire wall surface. From the observations using a high-speed video camera, the initial bubble is found to be a coalesced bubble into which small bubbles activated in succession along the heated surface are rolled. The growth rate of the initial bubble along the heated surface is not greatly affected by the thermal conductivity of the wall material but is affected markedly by the surface roughness. © 1997 Scripta Technica, Inc. Heat Trans Jpn Res. 25(1): 51–63, 1996     

Experimental identification of the phenomenon triggering the net vapor generation in upward subcooled flow boiling of water at low pressure

Observation of the bubble behavior was made using a high-speed camera to investigate the mechanisms to cause the net vapor generation in subcooled flow boiling. In the experiments, water was used as the test fluid, the flow direction was vertical upward, and the pressure was kept close to the atmospheric pressure. At high liquid subcooling close to the condition of the onset of nucleate boiling, all the bubbles were lifted off the heated surface immediately after the nucleation to disappear quickly in the subcooled bulk liquid due to condensation. It was found that the void fraction did not increase significantly unless the liquid subcooling became low enough for some bubbles to be reattached to the heated surface after the lift-off. When the reattachment took place, the bubble lifetime was substantially elongated since the bubbles slid up the vertical heated surface for a long distance after the reattachment. The reattachment therefore contributed to an increase in the void fraction. It was concluded that in the experimental conditions tested in this work, the bubble reattachment to the heated surface was a key phenomenon to cause the sharp increase of the void fraction at the point of net vapor generation.