超声驱动下激励参数对单泡空化振动的影响

根据考虑了液体可压缩性的改进的微气泡动力学方程,采用改进的初始半径对单泡超声空化现象进行了数值计算研究.结果表明,微气泡振动对一些参量很敏感:微气泡振动半径与初始半径的比值随振动频率的增大而减小;提高声场声压会加剧气泡崩塌程度,但过高的声压又不能使微气泡崩塌;微气泡崩塌速率随气泡初始半径的增加而增大,在一定范嘲内能保证空化泡稳定振动,在初始半径为1.6μm处空化程度最强,如果继续增大初始半径则空化程度减弱、甚至消失;微气泡崩塌程度随黏滞系数和表面张力的增大而减弱,过大的黏滞系数和表面张力会使微气泡崩塌难以发生.计算结果与他人的实验数据相比,发现液体的可压缩性使单泡空化强度增强,对最佳空化区域范围的确定有较大的影响.     

DYNAMICAL BEHAVIORS OF DOUBLE CAVITATION BUBBLES UNDER ULTRASONIC HONING

为了探讨超声珩磨作用下磨削区的空化机理,基于速度势叠加原理,考虑超声珩磨速度和珩磨压力,建立了磨削区两空化泡的动力学模型. 数值模拟了磨削区空化泡初始半径、两空化泡间距、超声声压幅值、珩磨压力、珩磨头转速对磨削区两空化泡动力学特性的影响. 研究表明,考虑两空化泡之间的相互作用时,要想获得良好的空化效果,可将两空化泡初始半径之比控制在3 倍以内;选择较高的超声波声压幅值与较低的珩磨压力,并且使超声波声压幅值与珩磨压力和液体静压力之差介于0.66~1.89MPa 之间;增大珩磨头转速空化泡溃灭也略有加速;通过试验测量材料表面粗糙度的方法间接验证了理论分析的合理性.     

超声珩磨作用下两空化泡动力学特性

为了探讨超声珩磨作用下磨削区的空化机理,基于速度势叠加原理,考虑超声珩磨速度和珩磨压力,建立了磨削区两空化泡的动力学模型. 数值模拟了磨削区空化泡初始半径、两空化泡间距、超声声压幅值、珩磨压力、珩磨头转速对磨削区两空化泡动力学特性的影响. 研究表明,考虑两空化泡之间的相互作用时,要想获得良好的空化效果,可将两空化泡初始半径之比控制在3 倍以内;选择较高的超声波声压幅值与较低的珩磨压力,并且使超声波声压幅值与珩磨压力和液体静压力之差介于0.66~1.89MPa 之间;增大珩磨头转速空化泡溃灭也略有加速;通过试验测量材料表面粗糙度的方法间接验证了理论分析的合理性.      

Parametric Analysis for a Single Collapsing Bubble

This work presents a sensitivity analysis for cavitation processes, studying in detail the effect of various model parameters on the bubble collapse. A complete model (Hauke et al. Phys Rev E 75:1–14, 2007) is used to obtain how different parameters influence the collapse in SBSL experiments, providing some clues on how to enhance the bubble implosion in real systems. The initial bubble radius, the frequency and the amplitude of the pressure wave are the most important parameters determining under which conditions cavitation occurs. The range of bubble sizes inducing strong implosions for different frequencies is computed; the initial radius is the most important parameter characterized the intensity of the cavitation processes. However, other parameters like the gas and liquid conductivity or the liquid viscosity can have an important effect under certain conditions. It is shown that mass transfer processes play an important role in order to correctly predict the trends related with the effect of the liquid temperature, which translates into the bubble dynamics. Moreover, under some particular circumstances, evaporation can be encountered during the bubble collapse; this can be profitably exploited in order to feed reactants when the most extreme conditions inside the bubbles are reached. Thus, this paper aims at providing a global assessment of the effect of the different parameters on the entire cycle of a single cavitating spherical bubble immersed in an ultrasonic field. This work has been partially supported by Ministerio de Ciencia y Tecnologia, under grant number CTM2004-06184-C02-02.     

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

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

空泡溃灭及空化噪声研究综述

本文从理论分析,数值计算,试验研究三方面对国内外空泡溃灭,空化噪声研究的最新进展及动态作了一简单的回顾,针对自己做的水翼及轴对称头体空化噪声试验结果,提出了一些尝试性解释。     

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

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

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.     

固壁空蚀数值研究

空蚀是空泡在固壁附近溃灭对固壁材料产生破坏的现象。本文将空泡界面假设为自由面，并由VOF(Volume of Fluid)中界面构造精度较高的Youngs方法求解，通过直接计算原始变量的Navier-Stokes方程，数值模拟了空泡距固壁不同位置时溃灭对固壁造成的空蚀破坏。计算发现空泡溃灭产生高压脉冲相对于高速射流对空蚀形成起主导作用；空泡在流场中位置不同，高压脉冲对固壁上的空蚀破坏结果不同，并给出了距离界限。     

