Robust numerical analysis of the dynamic bubble formation process in a viscous liquid

The dynamics of bubble formation from a submerged nozzle in a highly viscous liquid with relatively fast inflow gas velocity is studied numerically. The numerical simulations are carried out using a sharp interface coupled level set/volume-of-fluid (CLSVOF) method and the governing equations are solved through a hydrodynamic scheme with formal second-order accuracy. Numerical results agree well with experimental results and it is shown that the sharp interface CLSVOF method enables one to reproduce the bubble formation process for a wide range of inflow gas velocities. From numerical results, one can improve their understanding of the mechanisms regarding the dynamics of bubble formation. For example, it is found that for some sets of parameters that the bubble formation process reaches steady state after several bubbles are released from the nozzle. At steady state, bubbles uniformly rise freely in the viscous liquid. It is observed that the fluid flow around a formed bubble has a significant role in determining the overall dynamic process of bubble formation; e.g. the effect of the fluid flow from the preceding bubble can be seen on newly formed bubbles.     

Interactions between a central bubble and a surrounding bubble cluster

The interaction of multiple bubbles is a complex physical problem. A simplified case of multiple bubbles is studied theoretically with a bubble located at the center of a circular bubble cluster. All bubbles in the cluster are equally spaced and own the same initial conditions as the central bubble. The unified theory for bubble dynamics [35] is applied to model the interaction between the central bubble and the circular bubble cluster. To account for the effect of the propagation time of pressure waves, the emission source of the wave is obtained by interpolating the physical information on the time axis. An underwater explosion experiment with two bubbles of different scales is used to validate the theoretical model. The effect of the bubble cluster with a variation in scale on the pulsation characteristics of the central bubble is studied.     

Observations on braided thin wire nucleate boiling in microgravity

A microgravity experiment was conducted on the Space Shuttle Endeavor (STS-108) to observe sustained nucleate boiling of water. Subcooled water was boiled with a single strand and a braid of three 0.16 mm diameter and 80 mm long Nichrome resistive wires. A CCD video camera recorded the experiment while six thermistors recorded the temperature of the fluid at various distances from the heating element. This paper reports experimental results in observations, measurements, and data analysis. Bubble explosions were found to take place shortly after the onset of boiling for both the single and braid of wires. The explosion may produce a high heat transfer rate, as it generates a cloud of microbubbles. The number, size, and departure rate of the bubbles from the heater wire were measured and compared with theoretical models as a function of time. The temperature measurements revealed a complex temperature distribution in the fluid chamber due to bubbles ejected from the wire that carried thermal energy close to the temperature sensors. Drag forces on departing bubbles were calculated based on bubble movement and used to predict bubble propagation. Results from this experiment provided further understanding of nucleate boiling dynamics in microgravity for the eventual design and implementation of two-phase heat transfer systems in space applications.     

A boundary element method for the simulation of non-spherical bubbles and their interactions near a free surface

The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4 π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demonstrated in our calculation results.     

Bubble characteristics in a bubbling fluidized bed with a rotating distributor

In this paper, the effect of a novel rotating distributor for fluidized beds on the bubble size is studied. The distributor is a perforated plate that rotates around the vertical axis of the column.The formation of the bubbles on the rotating distributor is theoretically analyzed. The pierced length of the bubbles ascending in the bed were measured using optical probes. The probability distribution of bubble diameter was inferred from these experimental measurements using the maximum entropy method. The radial profile of the bubble diameter is presented for the static and rotating configurations at different gas velocities. The frequency of bubble passage and the distribution of bubbles in the cross section of the bed are also reported. Results were finally shown for different heights above the distributor.A radial decrease in the bubble size when the distributor rotates is found. The bubble growth with the bed height is also lower in the rotating case.     

Diffusion-induced growth of a gas bubble in a viscoelastic fluid

The diffusion-induced growth of a spherical gas bubble surrounded by a thin shell of viscoelastic fluid containing a limited amount of dissolved gas is analyzed. This is representative of a situation when a large number of bubbles grows in close proximity in a viscoelastic medium. The upper-convected Maxwell model is employed to describe the rheology of the fluid. Limited quantities of the dissolved gas available in the liquid shell mandates solution of the convection-diffusion equation, as opposed to using similarity solutions or polynomial profiles to describe the mass transport across the interface. Utilizing the properties of a potential field and a Lagrangian transformation, a new approach is introduced to solve the coupled system of integro-differential equations governing the bubble growth. The results indicate that, at the early stages of the growth, bubbles in a viscoelastic fluid grow faster than in a Newtonian fluid. However, eventually they attain the same steady-state configuration.     

On the collapse of cavities

