Comprehensive experimental investigation of the hydrodynamics of large-scale, 3D,oscillating bubble plumes

An extensive study of the most important hydrodynamic characteristics of fairly large-scale bubble plumes was conducted using several measurement techniques and a variety of tools to analyze the data. Particle image velocimetry (PIV), double-tip optical probes (OP) and photographic techniques were extensively applied to measure bubble and liquid velocities, void-fraction and bubble sizes. PIV measurements in a vertical plane crossing the centre of the injector provided the instantaneous velocity fields for both phases, as well as hydrodynamic parameters, such as the movement of the axis of the plume and its instantaneous shape. Statistical studies were performed using image processing to determine the distribution of the apparent instantaneous plume diameter and centreline position. An important finding was that there is little instantaneous spreading of the bubble plume core; the spreading of the time-averaged plume width (as measured from the time-averaged void-fraction and time-averaged liquid velocity fields) is largely due to plume meandering and oscillations. The liquid-phase stress tensor distributions obtained from the instantaneous velocity data indicate that, for the continuous phase, these stresses scale linearly with the local void-fraction in the range of 0.5% < α < 2.5%. The bubbles were found to be ellipsoidal, with shape factor e ≈ 0.5.     

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.     

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.     

Gas dispersion in horizontal pulp-fibre-suspension flow

Dispersion of gas into pulp-suspension horizontal flow was investigated downstream of 90° tees for ranges of fibre mass concentrations (0–3.0%), superficial liquid/pulp velocities (0.5–5.0 m/s) and superficial gas velocities (0.11–0.44 m/s) based on a gas mixing index, derived from the standard deviation of cross-sectional local gas holdup obtained from electrical resistance tomography. Mixing for dilute suspensions was similar to that for water, but differed significantly for higher suspension concentrations. Mixing worsened with increasing fibre mass concentration for the bubble flow regime, likely due to dense fibre networks in the core of the pipe causing bubbles to congregate near the wall. When buoyancy was significant, gas uniformity improved with increasing pulp concentration, since robust fibre networks caused liquid/pulp slugs to flow at the top of the pipe, whereas stratified flow was approached at lower concentrations. Mixing was less dependent on superficial liquid/pulp velocity at higher pulp concentrations, due to less variation in flow regimes.     

Dynamics of a Taylor bubble in steady and pulsatile co-current flow of Newtonian and shear-thinning liquids in a vertical tube

A computational analysis is carried out to ascertain the effects of steady and pulsatile co-current flow, on the dynamics of an air bubble rising in a vertical tube containing water or a solution of Carboxymethylcellulose (CMC) in water. The mass fraction (m_f) of CMC in the solution is varied in the range 0.1% ⩽ m_f ⩽ 1% to accommodate zero-shear dynamic viscosities in the range 0.009–2.99 Pa-s. It was found that the transient and time-averaged velocities of Taylor bubbles are independent of the bubble size under both steady as well as pulsatile co-current flows. The lengths of the Taylor bubbles under the Newtonian conditions are found to be consistently greater than the corresponding shear-thinning non-Newtonian conditions for any given zero-shear dynamic viscosity of the liquid. In contrast to observations in stagnant liquid columns, an increase in the dynamic viscosity of the liquid (under Newtonian conditions) results in a concomitant increase in the bubble velocity, for any given co-current liquid velocity. In shear-thinning liquids, the change in the bubble velocity with an increase in m_f is found to be relatively greater at higher co-current liquid velocities. During pulsatile shear-thinning flows, distinct ripples are observed to occur on the bubble surface at higher values of m_f, the locations of which remain stationary with reference to the tube for any given pulsatile flow frequency, while the bubble propagated upwards. In such a pulsatile shear-thinning flow, a localised increase in dynamic viscosity is accompanied near each ripple, which results in a localised re-circulation region inside the bubble, unlike a single re-circulation region that occurs in Newtonian liquids, or shear-thinning liquids with low values of m_f. It is also seen that as compared to frequency, the amplitude of pulsatile flow has a greater influence on the oscillating characteristics of the rising Taylor bubble. The amplitude of oscillation in the bubble velocity increases with an increase in the CMC mass fraction, for any given value of pulsatile flow amplitude.     

