A numerical and experimental study of a collapsing bubble-induced droplet ejector

This paper aims to study a novel drop-on-demand droplet generation mechanism in which the oscillation and deformation of a non-equilibrium bubble in close proximity to a free surface induce an axisymmetric liquid spike on the free surface. The evolution of the liquid spike and its deformation due to the effect of surface tension force lead to the formation of a droplet. The free surface can be accorded by either a circular hole on a horizontal flat plate or by the top opening/nozzle of a vertical cylinder. A high-speed camera capable of obtaining images at a frame rate of 15,000 fps is utilized to observe the droplet formation process. Numerical simulations corresponding to the experiments are performed using the boundary integral spatial solution coupled with the time integration, i.e., a mixed Eulerian–Lagrangian method. In the experiments the bubble is generated using a very low voltage (only 55 V) in contrast to the relatively much higher voltages usually employed in reported works. This is very attractive from a safety viewpoint and accords great simplification of the setup. A comparison is made between the numerical and experimental results. A reasonable agreement has been found. The influences of the main design parameters, namely, the bubble-free surface distance and the dimension of the hole/nozzle on the bubble dynamics and on the droplet formation process are discussed and the conditions of the bubble dynamics under which a satellite-free droplet can be generated are sought. Furthermore, the effects of different geometries, namely, the horizontal flat plate and the vertical cylinder on the bubble dynamics and on the droplet features are examined. One important feature of the proposed actuation mechanism is the capability of producing droplets much smaller than the nozzle size. The possible applications of this mechanism are those where the accurate direction of the ejected droplet is of great importance such as inkjet printing.      

Experimental investigation of the interaction between a pulsating bubble and a rigid cylinder

In this paper, the behavior of a bubble near a rigid cylinder is studied experimentally as the positions of bubble induction change, and several cylinders with different diameters are used in the experiment. The main results are as follows. The behavior of a bubble near a rigid cylinder is distinct from that near a rigid plate. When the cylinders are laid in deep water, there will occur three kinds of typical bubble shapes as the distance between bubble and cylinder increases. And the bubble shapes are different as the diameter of cylinder varies. When the cylinders are laid near a free surface, the behaviors of bubble near cylinders with different diameters are similar. For a certain distance between bubble and free surface, as the distance between bubble and cylinder increases, "double jet", "inclined jet" and "downward jet" will take place successively.     

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

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

Numerical study on the oscillation of a transient bubble near a confined free surface for droplet generation

In the present work, the oscillation of a spark-created bubble near a confined water–air interface and the ensuing droplet generation and ejection are studied numerically using the boundary element method. The interface is accorded by the top opening of either one of the following symmetrical configurations, which are distinguished by the value of angle between their vertical symmetry axis and lateral wall (i.e., θ): (i) a centrally perforated horizontal flat plate (θ = 90°) and (ii) vertically placed cylinder (θ = 0°), nozzle (θ > 0°) and diffuser (θ < 0°). Furthermore, the influences of the effective parameters such as the strength parameter (i.e., the intensity of local energy input), the bubble-free surface distance (standoff distance) and the nozzle size on the bubble dynamics and droplet formation and ejection processes are investigated. It was found that the moment at which the bubble attained its maximum volume was advanced as θ increased. In addition, by decreasing θ the attraction of the bubble toward the free surface during its expansion phase and its migration from the free surface during its contraction phase became stronger. Furthermore, for the nozzle case, by increasing θ, the volume of the droplet was increased. It was also found that by increasing the strength parameter, the volume of the droplet increased and its pinch-off happened earlier. Finally, as the standoff distance was increased, the volume of the droplet increased and its pinch-off was delayed.     

近场水下爆炸气泡脉动及水射流的实验与数值模拟研究

海上作战时,近场水下爆炸形成的水射流能造成水面舰船结构的严重局部毁伤。为了研究近场爆炸时舰船底部水射流的形成机理及规律,开展了TNT当量2.5 g的炸药在固支方板底部不同爆距下起爆的水下爆炸实验。结果表明,气泡坍塌形成水射流的过程随着爆距的增加由吸附式向非吸附式转化。接着,基于ABAQUS软件采用CEL方法开展了系列数值模拟,结果表明:爆距在0.821～0.867倍最大气泡半径时,存在吸附式射流向非吸附式射流转化的临界点;固支方板加快了气泡坍塌的进程,炸药与钢板间的距离越小则射流形成的时间越早;射流形成过程中最大速度和射流击中钢板时速度均随着爆距的增大先增大后减小,并在临界点附近达到最大值,射流速度最大可达621 m/s,射流击中钢板时速度最大可达269 m/s。最后,给出了射流开始形成时间、射流最大速度、射流最大速度出现时间、射流击中钢板速度和射流击中钢板时间与距离参数的函数关系式。     

