Enhanced drainage and coarsening in aqueous foams

Experiments are presented elucidating how the evolution of foam microstructure by gas diffusion from high to low pressure bubbles can significantly speed up the rate of gravitational drainage, and vice versa. This includes detailed data on the liquid-fraction dependence of the coarsening rate, and on the liquid-fraction and the bubble-size profiles across a sample. These results can be described by a "coarsening equation" for the increase of bubble growth rate for drier foams. Spatial variation of the average bubble size and liquid fraction can also affect the growth and drainage rates.     

A method for estimating time-dependent acoustic cross-sections of bubbles and bubble clouds prior to the steady state

Models for the acoustic cross-sections of gas bubbles undergoing steady-state pulsation in liquid have existed for some time. This article presents a theoretical scheme for estimating the cross-sections of single bubbles, and bubble clouds, from the start of insonation onward. In this period the presence of transients can significantly alter the cross-section from the steady-state value. The model combines numerical solutions of the Herring-Keller model with appropriate damping values to calculate the extinction cross-section of a bubble as a function of time in response to a continuous harmonic sound field (it is also shown how the model can be adapted to estimate the time-dependent scatter cross-section). The model is then extended to determine the extinction cross-section area of multiple bubbles of varying population distributions assuming no bubble-bubble interactions. The results have shown that the time taken to reach steady state is dependent on the closeness of the bubble to resonance, and on the driving pressure amplitude. In the response of the population as a whole, the time to reach steady state tends to decrease with increasing values of the driving pressure amplitude; and with the increasing values of the ratio of the numbers of bubbles having radii much larger than resonance to the number of resonant bubbles. The implications of these findings for the use of acoustic pulses are explored.     

II. Recirculation flow induced by a bubble stream rising in a viscous liquid

The recirculation flow induced by the rising motion of a bubble stream in a viscous fluid within an open-top rectangular enclosure is studied. The three-dimensional volume averaged conservation equations are solved by a control-volume method using a hybrid finite differencing scheme to describe the liquid phase hydrodynamics. The momentum exhange between the bubbles and the liquid phase is modeled with a source term equals to the volumetric buoyancy force acting on the gas in the bubble stream. The volumetric buoyancy force accounts for in line interactions between bubbles through the average gas volume fraction in the gas liquid column which depends on the size and the rising velocity of bubbles. The fluid flow within an open-top rectangular enclosure is further investigated by particle image velocimetry for a bubble stream rising in a water-glycerol solution. The measured fluid velocities in a vertical plane are compared with the predictions of the numerical model over a wide range of fluid viscosity (43 mPa s-800 mPa s) and gas flow rates. Finally, the recirculation flows resulting from the interaction of two neighbouring vertical bubble streams are studied. Received: 23 July 1997 / Revised: 19 December 1997 / Accepted: 11 May 1998     

Ultrasonic scattering cross sections of shell-encapsulated gas bubbles immersed in a viscoelastic liquid: first and second harmonics

The oscillations of gas bubbles, without shell, immersed in viscoelastic liquids and driven by an acoustic wave have been the subject of several investigations. They demonstrate that the viscosity coefficient and the spring constant of the liquid have significant influence on the scattering cross section of the gas bubble. For shell-encapsulated gas bubbles, the investigations have been concentrated to bubbles immersed in a pure viscous liquid. This present work computes the ultrasonic scattering cross section, first and second harmonics, of shell-encapsulated gas bubbles immersed in a viscoelastic liquid. The theoretical model of the bubble oscillation is based on the generalized Rayleigh-Plesset equation of motion of a spherical cavity immersed in a viscoelastic liquid represented by a three-parameter linear Oldroyd model. The scattering cross section is computed for Albunex type of bubble (shell thickness=15 nm, shell shear viscosity=1.77 Pas, shell modulus of rigidity=88.8 MPa) irradiated by a 3.5 MHz ultrasonic pressure wave with an amplitude of 30 kPa. The results demonstrate that encapsulated bubbles respond independently of the surrounding liquid being pure viscous or viscoelastic as long as the surrounding liquid shear viscosity is as low as 10(-3) Pas. Nevertheless, for higher shear viscosities, the bubble responds differently if the surrounding liquid is pure viscous or viscoelastic. In general, the scattering cross sections of first and second harmonics are larger for the viscoelastic liquid.     

