浮力气泡对水平壁面的回弹动力学特性

黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.      

浮力气泡对水平壁面的回弹动力学特性

黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid, VOF)来重构气泡的形状.从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数Ga和接触速度U_a对气泡回弹动力学特性的影响,分析了气泡碰撞过程中涡结构的变化.用回弹高度H、回弹周期T、长宽比A_r、浮升速度U、轴向位置z和回复系数Cr等参数来表征不同条件时气泡的运动和形状特性.研究结果表明,气泡的回弹运动特性对Ga十分敏感. Ga的增大可加剧气泡形变,促进气泡的回弹运动,增多回弹次数,增大回弹参数(T和H),提升回复系数.然而,接触速度并非决定气泡回弹动力学的控制参数, Ua的改变并不会改变回复系数.     

水中高压脉动气泡与浮体流固耦合特性研究

本文针对水中放电气泡与水面浮体流固耦合作用开展实验和数值研究, 采用边界积分法对气泡运动进行数值模拟, 利用辅助函数法提高非线性流固耦合问题的计算精度, 同时运用双节点法保证气-液-固三相交界线的计算稳定性. 实验中, 采用水下放电技术生成气泡, 使用高速摄影捕捉气泡动力学行为与浮体运动响应. 首先对比数值与实验结果, 二者吻合良好, 验证了数值计算模型的有效性和正确性. 然后通过对气泡与浮体的无量纲距离$\gamma_{s} $ (气泡最大半径为特征长度)进行系统研究发现: (1) $\gamma_{s} $从0.2增大至2时, 气泡在坍塌阶段分别形成了颈缩型环状射流(本文针对水中放电气泡与水面浮体流固耦合作用开展实验和数值研究,采用边界积分法对气泡运动进行数值模拟,利用辅助函数法提高非线性流固耦合问题的计算精度,同时运用双节点法保证气-液-固三相交界线的计算稳定性.实验中,采用水下放电技术生成气泡,使用高速摄影捕捉气泡动力学行为与浮体运动响应.首先对比数值与实验结果,二者吻合良好,验证了数值计算模型的有效性和正确性.然后通过对气泡与浮体的无量纲距离γ_s(气泡最大半径为特征长度)进行系统研究发现:(1)γ_s从0.2增大至2时,气泡在坍塌阶段分别形成了颈缩型环状射流(0.2≤γ_s≤0.3)、接触射流(0.4≤γ_s≤0.6)、非接触射流(0.7≤γ_s≤1)、对射流(1.1≤γ_s≤1.3)和反射流(1.4≤γ_s≤2)等5种典型射流模式;(2)正射流速度随γ_s先增大后减小再增大,并且当0.7≤γ_s≤0.9时,速度可达约1000 m/s;反射流速度随γ_s增大而增大;(3)在本文实验条件下,γ_s1.5时浮体对气泡的Bjerknes吸引力强于自由液面的Bjerknes排斥力导致气泡在坍塌阶段向浮体迁移;当γ_s≥1.5时自由液面对气泡的排斥作用更强,气泡在坍塌阶段远离自由液面.     

含气泡油滴撞击油膜壁面时气泡的变形与破裂

油--气润滑过程中润滑油液滴受高速气流扰动易形成含气泡油滴,微气泡将对油滴撞击壁面时的运动过程以及壁面油膜 层的形成质量产生重要影响. 基于耦合的水平集--体积分数 方法,对含气泡油滴撞击油膜壁面行为进行数值模拟研究, 考察含气泡油滴撞击油膜壁面时气泡的变形运动过程,探讨气泡破裂的动力学机制,分析气泡大小、碰撞速度和液体黏度等因素对含气 泡油滴撞壁过程中气泡变形特征参数的影响规律. 研究表明:含气泡油滴撞击油膜壁面后气泡会发生变形,并破裂形成膜液滴;气泡随同 液滴运动过程中,气泡内外压力和速度梯度变化是使气泡发生破裂的主要诱因. 气泡大小对气泡破裂方式影响较大,气泡较小时发生单 点破裂,而气泡较大时更容易发生多处破裂. 不同大小气泡受力差异较大,气泡大小与破裂发生时刻没有明显相关性. 碰撞速度和液体 黏度对气泡的变形、破裂和破裂发生时刻都具有一定的影响. 碰撞速度越大,油滴动能越大,更容易产生气泡变形和破裂现象. 液体黏 度增大,在油滴撞壁运动前期促进气泡变形,而在运动后期可以阻延气泡破裂行为发生.      

垂直壁面附近上升单气泡的弹跳动力学研究

采用高速摄影技术结合阴影法,对静止水中垂直壁面附近上升单气泡运动进行实验研究,对比气泡尺度及气泡喷嘴与壁面之间的初始无量纲距离(S~*)对气泡上升运动特性的影响,分析气泡与壁面碰撞前后,壁面效应与气泡动力学机制及能量变化规律.结果表明,对于雷诺数Re≈580～1100,无量纲距离S~*2～3时,气泡与壁面碰撞且气泡轨迹由无约束条件下的三维螺旋转变成二维之字形周期运动;当S~* 2～3时,壁面效应减弱,有壁面约束的气泡运动与无约束气泡运动特性趋于一致.气泡与壁面碰撞前后,壁面效应导致横向速度峰值下降为原峰值的70%,垂直速度下降50%;气泡与壁面碰撞前,通过气泡中心与壁面距离(x/R)和修正的斯托克斯数相关式可预测垂直速度的变化规律.上升气泡与壁面碰撞过程中,气泡表面变形能量单向传输给气泡横向动能,使得可变形气泡能够保持相对恒定的弹跳运动.提出了气泡在与壁面反复弹跳时的平均阻力系数的预测模型,能够很好地描述实验数据反映出的对雷诺数Re、韦伯数We和奥特沃斯数Eo等各无量纲参数的标度规律.     

垂直壁面附近上升单气泡的弹跳动力学研究

采用高速摄影技术结合阴影法,对静止水中垂直壁面附近上升单气泡运动进行实验研究,对比气泡尺度及气泡喷嘴与壁面之间的初始无量纲距离 ($S^{\ast}$)对气泡上升运动特性的影响,分析气泡与壁面碰撞前后,壁面效应与气泡动力学机制及能量变化规律.结果表明,对于雷诺数$Re \approx 580 \sim 1100$,无量纲距离$S^{\ast } <2 \sim3$时,气泡与壁面碰撞且气泡轨迹由无约束条件下的三维螺旋转变成二维之字形周期运动;当$S^{\ast } >2 \sim3$时,壁面效应减弱,有壁面约束的气泡运动与无约束气泡运动特性趋于一致.气泡与壁面碰撞前后,壁面效应导致横向速度峰值下降为原峰值的70%,垂直速度下降50%;气泡与壁面碰撞前,通过气泡中心与壁面距离($x/R$)和修正的斯托克斯数相关式可预测垂直速度的变化规律.上升气泡与壁面碰撞过程中,气泡表面变形能量单向传输给气泡横向动能,使得可变形气泡能够保持相对恒定的弹跳运动.提出了气泡在与壁面反复弹跳时的平均阻力系数的预测模型,能够很好地描述实验数据反映出的对雷诺数${Re}$、韦伯数${We}$和奥特沃斯数${Eo}$等各无量纲参数的标度规律.      

单个上升气泡与异质液膜相互作用的数值模拟

为了探讨在流体区域内气泡上升然后与异质液膜接触的过程中,气泡速度、形状随时间的变化,以及区域内液膜的运动特性。本文采用FTM方法模拟研究了在2D×D的空间范围内单个气泡从静止开始上升到与异质液膜接触过程中的运动特性。研究结果表明:气泡上升到与液膜一定距离后,会使液膜中间凸起变薄,并将液膜挤向两侧,气泡的竖直速度先增大后减小;液膜厚度不同时,气泡上升初期相同时刻的纵横比?基本相同;气泡上升后期相同时刻的纵横比随着液膜厚度的增厚而增大;不同表面张力下,相同时刻气泡纵横比随着表面张力的减小而减小;表面张力越小,气泡顶部和底端的压力差会越小,上升速度越小。     

壁面处气泡在静止流场和高速水流中溃灭过程的计算仿真

通过数值仿真计算,模拟近壁面以及附壁面气泡在静止流场和高速水流中的溃灭过程,研究气蚀作用机理.结果表明:气泡与壁面的距离和水流的速度影响其溃灭时间;附壁面气泡在高速水流中完全溃灭的时间最短,而在静止流场中最长,远离壁面将增加气泡的不稳定性;当气泡距离壁面一定距离溃灭时,射流不能直接作用于壁面,壁面承受冲击波的最大压力远小于气泡溃灭中心的压力;当气泡溃灭中心在壁面时,射流直接作用于壁面产生微小而严重的点破坏,而冲击波则使材料产生交变应力,造成环形破坏;当气泡在高速水流中溃灭时将产生逆流斜向射流,这可能是水力机械过流部件产生鱼鳞坑和波纹状破坏的主要原因.     

含气泡油滴撞击油膜壁面时气泡的变形与破裂

油-气润滑过程中润滑油液滴受高速气流扰动易形成含气泡油滴,微气泡将对油滴撞击壁面时的运动过程以及壁面油膜层的形成质量产生重要影响.基于耦合的水平集-体积分数方法,对含气泡油滴撞击油膜壁面行为进行数值模拟研究,考察含气泡油滴撞击油膜壁面时气泡的变形运动过程,探讨气泡破裂的动力学机制,分析气泡大小、碰撞速度和液体黏度等因素对含气泡油滴撞壁过程中气泡变形特征参数的影响规律.研究表明:含气泡油滴撞击油膜壁面后气泡会发生变形,并破裂形成膜液滴;气泡随同液滴运动过程中,气泡内外压力和速度梯度变化是使气泡发生破裂的主要诱因.气泡大小对气泡破裂方式影响较大,气泡较小时发生单点破裂,而气泡较大时更容易发生多处破裂.不同大小气泡受力差异较大,气泡大小与破裂发生时刻没有明显相关性.碰撞速度和液体黏度对气泡的变形、破裂和破裂发生时刻都具有一定的影响.碰撞速度越大,油滴动能越大,更容易产生气泡变形和破裂现象.液体黏度增大,在油滴撞壁运动前期促进气泡变形,而在运动后期可以阻延气泡破裂行为发生.     

电场作用下气泡上升行为特性的数值计算研究

利用电场控制气泡形态及运动,强化气液相间传热传质是电流体动力学的重要研究内容之一. 然而目前多数研究集中在非电场下的气泡动力学上,对于电场下的气泡行为特性及电场的作用机制仍需开展深入研究. 本研究对电场作用下单个气泡在流体中上升过程的动力学行为进行了数值模拟研究. 在建立二维模型的基础上求解电场方程与Navier-Stokes方程,并采用水平集方法捕捉了上升气泡的位置及形状. 模拟结果的准确性与有效性通过与前人实验和数值结果进行对比得到了验证. 通过改变雷诺数、邦德数和电邦德数等不同参数研究了电场下液体黏度、表面张力和电场力对气泡运动变形的影响. 计算结果表明,电场对气泡的动态特性有显著影响. 非电场情况下液体黏度和表面张力较大时气泡基本维持球状,反之气泡发生变形并逐步达到稳定状态. 此外,电场作用使气泡在初始上升阶段发生剧烈形变,随着不断上升,气泡形变程度不断减小,且气泡的上升速度和长径比均出现振荡. 垂直电场使气泡的上升速度有较大的提高,且随着电邦德数的增大,难以达到相对稳定的状态.      

Numerical study of bubble rise and interaction in a viscous liquid

In this article, the flow instabilities during the rise of a single bubble in a narrow vertical tube are studied using a transient two-dimensional/axisymmetric model. To predict the shape of the bubble deformation, the Navier-Stokes equations in addition to an advection equation for liquid volume fraction are solved. A modified volume-of-fluid technique based on Youngs' algorithm is used to track the bubble deformation. To validate the model, the results of simulations for terminal rise velocity and bubble shape are compared with those of the experiments. The effect of different parameters such as initial bubble radius, channel height, liquid viscosity and surface tension on the shape and rise velocity of the bubble is investigated.     

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.     

Co-current flow effects on a rising Taylor bubble

The effects of co-current flows on a rising Taylor bubble are systematically investigated by a front tracking method coupled with a finite difference scheme based on a projection approach. Both the upward (the co-current flows the same direction as the buoyancy force) and the downward (the co-current moves in the opposite direction of the buoyancy force) co-currents are examined. It is found that the upward co-current tends to elongate the bubble, while the downward co-current makes the bubble fatter and shorter. For large N_f (the inverse viscosity number), the upward co-current also elongates the skirted tail and makes the tail oscillate, while the downward co-current shortens the tail and even changes a dimpled bottom to a round shape. The upward co-current promotes the separation at the tail, while the downward co-current suppresses the separation. The terminal velocity of the Taylor bubble rising in a moving flow is a linear combination of the mean velocity (U_C) of the co-current and the terminal velocity (U₀) of the bubble rising in the stagnant liquid, and the constant is around 2 which agrees with the literature. The wake length is linearly proportional to the velocity ratio (U_C/U₀). The co-currents affect the distribution of the wall shear stresses near the bubble, but not the maximum.     

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

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

海洋大气间气泡的气体输运

研究球形小气泡在理想流体的波浪场中的气体扩散过程,把小雷诺数下均匀来流绕流球形气泡的气体交换结果与气泡运动方程耦合在一起进行求解．讨论了溶解于水中的气体浓度、波浪、气泡半径、气泡初始深度对单个气泡气体扩散量的影响．由于气泡云对气体的输运,溶解于水中的气体可出现过饱和状态．对１０ｍ／ｓ风速下气泡云的气体输运量进行了计算,得到水中Ｏ２的过饱和度可达１．８９％～３９２％,与实际观测值一致．     

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     

The motion of a small bubble in waves

The motion of a single spherical small bubble due to buoyancy in the ideal fluid with waves is investigated theoretically and experimentally in this article. Assuming that the bubble has no effect on the wave field, equations of a bubble motion are obtained and solved. It is found that the nonlinear effect increases with the increase of the bubble radius and the rising time. The rising time and the motion orbit are given by calculations and experiments. When the radius of a bubble is smaller than 0.5mm and the distance from the free surface is greater than the wave height, the results of the present theory are in close agreement with measurements. The project supported by the National Natural Science Foundation of China     

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.