黏性对内空泡诱导柱状液滴界面不稳定性的影响
INFLUENCE OF LIQUID VISCOSITY ON THE INTERFACIAL INSTABILITY OF CYLINDRICAL DROPLET INDUCED BY A CAVITATION BUBBLE
-
摘要: Rayleigh-Taylor不稳定性存在于爆炸、液滴形成和液体喷雾等工程应用过程中, 是流体力学关注的经典问题之一. 内空泡振荡诱导液滴界面演化问题是其研究中基本模型之一, 空泡振荡作用下液滴界面发生扰动并发展, 其特征形态主要表现为破碎、通气和稳定. 液体黏性是影响界面不稳定性发展的重要因素, 文章通过建立高精度的数值模拟方法, 开展液体黏性对内空泡诱导柱状液滴界面不稳定性的影响研究. 在数值模拟中, 基于开源OpenFOAM框架的多相可压缩求解器直接求解Navier-Stokes方程, 采用isoAdvector的几何流体体积法捕捉界面演化特征. 结果表明, 液体黏性的增加会减缓空泡的收缩, 进而减缓液滴界面扰动的发展, 该影响下通气工况液滴通气发生时间增加, 而稳定工况最大扰动幅值减小. 最大扰动幅值的减小直接影响了液滴的特征形态, 基于一系列数值模拟结果归纳得到液滴不稳定性相图. 在文章讨论的参数范围内, 随着黏性增加, 小液滴(Rd0 < 2 mm)的形态从破碎转变为通气进而变成稳定; 中液滴(2 mm < Rd0 < 3 mm)的形态从通气转变为稳定, 不出现破碎形态; 而大液滴(Rd0 > 3 mm)的形态不随液体黏性改变, 液滴形态会随着液体黏性增加趋于稳定.Abstract: Rayleigh-Taylor instability is one of the classical problems in fluid mechanics, which is widely used in underwater explosions, droplet formation, liquid sprays, and other engineering applications. The interfacial instability of cylindrical droplet induced by a cavitation bubble is one of the classical models in its study. As the perturbation of the interfacial instability develops, three distinct characteristics of the droplet deformation are obtained: (i) a splashing state; (ii) a ventilating state; and (iii) a stable state. Considering that liquid viscosity is an important factor affecting the development of interfacial instability, the influence of liquid viscosity on the interfacial instability of a cylindrical droplet induced by a cavitation bubble is investigated through a high-precision numerical simulation method in this paper. A direct numerical simulation is set up based on the compressible multiphase solver in the OpenFOAM framework. The interfaces are captured using a geometric fluid volume method named isoAdvector. The results show that within the parameters discussed in the paper, an increase in liquid viscosity slows down the oscillation of the bubble, thus reducing the perturbation growth rate during the bubble collapse. In these cases, the onset time of ventilation increases for a ventilating state, and the maximum perturbation reduces for a stable state. This reduction directly affects the distinct characteristics of the droplets. A phase diagram of the interfacial instability of droplets with different initial radius and liquid viscosity is obtained based on the results of numerical simulations. With the increase of liquid viscosity, the small droplets (Rd0 < 2 mm) change from a splashing state to a ventilating state and then a stable state; the medium droplets (2 mm < Rd0 < 3 mm) change from a ventilating state to a stable state; whereas, the large droplets (Rd0 > 3 mm) does not change and stay only in a stable state. The droplet tends to stabilize as the liquid viscosity increases.