两空泡运动特性研究

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

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

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

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

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

收缩扩张管内液氮空化流动演化过程试验研究

本文基于低温空化试验平台研究了收缩扩张流道内液氮非定常空化流动的演化过程. 试验采用高时空分辨率的高速摄像机对77 K液氮在不同空化数σ下空穴结构的演变进行了精细化的分析和研究. 利用试验得到的空穴长度和面积等数据, 定量分析了液氮空化流动的非定常特性与时空演变规律. 研究结果表明: (1)在相似来流速度和温度条件下, 随着空化数的减小, 液氮空化流动呈现四种典型流型, 空穴长度在2.5 h以内为初生空化、空穴长度在2.5 h ~ 7.5 h之间为片状空化、空穴长度在7.5 h ~ 15 h之间为大尺度云状空化, 空穴长度超过15 h为双云状空化, 且在大尺度云状空化和双云状空化阶段均捕捉到了回射流现象; (2)液氮空化流动从初生空化到双云状空化, 脱落空穴的尺度逐渐增大, 空穴面积脉动的幅值和准周期均有所增加. 同时, 在大尺度云状空化与双云状空化阶段, 喉口处堵塞效应对空化流动的影响显著增强; (3)相比于初生空化, 片状空化、大尺度云状空化以及双云状空化中脱落空穴的移动距离依次增加了0.97倍、2.65倍与2.68倍, 溃灭时间依次增加了1.18倍、3.59倍与4.47倍, 但溃灭速度依次减小了0.10倍、0.20倍与0.30倍. 除此之外, 对于双云状空化阶段, 存在两种显著不同的脱落空穴演化过程.      

垂直发射条件下水下航行体头型对通气空泡流动及压力特性的影响分析

为了揭示垂直发射条件下水下航行体头型对通气空泡演化过程的力学影响机理,首先基于有限体积法,结合改进型延迟分离涡模型、流体体积多相流模型及重叠网格技术建立了垂直发射条件下通气空泡的数值计算模型.其次,将计算结果与垂直发射实验进行对比,验证了所提出的数值方法对通气云空泡的预测具有较高精度,说明了该方法在通气空泡复杂非定常计算中的适用性.最后,对比研究了相同工况下流线头型和钝头头型航行体通气空泡流动特性和压力特性的差异,从涡量动力学的角度分析了差异产生的原因,结果表明:相比于流线头型航行体,钝头航行体通气空泡气液交界面处速度梯度较小,受到重力和浮力的影响更大,在瑞利-泰勒不稳定性机制的作用下,通气空泡更早发生非线性失稳,空泡失稳区域呈现更为剧烈的浮动行为以及空泡脱落等非定常流动特性;较强的空泡非定常流动特性影响了钝头航行体通气空泡末端的流动分离,从而抑制了空泡末端滞止高压的高幅值特性.     

激波与转捩边界层干扰非定常特性数值分析

激波与边界层干扰的非定常问题是高速飞行器气动设计中基础研究内容之一.以往研究主要针对层流和湍流干扰,在分离激波低频振荡及其内在机理方面存在着上游机制和下游机制两类截然不同的理论解释.分析激波与转捩边界层干扰下非定常运动现象有助于进一步加深理解边界层状态以及分离泡结构对低频振荡特性的影响规律,为揭示其产生机理指出新的方向.采用直接数值模拟方法对来流马赫数2.9,24?压缩拐角内激波与转捩边界层干扰下激波的非定常运动特性进行了数值分析.通过在拐角上游平板特定的流向位置添加吹吸扰动激发流动转捩,使得进入拐角的边界层处于转捩初期阶段.在验证了计算程序可靠性的基础上,详细分析了转捩干扰下激波运动的间歇性和振荡特征,着重研究了分离泡展向三维结构对激波振荡特性的影响规律,最后还初步探索了转捩干扰下激波低频振荡产生的物理机制.研究结果表明:分离激波的非定常运动仍存在强间歇性和低频振荡特征,其时间尺度约为上游无干扰区内脉动信号特征尺度的10倍量级;分离泡展向三维结构不会对分离激波的低频振荡特征产生实质影响.依据瞬态脉动流场的低通滤波结果,转捩干扰下激波低频振荡的诱因来源于拐角干扰区下游,与流场中分离泡的收缩/膨胀运动存在一定的关联.     

黏性对内空泡诱导柱状液滴界面不稳定性的影响

