垂直壁面附近上升单气泡的弹跳动力学研究

采用高速摄影技术结合阴影法,对静止水中垂直壁面附近上升单气泡运动进行实验研究,对比气泡尺度及气泡喷嘴与壁面之间的初始无量纲距离 ($S^{\ast}$)对气泡上升运动特性的影响,分析气泡与壁面碰撞前后,壁面效应与气泡动力学机制及能量变化规律.结果表明,对于雷诺数$Re \approx 580 \sim 1100$,无量纲距离$S^{\ast } <2 \sim3$时,气泡与壁面碰撞且气泡轨迹由无约束条件下的三维螺旋转变成二维之字形周期运动;当$S^{\ast } >2 \sim3$时,壁面效应减弱,有壁面约束的气泡运动与无约束气泡运动特性趋于一致.气泡与壁面碰撞前后,壁面效应导致横向速度峰值下降为原峰值的70%,垂直速度下降50%;气泡与壁面碰撞前,通过气泡中心与壁面距离($x/R$)和修正的斯托克斯数相关式可预测垂直速度的变化规律.上升气泡与壁面碰撞过程中,气泡表面变形能量单向传输给气泡横向动能,使得可变形气泡能够保持相对恒定的弹跳运动.提出了气泡在与壁面反复弹跳时的平均阻力系数的预测模型,能够很好地描述实验数据反映出的对雷诺数${Re}$、韦伯数${We}$和奥特沃斯数${Eo}$等各无量纲参数的标度规律.      

Numerical Simulation of Three-Dimensional Bubble Oscillations by a Generalized Vortex Method

A numerical method is implemented for simulating the simultaneous three-dimensional volume and shape oscillations of a compressible vapor or gas bubble suspended in an inviscid ambient fluid in the presence of interfacial tension. The flow generated by the bubble expansion, contraction, and deformation is represented by an interfacial distribution of potential dipoles supplemented by a point source situated inside the bubble, accounting for changes in the bubble volume. The mathematical formulation is completed by setting the strength of the point source proportional to the integral of the density of the double-layer potential over the interface. The motion of marker points distributed over the interface is computed using a boundary-element implementation of Baker's generalized vortex method in which the normal component of the interfacial velocity is computed in terms of tangential derivatives of the vector potential associated with the dipoles, whereas the tangential component of the interfacial velocity is computed in terms of the surface gradient of the scalar harmonic potential. The density of the double-layer distribution is computed by solving an integral equation of the second kind using an iterative method, while the evolution of the interfacial distribution of the harmonic potential is computed using Bernoulli's equation for irrotational flow. The onset of interfacial irregularities due to numerical instabilities is prevented by truncating the Fourier–Legendre spectrum of the interfacial distribution of the harmonic potential. With smoothing implemented, the numerical method is capable of describing simultaneous volume and shape oscillations for an indefinite period of time. Received 7 September 2001 and accepted 30 April 2002 Published online 30 October 2002 RID="*" ID="*" This research was supported by a grant provided by NASA. Communicated by J.R. Blake     

Marangoni效应下气泡上升过程的FTM模拟

采用界面跟踪法FTM(Front-Tracking Method),并结合CSF(continuum surface force)模型,研究了在垂直方向上温度分布不均匀的对称流场中由Marangoni效应引起的气泡上升运动问题。模拟了在不同的M a数和Pr数下单气泡及同轴双气泡的运动。研究结果表明,在不同的M a数下气泡的运动差异较大,M a数越大,气泡运动至稳态时的速度越大,且气泡运动的最大速度值与M a数呈正相关关系;增大Pr数所造成的粘度增大或热扩散率减小将削弱气泡的迁移运动;Marangoni对流中双气泡的局部运动证实了温度梯度和气泡运动速度紧密相关。     

On the interaction between two consecutive elongated bubbles in a vertical pipe

