气泡引起的皇冠型水冢实验与数值研究

气泡在近自由液面运动过程中与自由液面发生强烈的耦合作用, 自由液面会出现极其复杂的物理现象, 形成的水冢类型繁多. 本文将针对皇冠型水冢, 基于势流理论, 建立气泡与自由液面耦合作用的数值模型, 并在气泡完成射流向下运动时, 忽略气泡对自由液面的影响, 继续模拟自由液面的运动过程. 同时利用高速摄影对近自由面的电火花气泡进行实验研究, 数值结果与实验符合良好, 相对误差在10%以内. 通过数值计算, 发现了围裙卷缩和主峰珠化等特殊的物理现象, 研究了气泡初始条件与韦伯数对皇冠型水冢动态特性的影响, 旨为皇冠型水冢的研究提供参考. 关键词： 气泡 自由液面 皇冠型水冢 韦伯数     

气泡与弹性膜的耦合效应数值模拟

本文在前人研究的基础上,计入浮力、表面张力、不同流体密度比等因素对弹性膜附近气泡运动的影响,结合不可压缩理想流体理论,建立气泡与弹性膜耦合动力学数值模型,采用边界元方法进行求解,计算值与Turangan等的实验结果符合良好.通过对弹性膜附近气泡运动的数值模拟,详细分析了弹性膜两侧为同种密度液体以及不同密度液体时气泡的运动,随后又分析气泡在弹性膜和浮力的共同作用下气泡的射流特性.旨在为相关气泡与弹性膜相互作用特性的研究提供参考. 关键词： 弹性膜 边界元法 气泡 浮力     

基于边界元法的近平板圆孔气泡动力学行为研究

研究了带有圆孔的平板附近气泡动力学特性. 基于不可压缩势流理论, 建立了平板圆形破口附近气泡运动数值模型, 并针对气泡初始位置距离破口很近而导致计算结果发散的数值缺陷, 采用气泡壁和壁面融合的方法, 将流场分离为两个半无限域问题进行求解, 实现了在不同无量纲参数范围内的数值模拟, 数值结果与实验结果符合良好. 通过对圆孔附近气泡运动特性的研究发现, 圆孔对气泡的影响基本与壁面相反, 在膨胀阶段对气泡产生腔吸作用, 收缩阶段产生排斥, 在特定的工况下会产生对射流现象. 最后分析了气泡壁与壁面融合, 流场分离后的气泡动态特性以及各工况参数对其影响规律. 关键词： 气泡 边界元 射流 圆孔     

超声场下刚性界面附近溃灭空化气泡的速度分析

为了揭示刚性界面附近气泡空化参数与微射流的相互关系, 从两气泡控制方程出发, 利用镜像原理, 建立了考虑刚性壁面作用的空化泡动力学模型. 数值对比了刚性界面与自由界面下气泡的运动特性, 并分析了气泡初始半径、气泡到固壁面的距离、声压幅值和超声频率对气泡溃灭的影响. 在此基础上, 建立了气泡溃灭速度和微射流的相互关系. 结果表明: 刚性界面对气泡振动主要起到抑制作用; 气泡溃灭的剧烈程度随气泡初始半径和超声频率的增加而降低, 随着气泡到固壁面距离的增加而增加; 声压幅值存在最优值, 固壁面附近的气泡在该最优值下气泡溃灭最为剧烈; 通过研究气泡溃灭速度和微射流的关系发现, 调节气泡溃灭速度可以达到间接控制微射流的目的.     

水中上浮气泡动态特性研究

针对水中的上浮气泡,计入表面张力的作用,然后将气泡边界层内黏性效应分为法向和切向两部分,其中附加法向黏性应力通过Young-Laplace关系考虑;附加切向黏性力是基于黏性耗散能量等效原理,引入黏性修正压力代替.首先建立了轴对称和三维上浮气泡边界元模型,将数值结果与理论值和实验值进行对比分析,有良好的符合度,验证了数值模型的有效性;然后针对毫米量级上浮气泡的平衡速度与形态,讨论了气泡初始条件、表面张力和黏性对气泡上浮过程中动力学行为的影响;最后,提出了一种处理三维上浮气泡融合的数值方法,计算结果与实验现象符合良好,并且能够反映气泡融合后的复杂现象细节.     