Effect of gas expansion on the velocity of a Taylor bubble: PIV measurements

The effect of gas expansion on the velocity of a Taylor bubble was studied experimentally. The velocity field in the liquid ahead of a Taylor bubble was measured by particle image velocimetry (PIV), and the bubble velocity was measured with two pairs of laser diodes and photocells. The experiments were done in a 7.0 m long vertical tube with a 32 mm internal diameter. Solutions of carboxymethylcellulose (CMC) polymer with weight percentages between 0.01% and 0.1% were used. The expansion of slug gas induces an increase in the bubble velocity and a corresponding displacement of the liquid ahead of the bubble. The velocity of the bubble increases by an amount equal to the maximum velocity in the liquid displaced. For the solutions studied, the induced velocity profile was parabolic and the bubble velocity increase was equal to the liquid velocity at the tube axis, i.e., twice the mean velocity in the liquid displaced. The corrected velocity obtained by subtracting the velocity increase from the value of the bubble velocity is independent of the bubble length.     

One-dimensional interfacial area transport of vertical upward bubbly flow in narrow rectangular channel

The design and safety analysis for miniature heat exchangers, the cooling system of high performance microelectronics, research nuclear reactors, fusion reactors and the cooling system of the spallation neutron source targets requires the knowledge of the gas–liquid two-phase flow in a narrow rectangular channel. In this study, flow measurements of vertical upward air–water flows in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm were performed at seven axial locations by using the imaging processing technique. The local frictional pressure loss gradients were also measured by a differential pressure cell. In the experiment, the superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. The developing two-phase flow was characterized by the significant axial changes of the local flow parameters due to the bubble coalescence and breakup in the tested flow conditions. The existing two-phase frictional multiplier correlations such as Chisholm, 1967, Mishima et al., 1993 and Lee and Lee (2001) were verified to give a good prediction for the measured two-phase frictional multiplier. The predictions of the drift-flux model with the rectangular channel distribution parameter correlation of Ishii (1977) and several existing drift velocity correlations of Ishii, 1977, Hibiki and Ishii, 2003 and Jones and Zuber (1979) agreed well with the measured void fractions and gas velocities. The interfacial area concentration (IAC) model of Hibiki and Ishii (2002) was modified by taking the channel width as the system length scale and the modified IAC model could predict the IAC and Sauter mean diameter acceptably.     

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.     

Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column

A new approach for simulating the formation of a froth layer in a slurry bubble column is proposed. Froth is considered a separate phase, comprised of a mixture of gas, liquid, and solid. The simulation was carried out using commercial flow simulation software (FIRE v2014) for particle sizes of 60–150 μm at solid concentrations of 0–40 vol%, and superficial gas velocities of 0.02–0.034 m/s in a slurry bubble column with a hydraulic diameter of 0.2 m and height of 1.2 m. Modelling calculations were conducted using a Eulerian–Eulerian multiphase approach with k–ε turbulence. The population balance equations for bubble breakup, bubble coalescence rate, and the interfacial exchange of mass and momentum were included in the computational fluid dynamics code by writing subroutines in Fortran to track the number density of different bubble sizes. Flow structure, radial gas holdup, and Sauter mean bubble diameter distributions at different column heights were predicted in the pulp zone, while froth volume fraction and density were predicted in the froth zone. The model was validated using available experimental data, and the predicted and experimental results showed reasonable agreement. To demonstrate the effect of increasing solid concentration on the coalescence rate, a solid-effect multiplier in the coalescence efficiency equation was used. The solid-effect multiplier decreased with increasing slurry concentration, causing an increase in bubble coalescence efficiency. A slight decrease in the coalescence efficiency was also observed owing to increasing particle size, which led to a decrease in Sauter mean bubble diameter. The froth volume fraction increased with solid concentration. These results provide an improved understanding of the dynamics of slurry bubble reactors in the presence of hydrophilic particles.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

Evaluation of flow patterns and elongated bubble characteristics during the flow boiling of halocarbon refrigerants in a micro-scale channel