The motion of elongated air bubbles in a vertical pipe filled with water is studied quantitatively using video imaging of the flow and subsequent digital image processing of the recorded sequence of images. Experiments are carried out to determine the influence of the separation distance between two consecutive bubbles (liquid slug length) upon the behavior of the trailing bubble in vertical slug flow. The details of the trailing bubble acceleration and merging process are observed and the instantaneous parameters of the trailing bubble, such as its shape, velocity, acceleration, etc., are measured as a function of the separation distance. The leading bubble is found to be unaffected by the trailing elongated bubble.     

竖直流道宽度对气泡运动行为影响的数值模拟

用数值方法模拟了竖直通道宽度对气泡在液体中的非定常运动、变形以及传热特性的影响。在这个模拟中,界面跟踪采用了VOF方法,并采用PL IC进行界面重构。主流场计算采用有限容积方法将控制方程离散,其中扩散项采用中心差分格式,对流项采用一阶迎风格式。用成熟的S IM PLE算法求解N-S方程的速度与压力的耦合问题。引入CSF模型处理运动界面的表面张力。利用所编制的程序计算了竖直流道中的单个气泡的形状、运动特性以及气泡内外流场与传热特性,并对竖直通道宽度在不同情况下,对气泡的形状、运动特性以及传热特性进行了进一步的研究。得到了一系列有价值的结果,并与实验结果比较。表明数值模拟结果与实验结果吻合的较好。     

The local measurement of the size and velocity of bubbles rising in liquids

An experimental study of the rise of small air bubbles (0.1 to 1.0 mm in diameter) in a quiescent pool is described. Local measurements of rise velocity were obtained as a function of height above the source nozzle using a laser-Doppler method. In addition, the bubble diameter was determined simultaneously from the same optical signals. Data are presented for various bubble diameters and spacings in bubble columns for both distilled water and a dilute polymer solution.It was found that for distilled water the rise velocity near the nozzle reaches the maximum observed in other studies before decelerating to its terminal velocity due to surfactant accumulation at its interface. The maximum rise velocity in dilute polymer was much lower for the same bubble diameter and reached its terminal velocity much faster. The results are shown to be in closer agreement with predictions for a solid sphere in this case.     

Bubble and liquid flow characteristics in a cylindrical bath during swirl motion of bubbling jet

Gas injection into a cylindrical bath through a centric bottom nozzle causes a swirl motion like rotary sloshing. Conditions indicating the initiation and cessation of the swirl motion have been made clear by many researchers. So far, the effect of the swirl motion on transport phenomena in the bath is not clear yet. The present study was made to clarify the bubble characteristics (void fraction, bubble frequency) and liquid flow characteristics (mean velocity, turbulence intensity, Reynolds shear stress) during swirl motion of bubbling jet. These two characteristics were investigated using an electro-resistivity probe and a two-dimensional LDV, respectively.     

Associated disturbances around a vortex ring in a stratified fluid

A motion of a vortex ring in a stratified fluid is accompanied by associated disturbances which, in the schlieren visualization in the field of a horizontal density gradient, have the shape of a symmetric four-petal configuration. The criterion of the existence of the disturbances is the Froude number Fr based on the motion velocity and the vertical vortex size. On the range Fr > 1, the disturbances are stable with respect to the variation of themotion regime and the distortion of the vortex shape. For Fr < 1 the disturbances disappear. Computer processing of the schlieren photographs showed that the experimental spatial dependences of the disturbance amplitude are close to the functions describing the distribution of the vertical velocity component in the inviscid flow past a sphere.     

A finite element method for free surface flows of incompressible fluids in three dimensions. Part II. Dynamic wetting lines