刚性柱附近浅水爆炸荷载特性研究

刚性柱附近浅水爆炸时冲击波传播、气泡射流受多种因素影响。考虑水面、水底、刚性柱与水下爆炸冲击波及气泡的耦合作用,基于LS-DYNA有限元软件,建立浅水爆炸全耦合模型,通过经验公式验证有限元模型的正确性。研究表明:采用炸药直径1/3～1/2中心渐变网格能够较好地保证数值模拟精度。在冲击波传播阶段,刚性柱迎爆区冲击波峰值上升并产生切断现象,冲击波下降段被"截断",而背爆区冲击波峰值衰减约50%,同时正压作用时间增加;在气泡脉动阶段,气泡在收缩阶段产生指向刚性柱的气泡射流,当刚性柱与炸药之间的距离约为一个气泡半径时,刚性柱附近的脉冲荷载增幅最大,脉冲荷载最大测点水深较爆心上移。     

近壁面气泡的运动规律研究

基于势流假设,建立气泡与壁面耦合数值模型,运用边界积分法求解,并开发三维计算程序,计算值与实验值符合很好.从气泡与壁面相互作用的基本现象入手,基于开发的程序系统地研究了刚性壁面附近气泡的动力学特性,其中包括水平壁面及倾斜壁面,研究壁面的Bjerknes效应与各特征参数之间的关系,并将各种工况的计算结果与基于Kelvin-impulse理论的Blake准则进行对比分析讨论,得出偏射流方向及壁面压力与气泡的特征参数有密切的关系,同时给出了Blake准则的适用范围.旨在为相关的近壁面气泡动态特性研究提供参考. 关键词： 气泡 壁面 边界积分 Bjerknes效应     

气泡群的动态物理特性研究

假设气泡周围流场为无黏、无旋、不可压缩的理想流体,建立气泡群相互作用的三维数值模型.将多极快速傅里叶变换方法(FFTM)与高阶边界元法(HOBEM)相结合求解气泡群的运动,在达到同样计算精度时显著加快了边界积分方程的求解速度,可以在合理的时间内模拟气泡群的动态物理特性.同时为维持气泡群模拟过程中的数值稳定性,引入了弹性网格技术(EMT),并用算例验证了数值模型及算法的有效性.基于建立的数值模型,研究了不同组合的气泡群之间的相互作用,模拟和解释了各类气泡运动的物理现象,讨论了影响气泡群膨胀、坍塌、迁移及射流 关键词： 气泡群 FFTM 射流 三维     

有黏条件气泡声散射特性和衰减谱数值研究*

通过研究有黏条件下的气泡散射模型,数值分析水中单气泡声散射特性,进一步结合Beer-Lambert定律将其扩展到多气泡体系的声衰减预测。数值结果表明,随着谐波阶数递增,散射强度分布数值结果趋于稳定且前向散射增强。同时发现,无因次尺寸参量ka=0.1为过渡区与纯散射区的分界线,且在共振区间具有明显的吸收效应。对多气泡体系的声衰减预测也表明,ka 0.1时,该文气泡散射模型声衰减计算与经典ECAH模型结果吻合,在低浓度条件下声衰减谱随着剪切黏度的增加呈增宽趋势,且与体积浓度成正比例递增。模型预测的声衰减随粒径、声波频率、体积浓度分布数值特征能够为颗粒两相体系粒径及浓度表征提供理论依据。     

高粘度流体中上升气泡的直接数值模拟

本文将基于Front Tracking方法求解气液两相界面迁移特性的直接数值模拟方法扩展到中低Re数情形 (1．89≤Re≤55．6),针对无边界以及垂直壁面附近高粘度流体中单个气泡(1．2 mm≤D≤6 mm)的上升过程进行模拟,研究气泡运动的机理以及气泡与壁面的相互作用。数值模拟准确再现了气泡的上升过程和变形,不同Re数下气泡的上升速度计算结果同经验关联式非常吻合,验证了该方法在中低Re数范围内的有效性。壁面附近流场的不对称引发壁面对气泡的升力是导致气泡逐渐偏离垂直壁面的原因,计算得到的壁面作用力系数同实验结果的对比也确定了本文参数范围内计算壁面对气泡作用力的模型。     

