Single cavitation bubble generation and observation of the bubble collapse flow induced by a pressure wave

This study utilizes a U-shape platform device to generate a single cavitation bubble for a detailed analysis of the flow field characteristics and the cause of the counter jet during the process of bubble collapse caused by sending a pressure wave. A high speed camera is used to record the flow field of the bubble collapse at different distances from a solid boundary. It is found that a Kelvin–Helmholtz vortex is formed when a liquid jet penetrates the bubble surface after the bubble is compressed and deformed. If the bubble center to the solid boundary is within one to three times the bubble’s radius, a stagnation ring will form on the boundary when impinged by the liquid jet. The fluid inside the stagnation ring will be squeezed toward the center of the ring to form a counter jet after the bubble collapses. At the critical position, where the bubble center from the solid boundary is about three times the bubble’s radius, the bubble collapse flow will vary. Depending on the strengths of the pressure waves applied, the collapse can produce a Kelvin–Helmholtz vortex, the Richtmyer–Meshkov instability, or the generation of a counter jet flow. If the bubble surface is in contact with the solid boundary, the liquid jet can only move inside-out without producing the stagnation ring and the counter jet; thus, the bubble collapses along the radial direction. The complex phenomenon of cavitation bubble collapse flows is clearly manifested in this study.     

水中高压脉动气泡水射流形成机理及载荷特性研究

具有脉动特性的气泡(如水下爆炸气泡、螺旋桨空泡和气枪气泡)动力学行为很大程度上取决于其边界条件. 实验已证实,近自由液面气泡在坍塌过程中常常产生背离自由液面的水射流现象,而近刚性边界气泡在坍塌阶段产生朝向壁面的高速水射流,严重威胁水中结构的局部强度. 前人基于 Rayleigh-Plesset 气泡理论和 “Bjerknes” 力来预测气泡射流方向,然而理论方法难以透彻的揭示气泡射流的初生、发展和砰击过程中丰富的力学机理. 本文首先采用水下高压放电技术产生气泡,并通过高速摄影对不同边界条件下气泡的运动特性进行实验研究. 然后,采用边界积分法模拟气泡非球状坍塌过程. 研究表明,边界条件改变了气泡周围的流场压力梯度方向,进而影响气泡射流初生位置;射流在发展阶段,气泡附近流场的局部高压区和射流之间存在“正反馈效应”,从而揭示了气泡射流速度在短时间内即可增加到百米每秒的力学机理. 射流砰击会在流场中造成局部高压区,随着气泡回弹,射流速度和砰击压力逐渐减小. 本文还探讨了无量纲距离参数对气泡运动及射流砰击载荷的影响,旨为近场水下爆炸等相关领域提供参考.      

Level set method for numerical simulation of a cavitation bubble,its growth,collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier–Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving interface is captured by the level set function, and the interface velocity is resolved by “one-side” velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the “shock wave”-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating. The project supported by the National Natural Science Foundation of China (10272032 and 10672043). The English text was polished by Keren Wang.     

The Art, Craft and Science of Modelling Jet Impact in a Collapsing Cavitation Bubble

One of the key characteristics of the asymmetric collapse of a cavitation bubble near a rigid boundary is the development of a high speed liquid jet that penetrates the interior of the bubble, impacting on the other side to yield a toroidal bubble. After the formation of the toroidal bubble, a vigorous splash may occur that can lead to pressures on the boundary an order of magnitude greater than the impact pressures associated with the jet. Qualitative agreement with available experimental data is found although, as the bubble approaches minimum volume, shock waves are also observed which further complicate our full understanding of the mechanisms for damage.     

Experimental and numerical study on the growth and collapse of a bubble in a narrow tube

The growth, expansion and collapse of a bubble in a narrow tube are studied using both experiments and numerical simulations. In experiment, the bubble is generated by an electric spark in a water tank and recorded by a high-speed camera system. In numerical simulation, the evolution of the bubble is solved by adopting axisymmetric boundary integral equation, considering the surface tension effect. The results of experiments and numerical simulations are compared and good agreements are achieved. Both of them show that a counter-jet forms and penetrates the bubble at the end of the collapse stage, before a ring type bubble forms. Under the attraction of the tube wall due to Bjerknes force, a ring jet is generated, pointing towards the tube. On the basis of this, some physical quantities like the pressure on the tube wall and kinetic energy are calculated in a case study. The effects of tube diameters and tube lengths on the bubble’s behaviors are also investigated.     

Spark-generated bubble collapse near or inside a circular aperture and the ensuing vortex ring and droplet formation