To date, few researchers have solved three‐dimensional free surface problems with dynamic wetting lines. This paper extends the free surface finite element method (FEM) described in a companion paper [Cairncross RA, Schunk PR, Baer TA, Sackinger PA, Rao RR. A finite element method for free surface flows of incompressible fluid in three dimensions. Part I. Boundary fitted mesh motion. International Journal for Numerical Methods in Fluids 2000; 33 : 375–403] to handle dynamic wetting. A generalization of the technique used in two‐dimensional modeling to circumvent double‐valued velocities at the wetting line, the so‐called kinematic paradox, is presented for a wetting line in three dimensions. This approach requires the fluid velocity normal to the contact line to be zero, the fluid velocity tangent to the contact line to be equal to the tangential component of web velocity, and the fluid velocity into the web to be zero. In addition, slip is allowed in a narrow strip along the substrate surface near the dynamic contact line. For realistic wetting line motion, a contact angle that varies with wetting speed is required because contact lines in three dimensions typically advance or recede at different rates depending upon location and/or have both advancing and receding portions. The theory is applied to capillary rise of static fluid in a corner, the initial motion of a Newtonian droplet down an inclined plane, and extrusion of a Newtonian fluid from a nozzle onto a moving substrate. The extrusion results are compared with experimental visualization. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

Characteristics of normal and tangential forces acting on a single lug during translational motion in sandy soil

Lugs (i.e., grousers) are routinely attached to the surfaces of wheels/tracks of mobile robots to enhance their ability to traverse loose sandy terrain. Much previous work has focused on how lug shape, e.g., height, affects performance; however, the goal of this study is to experimentally confirm the effects of lug motion on lug–soil forces. We measured normal and tangential forces acting on a single lug as functions of inclination angle, moving direction angle, sinkage length, horizontal displacement, and traveling speed. The experimental results were mathematically fitted by using least square method to facilitate quantitative analyses on effects of changes in these motion parameters. Moreover, we compared the measured tangential forces to values calculated from a conventional tangential force model to evaluate the effects of the lug-tip surface, which is generally ignored in existing terramechanics models. The conclusions from this study would be useful for estimating the traveling performance of locomotive mechanisms equipped with lugs, modeling interaction mechanics between lugged wheels and soil, etc.     

喷管对水下爆轰气泡形态和压力特征影响的实验研究

针对具有不同类型喷管的爆轰管在水下爆轰中形成的燃气射流问题,搭建了爆轰实验平台,研究了单次爆轰过程中尾部喷管对水下气泡形态与压力特征的影响。采用数字粒子图像测速技术对高速摄影机拍摄得到的气泡脉动图片进行流场可视化分析,得到各喷管工况下的气泡速度场。为了确认爆轰管内是否形成稳定爆轰波,并观察爆轰波在气液两相界面上的透反射特性,爆轰管尾部安装有2个动态压力传感器,同时在距离喷管一定距离处设置一个水下爆炸传感器,以监测水中传播的压力波。结果表明：扩张喷管工况下的气泡脉动过程与直喷管工况基本一致,但扩张喷管提高了燃气射流速度,气泡膨胀体积更大;因为燃气射流的持续性,收敛喷管工况下的气泡脉动过程具有明显差异,气泡膨胀体积较小,但气泡二次脉动时长相较于一次脉动时长衰减更小;扩张喷管提高了气泡脉动强度,扩张喷管工况下的气泡脉动压力与透射冲击波压力远大于直喷管工况下的气泡脉动压力与透射冲击波压力;收敛喷管工况下的气泡脉动压力与透射冲击波压力都较小,但收敛喷管燃气射流的持续性减缓了气泡脉动压力的衰减速度。相比于直喷管,扩张喷管工况下的气泡脉动时间、气泡脉动压力与透射冲击波压力都更大。收敛喷管工况下的气泡...     

Bubbly jets in stagnant water