The collapse of a single cavity, or a cloud of bubbles has several physical consequences when in proximity to a structure or resident within a material during deformation. The earliest recognized of these was cavitation erosion of the propellers of steam ships. However, other processes include the rapid collapse of cavities leading to hot spots in explosives from which reaction ensues, or the more recent phenomenon of light generation by oscillating single bubbles or clouds. In the collapse of a cavity, the least considered but the most important mechanism is asymmetric closure. One of the consequences of this is the formation of jets leading to local high pressures and shears that result in the damage or reaction mechanisms observed. The challenge for the future remains in understanding the effects of cloud cavitation since it is likely that only one bubble in perhaps millions in a cloud catalyses an event. The review follows the author's work in the understanding of shock-induced cavity collapse and highlights several results which indicate the importance of this problem in a variety of fields. Received 27 July 2001/ Accepted 25 January 2002     

Sound generation on bubble coalescence following detachment

A system in which bubbles coalesced on formation was used to probe one mechanism by which bubbles create sound. The aim was to determine in which situations sound is produced and to predict its amplitude. A set of carefully co-ordinated high-speed video and acoustic timeseries showed that needle-formed bubbles generated loud bubble-acoustic emissions at the instant of coalescence of secondary bubbles with the primary bubble. As the air flow rate increased, the size and number of secondary bubbles increased, and the sound amplitude also increased. On coalescence, the sound pressure always rose initially. A dimensionless scaling found that the sound amplitude emitted scaled with the volume of the secondary bubble. This scaling was shown to be consistent with the sound-emission mechanism being the equalization of pressures in the coalescing bubbles. The trend in amplitude with bubble production rate was well predicted by the scaling.     

Numerical simulation of bubble formation on orifice plates with a moving contact line

The formation of bubbles on an orifice plate involves a moving contact line, especially in case of poor wetting conditions. The dynamics of the moving contact line and contact angle have a significant impact on the bubble departure size. Therefore, for the numerical simulation, an appropriate contact line boundary condition is essential for a correct prediction of the bubble formation. Numerical tests have been performed on two kinds of contact line models, one is a contact line velocity dependent model (Model-A, a commonly used model) and the other is a stick-slip model (Model-B). The calculation results using Model-A depend greatly on the prescribed maximum contact line velocity. With Model-B a parameter-independent prediction can be obtained provided that the mesh is sufficiently fine. The dynamic advancing and receding contact angles, which are two required inputs to both models, have a significant influence on the predicted bubble departure diameter, if the contact line moves beyond the inner rim of the orifice. The effect of wettability on the bubble departure size is realized via the variation of the maximum contact diameter. When the contact line sticks to a small region near the inner rim of the orifice, such as the bubble formation on a thin-walled nozzle, the effects of the contact angle and contact line models are negligible.     

Bubble–wall interaction and two-phase flow parameters on a full-scale boat boundary layer

Full scale bubbly flow experiments were performed on a 6 m flat bottom survey boat, measuring the void fraction, bubble velocity and size distributions as the bubbles naturally entrained at the bow of the boat interact with the boat’s boundary layer. Double-tip sapphire optical probes capable of measuring bubbles down to 50 μm in diameter were specifically designed and built for this experiment. The probes were positioned under the hull at the bow near the bubble entrainment region and at the stern at the exit of the bottom flat plate. Motorized positioners were used to vary the probe distance to the wall from 0 to 50 mm. The experiments were performed in fresh water (Coralville Lake, IA) and salt water (Panama City Beach, FL), at varying velocities with most data analysis performed at 10, 14 and 18 knots. The results indicate that the bubbles interact significantly with the boundary layer. At low velocity in fresh water, bubble accumulation under the hull and coalescence are evident by the presence of large bubbles at the stern. At high speeds bubble breakup dominates and very small bubbles are produced near the wall. It is also observed that salt water inhibits coalescence, even at low boat speeds. The void fraction increases with speed beyond 10 knots and peaks near the wall. Bubble velocities show slip with the wall at all speeds and exhibit large RMS fluctuations, increasing near the wall.     

On the dynamics and instability of bubbles formed during underwater explosions

Dynamics of explosion bubbles formed during underwater detonations are studied experimentally by exploding fuel (hydrogen and/or carbon monoxide)–oxygen mixture in a laboratory water tank. Sub-scale explosions are instrumented to provide detailed histories of bubble shape and pressure. Using geometric and dynamic scaling analyses it has been shown that these sub-scale bubbles are reasonable approximations of bubbles formed during deep sea underwater explosions. The explosion bubble undergoes pulsation and loses energy in each oscillation cycle. The observed energy loss, which cannot be fully explained by acoustic losses, is shown here to be partly due to the excitation of instability at the interface between the gaseous bubble and the surrounding water. Various possible mechanisms for the dissipation of bubble energy are addressed. The analysis of the experimental data gives quantitative evidence (confirmed by recent numerical studies) that the Rayleigh–Taylor instability is excited near the bubble minimum. The dynamics of the bubble oscillation observed in these experiments are in good agreement with experimental data obtained from deep sea explosions     

