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.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

Effects of hydrophilic particles on bubbly flow in slurry bubble column

To investigate the effects of hydrophilic particles on slurry bubble flows in a bubble column, distributions of the local gas holdup and the bubble frequency are measured using an electric conductivity probe. Particles are made of silica and their diameter is 100 μm. The particle volumetric concentration C_S is varied from 0 to 0.40. The measured data imply that the presence of particles promotes bubble coalescence. The film drainage time for two coalescing bubbles in a quasi two-dimensional bubble flow in a small vessel is also measured to quantitatively evaluate the particle effect on coalescence. A particle-effect multiplier is introduced into a coalescence efficiency model by taking into account the data of film drainage time and is implemented into a multi-fluid model. The main conclusions obtained are as follows: (1) the local gas holdup and bubble frequency in slurry bubble flows decrease with increasing the particle concentration, (2) the hydrophilic particles enhance bubble coalescence and the enhancement saturates at C_S ≃ 0.45, (3) the particle effect on coalescence is well accounted for by introducing the particle-effect multiplier to the film drainage time, and (4) the multi-fluid model can give good predictions for the distribution of the local gas holdup in the slurry bubble column.     

Particle and bubble dynamics in a creeping flow

Dynamics of a solid particle and non-deformable gaseous bubble in viscous fluid are studied analytically and numerically within the framework of creeping flow regime (flow at vanishingly small Reynolds numbers). Equations of motion for the particle and bubble include the consideration of the buoyancy force, Stokes drag force and memory-integral drag force. Exact analytical solutions are obtained and categorised in terms of inclusion (particle or bubble) density with respect to the density of a surrounding fluid. Through the analytical and numerical solutions, the dynamics of solid particle and air bubble in water have been found to behave differently especially at the early stages of motion, whereas some qualitative similarities exist in the long-term asymptotic.     

Visualization of a confined accelerated bubble

High-speed photography was used to study the collapse of a confined two-dimensional, air cavity in water, subjected to a propagating pressure disturbance. The 5–6 mm diameter cavity was confined in a rectangular duct. A sustained pressure disturbance was created by an accelerating piston in contact with the water 240 mm away from the bubble. The pressure increased from 0.1 MPa to about 0.12 MPa with a rise time of the order of 2 ms. The pressure pulse was not reflected until its arrival at the end of the duct, 320 mm from the piston. A microjet was produced at the proximal wall which penetrated the distal cavity wall, thereby producing a pair of bubbles which was thought to be regions of intense vorticity. The features of such confined bubble collapse were not found in previous investigations of unconfined bubble accelerations by weak pressure disturbances. Confinement apparently intensified the effect of the disturbance significantly. Received 18 August 1998 / Accepted 12 May 1999     

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.     

Numerical simulation of a flow around an impulsively started radially deforming circular cylinder

The development of a two‐dimensional viscous incompressible flow generated by a deformable circular cylinder impulsively started into rectilinear motion is studied numerically for the Reynolds numbers equal to 550 and 3000. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the streamfunction is solved by a Fourier–Galerkin method in the direction of the flow that is assumed to remain symmetrical and a second‐order finite difference for the radial direction. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The computed results are compared qualitatively with experimental and numerical results done before in the particular non‐deformable case. The comparison is found to be satisfactory. The influence of the deformation of the cylinder on the flow structure and the drag coefficient is then analyzed. Copyright © 1999 John Wiley & Sons, Ltd.     

Behavior of a small single bubble rising in a rotating flow field

