Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

Wake properties of a single gas bubble in three-dimensional liquid-solid fluidized bed

Experiments were performed to investigate the wake properties of a single gas bubble in a three-dimensional liquid-solid fluidized bed via a video camera moving at the same speed as the bubble. The solids holdup in the fluidized bed varied up to around 10%. The bubble size varied from 5 to 20 mm with corresponding bubble Reynolds numbers ranging from 1000 to 6500. The bubble was observed to have two types of wake configurations depending on the bubble size: the asymmetric/helical vortex wake for small bubbles and the symmetric wake for large bubbles. The bubble shape and relative rise velocity in the fluidized bed can be well-represented by correlations developed for single bubbles in liquid media, although the bubble shape in liquid-solid media is slightly more flattened compared to that in liquid media. The bubble rocking frequency was found to be independent of particle properties and to correspond in magnitude to the vortex shedding frequency in a two-dimensional liquid-solid fluidized bed. The average primary wake size in three dimensions is comparable to that in two dimensions.     

Dispersion of Gas Bubbles in Large-Scale Homogeneous Isotropic Turbulence

An approximate analysis is given of the dispersion of gas bubbles that rise at large Reynolds number through large-scale homogeneous, isotropic turbulence, characterized by the Kraichnan energy-spectrum function. A fairly well-established equation of motion of the bubbles, originally proposed by Thomas et al. [16], is used to derive a closed set of equations for the components of the dispersion tensor of the bubbles in a manner analogous to that used by Saffman [12] for fluid particles and by Pismen and Nir [10, 11] for solid particles. The equations are then solved to obtain the diffusivities and the intensities of bubble velocity fluctuations. Analytical solutions are compared with results from simulations of the bubble motion in a Gaussian random velocity field.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

Flow structures and dynamics in the wakes of sliding bubbles

An experimental investigation is reported for the flow structures in the wake of an air bubble sliding under an inclined surface in quiescent water. Time-resolved particle image velocimetry (PIV) is used to study the wakes of sliding bubbles for a range of measurement planes, bubble diameters and surface inclination angles. Additionally, key aspects of the bubble’s motion are measured simultaneously using a novel method that accounts for the motion of the bubble’s interface. Thus, vortex shedding may be linked to changes in the bubble shape and path.Analysis of the measured velocity and vorticity fields reveals a wake structure consisting of a near wake that moves in close proximity to the bubble, shedding vorticity at the inversion points of the bubble path. Downstream of the bubble in the far wake, these structures evolve into asymmetrical, oppositely-oriented hairpin vortices that are generated in the near wake. These hairpin vortices bear similarities to those observed behind freely rising bubbles and near-wall bluff bodies and are found to cause significant motion of the bulk fluid. This bulk fluid motion has the potential to offer significant convective cooling of adjacent heated surfaces, such as submerged electronics components.     

均匀电场中气泡上升特性的实验研究

电场中的气泡对于强化传热有显著的作用, 对于电场中气泡动力学特性的研究对增强换热器效率, 提高能源利用率有重要意义. 为了获得外加电场作用下气泡的动力学特性, 设计与搭建了可视化实验平台. 采用50 kV高压直流电源构建均匀电场, 高清摄像机拍摄实验图像. 引入电场强度、气泡体积与溶液介电常数作为变量, 探究其对于气泡动力学特性的影响. 观测了竖直与水平均匀电场中气泡的上升过程, 分析了不同变量下气泡变形状况与上升速度的变化. 引入气泡长宽比L/D用于表示气泡拉伸变形程度, 截取单个气泡上升过程分时段图像展示形态变化过程. 研究结果表明, 气泡沿电场方向伸长, 且电场强度越大, 变形越明显; 竖直电场中气泡伸长导致上升速度增大, 而水平电场中气泡上升速度减小. 气泡尺寸增大, 浮升力作用增强, 气泡上升速度增大. 溶液介电常数增加, 电场力作用明显增加, 气泡变形更加明显.      

Simulating air entrainment and vortex dynamics in a hydraulic jump