A numerical scheme for Euler–Lagrange simulation of bubbly flows in complex systems

An Eulerian–Lagrangian approach is developed for the simulation of turbulent bubbly flows in complex systems. The liquid phase is treated as a continuum and the Navier–Stokes equations are solved in an unstructured grid, finite volume framework for turbulent flows. The dynamics of the disperse phase is modeled in a Lagrangian frame and includes models for the motion of each individual bubble, bubble size variations due to the local pressure changes, and interactions among the bubbles and with boundaries. The bubble growth/collapse is modeled by the Rayleigh–Plesset (RP) equation. Three modeling approaches are considered: (a) one‐way coupling, where the influence of the bubble on the fluid flow is neglected, (b) two‐way coupling, where the momentum‐exchange between the fluid and the bubbles is modeled, and (c) volumetric coupling, where the volumetric displacement of the fluid by the bubble motion and the momentum‐exchange are modeled. A novel adaptive time‐stepping scheme based on stability‐analysis of the non‐linear bubble dynamics equations is developed. The numerical approach is verified for various single bubble test cases to show second‐order accuracy. Interactions of multiple bubbles with vortical flows are simulated to study the effectiveness of the volumetric coupling approach in predicting the flow features observed experimentally. Finally, the numerical approach is used to perform a large‐eddy simulation in two configurations: (i) flow over a cavity to predict small‐scale cavitation and inception and (ii) a rising dense bubble plume in a stationary water column. The results show good predictive capability of the numerical algorithm in capturing complex flow features. Copyright © 2010 John Wiley & Sons, Ltd.     

圆锥泡声致发光气泡动力学过程的理论分析

在绝热压缩模型的基础上, 详细讨论了圆锥泡声致发光中气泡运动的动力学过程，得到 了气泡塌陷速度方程、气泡内压强方程以及温度方程. 结果显示在气泡进入圆锥腔的初始阶 段，气泡的塌陷速度随着压缩半径的不断减小近似线性地增加；然后随着压缩半径的进一步 减小，气泡塌陷的加速度逐渐减小；当气泡塌陷速度达到最大值后，随着气泡压缩半径的 进一步减小, 塌陷速度迅速下降至零. 在假设初始气压为1000\,Pa的基础上，理论分析 得到气泡的最高塌陷速度可以达到5.8\,m/s; 气泡的最小压缩半径可以达 到1.37\,cm, 相应的气泡内极限压强超过$4.5\times10^5$\,Pa, 极限温度超 过3\,150\,K, 而液流能够提供给气泡的能量达到0.02\,J. 理论推导得到的结果 可以比较好地用来解释实验中的现象. 最后分析得到气泡内的初始气 压对气泡所能达到的极端条件有着重要的影响.     

泡内气体热力学性质对空泡溃灭的影响

数值研究固壁附近轴对称空泡溃灭问题. 忽略泡内气体与周围流体之间的质量和热交换, 假设气体瞬时处于热平衡状态, 通过引入不同的热力学模型, 考察泡内气体在空泡溃灭过程中的作用. 采用原始变量的Navier-Stokes方程作为流场的控制方程, 用流体体积方法跟踪运动空泡壁. 数值结果显示空泡溃灭过程中, 伴随空泡变形, 空泡发出多个高压脉冲和高速射流. 对于不同的热力学模型, 等温, 绝热和准绝热过程, 绝热过程能够最大程度抑制空泡溃灭, 从而减弱空泡溃灭对固壁造成的空蚀破坏. 在绝热及其类似过程中, 出现空泡回弹现象.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

A numerical model for cloud cavitation based on bubble cluster

The cavitation cloud of different internal structures results in different collapse pressures owing to the interaction among bubbles. The internal structure of cloud cavitation is required to accurately predict collapse pressure. A cavitation model was developed through dimensional analysis and direct numerical simulation of collapse of bubble cluster. Bubble number density was included in proposed model to characterize the internal structure of bubble cloud. Implemented on flows over a projectile, the proposed model predicts a higher collapse pressure compared with Singhal model. Results indicate that the collapse pressure of detached cavitation cloud is affected by bubble number density.     

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

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

固液两相流体中的空泡溃灭计算

本文导出了固液两相流体中球空泡溃灭的运动方程,计算并讨论了空泡溃灭过程中的颗粒运动和颗粒对空泡溃灭的影响,得到了固相浓度、颗粒尺寸等因素与空泡溃灭之间的定性关系。在分析过程中,考虑了液体与固体颗粒之间的阻力耦合作用。     

Bubble collapse in solid-liquid two-phase fluid

Equations of motion for bubble collapse in solid-liquid two-phase fluid have been derived, in which the resistance coupling effects between the liquid and solid particles have been considered. The motion of particles during the bubble collapse and the effects of particles on bubble collapse have been calculated and discussed. Qualitative relations between the concentration and the size of the particles and the rate of bubble collapse have been obtained. Institute of Water Conservancy and Hydroelectric Power Research