大小肥皂泡连通后的讨论与计算

从弯曲液面附加压强、气体状态方程以及"球缺"的有关公式等出发,对大小肥皂泡连通后的变化作了全面的阐释,导出了最后大泡半径的近似公式,还介绍了演示实验,拍得了3个泡的照片,测出了它们的半径,验证了这个近似公式.并对产生误差的原因作了分析.     

Quasi-one-dimensional foam drainage

Foam drainage is considered in a froth flotation cell. Air flow through the foam is described by a simple two-dimensional deceleration flow, modelling the foam spilling over a weir. Foam microstructure is given in terms of the number of channels (Plateau borders) per unit area, which scales as the inverse square of bubble size. The Plateau border number density decreases with height in the foam, and also decreases horizontally as the weir is approached. Foam drainage equations, applicable in the dry foam limit, are described. These can be used to determine the average cross-sectional area of a Plateau border, denoted A, as a function of position in the foam. Quasi-one-dimensional solutions are available in which A only varies vertically, in spite of the two-dimensional nature of the air flow and Plateau border number density fields. For such situations the liquid drainage relative to the air flow is purely vertical. The parametric behaviour of the system is investigated with respect to a number of dimensionless parameters: K (the strength of capillary suction relative to gravity), α (the deceleration of the air flow), and n and h (respectively, the horizontal and vertical variations of the Plateau border number density). The parameter K is small, implying the existence of boundary layer solutions: capillary suction is negligible except in thin layers near the bottom boundary. The boundary layer thickness (when converted back to dimensional variables) is independent of the height of the foam. The deceleration parameter α affects the Plateau border area on the top boundary: weaker decelerations give larger Plateau border areas at the surface. For weak decelerations, there is rapid convergence of the boundary layer solutions at the bottom onto ones with negligible capillary suction higher up. For strong decelerations, two branches of solutions for A are possible in the K = 0 limit: one is smooth, and the other has a distinct kink. The full system, with small but non-zero capillary suction, lies relatively close to the kinked solution branch, but convergence from the lower boundary layer onto this branch is distinctly slow. Variations in the Plateau border number density (non-zero n and h) increase individual Plateau border areas relative to the case of uniformly sized bubbles. For strong decelerations and negligible capillarity, solutions closely follow the kinked solution branch if bubble sizes are only slightly non-uniform. As the extent of non-uniformity increases, the Plateau border area reaches a maximum corresponding to no net upward velocity of foam liquid. In the case of vertical variation of number density, liquid content profiles and Plateau border area profiles cease to be simply proportional to one another. Plateau border areas match at the top of the foam independent of h, implying a considerable difference in liquid content for foams which exhibit different number density profiles. Received 3 July 2001     

Dry microfoams: formation and flow in a confined channel

We present an experimental investigation of the agglomeration of microbubbles into a 2D microfoam and its flow in a rectangular microchannel. Using a flow-focusing method, we produce the foam in situ on a microfluidic chip for a large range of liquid fractions, down to a few percent in liquid. We can monitor the transition from separated bubbles to the desired microfoam, in which bubbles are closely packed and separated by thin films. We find that bubble formation frequency is limited by the liquid flow rate, whatever the gas pressure. The formation frequency creates a modulation of the foam flow, rapidly damped along the channel. The average foam flow rate depends non-linearly on the applied gas pressure, displaying a threshold pressure due to capillarity. Strong discontinuities in the flow rate appear when the number of bubbles in the channel width changes, reflecting the discrete nature of the foam topology. We also produce an ultra flat foam, reducing the channel height from 250 μm to 8 μm, resulting in a height to diameter ratio of 0.02; we notice a marked change in bubble shape during the flow.     

I. Formation and rise of a bubble stream in a viscous liquid

The continuous emission of gas bubbles from a single ejection orifice immersed in a viscous fluid is considered. We first present a semi empirical model of spherical bubble growth under constant flow conditions to predict the bubble volume at the detachment stage. In a second part, we propose a physical model to describe the rise velocity of in-line interacting bubbles and we derive an expression for the net viscous force acting on the surrounding fluid. Experimental results for air/water-glycerol systems are presented for a wide range of fluid viscosity and compared with theoretical predictions. An imagery technique was used to determine the bubble size and rise velocity. The effects of fluid viscosity, gas flow rate, orifice diameter and liquid depth on the bubble stream dynamic were analyzed. We have further studied the effect of large scale recirculation flow and the influence of a neighbouring bubble stream on the bubble growth and rising velocity. Received: 23 July 1997 / Revised: 16 December 1997 / Accepted: 11 May 1998     