This paper aims to study the oscillation of a sparkgenerated submerged bubble located near or inside a circular aperture made in a flat plate using high-speed visualization technique. In the case of a bubble oscillating near an aperture the initial free surface of the water was set at the bottom surface of the plate. The effects of aperture size and bubblefree surface distance on the bubble behavior as well as on the ensuing droplet dynamics are investigated. It was found that the direction of the bubble reentrant jet was towards the aperture or away from it respectively when the normalized aperture size was smaller or greater than a certain critical value. In addition, a toroidal vortex ring was observed to form, which rotated inwards as it moved away from the aperture. It was also found that if the bubble was incepted at a distance sufficiently away from a supercritical size aperture a single droplet could be produced. In the case of a bub- ble initiated in the middle of a circular aperture submerged just beneath the water free surface, the bubble was found to take the shape of an ellipsoid during its expansion. Then a reentrant jet was initiated and pierced the bubble from its top side.     

A calculation of the parameters of the high-speed jet formed in the collapse of a bubble

As is known, the collapse of vapor bubbles in a liquid can cause the intensive destruction of solid boundary surfaces. Experimental and theoretical investigations of bubble collapse have led to the conclusion that the surface of a bubble can deform and a liquid jet directed toward the solid surface can form in the process [1, 2]. In the theoretical reports [3, 4] too low jet velocities were obtained, inadequate to explain the destruction of the surface in a single impact. In [5] it was found as a result of numerical calculations that the formation of jets possessing enormous velocities is possible. It was also found that two fundamentally different schemes of jet formation are possible in the collapse of a bubble near a wall. The transition from one scheme to the other occurs upon a relatively small change in the initial shape of the bubble. In the present report we investigate the case of sufficiently small initial deformations of a bubble when the region occupied by the bubble remains simply connected during the formation of the jet; i.e., the separation of a small bubble from the bubble does not occur. In the case of the second scheme of bubble collapse near a wall the connectedness of the free boundary is disrupted and a small bubble separates off during the formation of the jet.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 94–99, May–June, 1979.     

气泡多周期运动时引起的流场压力与速度

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.      

NUMERICAL ANALYSIS ON THE VELOCITY AND PRESSURE FIELDS INDUCED BYMULTI-OSCILLATIONS OF AN UNDERWATER EXPLOSION BUBBLE

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.     

Shock-induced collapse of a vapor nanobubble near solid boundaries

The collapse of a nano-bubble near a solid wall is addressed here exploiting a phase field model recently used to describe the process in free space. Bubble collapse is triggered by a normal shock wave in the liquid. The dynamics is explored for different bubble wall normal distances and triggering shock intensities. Overall the dynamics is characterized by a sequence of collapses and rebounds of the pure vapor bubble accompanied by the emission of shock waves in the liquid. The shocks are reflected by the wall to impinge back on the re-expanding bubble. The presence of the wall and the impinging shock wave break the symmetry of the system, leading, for sufficiently strong intensity of the incoming shock wave, to the poration of the bubble and the formation of an annular structure and a liquid jet. Intense peaks of pressure and temperatures are found also at the wall, confirming that the strong localized loading combined with the jet impinging the wall is a potential source of substrate damage induced by the cavitation.     

Cavitation bubble dynamics — new tools for an intricate problem

With the help of laser produced bubbles in water and high speed photography and holography sophisticated experiments on cavitation bubble dynamics can be conducted. The observation of a bubble vortex ring after jet formation upon collapse of a spherical bubble in front of a plane solid boundary is reported. The vortex ring may expand and contract several times until it disintegrates into a ring of bubbles by some instability finally taking over. A critical discussion of our qualitative understanding of jet formation is included. In a second part the problem of the acoustic cavitation noise spectrum is discussed. Numerically obtained ‘visible cavitation noise’ plots from a single bubble already resemble those obtained experimentally from acoustic cavitation. A discussion shows that the theory should be extended to self-consistency.     

单液滴内空化气泡的生长及溃灭研究

超空化燃油射流使得喷雾中部分燃油分裂液滴内含有空化气泡;空化气泡的生长及溃灭对液滴的分裂与雾化具有重要影响. 基于VOF 方法首次对超空化条件下燃油液滴内空化气泡的生长及溃灭过程进行了数值模拟. 通过研究发现,单液滴内空化气泡的生长过程可以按控制机理划分为表面张力控制阶段、综合竞争阶段和惯性力控制阶段;在第I 阶段,空泡的生长主要受表面张力的控制作用,惯性力对空泡生长的促进作用及黏性力对空泡生长的抑制作用可以忽略;在第II 阶段,空泡的生长受表面张力、惯性力及黏性力三者的综合作用,空泡的生长速率是促进空泡生长的惯性力和抑制空泡生长的表面张力及黏性力相互竞争、共同作用的结果;在第III 阶段,空泡的生长主要受惯性力的控制作用,抑制空泡生长的表面张力及黏性力的作用基本可以忽略. 单液滴内空化气泡的溃灭过程由多个溃灭阶段和反弹阶段构成,类似于有阻尼弹簧振子的振动过程;根据每个溃灭周期结束时空泡半径随时间的变化历程,可以将空泡的溃灭分为快速溃灭期、缓慢溃灭期以及稳定期;溃灭初期空泡溃灭压力的变化非常剧烈,但空泡溃灭体积的变化则要相对平缓得多;空泡反弹压力随时间的变化与空泡反弹体积随时间的变化基本对应.      