Air–water bubbly jets are studied experimentally in a relatively large water tank with a gas volume fraction, C_o, of up to 80% and nozzle Reynolds number, Re, ranging from 3500 to 17,700. Measurements of bubble properties and mean axial water velocity are obtained and two groups of experiments are identified, one with relatively uniform bubble sizes and another with large and irregular bubbles. For the first group, dimensionless relationships are obtained to describe bubble properties and mean liquid flow structure as functions of C_o and Re. Measurements of bubble slip velocity and estimates of the drag coefficient are also provided and compared to those for isolated bubbles from the literature. The study confirms the importance of bubble interactions to the dynamics of bubbly flows. Bubble breakup processes are also investigated for bubbly jets. It was found that a nozzle Reynolds number larger than 8000 is needed to cause breakup of larger bubbles into smaller bubbles and to produce a more uniform bubble size distribution. Moreover, the Weber number based on the mean water velocity appears to be a better criteria than the Weber number based on the bubble slip velocity to describe the onset of bubble breakup away from the nozzle, which occurs at a Weber number larger than 25.     

Computer simulation of the trajectories of large water jets

A three-dimensional computer simulation of the motion of a water jet is described which includes the effects of wind from any direction. The simulation is useful in the design of fire-fighting systems, particularly those used in offshore situations. The equations of motion are presented in vector form and the problem of the fluid dynamic drag variation is discussed. Semi-empirical approximations for the drag components along and across the jet are presented which involve four unknown constants. These are reduced to three by using previous data on the efficiency of vertical jets. To fix the remaining constants, information was available from a series of large jets tests carried out to prove an offshore fire fighting system. In these tests different nozzle shapes were tried and, using the best of these shapes, a large number of trajectories were measured photographically. These were used to fix the simulation drag constants and good agreement is shown between measurements and predictions. The simulation enables the effects of flow-rate, pressure, nozzle size, elevation and wind strength to be evaluated in the system design     

Flow characteristics of gas-liquid-particle mixing in a gas-stirred ladle system with throughflow

An experimental study is performed on air-liquid-particle mixing, resulting from an air-particle mixture injected into a liquid flowing through a slender ladle. Flow visualization combined with image processing is employed to investigate the bubble and particle behavior at the nozzle outlet. Effort is directed to particle discrimination in both the liquid and the bubbles to determine particle distribution, which affects the mixing performance of gas bubbles, solid particles and liquid. A real-time movement of bubble and particle behavior can be visualized by means of image processing with the use of a slow-motion video recording. It is disclosed that the particles injected through the nozzle may stick on the inner surface of the gas bubble, break through the bubble surface, or mingle with the gas stream to form a two-phase jet, depending on the particle-to-gas mass flow rate ratio. It is observed that when a solid-gas two-phase jet penetrates deeper in the horizontal direction, the particles and bubbles rise along the vertical sidewall and simultaneously spread in the transverse direction, thus promoting a better liquid-particle mixing. The application of the slow-motion video recording results in quantitative evaluations of both the penetration depth of particles or of gas-particles from the injection nozzle and the velocity distribution along the sidewall.List of symbols B Width of water vessel, m - B _n Nozzle location on bottom surface of water vessel, m - d _o Diameter of a gas-particle injection nozzle, m - H Height of water vessel, m - H _n Nozzle location on vertical surface of water vessel, m - L Penetration length of particles or of particles and gas from the nozzle, m - Q _g Volumetric flow rate of gas, m³/s - Q _l Volumetric flow rate of water, m³/s - Q _s Volumetric flow rate of particle, m³/s - Re _g Gas Reynolds number based on inner diameter of the air-particle injection nozzle - t Time, sec. - W Thickness of water vessel, m - x Transverse coordinate, m - y Longitudinal coordinate, m - Mass flow rate ratio of particles to gas Visiting scholar on leave from the Mechanical Engineering Department, Kagoshima University, Kagoshima, JapanThe work reported was supported by the National Science Foundation under the Grant No. CTS-8921584     

Turbulent flow past a bubble and an ellipsoid using shadow-image and PIV techniques