Rayleigh-Taylor不稳定性存在于爆炸、液滴形成和液体喷雾等工程应用过程中,是流体力学关注的经典问题之一.内空泡振荡诱导液滴界面演化问题是其研究中基本模型之一,空泡振荡作用下液滴界面发生扰动并发展,其特征形态主要表现为破碎、通气和稳定.液体黏性是影响界面不稳定性发展的重要因素,文章通过建立高精度的数值模拟方法,开展液体黏性对内空泡诱导柱状液滴界面不稳定性的影响研究.在数值模拟中,基于开源OpenFOAM框架的多相可压缩求解器直接求解Navier-Stokes方程,采用isoAdvector的几何流体体积法捕捉界面演化特征.结果表明,液体黏性的增加会减缓空泡的收缩,进而减缓液滴界面扰动的发展,该影响下通气工况液滴通气发生时间增加,而稳定工况最大扰动幅值减小.最大扰动幅值的减小直接影响了液滴的特征形态,基于一系列数值模拟结果归纳得到液滴不稳定性相图.在文章讨论的参数范围内,随着黏性增加,小液滴(R_d0 <2 mm)的形态从破碎转变为通气进而变成稳定;中液滴(2 mm d0 <3 mm)的形态从通气转变为稳定,不出现...     

可压缩空化流动空穴演化及压力脉动特性实验研究

空化流动具有高度的压缩性,空化流动非定常特性及其流体动力与压缩性密切相关.为研究可压缩空化流动空泡脱落的回射流和激波机制下周期性空穴结构演化及其诱导流体动力特性,本文采用多场同步测试方法对典型云状空化流动进行了实验研究,获得了文丘里管扩张段内部云状空化空穴形态演化及其诱导同步壁面压力脉动信号.并基于数字图像处理技术,对附着型片状空穴和脱落型云状空穴结构演化进行了精细化定量分析.结果表明:可压缩空化流动回射流机制下,空穴演化呈现附着型空穴生长-回射流产生及发展-附着型空穴失稳断裂及大尺度空泡云团产生脱落的非定常过程,激波机制下空穴演化具有附着型空穴生长-激波产生及传播-附着型空穴失稳断裂及大尺度空泡云团脱落的非定常特征,激波传播时间占空穴脱落周期小于回射流推进.激波与空穴相互作用导致空穴内部含气率瞬间大范围大幅度下降,诱导复杂流体动力.激波传播过程中,空泡内部压力脉动大幅增加,激波前缘诱导压力脉冲,而回射流推进过程中,壁面压力脉动相对平稳,回射流头部存在小幅增加.不同机制下空穴结构存在显著差异,具有不同的相间质量传输过程.     

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

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

自由场空泡溃灭过程能量转化机制研究

综合应用实验与数值模拟方法, 深入讨论了自由场空泡溃灭过程中的能量转化机制. 在实验研究中, 应用纹影法记录了空泡溃灭的演变过程, 提取了空泡在溃灭过程中的半径, 溃灭速度等数据, 结合空泡势能和动能方程, 描述了空泡能量的转化过程. 在开展数值模拟分析时, 运用弱可压缩流体质量守恒方程和动量方程, 建立了三维数值模型用以模拟空泡在自由场中的溃灭过程, 并且由结果中获取了空泡溃灭过程中的压力及速度变化规律, 揭示了空泡在溃灭过程中能量转化机制. 研究结果表明: (1) 自由场空泡在溃灭过程中, 空泡势能与空泡半径具有相同的演化趋势, 空泡动能与势能变化趋势相反; 当空泡达到最大半径处时, 空泡势能最大, 流场动能为零. (2) 溃灭后期在空泡周围会形成高压区域, 该区域的压力梯度与速度梯度较高, 随着空泡收缩, 高压区域面积逐渐减小. (3) 空泡在自由场中发生溃灭时, 空泡势能不断转化为流场动能, 在溃灭时刻可以明显观察到冲击波现象, 空泡的大部分能量会在此时转化为冲击波的波能.      

Bubbles in liquids with phase transition