The air entrainment characteristics of three separate Froude number hydraulic jumps are investigated numerically using an unsteady RANS, realizable k–ε turbulence model, with a Volume of Fluid treatment for the free surface. Mean velocity profiles, average void fraction, and Sauter mean diameter compare favorably with experimental data reported in literature. In all simulations, time-averaged void fraction profiles show good agreement with experimental values in the turbulent shear layer and an accurate representation of interfacial aeration at the free surface. Sauter mean diameter is well represented in the shear layer, and free surface entrainment results indicate bubble size remains relatively unchanged throughout the depth of the jump. Several different grid resolutions are tested in the simulations. Significant improvements in void fraction and bubble size comparison are seen when the diameter to grid size ratio of the largest bubbles in the shear layer surpasses eight. A three-dimensional simulation is carried out for one Froude number jump, showing an improvement in the prediction of entrained air and bubble size compared with two-dimensional results at a substantial increase in computation time. An analysis of three-dimensional vorticity shows a complex interaction between spanwise and streamwise vortical structures and entrained air bubbles. The jump is similar to a turbulent mixing layer, constrained by the free surface, with vortex pairing and subsequent fluctuations in free surface elevation. Downstream fluctuations of the toe are associated with a roll up of the primary spanwise vortex, fluctuations of the free surface, and counter-rotating streamwise vortex pairs. The action of these flow structures is likely responsible for the improvement in three-dimensional results.     

The transient buoyancy driven motion of bubbles across a two-dimensional quiescent domain

The transient buoyancy driven motion of two-dimensional bubbles across a domain bounded by two horizontal walls is studied by direct numerical simulations. The bubbles are initially released next to the lower wall and as they rise, they disperse. Eventually all the bubbles collect at the top wall. The goal of the study is to examine how a simple one-dimensional model for the averaged void fraction captures the unsteady bubble motion. By using void fraction dependent velocities, where the exact dependency is obtained from simulations of homogeneous bubbly flows, the overall dispersion of the bubbles is predicted. Significant differences remain, however. We suggest that bubble dispersion by the bubble induced liquid velocity must be included, and by using a simple model for the bubble dispersion we show improved agreement.     

Visualization of gas–liquid mass transfer and wake structure of rising bubbles using pH-sensitive PLIF

A planar laser-induced fluorescence (PLIF) technique for visualizing gas–liquid mass transfer and wake structure of rising gas bubbles is described. The method uses an aqueous solution of the pH-sensitive dye Naphthofluorescein and CO₂ as a tracer gas. It features a high spatial resolution and frame rates of up to 500 Hz, providing the ability to capture cinematographic image sequences. By steering the laser beam with a set of two programmable scanning mirrors, sequences of three-dimensional LIF images can be recorded. The technique is applied to freely rising bubbles with diameters between 0.5 and 5 mm, which perform rectilinear, oscillatory or irregular motions. The resulting PLIF image sequences reveal the evolution of characteristic patterns in the near and far wake of the bubbles and prove the potential of the technique to provide new and detailed insights into the spatio-temporal dynamics of mass transfer of rising gas bubbles. The image sequences further allow the estimation of bubble size and rise velocity. The analysis of bubble rise velocities in the Naphthofluorescein solution indicates that surfactant-contaminated conditions are encountered.     

Air transport in chute flows

The air bubble rise velocity in still water depends mainly on the bubble size and is basically influenced by buoyancy, viscosity and surface tension. In high-speed flows the number of forces acting on air bubbles increases with turbulence, non-hydrostatic pressure gradient, shear forces, bubble clouds and free-surface entrainment. Air bubbles in these flows are used for cavitation protection of hydraulic structures such as chutes, spillways and bottom outlets. Here, air is normally added by means of aerators upstream of regions where the cavitation number falls below a critical value mainly to reduce the sonic velocity of the fluid and cushion the cavitation bubble collapse process. The distance between successive aerators depends basically on the bubble rise velocity. Until today, the bubble rise velocity in high-speed flows was not thoroughly investigated because of limited laboratory instrumentation. The present project focused on the streamwise development of air concentrations in high-speed flows along a 14 m long model chute. The bubble rise velocity was indirectly derived from the air detrainment gradient of the air concentration contour lines downstream of an aeration device. It accounts for the main hydraulic parameters chute slope, Froude number and air concentration. It is demonstrated that the bubble rise velocity in high-speed flow and stagnant water differ significantly due to fracturing processes, turbulence, and the ambient air concentration.     

The effect of viscoelasticity on a rising gas bubble