Experimental Study on Bubble Movement Characteristics during Underwater Pyrotechnic Combustion

High-speed photography was used to study bubble movement characteristics during underwater pyrotechnic combustion. The results show that bubble behaviors include bubble formation at the nozzle, departure from the nozzle, bubble coalescence, and bubble breakup. Compared with cavitation bubbles and fluidization bubbles, the nozzle bubbles formed during underwater pyrotechnic combustion feature larger diameters, up to centimeters, and darker, and more irregular shapes. During large bubble coalescence, two bubbles approach each other, generate a channel for transfer of mass and heat, and finally coalesce. The bubbles contain high-temperature gases and solid residues generated during pyrotechnic combustion, which lead to non-uniform forces on the bubble surface and make the bubbles more prone to breakup. Because of the high-temperature solid grains, the surrounding liquid vaporizes to form bubbles.     

Visualization of capillary boiling

Capillary boiling has been experimentally analyzed using an optical technique, based on the extraction of gray levels at specific locations from video frames. Boiling is achieved by placing a small capillary in a large container filled with water at boiling temperature. A heating wire located inside the capillary provides the additional heat that triggers the evaporation. The phenomenon is videotaped and digitized. A small region of the filmed field is chosen and the gray level of all the digitized frames is stored in a file which is then analyzed. This constitutes in essence the possibility of using a large number of non-intrusive, fixed (Eulerian) virtual sensors. The information extracted concerns the global properties of the phenomenon such as the time between the departure of two subsequent bubbles, and more detailed properties such as the shape and volume of the departing bubbles or even the motion of liquid packets inside the capillary. Moreover, this technique permits the study of time and space correlations. It is observed that the bubbles depart from the capillary tip with a variable frequency that can be approximately grouped around two values. Clear evidence that this is due to two different mechanisms is given. The time between the departure of two bubbles is correlated with the volume of the departing bubbles. Received: 21 April 1999/Accepted: 18 July 2000     

Artificial boiling heat transfer in the free convection to carbonic acid solution

Free convection phenomenon has been experimentally investigated around a horizontal rod heater in carbonic acid solution. Because of the tendency of the solution to desorb carbon dioxide gas when temperature is increased, bubbles appear when cylinder surface is heated. The bubbles consists mainly carbon dioxide and also a negligible amount of water vapor. The results present that dissolved carbon dioxide in water significantly enhances the heat transfer coefficient in compare to pure free convection regime. This is mainly due to the microscale mixing on the heat transfer surface, which is induced by bubble formation. In this investigation, experiments are performed at different bulk temperatures between 288 and 333 K and heat fluxes up to 400 kW m⁻² at atmospheric pressure. Bubble departure diameter, nucleation site density and heat transfer coefficient have been experimentally measured. A model has been proposed to predict the heat transfer coefficient.     

Three-dimensional flow pattern visualization and bubble size distributions in stationary and transient upward flashing flow

For the first time, an experimental three-dimensional reconstruction and visualization of stationary and transient flashing flow in a vertical pipe (47 mm diameter) is presented. The measurements have been performed by means of wire-mesh sensors. This type of sensor delivers two-dimensional void-fraction distributions in the pipe cross-section where it is mounted with a maximum sampling rate of 10,000 frames per second. A sampling rate of 1200 frames per second has been used in this work. Steam bubbles have been identified from the wire-mesh data and their complete three-dimensional reconstruction has been performed by taking into account the steam bubble velocity. For the estimation of the bubble velocity, two wire-mesh sensors positioned at a small axial distance from each other have been used. The velocity has been determined by cross-correlation of the two wire-mesh signals, by direct identification of the traveling time of the steam bubbles between the two sensors and by means of a drift-flux model. A comparison between the three methods of bubbles velocity measurement is reported. Stationary and time-dependent bubble size distributions have been derived. The stationary radial void-fraction profiles have been decomposed according to bubble size classes and compared with the results obtained with an equilibrium model.     

2D lattice Boltzmann investigation of saturated pool boiling using a tunable surface tension model: Prandtl number effects on film boiling