In the forthcoming second part of this paper a system of balance laws for a multi-phase mixture with many dispersed bubbles in liquid is derived where phase transition is taken into account. The exchange terms for mass, momentum and energy explicitly depend on evolution laws for total mass, radius and temperature of single bubbles. Therefore in the current paper we consider a single bubble of vapor and inert gas surrounded by the corresponding liquid phase. The creation of bubbles, e.g. by nucleation is not taken into account. We study the behavior of this bubble due to condensation and evaporation at the interface. The aim is to find evolution laws for total mass, radius and temperature of the bubble, which should be as simple as possible but consider all relevant physical effects. Special attention is given to the effects of surface tension and heat production on the bubble dynamics as well as the propagation of acoustic elastic waves by including slight compressibility of the liquid phase. Separately we study the influence of the three phenomena heat conduction, elastic waves and phase transition on the evolution of the bubble. We find ordinary differential equations that describe the bubble dynamics. It turns out that the elastic waves in the liquid are of greatest importance to the dynamics of the bubble radius. The phase transition has a strong influence on the evolution of the temperature, in particular at the interface. Furthermore the phase transition leads to a drastic change of the water content in the bubble. It is shown that a rebounding bubble is only possible, if it contains in addition an inert gas. In Part 2 of the current paper the equations derived are sought in order to close the system of equations for multi-phase mixture balance laws for dispersed bubbles in liquids involving phase change.     

Evolution of a Small Distortion of the Spherical Shape of a Gas Bubble under Strong Expansion-Compression

The evolution of a small distortion of the spherical shape of a gas bubble which undergoes strong radial expansion-compression upon a single oscillation of the ambient liquid pressure under a harmonic law are analyzed by numerical experiments. It is assumed that the distortions of the spherical bubble shape are axisymmetric and have the form of individual spherical surface harmonics with numbers of 2–5. Bubble-shape oscillations prior to the beginning of expansion are taken into account. Generally, the distortion value during bubble expansion-compression depends on the phase of bubble-shape oscillation at the beginning of the expansion (initial phase). Emphasis is placed on the dependence of the maximum distortions in the initial phase at certain characteristic times of bubble expansion-compression on the amplitude of the external excitation, liquid viscosity, and distortion mode (harmonic number). The parameters of the problem are typical of the stable periodic sonolumiescence of an individual air bubble in water at room temperature. An exception is the liquid pressure oscillation amplitude, which is varied up to values that are five times the static pressure. That large excitation amplitudes are beyond the stability threshold of periodic oscillations of spherical bubbles. Their consideration is of interest from the point of view of increasing the compression ratio of the bubble gas, i. e., increasing the maximum temperature and density achievable in the final compression stage.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 4, pp. 17–28, July–August, 2005.     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

The effect of high viscosity on compressible and incompressible Rayleigh–Plesset-type bubble models

Free oscillations of a single spherical gas bubble in glycerol have been examined numerically and experimentally at different ambient temperatures and pressures. The bubble was generated using a Q-switched Nd:YAG laser and the unsteady radius measurement was based on a shadowing technique of a He–Ne laser beam. The measurements were compared to computations obtained from two models, first taking into consideration the liquid compressibility and then assuming an incompressible liquid domain, respectively. In both cases the temperature fields inside and outside the bubble were computed by solving the energy equation in both phases as the thermodynamic processes have great importance to the bubble behavior. For high amplitude oscillations the incompressible model provides poor agreement with the measurements and the modeling of the liquid compressibility becomes necessary. In contrast to the standard method, a practical region of applicability for the incompressible approach was determined as a function of the instantaneous Mach and Reynolds numbers, rather than specifying a simple threshold Mach number.     

The influence of thermodynamic gas parameters on laser-induced bubble dynamics in water

The oscillating properties of laser-induced cavitation bubbles in water are investigated by means of a fiber-optic sensor based on optical beam deflection. The experimental results show two important points. One is that the smaller the bubble radius the more quickly the bubble surface moves. Thus, the variations of the temperature and the pressure inside the bubble will be close to those of an adiabatic process. The other is that the high-energy vapor inside the newborn bubble diffuses and coagulates rapidly through violent expansion and thermal conduction. Thus, the gas content of the bubble reduces significantly in the first oscillation. Numerical simulation is made for the bubble model with consideration of liquid viscosity, surface tension, and gas content. Through modification of the polytropic index and the gas content parameter, two parameters of this model, the numerical results fit the experimental results well.     

Bubbles in liquids with phase transition—part 2: on balance laws for mixture theories of disperse vapor bubbles in liquid with phase change