Motion and filling of cavities in a boundless liquid and close to a plane

The motion of bubbles in liquids has been studied in many earlier papers [1–8]. In this paper methods of the projection type are applied to the problem of a cavity in an ideal, incompressible liquid in the absence of vortices. The collapse of a bubble having a finite initial velocity in a boundless liquid is considered; also considered is the collapse of a stationary bubble close to a solid wall. Using the small-parameter method the generation of a jet is examined analytically. A numerical computing method not involving small parameters is developed; it is based on calculating the projection by numerical computation of the corresponding integrals. The method combines economy and simplicity of application with a high accuracy in the region in which the representation of the velocity potential by a series of spherical functions remains effective.     

瞬态空化泡演变过程的数值模拟

采用边界积分方程方法，对无粘流体中三个空化泡以及自由面附近二个空化泡相互作用的演变过程进行了数值模拟。计算中边界用二阶有限元离散，影响系数矩阵非对角线元素用六点高斯数值积分方法计算，对第一类、第二类完全椭圆积分用高次多项式近似，对计算系数矩阵对角线元素中遇到的奇异积分进行了特殊处理。结果表明，在不同的给定参数下，空化泡的溃灭形态各异，柱状射流和环形射流都有可能发生，使空化泡演变成双泡或环形泡。     

Flow of a shear-thinning liquid in a cylindrical container with a rotating end wall

The flow patterns produced by rotating one end wall of a circular cylinder completely filled with a strongly shear-thinning viscoelastic liquid have been investigated using the laser-induced fluorescence flow visualization technique. An intense toroidal vortex is produced in the vicinity of the rotating end wall with outward spiraling flow over the end wall itself. This vortex drives a second countercirculating vortex of low intensity in the region of the stationary end wall. Under some circumstances an axial jet of fluid is observed moving away from the rotating end wall. This jet showed evidence of instability, whereas all flows were otherwise completely steady. The double-vortex structure is different from those recently observed in either a Newtonian or slightly shear-thinning liquid or in the low Reynolds number flow of an elastic liquid. There are, however, similarities with older work for a viscoelastic liquid at relatively high Reynolds numbers. The observations highlight the suitability of the cylinder/rotating end wall configuration as a sensitive test case for computational work.     

Wake structures behind an oscillating bubble rising close to a vertical wall

In the present study, the wake structures behind an oscillating (zigzagging in a plane) air bubble, rising in a close vicinity of a vertical wall are experimentally investigated using a high-speed two-phase particle image velocimetry. While varying the distance between the rising bubble and the wall, the spatial and temporal variations in the spanwise and streamwise vorticity components contained in the wake vortices, in addition to the bubble trajectory, are measured in a tank filled with water. In particular, the Lagrangian streamwise vorticity fields in the bubble wake have been reconstructed and investigated in detail with different conditions. Without the wall, it is confirmed that there exist counter-rotating streamwise vortex tubes in the bubble wake, agreeing with the case of a two-dimensional zigzagging bubble, as reported in the literatures. It is also found that the hairpin vortex chain structures, initially attached to the bubble rear, evolve to detached vortex ring structures as the bubble rises in an oscillating path. While the detailed vortex structures show up quite differently from the reference case depending on the distance to the wall (e.g., actual bubble-wall collision), in general the wake behind the bubble as it moves toward and away from the wall can be summarized as: (i) transition to the detached vortex ring structures is accelerated; (ii) streamwise length of vortex tubes is shortened (evolution is interfered); (iii) counter-rotating vortex tubes approaching the wall tend to slightly bounce off and slide away (being dissipated fast) from each other on the wall; and (iv) boundary-layer like secondary flow structures are induced on the wall due to additional viscous effects. These wall-induced wake modification indicates that more fluid energy is wasted due to the wall interference, rather than being used to force the lateral movement of the bubble, which agrees with the reduced amplitude and wavelength of the oscillating bubble path on the wall. Finally, this explanation has been further confirmed by estimating the vortex-induced lateral forces acting on the bubble for each case.     

近距离下射流冲击平板PIV实验研究

