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

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

气泡夹带对壁面油膜流动特性的影响

油-气润滑系统工作过程中,润滑油膜受微油滴冲击和压缩空气扰动影响易形成气泡夹带现象,气泡夹带行为将对壁面润滑油膜层的形成及流动过程产生重要影响。基于VOF数值模拟方法,对含气泡油膜沿倾斜壁面的流动行为进行研究,考察了气泡的存在对油膜形态和流动速度的影响规律,以及气泡破裂阶段空腔邻域内流体压力变化特性。研究表明,油膜夹带气泡的形变和迁移诱发气泡周围微流场的速度扰动现象,导致气液界面处产生非均匀速度梯度分布,进而引发油膜表面的形态波动。气泡发生破裂时,油膜空穴部位发生明显的正负压力波动现象,气泡附近壁面将承受一定的交变载荷作用。     

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

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

液滴撞击超疏水壁面反弹及破碎行为研究

液滴撞击壁面时,壁面亲水性对液滴撞击壁面后的变化历程具有重要的影响。利用相界面追踪的复合Level Set-VOF方法对液滴撞击超疏水壁面的运动进行了研究。研究结果表明,撞击速度较小时,液滴撞壁后发生反弹;撞击速度较大时,液滴撞壁后会发生破碎现象;初始粒径的增大和表面张力的减小,有利于液滴撞壁后产生铺展破碎现象;撞击角度对撞壁后的液滴行为具有较大的影响。通过数值模拟,给出了一定条件下液滴垂直及倾斜撞击超疏水壁面反弹及破碎的临界条件。     

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

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

MOVEMENT AND SPLASHING OF A DROPLET IMPACTING ON A WET WALL

壁面液体层的存在对液滴撞击壁面的运动具有重要的影响。采用气液两相流动相界面追踪的水平集和流体体积复合方法和壁面润湿模型，实现了液滴撞击湿润壁面运动的数值求解；在此基础上，开展了液滴撞击湿润壁面运动的研究。研究结果表明：液滴以不同速度撞击湿润壁面时，会呈现出黏附铺展、波动运动、皇冠几何体运动以及飞溅运动等几种不同的运动形态，液滴撞击湿润壁面后的压力分布是不同运动形态形成的主要原因；飞溅运动是一定条件下皇冠几何体运动的一种特殊形态，液滴从皇冠几何体侧壁顶端的飞溅分离满足毛细破碎理论；撞击速度对分离液滴的运动方向影响较小，而对壁面液体层厚度的影响则较大；撞击速度和壁面液体层厚度对分离液滴形态、飞溅分离位置、飞溅速度以及飞溅发生时刻等都具有一定的影响。     

气泡碰撞亲疏水曲壁的行为特性研究

气泡碰撞固壁行为和影响因素的研究一直以来都是科学界关注的重点之一, 其在矿物浮选、气膜减阻等工业领域中的应用也极具科研价值. 论文聚焦曲壁对于气泡撞击行为特性的影响研究. 采用高速摄像技术记录气泡碰撞不同曲率半径下亲疏水曲壁的撞击过程, 分析了曲壁润湿性、曲率半径对气泡碰撞固体曲壁的影响规律. 结果表明, 气泡碰撞亲水曲壁时会发生多次弹跳直至离开曲壁; 曲率半径越大, 弹跳次数越少, 且第一次反弹的最远距离越近, 再次发生碰壁时的速度越小. 而碰撞疏水曲壁时会出现碰撞?滑移?附着的现象, 此外针对液膜挤压破裂的现象, 建立理论模型推导出液膜诱导时间的预测公式, 其主要与液膜厚度、液膜临界破裂厚度和液膜被压缩速度有关, 预测误差小于5.0%.      

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

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

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

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

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

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

Low-Reynolds-number droplet motion in a square microfluidic channel