We study averaging methods for the derivation of mixture equations for disperse vapor bubbles in liquids. The carrier liquid is modeled as a continuum, whereas simplified assumptions are made for the disperse bubble phase. An approach due to Petrov and Voinov is extended to derive mixture equations for the case that there is a phase transition between the carrier liquid and the vapor bubbles in water. We end up with a system of balance laws for a multi-phase mixture, which is completely in divergence form. Additional non-differential source terms describe the exchange of mass, momentum and energy between the phases. The sources depend explicitly on evolution laws for the total mass, the radius and the temperature of single bubbles. These evolution laws are derived in a prior article (Dreyer et al. in Cont Mech Thermodyn. doi:10.1007/s00161-0225-6, 2011) and are used to close the system. Finally, numerical examples are presented.     

Seed bubbles stabilize flow and heat transfer in parallel microchannels

Seed bubbles are generated on microheaters located at the microchannel upstream and driven by a pulse voltage signal, to improve flow and heat transfer performance in microchannels. The present study investigates how seed bubbles stabilize flow and heat transfer in micro-boiling systems. For the forced convection flow, when heat flux at the wall surface is continuously increased, flow instability is self-sustained in microchannels with large oscillation amplitudes and long periods. Introduction of seed bubbles in time sequence improves flow and heat transfer performance significantly. Low frequency (∼10 Hz) seed bubbles not only decrease oscillation amplitudes of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures, but also shorten oscillation cycle periods. High frequency (∼100 Hz or high) seed bubbles completely suppress the flow instability and the heat transfer system displays stable parameters of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures. Flow visualizations show that a quasi-stable boundary interface from spheric bubble to elongated bubble is maintained in a very narrow distance range at any time. The seed bubble technique almost does not increase the pressure drop across microsystems, which is thoroughly different from those reported in the literature. The higher the seed bubble frequency, the more decreased heating surface temperatures are. A saturation seed bubble frequency of 1000–2000 Hz can be reached, at which heat transfer enhancement attains the maximum degree, inferring a complete thermal equilibrium of vapor and liquid phases in microchannels. Benefits of the seed bubble technique are the stabilization of flow and heat transfer, decreasing heating surface temperatures and improving temperature uniformity of the heating surface.     

INVESTIGATION OF BUBBLE-BUBBLE INTERACTION EFFECT DURING THE COLLAPSE OF MULTI-BUBBLE SYSTEM

通过直接数值模拟方法对多泡在压力驱动下的溃灭过程进行了研究。气相满足理想气体正压模型,液相为不可压 流体,采用基于压力的方法求解多泡的两相流场。数值研究表明,在多泡流场中,中心气泡的溃灭过程明显不同于单泡,存在总体溃灭延迟现象和后期加速现象。随着周围气泡数的增多或气泡间距的减小,中心气泡的溃灭时间增长,溃灭时的压力峰值增大。结合理论定性分析发现,气泡运动不仅受远场压力的驱动,还受周围气泡诱导压力场的影响。周围气泡的诱导压力经历先减小后增大的过程,从而使受其影响的中心气泡产生先延迟后加速的特征。     

The transient buoyancy driven motion of bubbles across a two-dimensional quiescent domain

The transient buoyancy driven motion of two-dimensional bubbles across a domain bounded by two horizontal walls is studied by direct numerical simulations. The bubbles are initially released next to the lower wall and as they rise, they disperse. Eventually all the bubbles collect at the top wall. The goal of the study is to examine how a simple one-dimensional model for the averaged void fraction captures the unsteady bubble motion. By using void fraction dependent velocities, where the exact dependency is obtained from simulations of homogeneous bubbly flows, the overall dispersion of the bubbles is predicted. Significant differences remain, however. We suggest that bubble dispersion by the bubble induced liquid velocity must be included, and by using a simple model for the bubble dispersion we show improved agreement.     

Evolution of distortions of the spherical shape of a cavitation bubble in acoustic supercompression

The evolution of small perturbations of the spherical shape of a vapor bubble in the process of its single strong expansion and compression in deuterated acetone is studied. In the mathematical model used the motion of vapor and liquid is broken down into the spherical component and its small nonspherical perturbation. The spherical component is described by the fluid dynamics equations with account for time-dependent heat conduction and evaporation and condensation on the liquid-vapor interface using equations of state constructed from experimental data. In describing the nonspherical component the liquid viscosity and the surface tension are taken into account, while the effect of the bubble content is disregarded. Certain simple analytical formulas are presented which describe the bubble radius at the moment of maximum expansion, its variation in the compression stage, and the evolution of the bubble sphericity distortion in compression.