Symmetry-breaking transitions and dissociation of two-dimensional Plateau borders

We have experimentally studied the dissociation/coalescence of internal Plateau borders (PBs) in simple monolayer bubble clusters, as a result of changing the liquid fraction. At large liquid content, the clusters consist of n bubbles of the same size, symmetrically placed around an internal n-sided PB (n-PB). On decreasing the liquid fraction we observed symmetry-breaking transitions in the 4- and 5-bubble clusters (but not in the 3-bubble cluster), followed by dissociation of the PBs. We used the Surface Evolver to determine the various equilibrium configurations of the corresponding two-dimensional wet clusters and their surface energies. The sequence of 4-bubble cluster configurations observed on varying the liquid fraction correlates qualitatively with that predicted on the basis of Surface Evolver calculations. The same is not true of the 5-bubble cluster.     

球状泡群内气泡的耦合振动

振动气泡形成辐射场影响其他气泡的运动, 故多气泡体系中气泡处于耦合振动状态. 本文在气泡群振动模型的基础上, 考虑气泡间耦合振动的影响, 得到了均匀球状泡群内振动气泡的动力学方程, 以此为基础分析了气泡的非线性声响应特征. 气泡间的耦合振动增加了系统对每个气泡的约束, 降低了气泡的自然共振频率, 增强了气泡的非线性声响应. 随着气泡数密度的增加, 振动气泡受到的抑制增强; 增加液体静压力同样可抑制泡群内气泡的振动, 且存在静压力敏感区(1–2 atm, 1 atm=1.01325×10⁵ Pa); 驱动声波对气泡振动影响很大, 随着声波频率的增加, 能够形成空化影响的气泡尺度范围变窄. 在同样的声条件、泡群尺寸以及气泡内外环境下, 初始半径小于5 μm 的气泡具有较强的声响应. 气泡耦合振动会削弱单个气泡的空化影响, 但可延长多气泡系统空化泡崩溃发生的时间间隔和增大作用范围, 整体空化效应增强.     

Calculation of shock wave propagation in water containing reactive gas bubbles

The entry of a shock wave from air into water containing reactive gas (stoichiometric acetylene–oxygen mixture) bubbles uniformly distributed over the volume of the liquid has been numerically investigated using equations describing two-phase compressible viscous reactive flow. It has been demonstrated that a steady-state supersonic self-sustaining reaction front with rapid and complete fuel burnout in the leading shock wave can propagate in this bubbly medium. This reaction front can be treated as a detonation-like front or “bubble detonation.” The calculated and measured velocities of the bubble detonation wave have been compared at initial gas volume fraction of 2 to 6%. The observed and calculated data are in satisfactory qualitative and quantitative agreement. The structure of the bubble detonation wave has been numerically studied. In this wave, the gas volume fraction behind the leading front is approximately 3–4 times higher than in the pressure wave that propagates in water with air bubbles when the other initial conditions are the same. The bubble detonation wave can form after the penetration of the shock wave to a small depth (~300 mm) into the column of the bubbly medium. The model suggested here can be used to find optimum conditions for maximizing the efficiency of momentum transfer from the pressure wave to the bubbly medium in promising hydrojet pulse detonation engines.     

Observations of the collapses and rebounds of millimeter-sized lithotripsy bubbles

Bubbles excited by lithotripter shock waves undergo a prolonged growth followed by an inertial collapse and rebounds. In addition to the relevance for clinical lithotripsy treatments, such bubbles can be used to study the mechanics of inertial collapses. In particular, both phase change and diffusion among vapor and noncondensable gas molecules inside the bubble are known to alter the collapse dynamics of individual bubbles. Accordingly, the role of heat and mass transport during inertial collapses is explored by experimentally observing the collapses and rebounds of lithotripsy bubbles for water temperatures ranging from 20 to 60 °C and dissolved gas concentrations from 10 to 85% of saturation. Bubble responses were characterized through high-speed photography and acoustic measurements that identified the timing of individual bubble collapses. Maximum bubble diameters before and after collapse were estimated and the corresponding ratio of volumes was used to estimate the fraction of energy retained by the bubble through collapse. The rebounds demonstrated statistically significant dependencies on both dissolved gas concentration and temperature. In many observations, liquid jets indicating asymmetric bubble collapses were visible. Bubble rebounds were sensitive to these asymmetries primarily for water conditions corresponding to the most dissipative collapses.     