In this study, we investigate computationally the low-Reynolds-number droplet motion in a square micro-channel, a problem frequently encountered in microfluidic devices, enhanced oil recovery and coating processes. The droplet deformation and motion are determined via a three-dimensional spectral boundary element method for wall-bounded flows. The effects of the flow rate, viscosity ratio and droplet size on the interfacial dynamics are identified for droplets smaller and larger than the channel size and for a wide range of viscosity ratio. Owing to the stronger hydrodynamic forces in the thin lubrication film between the droplet interface and the solid walls, large droplets exhibit larger deformation and smaller velocity. Under the same average velocity, a droplet in a channel shows a significantly smaller deformation and higher velocity than in a cylindrical tube with the same size, owing to the existence of the corners’ area in the channel which permits flow of the surrounding fluid. A suitable periodic boundary implementation for our spectral element method is developed to study the dynamics of an array of identical droplets moving in the channel. In this case, the droplet deformation and velocity are reduced as their separation decreases; the reduction is influenced by the flow rate, viscosity ratio and more significantly the droplet size.     

Study on heat transfer performance of spray cooling: model and analysis

In consideration of droplet–film impaction, film formation, film motion, bubble boiling (both wall nucleation bubbles and secondary nucleation bubbles), droplet–bubble interaction, bulk air convection and radiation, a model to predict the heat and mass transfer in spray cooling was presented in this paper. The droplet–film impaction was modeled based on an empirical correlation related with droplet Weber number. The film formation, film motion, bubble growth, and bubble motion were modeled based on dynamics fundamentals. The model was validated by the experimental results provided in this paper, and a favorable comparison was demonstrated with a deviation below 10%. The film thickness, film velocity, and non-uniform surface temperature distribution were obtained numerically, and then analyzed. A parameters sensitivity analysis was made to obtain the influence of spray angle, surface heat flux density, and spray flow rate on the surface temperature distribution, respectively. It can be concluded that the heat transfer induced by droplet–film impaction and film-surface convection is dominant in spray cooling under conditions that the heated surface is not superheated. However, the effect of boiling bubbles increases rapidly while the heated surface becomes superheated.     

三维气泡与自由表面相互作用的直接数值模拟

采用VOF中的PLIC界面重构方法数值模拟了三维气泡在液流中上升并与自由表面相互作用的运动.分别考察了不同初始高度,有无来流及有无再生气泡对气泡上升高度、上升速度、压力及与自由表面相互作用等的影响.结果表明:气泡初始位置越低,顶端上升的高度越大,自由面隆起的范围更广.越靠近自由表面,底部射流横向发展越窄,而向上的压力梯度,气泡上升速度,底部射流上升高度越大,反之则反;但如果底部射流均在接近自由表面以前已横向发展充分,则差别不大.气泡外形、上升高度、破裂时间以及上升速度与来流无关.产生再生气泡后,原生气泡与再生气泡相吸,相互加速对方的上升;自由表面抬升的高度增幅较大,范围拓宽,上升速度也大大增加,且再生气泡越多,自由表面隆起的范围越大.     

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.     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

Hydrodynamic and mass transfer characteristics of slug flow in a vertical pipe with and without dispersed small bubbles

In this work, we present a numerical study to investigate the hydrodynamic characteristics of slug flow and the mechanism of slug flow induced CO₂ corrosion with and without dispersed small bubbles. The simulations are performed using the coupled model put forward by the authors in previous paper, which can deal with the multiphase flow with the gas–liquid interfaces of different length scales. A quasi slug flow, where two hypotheses are imposed, is built to approximate real slug flow. In the region ahead of the Taylor bubble and the liquid film region, the presence of dispersed small bubbles has less impacts on velocity field, because there are no non-regular intensive disturbance forces or centrifugal forces breaking the balance of the liquid and the dispersed small bubbles. In the liquid slug region, the strong centrifugal forces generated by the recirculation below the Taylor bubble lead to the effect of heterogeneity, which makes the profile of the radial liquid velocity component sharper with higher volume fraction of dispersed small bubbles. The volume fraction has a maximum value in the range of r/R = 0.5–0.6. Meanwhile, it is usually higher than 0.35, which means that larger dispersed bubbles can be formed by coalescences in this region. These calculated results are in good agreement with experimental results. The wall shear stress and the mass transfer coefficient with dispersed small bubbles are higher than those without dispersed small bubbles due to enhanced fluctuations. For short Taylor bubble length, the average mass transfer coefficient is increased when the gas or liquid superficial velocity is increased. However, there may be an inflection point at low mixture superficial velocities. For the slug with dispersed small bubbles, the product scales still cannot be damaged directly despite higher wall shear stress. In fact, the alternate wall shear stress and the pressure fluctuations perpendicular to the pipe wall with high frequency are the main cause for breaking the product scales.