A small single bubble was generated with a single-hole nozzle facing upward in a water bath contained in a rotating cylindrical vessel. The bubble size falls in the surface tension force dominant regime. The vertical, radial, and tangential migration velocities of the bubble were measured with two CCD cameras and a high-speed video camera. The tangential velocity component of water flow was measured with particle image velocimetry. A helical motion of the bubble was observed under every experimental condition. The direction of the helical motion was the same as that of the tangential velocity component. This helical motion is associated with the large initial shape deformation of the bubble near the nozzle exit and the subsequent regular shedding of vortices behind it. The period and amplitude of the helical motion were obtained by analyzing the trajectory of the bubble. These quantities were non-dimensionalized by the volume equivalent bubble diameter and the terminal bubble velocity in the vertical direction and correlated as functions of the Eotvos number. Empirical equations were proposed for the period and amplitude. Originally published in the Journal of JSEM, Vol. 4, No. 2, pp. 38–45 (2004).     

Spin-up and wrapping of cryogenic helium bubble around dewar well from rest in microgravity

Time dependent animation of spherically shaped bubble sitting at one-side of the dewar container, which is spin-up from rest to certain rotating speed of interests, has been investigated. Time dependent deformation of bubble with or without a completion of wrapping around the inner well of the rotating dewar, are numerically studied and simulated. Some similarity parameters are considered for a completion of bubble wrapping around the dewar well. Illustrative examples show that the degree of the completion of bubble wrapping around the dewar well increases with increasing container rotating speed and Weber number.     

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.     

A simultaneous PIV and heat transfer study of bubble interaction with free convection flow

Whole field velocity and point temperature and surface heat flux measurements were performed to characterise the interaction of a single rising ellipsoidal air bubble with the free convection flow from a heated flat surface immersed in water at different angles of inclination. Two thermocouples and a hot film sensor were used to characterise heat transfer from the surface, while a time-resolved digital particle image velocimetry technique was used to map the bubble induced flow in a plane parallel to the surface. Heat flux fluctuations, preceding and following the bubble passage, were shown to correlate with the variation in both local flow velocities and fluid temperatures. The largest increases in heat transfer were recorded when both flow and temperature effects combined to enhance the convective cooling simultaneously. Such conditions were shown to be most likely met when the block was inclined at 45°, thus forcing the bubble to slide closer to the heated surface and hence to the thermal boundary layer.     

Interfacial aeration and bubble count rate distributions in a supercritical flow past a backward-facing step

An example of high-velocity open channel flows is a supercritical flow past an abrupt drop. In such a geometry, the basic air–water flow properties were measured, including distributions of void fraction and bubble count rate, and local air and water chord size distributions, at and downstream of the backward-facing step. The bubble count rate distributions were compared with a conceptual model of streamwise distribution of air and water chords which yields a quasi-parabolic relationship between bubble count rate and void fraction. The proposed model was an attempt to explain the experimental relationship between bubble count rate and void fraction, rather a meticulous breakdown of the complex air–water structure.     

Particle tracking velocimetry of the flow field around a collapsing cavitation bubble

The velocity field in the vicinity of a laser-generated cavitation bubble in water is investigated by means of particle tracking velocimetry (PTV). Two situations are explored: a bubble collapsing spherically and a bubble collapsing aspherically near a rigid wall. In the first case, the accuracy of the PTV method is assessed by comparing the experimental data with the flow field around the bubble as obtained from numerical simulations of the radial bubble dynamics. The numerical results are matched to the experimental radius–time curve extracted from high-speed photographs by tuning the model parameters. Trajectories of tracer particles are calculated and used to model the experimental process of the PTV measurement. For the second case of a bubble collapsing near a rigid wall, both the bubble shape and the velocity distribution in the fluid around the bubble are measured for different standoff parameters γ at several instants in time. The results for γ > 1 are compared with the corresponding results of a boundary-integral simulation. For both cases, good agreement between simulation and experiment is found.     

Void fraction measurement of bubble and slug flow in a small channel using the multivision technique

