曲面固壁附近空化泡溃灭动力学的 流体体积数值研究*

流体体积法(VOF)可以便捷、高效地实现对多相流界面的捕捉和追踪。本文基于VOF方法,对单个空化泡在曲面固壁附近的运动进行了数值模拟,从实验对比、压力场、速度场、温度场演化、溃灭时间、射流速度、固壁温度等方面分析了空化泡溃灭过程的热动力学影响。结果表明,数值模拟得到的空化泡形态演化与实验观测到的现象一致,随着位置参数、泡内外压差及曲面固壁尺寸的改变,空化泡热动力学行为也将发生变化,受到流体运动及射流冲击的影响,溃灭瞬间产生的高温高压使得曲面固壁温度升高。本文研究的曲面固壁附近空化泡溃灭效应,揭示了空化泡与曲面固壁间的相互作用规律,对学术研究及工程应用都具有重要意义。     

Dynamic behavior of a single bubble between the free surface and rigid wall

The objective of this paper is to apply high-speed photography and schlieren method to investigate the bubble dynamics between the free surface and a rigid wall. The temporal evolution of the bubble shape and the free surface motion are recorded by two synchronous high-speed cameras. Experiments are carried out for a single bubble generated at various normalized stand-off distances from bubble center to the free surface and to the rigid wall. The results show that (1) three distinctive patterns are identified with the morphology of the bubble and free surface, namely single toroidal bubble without spike (STB), single toroidal bubble with a spike (STBS) and double toroidal bubbles with a spike (DTBS). (2) The dynamic characteristics of the bubble at collapse and rebound stage vary evidently at different patterns, including the bubble shape variations and free surface motion. In detail, the schlieren images show the formation and propagation of shock waves, which explains the radiative process of bubble collapse energy. (3) Qualitative comparisons are carried out for the bubble and free surface at the same pattern. And quantitative analyses are conducted for the jet velocity, bubble collapse position, bubble collapse time and spike height, etc. for different values of bubble-rigid wall distance.     

Evolution of the shape of a bubble deformed at the zero time in a viscous fluid

The boundary-value problem of the shape of a charged bubble in a viscous incompressible insulating fluid as a function of time is solved in an approximation linear in the initial deformation amplitude. Both centrosymmetric pulsations and vibrations at constant volume are considered. It is shown that the time evolution of the shape of the bubble, as well as of the velocity and pressure fields of the fluid in its neighborhood, can be represented by finite sums over the numbers of initially excited modes, which involve two terms. The first term is a sum over the roots of the dispersion relation; the second, an improper integral. In the low-and high-viscosity limits, the relevant analytical expressions simplify, i.e., become free of integrals. It turns out that the damping constant of bubble surface vibrations varies with viscosity nonmonotonically in the case of radial pulsations and this dependence is different in the different limits. The frequencies of radial pulsations and surface vibrations vary with viscosity monotonically.     

Liquid jet directionality and droplet behavior during emulsification of two liquids due to acoustic cavitation

The present study numerically investigates liquid-jet characteristics of acoustic cavitation during emulsification in water/gallium/air and water/silicone oil/air systems. It is found that a high-speed liquid jet is generated when acoustic cavitation occurs near a minute droplet of one liquid in another. The velocity of liquid jet significantly depends on the ultrasonic pressure monotonically increasing as the pressure amplitude increases. Also, the initial distance between cavitation bubble and liquid droplet affects the jet velocity significantly. The results revealed that the velocity takes maximum values when the initial distance between the droplet and cavitation bubble is moderate. Surprisingly, the liquid jet direction was found to depend on the droplet properties. Specifically, the direction of liquid jet is toward the droplet in the case of water/gallium/air system, and vice versa the jet is directed from the droplet in the case of water/silicone oil/air system. The jet directionality can be explained by location of the high-pressure spot generated during the bubble contraction.     

Numerical study on the splitting of a vapor bubble in the ultrasonic assisted EDM process with the curved tool and workpiece