曲面固壁附近空化泡溃灭动力学的 流体体积数值研究*

流体体积法(VOF)可以便捷、高效地实现对多相流界面的捕捉和追踪。本文基于VOF方法,对单个空化泡在曲面固壁附近的运动进行了数值模拟,从实验对比、压力场、速度场、温度场演化、溃灭时间、射流速度、固壁温度等方面分析了空化泡溃灭过程的热动力学影响。结果表明,数值模拟得到的空化泡形态演化与实验观测到的现象一致,随着位置参数、泡内外压差及曲面固壁尺寸的改变,空化泡热动力学行为也将发生变化,受到流体运动及射流冲击的影响,溃灭瞬间产生的高温高压使得曲面固壁温度升高。本文研究的曲面固壁附近空化泡溃灭效应,揭示了空化泡与曲面固壁间的相互作用规律,对学术研究及工程应用都具有重要意义。     

上浮气泡在壁面处的弹跳特性研究

本文针对毫米量级的上浮气泡在壁面处的弹跳现象进行数值研究.基于势流方法求解气泡的运动,同时考虑气泡的表面张力作用.在伯努利方程中,对气泡与壁面之间水膜中因黏性引起的压力梯度进行修正,开发相应的计算程序,计算值与实验值符合良好.从气泡弹跳的基本现象入手,研究了特征参数对气泡弹跳过程的动态特性以及最终平衡形态的影响.发现随着泡在撞击壁面之前上浮距离增大,气泡回弹距离和弹跳周期增加,但是当上浮距离增加到一定程度后将不会影响气泡的弹跳特性;表面张力是影响气泡弹跳特性的重要因素,气泡的弹跳周期随其增大逐渐减小,但回弹距离却呈现先增后减的规律;最后,影响气泡最终平衡形态的主要因素是气泡的浮力参数与韦伯数.     

Cavitation Bubble Collapse near a Curved Wall by the Multiple-Relaxation-Time Shan-Chen Lattice Boltzmann Model

The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received inereasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time Shan-Chen model is built to study the cavitation bubble collapse. Using this model, the cavitation phenomena induced by density perturbation are simulated to obtain the coexistence densities at certain temperature and to demonstrate the Young-Laplace equation. Then, the cavitation bubble collapse near a curved rigid wall and the consequent high-speed jet towards the wall are simulated. Moreover, the influences of initial pressure difference and bubble-wall distance on the cavitation bubble collapse are investigated.     

Thermodynamic effect of single bubble near a rigid wall

The objective of this paper is to numerically investigate the thermodynamic effect during bubble collapse near a rigid boundary. A compressible fluid model is introduced to accurately capture the transient process of bubble shapes and temperature, as well as corresponding pressure, and velocity. The accuracy of the numerical model is verified by the experimental data of bubble shapes, and Keller-Kolodner equation as well as its thermodynamic equation. The results show that a bubble near the rigid boundary presents high-speed jet in collapse stage and counter jet in rebound stage, respectively. In the collapse stage, the bubble margin will shrink rapidly and do the positive work on the compressible vapor inside the bubble, then a significant amount of heat will be generated, and finally the generation of high-speed jet drives the low-temperature liquid outside the bubble to occupy the position of high-temperature vapor inside the bubble. In the rebound stage, the counter jet moving away from the rigid boundary takes part of heat away from the sub-bubble, which avoids the external work of the expansion of the sub-bubble and the temperature reduction caused by the dissipation effect of the vortex structure. In addition, the initial standoff has a significant effect on the thermodynamics of bubble oscillation. The temperature keeps increasing with the increase of the initial standoff in the collapse stage, while it shows a downward trend with the increase of the initial standoff in the rebound stage. That’s because the high-speed jet and counter jet of bubble gradually disappear when the initial standoff increases, which is the important reason for the opposite evolution trend of temperature in collapse and rebound stage.     