In the present study, quasi-diabatic two-phase flow pattern visualizations and measurements of elongated bubble velocity, frequency and length were performed. The tests were run for R134a and R245fa evaporating in a stainless steel tube with diameter of 2.32 mm, mass velocities ranging from 50 to 600 kg/m² s and saturation temperatures of 22 °C, 31 °C and 41 °C. The tube was heated by applying a direct DC current to its surface. Images from a high-speed video-camera (8000 frames/s) obtained through a transparent tube just downstream the heated sections were used to identify the following flow patterns: bubbly, elongated bubbles, churn and annular flows. The visualized flow patterns were compared against the predictions provided by Barnea et al. (1983) [1], Felcar et al. (2007) [10], Revellin and Thome (2007) [3] and Ong and Thome (2009) [11]. From this comparison, it was found that the methods proposed by Felcar et al. (2007) [10] and Ong and Thome (2009) [1] predicted relatively well the present database. Additionally, elongated bubble velocities, frequencies and lengths were determined based on the analysis of high-speed videos. Results suggested that the elongated bubble velocity depends on mass velocity, vapor quality and saturation temperature. The bubble velocity increases with increasing mass velocity and vapor quality and decreases with increasing saturation temperature. Additionally, bubble velocity was correlated as linear functions of the two-phase superficial velocity.     

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.     

The effect of aspect ratio in counter-current gas-liquid bubble columns: Experimental results and gas holdup correlations

It is generally admitted that the gas holdup is independent of the column dimensions and gas sparger design if three criteria are satisfied: the diameter of the bubble column is larger than 0.15 m, gas sparger openings are larger than 1–2 mm and the aspect ratio is larger than 5. This paper contributes to the existing discussion; in particular, the effect of the aspect ratio (within the range 1–15) in a counter-current gas-liquid bubble column has been experimentally studied and a new gas holdup correlation to estimate the influence of aspect ratio, operation mode and working fluid on the gas holdup has been proposed. The bubble column, equipped with a spider gas sparger, is 5.3 m in height, has an inner diameter of 0.24 m; gas superficial velocities in the range of 0.004–0.23 m/s have been considered, and, for the runs with water moving counter-currently to the gas phase, the liquid has been recirculated at a superficial velocity of −0.0846 m/s. Filtered air has been used as the gaseous phase in all the experiments, while the liquid phase has included tap water and different aqueous solutions of sodium chloride as electrolyte. Gas holdup measurements have been used to investigate the flow regime transitions and the global bubble column hydrodynamics. The counter-current mode has turned out to increase the gas holdup and destabilize the homogeneous flow regime; the presence of electrolytes has resulted in increasing the gas holdup and stabilizing the homogeneous flow regime; the aspect ratio, up to a critical value, has turned out to decrease the gas holdup and destabilize the homogeneous flow regime. The critical value of the aspect ratio ranged between 5 and 10, depending on the bubble column operation (i.e., batch or counter-current modes) and liquid phase properties. Since no correlation has been found in the literature that can correctly predict the gas holdup under the investigated conditions, a new scheme of gas holdup correlation has been proposed. Starting from considerations concerning the flow regime transition, corrective parameters are included in the gas holdup correlation to account for the effect of the changes introduced by the aspect ratio, operation mode and working fluid. The proposed correlation has been found to predict fairly well the present experimental data as well as previously published gas holdup data.     

One-group interfacial area transport equation and its sink and source terms in narrow rectangular channel

The characteristics of two-phase flow in a narrow rectangular channel are expected to be different from those in other channel geometries, because of the significant restriction of the bubble shape which, consequently, may affect the heat removal by boiling under various operating conditions. The objective of this study is to develop an interfacial area transport equation with the sink and source terms being properly modeled for the gas–liquid two-phase flow in a narrow rectangular channel. By taking into account the crushed characteristics of the bubbles a new one-group interfacial area transport equation was derived for the two-phase flow in a narrow rectangular channel. The random collisions between bubbles and the impacts of turbulent eddies with bubbles were modeled for the bubble coalescence and breakup respectively in the two-phase flow in a narrow rectangular channel. The newly-developed one-group interfacial area transport equation with the derived sink and source terms was evaluated by using the area-averaged flow parameters of vertical upwardly-moving adiabatic air–water two-phase flows measured in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm. The flow conditions of the data set covered spherical bubbly, crushed pancake bubbly, crushed cap-bubbly and crushed slug flow regimes and their superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. Good agreement with the average relative deviation of 9.98% was obtained between the predicted and measured interfacial area concentrations in this study.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Experimental investigation of the influence of column scale,gas density and liquid properties on gas holdup in bubble columns