An experimental investigation on flow around an oscillating bubble and solid ellipsoid with a flat bottom was conducted. A single air bubble (equivalent diameter D_e=9.12 mm) was attached to a small disk (1 mm) at the end of a needle and suspended across a vertical square channel (100 mm) by wire wherein water flowed downward at a constant flowrate. The solid ellipsoid (D_e9.1 mm) was suspended across the square channel in the same manner. The equivalent diameter-based Reynolds and Eotvos number range, 1950<Re<2250 and 11<Eo<11.5, placed the bubble in the ‘wobbly’ regime while the flow in its wake was turbulent. A constant flowrate and one bubble size was used such that flow in the wake was turbulent. Velocity measurements of the flow field around the bubble or solid were made using a one CCD camera Digital Particle Image Velocimetry (DPIV) system enhanced by Laser Induced Fluorescence (LIF). The shape of the bubble or solid was simultaneously recorded along with the velocity using a second CCD camera and an Infrared Shadow Technique (IST). In this way both the flow-field and the boundary of the bubble (solid) were measured. The velocity vector plots of flow around and in the wake of a bubble/solid, supplemented by profiles and contours of the average and root-mean-square velocities, vorticity, Reynolds stress and turbulent kinetic energy, revealed differences in the wake flow structure behind a bubble and solid. One of the significant differences was in the inherent, oscillatory motion of the bubble which not only produced vorticity in the near-wake, but as a result of apparent vorticity stretching distributed the turbulent kinetic energy associated with this flow more uniformly on its wake, in contrast to the solid.     

Analytical and Experimental Study of Pressure Dynamics in a Pulsed Water Jet Device

Pulsed high-velocity water jets are of interest for breaking rocks and other materials. This paper describes a straightforward way of generating single water pulse with a hammer impacting a piston that rests on top of a chamber filled with water. This impacting action pressurises the water, expelling it at high velocity through a nozzle. A theoretical investigation is outlined aimed at gaining a better understanding of this system for generating water pulses. A computational model is developed to simulate the pressure dynamics in the chamber based on continuity and momentum equations for a compressible viscous flow. This model is used to optimise the relative sizes of the hammer and piston as well as the height of the water column to produce the highest velocity water pulse. The model was validated by building an experimental apparatus. In these experiments maximum pressures of about 200 MPa were measured inside the chamber over a time period of about 560???s. This produced a water pulse with maximum velocity of 600 m/s. Experiments were conducted with nozzle diameters between about 1 mm and 4 mm to study the effect of discharge volume on the pressure history. The results illustrate that although the peak attainable pressure decreases with an increase in nozzle diameter, the duration of the elevated pressure remains similar for all nozzles.     

利用格子Boltzmann方法模拟单个气泡在复杂流道内的运动特性

利用格子Boltzmann方法模拟了单个气泡在具有三个半圆形喉部的复杂流道内的上升过程.通过分析气泡运动过程中的形态、运动轨迹及运动速度的变化,研究复杂流道对气泡运动特性的影响.研究结果表明,在上升过程中,由于壁面的影响,气泡的形状发生严重的变形,运动轨迹也发生相应的偏转.通过实验结果的对比,验证了模拟结果的正确性.结果表明格子Boltzmann方法可以用于模拟具有复杂边界的两相流问题.     

An initially spherical bubble rising near a vertical wall

The dynamics of a single rising bubble in the vicinity of a vertical wall is explored via three-dimensional numerical simulations. A finite volume method is used to solve the incompressible Navier–Stokes equations. The gas–liquid interface is reconstructed by volume-of-fraction (VOF) method. The trajectory, velocity, shape and vorticity of the bubble are analyzed in detail. The numerical results show that the presence of the wall imposes a repulsion on the bubble and that the bubble migrates away from the wall upon release. The onset of bubble path oscillation is found to occur earlier than for a freely rise counterpart and also at a lower Galilei number. Interestingly, we find that the vertical wall serves as a destabilizing factor in the wall-normal direction but a stabilizing factor in the spanwise direction. The increase of bubble inertia is discovered to enhance the influence of the wall. Furthermore, the bubble oscillations seem insensitive to the variation of the initial bubble-wall distance.     

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.