Experimental studies of bubble dynamics inside a corner

The dynamics of a bubble near a corner formed by two flat rigid boundaries (walls), is studied experimentally using a spark-generated bubble. The expansion, collapse, rebound, re-collapse and migration of the bubble, along with jetting and protrusion, are captured using a high-speed camera. Our experimental observations reveal the behaviour of the bubble in terms of the corner angle and the dimensionless standoff distances to the near and far walls in terms of the maximum bubble radius. The bubble remains approximately spherical during expansion except for its surface becoming flattened when in close proximity to a wall. When a bubble is initiated at the bisector of the two walls, the bubble becomes oblate along the bisector during the late stages of collapse. A jet forms towards the end of collapse, pointing to the corner. The closer the bubble to the two walls, the more oblate along the bisector the bubble becomes, and the wider the jet. A bubble initiated near one of the two walls is mainly influenced by the nearer wall. The jet formed is pointing to the near wall but inclined towards the corner. After the jet penetrates through the bubble surface, the bubble becomes a bubble ring, and a bubble protrusion forms following the jet. The bubble ring collapses and subsequently disappears, while the protrusion firstly expands, and then collapses and migrates to the corner.     

Dynamics of a cavitation bubble near a solid wall

The cavitation bubble dynamics, the variation of pressure and velocity fields of the surrounding liquid in the process of the bubble axisymmetric compression near a planar solid wall are considered. It is assumed that the liquid is at rest at the initial moment of time, and the bubble has a spheroidal shape. The liquid is assumed inviscid and incompressible, its motion being potential. The bubble surface deformation and the liquid velocity on the surface are computed by the Euler scheme using the boundary element method until the moment of the collision of some parts of the bubble surface with one another. The influence of the distance of the bubble from the wall and its initial nonsphericity on the liquid pressure and velocity fields, the bubble shape, and the pressure inside the bubble at the end of the time interval under consideration are studied. The maximum pressure in liquid is shown to realize at the bottom of the cumulative jet arising at the bubble collapse with direction to the wall. In the upper part of this jet, the velocity and pressure are practically constant, and the pressure in the jet is approximately equal to the pressure in the bubble.     

Influence of rigid wall on the nonlinear pulsation of nearby bubble

This paper mainly focuses on the nonlinear pulsation of a bubble near the rigid wall. Dynamics of near-wall bubble and free bubble are discussed and compared in details. Investigation reveals as the driving acoustic pressure amplitude increases, nonlinear pulsation of bubble becomes intense gradually. Besides, decreasing the viscosity of host liquid is advantageous for the nonlinear pulsation of bubble. Bifurcation diagrams of bubble radius show acoustic reflection of the rigid wall makes the initial bifurcation appear at low driving acoustic amplitude and on bubble with small ambient radius, and makes the bifurcation still exist for bubble in high-viscosity liquids. That indicates the rigid wall will produce enhancement on the nonlinearity of nearby bubble. As the bubble approaches the wall, the enhancement becomes strong. Moreover, research on the influence of driving frequency shows the rigid wall makes the frequency band corresponding to chaos around the resonant frequency of free bubble shift downward.     

不同沙粒底面下气泡脉动特性实验研究

通过高速摄影系统对电火花气泡与不同沙粒底面间的相互作用进行了实验研究,并改变气泡与沙粒底面之间的距离.实验结果表明:气泡在与沙粒底面的相互作用中会产生两种明显不同的现象,即形成与近刚性壁面类似的气泡射流以及“蘑菇状”气泡,“蘑菇状”气泡撕裂形成两个气泡,随后产生两个反方向的沿轴线方向的射流.沙粒底面边界具有刚性与弹性两种特征.另外,随着气泡与沙粒底面之间的距离d的增大,气泡脉动周期先增大然后减小,存在气泡脉动周期峰值.对于不同的沙粒底面边界,出现气泡脉动周期峰值的距离d随着沙底粒径的增大而越小.