微小水滴撞击深水液池空腔运动的数值模拟及机理研究

Numerical investigation of cavity formation mechanism for micron-waterdrop impact on deep pool

为了探究微米级微小水滴撞击深水液池运动中空腔的成长过程与机理,采用自适应网格技术和流体体积方法对撞击速度为2.5—6.5 m/s的微小水滴撞击深水液池的运动进行数值模拟研究,考察不同撞击速度下水滴撞击深水液池后的水体混掺、毛细波传播、空腔变形规律以及气泡截留过程,并深入探究空腔运动的动力学机制.研究结果表明,不同撞击速度下,在忽略毛细波作用、空腔深度h∈(D, h_(max))的前提下,空腔深度随时间的成长仍满足t∝h^(5/2)的关系;液滴撞击产生的空腔形状有U形和半球形两种,前者一般向V形转变,后者空腔底部会变为圆柱形,产生细长射流,并有可能发生气泡截留现象;在撞击速度较低时,低压区首先在空腔侧壁与底部交界处产生,随后在靠近液面以及空腔底部靠近中心区域各产生一个较大的涡环;在撞击速度较高,产生细长射流时,涡环的生成被抑制,低压区首先在波浪底部与侧壁上交界处产生,随后空腔底部变为圆柱状,空腔侧壁首先坍塌形成气泡截留.

As one of the most fundamental and iconic fluid motion, droplet impact exists widely in scientific technologies and natural environment, and the phenomenon has been studied both for fundamental mechanism and for industrial applications in aerospace engineering, inkjet printing, agricultural irrigation and hydraulic structure erosion. Therefore,it is of great significance to study such basic movements for understanding the interfacial deformation of gas and liquid flow and improving the applications of droplet impact movement in engineering. Droplet impacting on a deep liquid pool has been extensively investigated for droplets with millimeter diameter. In this article, focusing on the cavity formation mechanism during a Micron-sized waterdrop impact on a deep pool, we perform systematic numerical simulations with adaptive mesh refinement technique and volume of fluid method to study the impact of a 290 μm water droplet on a deep water pool at velocities in a range of 2.5–6.5 m/s. The free surface motion, geometric variation of the cavity, local pressure field and vorticity field at selected times are presented to identify the pool-drop water mixing, capillary wave propagation, cavity formation, vortex ring generation and bubble entrapment phenomenon, and the dynamic mechanism of cavity motion is further explored. It is found that under the premise of neglecting the surface tension effects on the cavity whose depth is in a range of h ∈（D, h（max））, where D is the radius of initial droplet and hmax is the maximum depth,the cavity growth time to reach its maximum depth still scales as t ∝ h^5/2, where t is time, but in the end, the formation of the bottom of the cavity is driven by capillary waves. There are two types of the initial cavity shapes： one is U-shape and the other is hemispherical shape, the former one generally changes into V-shape, and in the latter case, the bottom of the cavity will gradually transform into cylindrical shape, resulting in a thin jet and possible bubble entrapment. Cavity collapse is closely related to capillary wave propagation. When the impact velocity is low（Fr = 567.1, Re = 1595,W e = 121.8）, the low-pressure zone is initially generated at the junction between the cavity sidewall and the bottom, a large vortex ring is then generated near the free surface and the bottom of the cavity, respectively. Under high impact velocities（Fr = 792.1, Re = 1885, We = 170.2）, the thin jet is observed, the generation of the vortex ring is suppressed.The low-pressure zone is first generated at the junction between the wave bottom and the cavity sidewall, after the cavity becomes cylindrical, the cavity collapses before the capillary wave arrives at the bottom, causing a bubble entrapment.