不同环境下气泡脉动特性实验研究

水下爆炸气泡对舰船结构造成严重的毁伤, 海上实船爆炸实验是考核水下爆炸气泡对舰船毁伤威力最直接、最有效的方法, 然而真实情况的舰船水下爆炸实验难以进行.为此, 在已有研究成果的基础上, 设计实验电路, 利用电容在相对较低的电压下放电打火产生的电火花气泡来模拟水下爆炸气泡, 设计多组实验工况来模拟气泡在不同环境下的运动特性. 研究气泡在不同环境下的脉动特性、射流特征, 揭示一些特殊的实验现象, 总结环境对气泡运动的影响规律.      

壁面处气泡在静止流场和高速水流中溃灭过程的计算仿真

通过数值仿真计算,模拟近壁面以及附壁面气泡在静止流场和高速水流中的溃灭过程,研究气蚀作用机理.结果表明:气泡与壁面的距离和水流的速度影响其溃灭时间;附壁面气泡在高速水流中完全溃灭的时间最短,而在静止流场中最长,远离壁面将增加气泡的不稳定性;当气泡距离壁面一定距离溃灭时,射流不能直接作用于壁面,壁面承受冲击波的最大压力远小于气泡溃灭中心的压力;当气泡溃灭中心在壁面时,射流直接作用于壁面产生微小而严重的点破坏,而冲击波则使材料产生交变应力,造成环形破坏;当气泡在高速水流中溃灭时将产生逆流斜向射流,这可能是水力机械过流部件产生鱼鳞坑和波纹状破坏的主要原因.     

水中爆炸气泡与水底边界相互作用的水射流现象

应用开尔文冲量原理,分析了水中爆炸产生的气泡与水底边界相互作用产生水射流的主要影响因 素,开展了气泡与水平放置的钢板相互作用水射流现象的实验研究工作。通过改变炸药到钢板的距离,得到 不同爆距条件下,气泡与钢板相互作用形成水射流的动态过程图像和气泡2次脉动压力。理论分析和实验结 果均表明:随着爆距的减小,Bjerknes力迅速增加,气泡与水平放置钢板相互作用产生水射流的临界条件为爆 距小于气泡最大半径。水射流具有非常明显的方向性,对目标毁伤作用明显。     

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.     

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

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

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

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

三维气泡与自由表面相互作用的直接数值模拟

采用VOF中的PLIC界面重构方法数值模拟了三维气泡在液流中上升并与自由表面相互作用的运动.分别考察了不同初始高度,有无来流及有无再生气泡对气泡上升高度、上升速度、压力及与自由表面相互作用等的影响.结果表明:气泡初始位置越低,顶端上升的高度越大,自由面隆起的范围更广.越靠近自由表面,底部射流横向发展越窄,而向上的压力梯度,气泡上升速度,底部射流上升高度越大,反之则反;但如果底部射流均在接近自由表面以前已横向发展充分,则差别不大.气泡外形、上升高度、破裂时间以及上升速度与来流无关.产生再生气泡后,原生气泡与再生气泡相吸,相互加速对方的上升;自由表面抬升的高度增幅较大,范围拓宽,上升速度也大大增加,且再生气泡越多,自由表面隆起的范围越大.     

狭长直管约束条件水下电爆炸所产生的气泡运动和界面射流

以高速摄影为主要手段,揭示直管中爆炸诱导气泡和射流的典型演变过程,并测试爆炸深度和爆炸能量对该现象的影响。研究发现直管中爆炸诱导的表面射流分为光滑和粗糙的两段,这区别于自由表面射流的形态;爆炸气泡的发展经历一个先膨胀再坍缩的过程,其中封闭坍缩以气泡顶部形成内向射流为特征。表面射流速度主要来自爆炸早期短时间内气泡膨胀赋予水体的冲量,且整体上与起爆能量成正相关,而与爆炸深度成反相关;用准一维的简化模型能够很好地描述它们之间的依赖关系,计算结果不仅在趋势上与实验结果一致,数值上也能很好吻合。     

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

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

Convective and diffusive surfactant transfer in multiphase liquid systems