运用时间分辨粒子成像测速系统(time-resolved particle image velocimetry, TR-PIV)对近距离下射流冲击平板时的流场进行了直接测量, 通过对两个正交的平面流场开展测量, 揭示了冲击距离和雷诺数对射流间隙内三维流动特征及涡系结构演化规律的影响. 结果表明: 射流间隙存在三种典型的涡系结构, 分别为双涡环模式、单涡环模式和卷吸模式, 但在大流量湍流状态下, 射流可能会冲破涡环, 形成随机的高速出流, 各流动模式的出现主要与射流流态及壁面约束作用有关. 运用涡量分析对三种典型涡系结构的能量传递和损失特性进行了比较. 结果表明: 近距离冲击时, 射流的能量通过涡环模式向外传递. 在双涡环模式下, 两个涡环的旋向相反, 端面的约束作用使得两个涡环都被严格约束在射流棒端面之内, 且一次涡环强度显著大于二次涡环强度. 最后, 运用本征正交分解方法对射流间隙内的流动模态及其能量分布进行了分析. 单涡和双涡模式前十阶模态分析结果表明: 能量脉动在较低阶时即以配对的模式出现, 这表明一次涡环与二次涡环均具有良好的对称性, 同时在双涡模式中, 一次涡环是占主导作用的大尺度流动结构. 卷吸模式的前三阶模态分析表明: 射流的能量集中在射流上游, 能量随紊动扩散急剧衰减.      

环形喷管喷口气泡演化的实验研究

水下气泡的发展演化及气泡动力学行为是气液两相动力学的基础理论与水下射流应用的重要基础. 环形喷管/喷口形成的气泡及气体射流具有其不同于圆孔实心射流的特殊表现与规律机制,随着同心筒破水发射等特殊应用的出现,环形喷口气体射流/泡流的基础现象观测和机制分析成为迫切的需求. 基于环形喷管的设计和水下射流条件的分析,设计建立了一套环形喷管水箱实验系统,对水下环形喷管喷口气泡发展演化过程进行了初步的实验研究. 为观测研究气体通过环形喷管气泡生长发展过程,在较低压力、较低流速下,采用高速摄影仪记录气泡生长及发展演化过程. 结合对气泡发展演化过程的图像处理与分析,研究分析了环形喷口气泡形成区制、气泡生长过程形态发展特点、以及气泡形成时间及气泡体积变化特点. 研究表明:在本实验气体流量范围内(50.8~237.3 dm3/min),环形喷口气泡发展演化过程呈现较为明显的三周期区制,前泡尾流影响是环形气泡呈三周期区制的主要原因;不同周期内的气泡形成时间具有较稳定规律,并受到流量影响;气泡生长过程中有较为明显的下沉、回升特征;气泡表面张力、液体惯性与流动的共同作用,造成了典型的气泡顶部坍塌现象.      

Particle removal from a surface by a bounded vortex flow

A bounded vortex flow consists of an axisymmetric vortex that is confined top and bottom between two plates (the “confinement plate” and “impingement plate”, respectively) and surrounded laterally by a swirling annular slot jet. The bottom of the vortex terminates on the boundary layer along the impingement plate and the top of the vortex is drawn into a suction port positioned at the center of the confinement plate. The circumferential flow within the annular jet is important for supplying circulation to the central wall-normal vortex. This flow field is proposed as a method for mitigation of dust build-up on a surface, where the vortex–jet combination supplements the more traditional vacuum port by enhancing the surface shear stress and related particle transport rate. The paper reports on a computational study of the velocity field and particle transport by a bounded vortex flow. Fluid flow computations are performed using a finite-volume approach for an incompressible fluid and particle transport is simulated using a discrete-element method. Computations are performed for different values of two dimensionless parameters – the ratio of the plate separation distance and the average radial location of the jet inlet (the dimensionless confinement height) and the ratio of flow rate withdrawn at the suction outlet and that injected by the jet (the flow rate ratio). For small values of the flow rate ratio, the impinging jet streamlines pass down to the boundary layer along the bottom surface and then travel up the vortex core. By contrast, for large values of flow rate ratio, the annular jet is quickly entrained into the suction outlet and no wall-normal vortex is formed. Particles are observed to roll along the impingement surface in a direction determined by the fluid shear stress lines. Particles roll outward when they lie beyond a separatrix curve of the surface shear stress lines, where particles within this separatrix curve roll inward, piling up at the center of the flow field. A toroidal vortex ring forms for the small confinement height case with flow rate ratio equal to unity, which yields double separatrix curves in the shear stress lines. The inward rolling particles intermittently lift up due to collision forces and burst away from the impingement surface, eventually to become entrained into the flow out the suction port or resettling back onto the impingement surface.