The dynamics of cavitation bubbles in a sealed vessel

A model of cavitation bubbles is derived in liquid confined in an elastic sealed vessel driven by ultrasound. In this model, an assumption that the pressure acting on the sealed vessel due to bubble pulsations is proportional to total volume change of bubbles is made. Numerical simulations are carried out for a single bubble and for bubbles. The results show that the pulsation of a single bubble can be suppressed to a large extent in sealed vessel, and that of two matched bubbles with same ambient radius can be further suppressed. However, when two mismatched bubbles have the same ambient radii, an interesting breathing phenomenon takes place, where one bubble pulsates inversely with the other one. Due to this breathing phenomenon the suppression effect becomes weak, so the maximum radii of two mismatched bubbles can be larger than that of a single bubble or that of two matched bubbles in sealed vessel. Besides that, for two mismatched bubbles with different ambient radii, the small one in sealed vessel under some certain parameters can pulsate as strong as or even stronger than that of a single bubble in an open vessel.     

Dynamics of a bubble cluster in an acoustic field

A mathematical model describing the dynamics of clustered gas bubbles under the effect of an acoustic field is presented. The proposed model is used as the basis for an analytical study of small bubble oscillations in monodisperse and polydisperse clusters and for a numerical study of nonlinear bubble oscillations under high-amplitude external pressures. The following effects are found to occur in a polydisperse cluster: a synchronization of the collapse phases of bubbles with different radii and a collapse intensification for bubbles of one size in the presence of bubbles of another size. These effects are explained by the interaction between bubbles of different radii in the cluster.     

Bubble Velocity and Size Measurement with a Four‐Point Optical Fiber Probe

The possibility to measure the velocity and size of individual bubbles in a high‐void fraction bubbly flow is investigated by using a four‐point optical fiber probe. The air bubbles have an initial spherical equivalent diameter ranging from 4 to 10 mm and the void fraction is up to 0.3. Firstly, single bubble experiments show that intrusiveness effects, i.e. bubble deformations due to the probe, are negligible provided that the bubble approaches the probe at the axis of the central fiber. A selection criterion is utilized for multiple bubble experiments. A good compromise can be found between the required accuracy, the duration of the measurements and the number of validated bubbles required for reliable statistical averaging. In an air‐water high‐void fraction vertical bubbly pipe flow, the void fraction obtained with the instrument is found to be in good agreement with both local single‐fiber probe measurements, and with the volume average void fraction obtained from pressure gradient measurements. The area average volumetric gas flow rate, based on the bubble velocity and void fraction as measured with the four‐point probe, agree with the measured gas flow rate. Also, the liquid velocity is measured by means of a laser‐Doppler anemometer, to investigate the slip velocity. The results show that reliable and interesting measurements can be obtained by using a four‐point optical fiber probe in high void fraction flows.     

Interaction between phases in the liquid–gas system

This work analyzes the equilibrium between a liquid and a gas over this liquid separated by an interface. Various gas forms exist inside the liquid: dissolved gas molecules attached to solvent molecules, free gas molecules, and gaseous bubbles. Thermodynamic equilibrium is maintained between two phases; the first phase is the liquid containing dissolved and free molecules, and the second phase is the gas over the liquid and bubbles inside it. Kinetics of gas transition between the internal and external gas proceeds through bubbles and includes the processes of bubbles floating up and bubble growth as a result of association due to the Smoluchowski mechanism. Evolution of a gas in the liquid is considered using the example of oxygen in water, and numerical parameters of this system are given. In the regime under consideration for an oxygen–water system, transport of oxygen into the surrounding air proceeds through micron-size bubbles with lifetimes of hours. This regime is realized if the total number of oxygen molecules in water is small compared with the numbers of solvated and free molecules in the liquid.     

Interaction of a bubble and a bubble cluster in an ultrasonic field