Measurements of gas holdups in bubble columns of 0.16, 0.30 and 0.33 m diameter were carried out. These columns were operated in co-current flow of gas and liquid phases and in semibatch mode. The column of 0.33 m diameter was operated at elevated pressures of up to 3.6 MPa. Nitrogen was employed as the gas phase and deionized water, aqueous solutions of ethanol and acetone and pure acetone and cumene as the liquid phase. The effects of differing liquid properties, gas density (due to elevated pressure), temperature, column diameter and superficial liquid velocity on gas holdup were studied. The gas holdup measurements were utilized by differential pressure measurements at different positions along the height of the bubble columns which allowed for the identification of axial gas holdup profiles. A decrease of gas holdup with increasing column diameter and an increase of gas holdup with increasing pressure was observed. The effect of a slightly decreasing gas holdup with increasing liquid velocity was found to exist at smaller column diameters. The use of organic solvents as the liquid phase resulted in a significant increase in gas holdup compared to deionized water. It is found that published gas holdup models are mostly unable to predict the results obtained in this study.     

Study of material agglomeration during nozzle injection of multiviscosity liquid into a fluidized bed reactor by the electric conductance method

Fluidized bed agglomeration is an important and challenging problem for thermal cracking in fluid cokers. A low coker temperature can be problematic because the bitumen is injected into the fluidized bed with a different viscosity, resulting in formation of agglomerates of varying sizes, which slows the cracking reactions. In the present study, the bed material agglomeration process during nozzle injection of multiviscosity liquid was investigated in a fluidized bed operated at different mass ratios of the atomization gas to the liquid jets (GLR = 1%–3.5%) and gas velocities (3.9U_mf and 5.9U_mf) based on a conductance method using a water–sand system to simulate the hot bitumen–coke system at room temperature. During the tests of liquid-jet dispersion throughout the bed, different agglomeration stages are observed at both gas velocities. The critical amount of tert-butanol in the liquid jets that could lead to severe agglomeration of the bed materials (poor fluidization) at GLR = 1% is about 10 wt% at the low fluidizing gas velocity (3.9U_mf) and 18 wt% at the high gas velocity (5.9U_mf). This study provides a new approach for on-line monitoring of bed agglomeration during liquid injection to guarantee perfect contact between the atomized liquid and the bed particles.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Hot-film anemometry in bubbly flow I: bubble–probe interaction

The interaction between a bubble, which is rising in a descending water flow, and a hot-film anemometer was experimentally investigated using stereoscopic high-speed imaging. The mean downward water velocity varied from 0 up to 0.15 m/s, i.e., relatively low, allowing for an extended bubble–probe interaction. Moreover, the direction of the water causes the wake of the probe to play a role before the bubble touches the probe. The equivalent bubble radii were 0.4–2.8 mm and the bubble velocities relative to the probe ranged from 0.04 to 0.38 m/s. Image processing techniques were applied to reconstruct the bubbles’ path, shape, and orientation during the interaction process. As a result, three types of interactions were found, namely penetrating, bouncing, and splitting interactions. The image sequences were compared with the corresponding time series of the hot-film anemometer. From the time series the type of interaction cannot be deduced, at least not for the analyzed flow situation. Furthermore, we demonstrated that the residence time estimate from the hot-film data is systematically biased in our type of experiments. Finally, it was found that the velocity of a bubble may be altered considerably due to the interaction.     

Computational study of steady streaming from oscillating microbubbles with uniform and wavy wall motions

Steady streaming flow fields of a 5 μm bubble oscillating with uniform radial wall motion and a 500 μm bubble oscillating with wavy wall motion were simulated using a computational fluid dynamics method that incorporated fluid–structure interactions. The steady streaming flow fields for both bubbles were calculated, and they exhibited upward jet flow with two symmetrical counter-rotating vortices. The maximum streaming velocity ranged from a few to tens of millimeters per second. The simulated flow fields were compared with the theory and experimental measurements using particle image velocimetry. The simulation results agreed well with the theoretical and experimental data. Therefore, the proposed computational method would provide a useful tool to predict steady streaming flow fields of oscillating bubbles.