单空泡与自由液面相互作用规律研究进展

空化与空泡溃灭现象普遍存在于自然界、水力机械和生物医学等领域.空泡与自由面相互作用会产生瞬态强烈耦合,涉及到空泡非球形溃灭、自由面非线性变形及失稳等现象,是流体力学领域重要的前沿与基础问题.本文围绕这一热点,从空泡非球形演化和自由面变形规律角度出发,概述和归纳近年该领域的研究进展与成果.对于近自由面空泡的非球形演化,基于表征开尔文冲量的无量纲参数,重点关注了体积振荡、射流生成、水锤效应及溃灭冲击波生成等关键过程,介绍了关键参数的理论建模方法,获得了空泡溃灭过程中能量分配机制.针对自由液面变形演化,根据细射流和粗射流生成和发展,归纳了4种典型现象及特点:透明水层及水柱生成、不稳定与稳定冠状水裙结构.进一步总结了开尔文冲量理论、界面凹陷奇点概念和泰勒不稳定性等理论模型的建立和应用,讨论了气泡溃灭过程、液面射流和界面稳定性等主要机制.此外,本文也概述了空泡脉动对球状、圆柱状等非平面液面变形行为的影响,归纳了曲率对于液面变形的影响机制.最后,针对目前研究状况提出该领域研究中尚未解决的问题,期望对将来的空泡及空泡群与自由液面相互作用深入研究提供借鉴.     

淹没磨料射流的空泡运动分析

通过对淹没条件下磨料射流的空泡运动分析研究,建立了淹没磨料射流的空泡运动方程,揭示了淹没磨料射流中空泡的溃灭特性,数值模拟了淹没磨料射流的磨料体积浓度以及空泡所处流场压力对空泡运动及溃灭的影响规律.分析表明:淹没磨料射流中磨料的存在增大流体的粘性系数,增大空泡溃灭历时,减弱射流的空蚀破坏能力;流场压力的改变对空泡溃灭过程影响显著,压力越高,空泡溃灭历时越短.     

两空泡运动特性研究

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

Fluent环境中近壁面微空泡溃灭的仿真计算

基于FLUENT软件环境,采用VOF模型和非稳态方法求解Navier-Stokes方程,模拟了近壁面的空泡溃灭过程,同时计算了空泡溃灭处与壁面的距离对射流强度的影响.结果表明：在近壁面,空泡将形成非对称溃灭,因水锤作用,引发高速水射流在壁面产生高压而形成空蚀破坏;基于FLUENT环境的计算结果与已有的实验和计算结果相符,为研究空泡溃灭和空蚀机制提供了类比的数值计算方法.     

近固壁微米尺度空泡溃灭的数值研究

考虑空泡表面张力、液体黏性和气体可压缩性,采用VOF多相流模型对近固壁微米尺度空泡在静止流场溃灭过程进行了数值研究.获得了近固壁空泡溃灭过程的流场细节,分析了空泡与固壁的无量纲距离γ对空泡溃灭过程动力特性的影响,并揭示了不同γ条件下的固壁空蚀破坏机理.计算结果表明:随着γ的减小,泡心向固壁移动的趋势明显,射流形成前空泡上部高压区内压力减小,空泡溃灭时间延长,最大射流速度减小.模拟结果验证了空泡溃灭将产生冲击波和高速微射流,二者均会在固壁面产生脉冲压力,其是造成壁面损伤的两种主要原因.参数γ对固壁的空蚀破坏机理有重要影响.与微射流机制相比,以冲击波机制为主的空蚀破坏更显著.微射流冲击固壁的作用半径为10μm左右,将引起固壁"点"蚀坑的出现.当γ=2.0时,冲击波扫掠壁面的范围相对较广,有效作用半径约为1 mm,其导致固壁产生较大圆形蚀坑,且中心空蚀严重.     

绕水翼空化流动多尺度数值研究

空化的多尺度效应是一种涉及连续介质尺度、微尺度空化泡以及不同尺度间相互转化的复杂水动力学现象, 跨尺度模型的构建是解析该多尺度现象的关键. 本文基于欧拉-拉格朗日联合算法, 通过界面捕捉法求解欧拉体系下大尺度空穴演化, 通过拉格朗日体系下离散空泡模型求解亚网格尺度离散空泡的运动及生长溃灭. 同时, 通过判断空泡与网格尺度间的关系判定不同尺度空化泡的求解模型. 基于建立的多尺度算法对绕NACA66水翼空化流动进行模拟, 将数值结果与实验进行对比, 验证了数值计算方法的准确性. 研究结果表明, 离散空泡数量与空化发展阶段密切相关, 在附着型片状空穴生长阶段, 离散空泡数量波动较小, 离散空泡主要分布在气液交界面位置; 在回射流发展阶段, 离散空泡逐渐增加并分布在回射流扰动区; 在云状空穴溃灭阶段, 离散空泡数量增多且主要分布在气液掺混剧烈的空化云团溃灭区. 在各空化发展阶段, 离散空泡直径概率密度函数均符合伽玛分布. 空化湍流流场特性对拉格朗日空泡空间分布具有重要影响, 离散空泡主要分布在强湍脉动区、旋涡及回射流发展区域.      