This paper investigates the role of viscoelasticity on the dynamics of rising gas bubbles. The dynamics of bubbles rising in a viscoelastic liquid are characterised by three phenomena: the trailing edge cusp, negative wake, and the rise velocity jump discontinuity. There is much debate in the literature over the cause of the jump discontinuity, which is observed once the bubble exceeds a certain critical volume. In this paper, the employment of some choice modelling assumptions allows insights into the mechanisms of the jump discontinuity which cannot be ascertained experimentally. The ambient fluid is assumed incompressible and the flow irrotational, with viscoelastic effects included through the stress balance on the bubble surface. The governing equations are solved using the boundary element method. Some Newtonian predictions are discussed before investigating the role of viscoelasticity. The model predicts the trademark cusp at the trailing end of a rising bubble to a high resolution. However, the irrotational assumption precludes the prediction of the negative wake. The corresponding absence of the jump discontinuity supports the hypothesis that the negative wake is primarily responsible for the jump discontinuity, as mooted in previous studies.     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

Numerical simulation of dielectric bubbles coalescence under the effects of uniform magnetic field

In this research, the co-axial coalescence of a pair of gas bubbles rising in a viscous liquid column under the effects of an external uniform magnetic field is simulated numerically. Considered fluids are dielectric, and applied magnetic field is uniform. Effects of different strengths of magnetic field on the interaction of in-line rising bubbles and coalescence between them were investigated. For numerical modeling of the problem, a computer code was developed to solve the governing equations which are continuity, Navier–Stokes equation, magnetic field equation and level set and reinitialization of level set equations. The finite volume method is used for the discretization of the continuity and momentum equations using SIMPLE scheme where the finite difference method is used to discretization of the magnetic field equations. Also a level set method is used to capture the interface of two phases. The results are compared with available numerical and experimental results in the case of no-magnetic field effect which show a good agreement. It is found that uniform magnetic field accelerates the coalescence of the bubbles in dielectric fluids and enhances the rise velocity of the coalesced bubble.     

An experimental study of plane bubbles rising at inclination

The characteristics of bubbles rising in two-dimensional tanks is investigated. It is found that the influence of fluid properties, while negligible for vertically rising bubbles, measurably affects the inclined bubble rise velocity. The measured increase in the rise velocity at inclination relative to the vertical value is explained qualitatively. In addition, a comparison is made between the wave analogy correlation of Maneri & Mendelson and vertical rise velocity data. The comparisons show that the correlation predicts the three-dimensional effects of tank spacing when the data are reduced on an equivalent bubble radius basis.     

Stabilization of cylinder wakes in shallow water flows by means of roughness elements: an experimental study

Shallow wakes that occur in a wide range of natural flows tend to generate instabilities that develop into large, 2D coherent structures (2DCS). We present the results of an experimental study to stabilize shallow wakes by local, enhanced bottom roughness. Two successful stabilization strategies are compared to a base case of an unsteady bubble wake. First, localized bed roughness is placed in the lateral shear layers near the shoulders of the cylinder. Second, a local roughness element is placed at the end of the recirculation bubble, in the downstream region where large-scale vortices would normally shed. Dye visualization is used to assess the qualitative behavior of the wake, and two-component laser Doppler velocimetry (LDV) measurements are made to measure the Reynolds stress distributions and time-averaged velocity profiles. In both stabilization cases, a minimum patch size of the enhanced roughness elements is required for stabilization, which depends on the momentum thickness of the shear layers and the locations of enhanced Reynolds shear stresses. The main effect of the wake stabilization is a reduction in momentum exchange with the ambient flow due to damping of the large 2DCS. This reduction in eddy diffusivity results in a narrower wake and a slower decay of the centerline velocity deficit with downstream distance compared to the base case of an unsteady bubble wake.     

UDV measurements of single bubble rising in a liquid metal Galinstan with a transverse magnetic field

Liquid metal is an important type of energy transport carrier in nuclear reactors, such as in accelerator-driven sub-critical systems, fusion reactors and spallation neutron source devices. It is necessary to conduct research for bubbles rising in a liquid metal under different magnetic field intensities. The Perspex container is positioned concentrically inside a transverse magnetic field, which provides a homogeneous DC longitudinal magnetic field that passes through the fluid district. The coils are supplied with maximum field strength of 1.97 T. The equivalent diameter of the bubble is 3.1–5.6 mm. The Ultrasonic Doppler Velocimetry (UDV) method is used to evaluate the internal flow velocity of opaque liquid metals. Research shows that the influence of the Lorenz force on the bubble ascension velocity is not simply positive or negative. The magnetic field inhibits the ascension velocity of small bubbles with diameters of 3.1 mm and 3.4 mm. The terminal velocity for large bubbles with diameters of 4.57 mm, 5.15 mm and 5.6 mm is higher under a weak magnetic field than without a magnetic field. The positive effect happens under strong magnetic intensity. The target is to obtain the hydro-dynamical relationships between the terminal velocity, drag coefficient, the Eötvös number, Reynolds number, and Stuart number in a strong magnetic field using a multiple regression method to reveal that the mechanism of the induced current's restraining influence determines the ascension velocity of the bubble in viscous electric liquids with a strong magnetic field.     