In this paper, the saturated pool boiling is investigated using lattice Boltzmann method. The written FORTRAN code is validated in two aspects: For flow, the thermodynamic consistency test and Laplace law are applied and for heat transfer, the space- and time- averaged Nusselt number is compared with Berenson analytical solution in film boiling regime. In addition, the results of bubble generation and departure are compared with some well-known analytical solutions to show the accuracy of the code. It is confirmed that bubble departure diameter and the departure frequency are related to the gravity acceleration with powers of ? 0.505 and 0.709, respectively, which is in a very good agreement with the existing analytical expressions. The present model has the ability to tune different surface tensions independent of liquid/vapor density ratio, which was unreachable using other existing numerical models of boiling. Thus, the sole effects of surface tension on boiling can also be taken into consideration using the present model. It is also shown that the departure diameter is related to the surface tension with a power of 0.485, which is in good agreement with the analytical expressions. Temperature contours are shown together with flow lines to have a better viewpoint for studying the bubble’s behavior. An intensive temperature gradient is observed in the necking area at the departure time. All the four boiling regimes in the boiling curve are simulated under constant temperature boundary condition. The Prandtl number effects on vapor bubble dynamics in the film boiling regime are investigated using the improved Shan and Chen model for the first time. Results revealed that bubbles are more resistant to depart from the vapor blanket with increasing the Prandtl number.     

Observation of multibubble phenomena in an ultrasonic reactor

A cylindrical ultrasonic reactor radiated by a BLT (bolted Langevin transducer) with a large horn operated at 27.3 kHz was employed to investigate the behaviors of bubbles in the present research. The characteristics of this reactor were analyzed using experimental methods. It was found that traveling transient bubbles can eliminate the effect of standing wave on cleaning efficiency. Fresh tap water and degassed water were employed as test liquids. In tap water, a conical cavitation structure was observed as the input voltage of transducer was increased over 100 V. The detail of the conical cavitation structure was anatomized and its shape varied with the change of input voltage. In degassed water, ultrasonic fireworks were observed and special attention was paid to their subsequences, snowflake bubble clouds and water jet.     

Injection de bulles en paroi d'un écoulement cisaillé visqueuxBubble injection at the wall of a viscous shear flow

Air bubble injection at the wall of a viscous shear flow is filmed using a high speed video camera. The temporal evolution of the bubbles equivalent radius and the position of their centre of mass throughout their growth are determined from image processing. The experimental results are then used to validate a model of forces acting on the bubbles during their growth and after their detachment within the limit of small bubble Reynolds numbers. To cite this article: G. Duhar, C. Colin, C. R. Mecanique 331 (2003).     

Liquid flow patterns about barbotage bubbles

This paper summarizes the results of a flow visualization study on the liquid motion around barbotage bubbles during growth and departure. Flow patterns, as well as for the first time, instantaneous velocities, are reported as a function of time and location about the bubbles. The experiments, employing the hydrogen-bubble technique and high-speed cine photography, were with: water as the liquid, air as the bubbled gas, orifice diameters of 0.116 and 0.252 cm, and different air flow rates; the two limiting cases of constant supply pressure and constant volumetric flow rate were covered. It was found that the liquid around a barbotage bubble assumes two velocity maxima, the first an outward maximum during bubble growth and the second in the opposite direction approximately at the time of bubble departure; further, liquid velocities were found to be higher close to the bubbling site. Certain differences in liquid velocities between the constant pressure and constant flow cases are explained in terms of available theoretical solutions to the bubble growth rate. Qualitative comparisons of the barbotage liquid flow patterns and those recently reported for boiling flow patterns are also presented.     

Opposed bubbly jets at different impact angles: Jet structure and bubble properties

The structure of two colliding water jets containing small gas bubbles is studied experimentally. The effects of the separation distance between jets, as well as the orientation angle, on the spatial distribution of bubbles have been considered. Results on the global structure of the final jet and bubble properties have been obtained using a high-speed video camera, and measurements of the positions of coalescence events are presented. Jets are introduced through inclined pipes (with a diameter of 0.7 mm) into a large water tank to avoid wall effects. Inclination angle has been changed from 0° to 45° with respect to the horizontal, resulting in a 0° up to 90° impact angle between jets. Generation of bubbles is controlled by a T-junction device where a regular slug-flow is created prior to injection. Bubble sizes have been measured, and a mean diameter of around 1 mm has been obtained using high values of the liquid flow rate. In the studied range of separation distances between the bubbly jets, a more homogeneous dispersion of bubbles is created as the distance between jets is decreased and the momentum flux of each jet is increased. Higher numbers of coalescences are observed when using smaller distance between jets, and the obtained measurements revealed that the number of bubble coalescence events is reduced significantly using high values of liquid flow rates.