Electrical discharge machining (EDM) is a powerful and modern method of machining. In the EDM process, a vapor bubble is generated between the tool and the workpiece in the dielectric liquid due to an electrical discharge. In this process dynamic behavior of the vapor bubble affects machining process. Vibration of the tool surface affects bubble behavior and consequently affects material removal rate (MRR). In this paper, dynamic behavior of the vapor bubble in an ultrasonic assisted EDM process after the appearance of the necking phenomenon is investigated. It is noteworthy that necking phenomenon occurs when the bubble takes the shape of an hour-glass. After the appearance of the necking phenomenon, the vapor bubble splits into two parts and two liquid jets are developed on the boundaries of the upper and lower parts of the vapor bubble. The liquid jet developed on the upper part of the bubble impinges to the tool and the liquid jet developed on the lower part of the bubble impinges to the workpiece. These liquid jets cause evacuation of debris from the gap between the tool and the workpiece and also cause erosion of the workpiece and the tool. Curved tool and workpiece affect the shape and the velocity of the liquid jets during splitting of the vapor bubble. In this paper dynamics of the vapor bubble after its splitting near the curved tool and workpiece is investigated in three cases. In the first case surfaces of the tool and the workpiece are flat, in the second case surfaces of the tool and the workpiece are convex and in the third case surfaces of the tool and workpiece are concave. Numerical results show that in the third case, the velocity of liquid jets which are developed on the boundaries of the upper and lower parts of the vapor bubble after its splitting have the highest magnitude and their shape are broader than the other cases.     

超声场下刚性界面附近溃灭空化气泡的速度分析

为了揭示刚性界面附近气泡空化参数与微射流的相互关系, 从两气泡控制方程出发, 利用镜像原理, 建立了考虑刚性壁面作用的空化泡动力学模型. 数值对比了刚性界面与自由界面下气泡的运动特性, 并分析了气泡初始半径、气泡到固壁面的距离、声压幅值和超声频率对气泡溃灭的影响. 在此基础上, 建立了气泡溃灭速度和微射流的相互关系. 结果表明: 刚性界面对气泡振动主要起到抑制作用; 气泡溃灭的剧烈程度随气泡初始半径和超声频率的增加而降低, 随着气泡到固壁面距离的增加而增加; 声压幅值存在最优值, 固壁面附近的气泡在该最优值下气泡溃灭最为剧烈; 通过研究气泡溃灭速度和微射流的关系发现, 调节气泡溃灭速度可以达到间接控制微射流的目的.     

Interaction of microbubbles with high intensity pulsed ultrasound

High intensity pulsed ultrasound, interacting with microbubble contrast agents, is potentially useful for drug delivery, cancer treatment, and tissue ablation, among other applications. To establish the fundamental understanding on the interaction of a microbubble (in an infinite volume of water) with an ultrasound pressure field, a numerical study is performed using the boundary element method. The response of the bubble, in terms of its shape at different times, the maximum bubble radius obtained, the oscillation time, the jet velocity, and its translational movement, is studied. The effect of ultrasound intensity and initial bubble size is examined as well. One important outcome is the determination of the conditions under which a clear jet will be formed in a microbubble in its interaction with a specific sound wave. The high speed jet is crucial for the aforementioned intended applications.     

High-speed imaging of ultrasonic emulsification using a water-gallium system

Aiming at elucidating ultrasonic emulsification mechanisms, the interaction between a single or multiple acoustic cavitation bubbles and gallium droplet interface was investigated using an high-speed imaging technique. To our best knowledge, the moment of emulsification and formation of fine droplets during ultrasound irradiation were observed for the first time. It was found that the detachment of fine gallium droplets occurs from the water-gallium interface during collapse of big cavitation bubbles. The results suggest that the maximum size of cavitation bubble before collapsing is of prime importance for emulsification phenomena. Previous numerical simulation revealed that the collapse of big cavitation bubble is followed by generation of high-velocity liquid jet directed toward the water-gallium interface. Such a jet is assumed to be the prime cause of liquid emulsification. The distance between cavitation bubbles and water-gallium interface was found to slightly affect the emulsification onset. The droplet fragmentation conditions are also discussed in terms of the balance between (1) interfacial and kinetic energies and (2) dynamic and Laplace pressure during droplet formation.     

Nonlinear corrections to the oscillation mode frequencies of an ideal liquid charged jet

An analytical asymptotic expression is derived for the time evolution of an ideal incompressible conducting liquid jet with a uniformly charged surface that moves with a constant velocity along the symmetry axis of an undisturbed cylindrical surface. The expression is obtained in the third order of smallness in jet oscillation amplitude under the conditions when the initial deformation of the equilibrium jet surface is specified by one (axisymmetric or nonaxisymmetric) mode. Analytical expressions are also found for the positions of internal nonlinear degenerate three-and four-mode resonances, which are typical of nonlinear corrections to the equilibrium shape of the jet, liquid velocity field potential in the jet, and electrostatic field near the jet, and also for a nonlinear frequency correction. It is established that the nonlinear oscillations of the jet take place about a surface depending on the type of initial (generally asymmetric) deformation rather than about the equilibrium cylindrical shape.     

Comparative visualized investigation of impact-driven high-speed liquid jets injected in submerged water and in ambient air