Non-intrusive measurements of bubble size and velocity

A non-intrusive measuring technique based on video-imaging has been developed for the measurement of bubble size, velocity and frequency. Measurements carried out with this method have been compared to those obtained by an optimized phase-Doppler system in standard configuration, for a wide range of bubble sizes produced from single injectors in a quiescent environment. The two measuring techniques have yielded velocities and frequencies that are in very good agreement while the size of spherical bubbles was consistently measured by both methods. The phase-Doppler system was also used to size oblate-spheroidal bubbles moving with their equatorial plane parallel to the scattering plane, yielding measurements reasonably close to the average radius of curvature of the bubbles in the neighborhood of the equatorial plane, as calculated from the video-imaging data. Both methods were used for detailed velocity measurements of the bubble-stream in the neighborhood of the injector tip. The observed bubble-velocity variation with the distance from the injector tip does not always display the usual increasing trend leading into the terminal velocity. When injection conditions are near the transition from discrete to jet injection mode and the bubbles are small, the latter decelerate into a terminal velocity due to direct interaction of successive bubbles at the injector tip. The measured terminal velocities of bubble-chains for a variety of bubble sizes and injection frequencies, are successfully predicted by using a far-field wake approximation to account for the drafting effect which is responsible for bubble-chain velocities higher than those of single bubbles.The authors gratefully acknowledge the support from the Louisiana Board of Reagents through grants LEQSF (1991-1992)-ENH-31 and LEQSF (1990-1992)-RD-A-07, ARCO Chemical Company through it's University Grants Program and the National Science Foundation through its Graduate Fellowship Program.     

Physical and transport properties governing bubble column operations

In this work, physical and transport properties governing fluid dynamics in bubble columns were studied in the heterogeneous liquid circulation regime of operation based on several published experimental data. It was confirmed that the axial dispersion phenomenon of the liquid phase could be defined by the constant Peclet number predicated on the two-phase slip velocity, which comprised the time-averaged bubbly flow velocity superimposed by the velocity of large bubbles rising in the fluid. In addition, the slip velocity of mean bubbles in the bubbly flow not only indicated a constant rising velocity of a particular liquid but also apparently reduced to that for the Higbie model with respect to the bubble diameter and its terminal rising velocity, when applied to Akita and Yoshida's dimensionless correlation for mass transfer operation (Akita and Yoshida, 1974). In the course of demonstration, it was confirmed that a semi-theoretical correlation for the gas holdup proposed by Mersmann (1978) is similar to the empirical correlation proposed by Akita and Yoshida (1973) with respect to the expression and the physical perspective.     

The interaction of background ocean air bubbles with a surface ship

A two-fluid model suitable for the calculation of the two-phase flow field around a naval surface ship is presented. This model couples the Reynolds-averaged Navier–Stokes (RANS) equations with equations for the evolution of the gas-phase momentum, volume fraction and bubble number density, thereby allowing the multidimensional calculation of the two-phase flow for monodisperse variable size bubbles. The bubble field modifies the liquid solution through changes in the liquid mass and momentum conservation equations. The model is applied to the case of the scavenging of wind-induced sea-background bubbles by an unpropelled US Navy frigate under non-zero Froude number boundary conditions at the free surface. This is an important test case, because it can be simulated experimentally with a model-scale ship in a towing tank. A significant modification of the background bubble field is predicted in the wake of the ship, where bubble depletion occurs along with a reduction in the bubble size due to dissolution. This effect is due to lateral phase distribution phenomena and the generation of an upwelling plume in the near wake that brings smaller bubbles up to the surface. © 1998 John Wiley & Sons, Ltd.     

Ellipsoidal model of the rise of a Taylor bubble in a round tube

The rise velocity of long gas bubbles (Taylor bubbles) in round tubes is modeled by an ovary ellipsoidal cap bubble rising in an irrotational flow of a viscous liquid. The analysis leads to an expression for the rise velocity which depends on the aspect ratio of the model ellipsoid and the Reynolds and Eötvös numbers. The aspect ratio of the best ellipsoid is selected to give the same rise velocity as the Taylor bubble at given values of the Eötvös and Reynolds numbers. The analysis leads to a prediction of the shape of the ovary ellipsoid which rises with same velocity as the Taylor bubble.