Laser interferometry was used to investigate diffusive and convective mass transfer in a multicomponent fluid mixture with a liquid–liquid or liquid–gas interface. For this purpose, an immobile gas bubble or insoluble fluid droplet, having the shape of a short cylinder with a free lateral surface, was inserted into a thin liquid layer. In the case of non-uniform distribution of the dissolved surfactant component, the Marangoni convection near the drop/bubble was initiated by the surface tension inhomogeneities, depending on the surfactant concentration. The applied experimental techniques allowed us to study the structure and evolution of the convective flows and concentration fields in a liquid layer, which due to its small thickness were nearly two-dimensional. Making use of both the vertical and horizontal orientation of the liquid layer, we investigated the mass transfer process at different levels of the interaction between gravity and capillary forces. During the experiments, we detected new solutocapillary phenomena, which were found to be caused by oscillatory regimes of solutal convection occurring around air bubbles and chlorobenzene drops in heterogeneous aqueous solutions of alcohol with a vertical surfactant concentration gradient. The role of the oscillatory instability in the processes of drop saturation by the surfactant from its water solution and an inverse process of surfactant extraction from the drop into the surrounding homogeneous fluid (water) was determined. A reasonable explanation for the driving mechanisms of the discovered effects has been proposed.     

两空泡运动特性研究

本文应用边界方法研究了两个相邻空泡的运动特性,得到了空泡的演化规律,以及空泡溃灭时的射流速度与溃灭时间的变化趋势,对于两个空泡之间的距离和半径比的影响进行了讨论。计算结果表明：不同大小的空泡在一起时则小泡会先溃灭,且人泡的存在时间与两泡的半径比成正比;大泡对小泡来说其作用相当于－固壁面,小泡会形成－指向大泡的溃灭射流。相同大小的空泡在一起溃灭时,会同时形成指向中间的射 流,与单空泡在固壁面附近的溃     

Vortex shedding induced by a solitary wave propagating over a submerged vertical plate

Experimental study was conducted on the vortex shedding process induced by the interaction between a solitary wave and a submerged vertical plate. Particle image velocimetry (PIV) was used for quantitative velocity measurement while a particle tracing technique was used for qualitative flow visualization. Vortices are generated at the tip of each side of the plate. The largest vortices at each side of the plate eventually grow to the size of the water depth. Although the fluid motion under the solitary wave is only translatory, vortices are shed in both the upstream and downstream directions due to the interaction of the generated vortices as well as the vortices with the plate and the bottom. The process can be divided into four phases: the formation of a separated shear layer, the generation and shedding of vortices, the formation of a vertical jet, and the impingement of the jet onto the free surface. Similarity velocity profiles were found both in the separated shear layer and in the vertical jet.     

近自由面的多个水下爆炸气泡相互作用研究

将气泡运动阶段周围的流场假设为无黏、无旋、不可压缩的理想流体,运用边界积分法模拟流场中气泡的运动,并开发了三维计算程序,计算值与实验值吻合较好.用该方法模拟了近自由面多气泡之间的相互作用,包括同相气泡和异相气泡. 通过计算发现,气泡的周期随两气泡中心的距离减小而增大,这是由于多气泡之间存在抑制作用,特别是对异相气泡,这种抑制作用更加明显,称之为多气泡之间的抑制效应. 无论有、无自由面存在,多气泡之间均存在抑制效应,由于抑制效应导致同相与异相气泡相互耦合作用的动态特性存在巨大的差异,这些现象可为将来研究多个同时或延时产生的水下爆炸气泡的威力提供参考.      

Characteristics of impact-driven high-speed liquid jets in water

This paper describes a preliminary investigation of the characteristics of high-speed water jets injected into water from an orifice. The high-speed jets were generated by the impact of a projectile launched by a horizontal single-stage powder gun and submerged in a water test chamber. The ensuing impact-driven high-speed water jets in the water were visualized by the shadowgraph technique, and the images were recorded by a high-speed digital video camera. The processes following such jet injection into water, the jet-induced shock waves, shock wave propagation, the bubble behavior, bubble collapse-induced rebound shock waves and bubble cloud re-generation were observed. Peak over-pressures of about 24 and 35 GPa measured by a Polyvinylidence difluoride (PVDF) piezoelectric film pressure sensor were generated by the jet impingement and the bubble impingement, respectively. The peak over-pressure was found to decrease exponentially as the stand-off distance between the PVDF pressure sensor and the nozzle exit increases.     

水下爆炸气泡射流载荷计算的新方法研究

由于存在强冲击、多相界面复杂运动等强非线性效应,目前对于水下爆炸气泡运动的计算方法无法给出较为可信的水射流运动特征及其载荷形式.本文基于多相可压缩流体的five-equation计算模型,引入界面函数和界面密度压缩技术,采用5阶WENO重构与HLLC近似Riemann求解器进行空间离散,时间离散采用3阶TVD Runge-Kutta法.通过上述方法来提高水下爆炸气泡溃灭过程中气-液界面的计算精度,捕捉水射流的产生、发展以及水锤冲击等典型演化过程.计算结果表明,新的计算方法能够对水下近壁面爆炸气泡的溃灭过程进行有效的模拟,初步揭示了水射流直接载荷的特征,为水射流的产生机理及其损伤效应的研究提供技术支撑.