This paper is a comparative study on the characteristics of high-speed liquid jets injected in surrounding water and air using shadowgraph technique. One of the main objectives is to investigate the effects of liquid’s physical properties, used to generate the high-speed liquid jets, on jet generation’s characteristics. Moreover, comparative investigations on effects of those liquid jets after injected in water and air are reported. The high-speed liquid jets were generated by the impact of a projectile launched by a horizontal single-stage power gun. The impact-driven high-speed liquid jets were visualized by shadowgraph technique and images were recorded by a high-speed digital video camera. The process of impact-driven high-speed liquid jet injection in air and water, oblique shock waves, jet-induced shock waves, shock waves propagation, the bubble behavior, bubble collapse-induced rebound shock waves and bubble cloud regeneration were clearly observed. It was found that different properties of liquid (surface tension and kinematic viscosity) affect the jet maximum velocity and shape of the jet. Bubble behaviors were only found for the jet injected in water. From the shadowgraph images, it is found that the maximum average jet velocity, expansion and contraction velocities of bubble in axial direction increase when the value of the multiplied result of surface tension by kinematic viscosity increases. Therefore, surface tension and kinematic viscosity are the significant physical properties that affect characteristics of high-speed liquid jets.     

空心旋转液体射流初始阶段特性的实验研究

以空心旋转射流的理论模型为基础［１］,通过实验对旋流喷嘴射流特性做进一步深入的研究,特别是实际射流形状的变化规律及喷嘴初始参数的变化规律．结果表明,旋流喷嘴的初始参数和空心射流内外表面压力差是决定射流形状的主要因素。     

自由表面磁流体射流不稳定性分析(II)

用微分解析方法分析了自由表面磁流体射流在两个非均匀、相互垂直的横向磁场中的稳定性行为。获得了在给定射流初始速度分布和磁场分布条件下射流的速度空间分布表达式;分析了射流扭转及向外溅射的位置。     

Numerical study on the desorption processes of oil droplets inside oil-contaminated sand under cavitation micro-jets

The removal of the adsorbed oil droplet is critical to deoiling treatment of oil-bearing solid waste. Ultrasonic cavitation is regarded as an extremely useful method to assist the oil droplets desorption in the deoiling treatment. In this paper, the effects of cavitation micro-jets on the oil droplets desorption were studied. The adsorbed states of oil droplets in the oil-contaminated sand were investigated using a microscope. Three representative absorbed states of the oil droplets can be summarized as: (1) the individual oil droplet adsorbed on the particle surface (2) the clustered oil droplets adsorbed on the particle surface; (3) the oil droplet adsorbed in a gap between particles. The micro-jet generation during the bubble collapse near a rigid wall under different acoustic pressure amplitudes at an ultrasonic frequency of 20 kHz was investigated numerically. The desorption processes of the oil droplets at the three representative absorbed states under micro-jets were also simulated subsequently. The results showed that the acoustic pressure has a great influence on the velocity of micro-jet, and the initial diameter of cavitation bubbles is significant for the cross-sectional area of micro-jets. The wall jet caused by a micro-jet impacting on the solid wall is the most important factor for the removal of the absorbed oil droplets. The oil droplet is broken by the jet impinging, and then it breaks away from the solid wall due to the shear force generated by the wall jet. In addition to a higher sound pressure, the cavitation bubble at a larger initial diameter is more important for the desorption of the clustered oil droplets. Conversely, the micro-jet generated by the cavitation bubble at a smaller initial diameter (0.1 mm) is more appropriate for the desorption of the oil droplet in a narrow or sharp-angled gap.     

Thermodynamic effect of single bubble near a rigid wall

The objective of this paper is to numerically investigate the thermodynamic effect during bubble collapse near a rigid boundary. A compressible fluid model is introduced to accurately capture the transient process of bubble shapes and temperature, as well as corresponding pressure, and velocity. The accuracy of the numerical model is verified by the experimental data of bubble shapes, and Keller-Kolodner equation as well as its thermodynamic equation. The results show that a bubble near the rigid boundary presents high-speed jet in collapse stage and counter jet in rebound stage, respectively. In the collapse stage, the bubble margin will shrink rapidly and do the positive work on the compressible vapor inside the bubble, then a significant amount of heat will be generated, and finally the generation of high-speed jet drives the low-temperature liquid outside the bubble to occupy the position of high-temperature vapor inside the bubble. In the rebound stage, the counter jet moving away from the rigid boundary takes part of heat away from the sub-bubble, which avoids the external work of the expansion of the sub-bubble and the temperature reduction caused by the dissipation effect of the vortex structure. In addition, the initial standoff has a significant effect on the thermodynamics of bubble oscillation. The temperature keeps increasing with the increase of the initial standoff in the collapse stage, while it shows a downward trend with the increase of the initial standoff in the rebound stage. That’s because the high-speed jet and counter jet of bubble gradually disappear when the initial standoff increases, which is the important reason for the opposite evolution trend of temperature in collapse and rebound stage.     