Using an appropriate approximation, we have formulated the interacting equation of multi-bubble motion for a system of a single bubble and a spherical bubble cluster. The behavior of the bubbles is observed in coupled and uncoupled states. The oscillation of bubbles inside the cluster is in a coupled state. The numerical simulation demonstrates that the secondary Bjerknes force can be influenced by the number density, initial radius, distance, driving frequency, and amplitude of ultrasound. However, if a bubble approaches a bubble cluster of the same initial radii, coupled oscillation would be induced and a repulsive force is evoked, which may be the reason why the bubble cluster can exist steadily. With the increment of the number density of the bubble cluster, a secondary Bjerknes force acting on the bubbles inside the cluster decreases due to the strong suppression of the coupled bubbles. It is shown that there may be an optimal number density for a bubble cluster which can generate an optimal cavitation effect in liquid for a stable driving ultrasound.     

Nonlinear oscillation and acoustic scattering of bubbles

The scattered acoustic pressure and scattered cross section of bubbles is studied using the scattered theory of bubbles. The nonlinear oscillations of bubbles and the scattering acoustic fields of a spherical bubble cluster are numerically simulated based on the bubble dynamic and fluid dynamic. The influences of the interaction between bubbles on scattering acoustic field of bubbles are researched. The results of numerical simulation show that the oscillation phases of bubbles are delayed to a certain extent at different positions in the bubble cluster, but the radii of bubbles during oscillation do not differ too much at different positions. Furthermore, directivity of the acoustic scattering of bubbles is obvious. The scattered acoustic pressures of bubbles are different at the different positions inside and outside of the bubble cluster. The scattering acoustic fields of a spherical bubble cluster depend on the driving pressure amplitude, driving frequency, the equilibrium radii of bubbles, bubble number and the radius of the spherical bubble cluster. These theoretical predictions provide a further understanding of physics behind ultrasonic technique and should be useful for guiding ultrasonic application.     

气泡在超声场中绕圈运动的高速摄影及其图像分析

借助高速摄影和图像分析技术对首次发现的附壁气泡的绕圈现象进行了实验研究,重点研究游移气泡的运动轨迹、附壁气泡的布阵过程、气泡的来源以及气泡的振动细节.研究发现游移绕圈气泡的运动轨迹呈现出不稳定、不规则、不光滑的特点.阵列气泡源于游移气泡,而游移气泡变成阵列气泡的方式主要是通过合并增大体积,从而减小所受的Bjerknes力,降低活性的方式实现的.游移气泡源于ALF(acoustic lichtenberg figure)空化云中大量空泡的合并,使以径向振动为主的空泡逐渐过渡到以表面波动为主的气泡.阵列气泡在Bjerknes力的作用下呈现出规则的表面波动,而体积更小受力更大的游移空泡的表面完全失稳,呈现极不规则的形貌,并对附近阵列气泡的表面波动产生影响.阵列气泡呈现出十分规则的排布,相邻阵列气泡之间的振动相位是相反的,表现为相互排斥.     

Rapid Plateau border size variations expected in three simple experiments on 2D liquid foams<Superscript>⋆</Superscript>

Up to a global scaling, the geometry of foams squeezed between two solid plates (2D GG foams) essentially depends on two independent parameters: the liquid volume fraction and the degree of squeezing (bubble thickness to diameter ratio). We describe it in two main asymptotic regimes: fully dry floor tiles, where the Plateau border radius is smaller than the distance between the solid plates, and dry pancakes, where it is larger. We predict a rapid variation of the Plateau border radius in one part of the pancake regime, namely when the Plateau border radius is larger than the inter-plate distance but smaller than the geometric mean of that distance and the bubble perimeter. This rapid variation is not related to any topological change in the foam: in all the regimes we consider, the bubbles remain in mutual lateral contact through films located at mid-height between both plates. We provide asymptotic predictions in different types of experiments on such 2D GG foams: when foam is being progressively dried or wetted, when it is being squeezed further or stretched, when it coarsens through film breakage or through inter-bubble gas diffusion. Our analysis is restricted to configurations close to equilibrium, as we do not include stresses resulting from bulk viscous flow or from non-homogeneous surfactant concentrations. We also assume that the inter-plate distance is sufficiently small for gravity to be negligible. The present work does not provide a method for measuring small Plateau border radii experimentally, but it indicates that large (and easily observable) Plateau borders should appear or disappear rather suddenly in some types of experiments with small inter-plate gaps. It also gives expected orders of magnitude that should be helpful for designing experiments on 2D GG foams.