Using the multivision technique, a new void fraction measurement method was developed for bubble and slug flow in a small channel. The multivision system was developed to obtain images of the two-phase flow in two perpendicular directions. The obtained images were processed—using image segmentation, image subtraction, Canny edge detection, binarization, and hole filling—to extract the phase boundaries and information about the bubble or slug parameters. With the extracted information, a new void fraction measurement model was developed and used to determine the void fraction of the two-phase flow. The proposed method was validated experimentally in horizontal and vertical channels with different inner diameters of 2.1, 2.9, and 4.0 mm. The proposed method of measuring the void fraction has better performance than the methods that use images acquired in only one direction, with a maximum absolute difference between the measured and reference values of less than 6%.     

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.     

A new approach to potential flow around axisymmetric bodies

A new method is introduced to solve potential flow problems around axisymmetric bodies. The approach relies on expressing the infinite series expansion of the Laplace equation solution in terms of a finite sum which preserves the Laplace solution for the potential function under a Neumann-type boundary condition. Then the coefficients of the finite sum are calculated in a least squares approximation sense using the Gram-Schmidt orthonormalization method. Sample benchmark problems are presented and discussed in some detail. The solutions are accurate and converged faster when a rather small number of terms were used. The method is simple and can be easily programmed.     

A rising bubble in a polymer solution

We propose a boundary integral method to study the shape of a bubble rising under gravity in a dilute polymer solution. Constitutive properties are modelled using a FENE model [M.D. Chilcott, J.M. Rallison, J. Non-Newtonian Fluid Mech. 29 (1988) 381] with a pure surface tension interface. We employ a birefringent strand representation [O.G. Harlen, J.M. Rallison, M.D. Chilcott, High-Deborah-number flows of dilute polymer, J. Non-Newtonian Fluid Mech.34 (1990) 319–349] of the wake to simulate the shape and the time-dependent motion of the bubble. Steady and non-steady solutions reproduce qualitatively the bubble deformation seen in experiment with a small region of very high curvature near the rear stagnation point of the bubble. We find a limit point for steady axisymmetric solutions if the polymer concentration is increased or the surface tension is decreased. Rise speed jump discontinuities were not found.     

Birefringence measurements in extensional flow of polyisobutylene around an expanding bubble

Birefringence in liquid polymers offers the possibility of obtaining information about stress in complex flows. In this work, this is done for extensional flows of polyisobutylene in a “breathing bubble” rheometer. In this type of rheometer, a bubble consisting of an incompressible, low-viscosity fluid (usually water) is injected into the sample with a nozzle. Expanding or collapsing the bubble by adding or removing water induces biaxial or uniaxial extension in the surrounding sample. The pressure difference between the bubble and the surroundings can be measured and compared to the predictions of constitutive equations. This measurement only gives one integral value for a complex flow history. In this paper, the birefringence around the bubble is measured in order to learn more about the flow. This is done by comparing pressure and birefringence results to those of standard constitutive equations for a polyisobutylene sample. A good agreement between the pressure and optical measurements and the theory is found with a single value of the stress-optical constant. Received: 25 June 1997 Accepted: 12 November 1997     

Simulation of flow around a thin,flexible obstruction by means of a deforming grid overlapping a fixed grid

This paper presents a numerical method for simulation of coupled flows, in which the fluid interacts with a thin deformable solid, such as flows in cardiovascular valves. The proposed method employs an arbitrary Lagrangian–Eulerian (ALE) method for flow near the solid, embodied in the outflow domain in which the mesh is fixed. The method was tested by modelling a two‐dimensional channel flow with a neo‐Hookean obstacle, an idealization of the coupled flow near a cardiovascular valve. The effects of the Reynolds number and the dimensionless elastic modulus of the material on the wall shear stress, the size of the downstream reverse flows, and the velocity and pressure profiles were investigated. The deformation of the obstacle, the pressure drop across the obstacle, and the size of the top reverse flow increased as the Reynolds number increased. Conversely, increasing the elastic modulus of the obstacle decreased the deformation of the obstacle and the size of the top reverse flows, but did not affect the pressure drop across the obstacle over the range studied. Copyright © 2007 John Wiley & Sons, Ltd.