脉冲磨料射流中球泡溃灭特性研究及数值分析

从理论上研究了脉冲磨料射流中空化球泡的溃灭特性,建立了脉冲磨料射流中空化球泡的动力学方程,数值计算分析了磨料浓度、磨料粒径、射流振荡频率和幅值以及流场压力对脉冲磨料射流中空泡溃灭过程的影响规律。研究表明,脉冲磨料射流中磨料的存在将增大空泡溃灭历时,从而减弱射流的空蚀破坏能力,其中磨料的粒径、浓度和流场压力影响幅度较大。     

高速航行体水下发射水动力学研究进展

高速航行体水下发射水动力学研究, 是具有重大工程应用背景的前沿基础问题.与之紧密相关的非定常空化流动, 特别是空泡稳定性、溃灭等问题, 是影响发射载荷及安全性的关键.本文首先简述了这一领域的主要科学问题, 归纳了主要控制参数和影响方式; 之后针对非定常空化流动问题, 综述了已有的实验观测手段及数值模拟方法; 总结了空泡发展、稳定性、溃灭及流动控制等重要物理机制、模型及各因素相互作用规律; 最后展望了该领域仍存在的主要科学问题与未来发展趋势.      

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

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

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

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

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

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

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.     

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.     

The influence of viscoelasticity on the collapse of cavitation bubbles near a rigid boundary

An understanding of the phenomena associated with cavitation is important in many areas of science and engineering. This paper is concerned with the influence of viscoelasticity on the dynamics of cavitation bubbles near rigid boundaries. Viscoelastic effects are modelled using a Maxwell constitutive equation, and a generalized Bernoulli equation is derived. The governing equations are solved using the boundary element method in which both the bubble surface and the potential are represented by cubic splines. The numerical scheme is validated through comparisons with results in the literature for the inviscid case. The introduction of viscoelasticity introduces some interesting bubble dynamics including the occurrence of oscillations during collapse. Most importantly, it is shown that viscoelasticity can serve to suppress the formation of a liquid jet. The subsequent reduced pressures compared with the inviscid case suggest that viscoelasticity has a mitigating effect on cavitation damage.     

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.     

The Peculiar Dynamics of Cavitation Bubbles

Details from cavitation bubble dynamics are reported: jet formation, counterjet formation, shock wave radiation and light emission. Multiple shock wave radiation from single bubble collapse with jet formation could be time resolved by high speed photography with 20 million frames per second. An explanation of counterjet formation is given. Pictures of the light emission (sonoluminescence) in acoustic cavitation are presented.     

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

由于存在强冲击、多相界面复杂运动等强非线性效应,目前对于水下爆炸气泡运动的计算方法无法给出较为可信的水射流运动特征及其载荷形式.本文基于多相可压缩流体的five-equation计算模型,引入界面函数和界面密度压缩技术,采用5阶WENO重构与HLLC近似Riemann求解器进行空间离散,时间离散采用3阶TVD Rung...     

Numerical simulation of trajectory and deformation of bubble in tip vortex

According to the behaviors of a bubble in the ship wake flow,the numerical simulation is divided into two stages,quasi-spherical motion and non-spherical motion,based on whether the bubble is captured by the vortex or not.The one-way coupled particle tracking method(PTM) and the boundary element method(BEM) are adopted to simulate these two stages,respectively.Meanwhile,the initial condition of the second stage is taken as the output of the first one,and the entire simulation is connected and completed.Based on the numerical results and the published experimental data,the cavitation inception is studied,and the wake bubble is tracked.Besides,the split of the bubble captured by the vortex and the following sub-bubbles are simulated,including motion,deformation,and collapse.The results provide some insights into the control on wake bubbles and optimization of the wake flow.     

Air entrapment during a high-speed liquid jet entry in a deep water surface

The fluid dynamic phenomena of a high speed liquid jet impact on a deep water surface have been studied using Imacon high-speed photography. Both framing and streak techniques are applied to investigate the initial impact stage and penetration velocity. The cavitation caused by air entrapment between two colliding liquid surfaces has been found. The bubble collapse experiences different stages in relation to the contact area, liquid shock wave, release wave and fluid convection.