Nonlinear oscillations of a charged jet of a conducting liquid under multimode initial deformation of its surface

The subject of consideration is a uniformly charged jet of an ideal incompressible conducting liquid moving with a constant velocity along the symmetry axis of an undisturbed cylindrical surface. An evolutionary expression for the jet shape is derived accurate to the second order of smallness in oscillation amplitude for the case when the initial deformation of the equilibrium surface is a superposition of a finite number of both axisymmetric and nonaxisymmetric modes. The flow velocity field in the jet and the electric field distribution near it are determined. The positions of internal nonlinear secondary combined three-mode resonances are found, which are typical of nonlinear corrections to the analytical expressions for the jet shape, flow velocity field potentials, and electrostatic field in the vicinity of the jet.     

圆锥边界附近激光空泡溃灭行为的研究

为了研究刚性圆锥边界锥角对激光空泡溃灭行为的影响,文章建立了虚拟平面边界模型, 同时采用阴影摄影术、光偏转法以及数值计算的手段对边界附近空泡溃灭过程进行了研究. 结果表明边界的锥角对空泡的形状、溃灭时间以及液体射流形成均有明显影响. 空泡形状偏离球形的程度和溃灭时间均随锥角的增大而增大,且增大锥角度可以促使射流的形成. 空泡溃灭时间的实验值同理论值具有较高的一致性,验证了虚拟平面边界假设及无量纲距离修正的有效性.     

影响空心旋转液体射流初始阶段运动的无量纲参数的分析

本文以描述空心旋转时流初始阶段运动规律的理论模型[1]为依据,重点讨论了影响射流运动的无量纲参数：韦伯数,液膜内外压力差,液膜初始内外环半径比,初始旋流度;初始径向速度等。从而进一步加深了对旋流喷嘴所产生的空心旋转射流初始阶段运动及发展规律的认识。     

对流梯形微通道内气液两相流动的数值研究

采用VOF方法,对梯形微通道内不可压缩气液两相流动进行了数值模拟研究,详细分析了气泡形成过程,以及当量直径、截面形状、液体表面张力和粘度等对气泡液柱形成过程和长度的影响,拟合出微通道气泡液柱长度计算公式。结果表明:气泡液柱的长度受表观气速和表观液速的影响较大;表面张力对气泡尺寸的影响较小,当液体粘度增加为水粘度的10倍时,形成的气泡形状不规则。增大表面张力,形成气泡的时间增加;增大粘度,形成气泡的时间减小。     

过热液体中蒸汽泡上升过程的格子Boltzmann三维数值模拟

应用格子Boltzmann相变模型,在三维空间研究蒸汽泡在过热液体中生长、上升和变形等动力学行为.为研究传热传质对蒸汽泡运动的影响,对比模拟相同条件下气泡在等温环境中上升的物理过程.结果表明：蒸汽泡在过热液体中上升发生的变形程度较小,意味着相变对蒸汽泡的影响和表面张力一样使汽泡保持初始的形状.蒸汽泡在过热液体中的上升速度较小,说明随着汽泡生长拖拽力的影响比浮力大.蒸汽泡生长率在初始阶段达到最大值,随后会趋于一个恒定的值.随着汽泡体积增大和上升速度的增加,其对流场的扰动也越来越剧烈.蒸汽泡生长和上升引起的对流运动对温度场的演化造成很大的影响.     

Experimental investigation on the effects of the standoff distance and the initial radius on the dynamics of a single bubble near a rigid wall in an ultrasonic field

Bubble behaviors near a boundary in an ultrasonic field are the fundamental forms of acoustic cavitation and of substantial importance in various applications, such as industry cleaning, chemical engineering and food processing. The effects of two important factors that strongly affect the dynamics of a single acoustic cavitation bubble, namely, the initial bubble radius and the standoff distance, were investigated in this work. The temporal evolution of the bubble was recorded using high speed microphotography. Meanwhile, the time of bubble collapse and the characteristics of the liquid jets were analyzed. The results demonstrate that the intensity of the acoustic cavitation, which is characterized by the time of bubble collapse and the liquid jet speed, reaches the optimum level under suitable values of the initial bubble radius and the normalized standoff distance. As the initial bubble radius and the normalized standoff distance increase or decrease from the optimal values, the time of the bubble collapse increases, and the first liquid jet’s speed decreases substantially, whereas the speeds of the second and third liquid jets exhibit no substantial changes. These results on bubble dynamics in an ultrasonic field are important for identifying or correcting the mechanisms of acoustic cavitation and for facilitating its optimization and application.