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

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

本站小编 Free考研考试/2022-01-01

王悦柔, 王军锋^,²⁾, 刘海龙江苏大学能源与动力工程学院,江苏镇江 212013

NUMERICAL SIMULATION ON BUBBLE RINSING BEHAVIORS UNDER ELECTRIC FIELD ¹⁾

Wang Yuerou, Wang Junfeng^,²⁾, Liu HailongSchool of Energy and Power Engineering,Jiangsu University, Zhenjiang 212013,Jinagsu,China

通讯作者: 2) 王军锋,教授,主要研究方向：荷电多相流理论及工程应用研究. E-mail:wangjunfeng@ujs.edu.cn

收稿日期:2019-07-18接受日期:2019-12-9网络出版日期:2020-01-18

基金资助:

1) 国家自然科学基金资助项目.51761145011
国家自然科学基金资助项目.51876086）

Received:2019-07-18Accepted:2019-12-9Online:2020-01-18
作者简介 About authors

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

Abstract
Employing an electric field to control bubble morphology and motion, and enhance heat and mass transfer between gas and liquid are one of the important research contents of electrohydrodynamics (EHD). However, most of the current researches focus on the dynamics of bubbles under non-electric fields. Further research is needed on the behaviors of bubbles and the mechanism of electric fields under electric fields. In the present study, the dynamic behaviors of a single bubble rising in a fluid under an external electric field are numerically simulated. Based on the model established on two-dimensional, the equation of electric field and Navier-Stokes are solved, and the level set method is used to accurately capture the position and morphology of the rising bubble. The accuracy and validity of the present simulation results are verified by comparing with previous experiments and numerical results. The effects of liquid viscosity, surface tension and electric field force on bubble motion and deformation under electric field are examined by employing the Reynolds ($Re$), Bond ($Bo$) and electrical bond ($Bo_{\rm e}$) numbers. The calculation results show that the electric field has a significant influence on the dynamic characteristics of the bubble. In the case of a non-electric field, the bubble basically maintains a spherical shape when the viscosity of the liquid and the surface tension are large. Instead, the bubbles deform and gradually reach a steady state. Additionally, the bubble is strongly deformed at the initial rising stage by attributing to the electric field. With the bubble further rising, the deformation weakened, and both of the bubble velocity and aspect ratio are oscillated. The effect of vertical electric field causes a greater increase in the velocity of rising bubble. As the number of $Bo_{\rm e}$ increases, the bubble has an intense oscillation, making it more difficult to achieve a relatively stable state.
Keywords：electric field;bubble;deformation;numerical simulation;dimensionless number

PDF (8884KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex 收藏本文
本文引用格式
王悦柔, 王军锋, 刘海龙. 电场作用下气泡上升行为特性的数值计算研究 ¹⁾. 力学学报[J], 2020, 52(1): 31-39 DOI:10.6052/0459-1879-19-193
Wang Yuerou, Wang Junfeng, Liu Hailong. NUMERICAL SIMULATION ON BUBBLE RINSING BEHAVIORS UNDER ELECTRIC FIELD ¹⁾. Chinese Journal of Theoretical and Applied Mechanics[J], 2020, 52(1): 31-39 DOI:10.6052/0459-1879-19-193

引言

气液多相流不仅存在于自然现象中,在工业生产中也有广泛的应用^[1-7],如核池沸腾、液体分离、废水处理和化学反应等. 对流体中上升气泡的运动特性及规律的认识对相关工业装置的设计和运行十分重要.为此研究人员已经对黏性流体中的气泡动力学^[8-20] 进行了的大量研究. Clift等^[8]最初对气泡上升过程的研究进行了回顾,表明气泡的运动可以用无量纲数描述,即雷诺数$Re$、莫顿数${Mo}$和欧特沃斯数$Eo$,后者也可以称为邦德数$Bo$.并且他认为在低$Re$和$Bo$下气泡能维持球状,随着$Re$和$Bo$的增加会产生不同的气泡形态,如椭圆帽状、球帽状、椭球状以及摆动形状.之后的研究表明^[9-11],在$Re$和$Bo$足够大的情况下,气泡可能会变为环形.Tripathi等^[12] 在$Mo$较低的情况下发现了一种新的气泡破裂模式,气泡变为球状并产生卫星气泡.王焕然等^[13] 发现在$Mo$较低,$Re$较高的情况下,气泡形状表现出周期振荡变化.此外,其他参数也可能影响气泡上升机制. Chen等^[14] 研究了雷诺数、邦德数、密度和黏度比对上升气泡的影响.结果表明,$Re$和$Bo$以及密度比的增大会增强气泡下方液体射流的影响,从而导致气泡中心穿孔.Amaya-Bower^[15]分析表明,黏度比对气泡的终端速度和形状的影响更大,随着黏度比的降低,剪切应力增大,气泡变形程度和终端速度减小.Alizadeh等^[16] 在二维和三维数值研究的基础上验证了表面张力对气泡形变行为的影响符合拉普拉斯定律.张洋等^[17] 通过数值研究表明当域约束比$Cr \geqslant 8$时可忽略壁面对裙带气泡特性的影响.

目前,多数研究集中在对非电场作用下的气泡动力学特性上.将电场应用于不同介质的多相体系中,能够有效地控制气泡的运动及变形^[21-27],使气泡获得较好的流动状态和上升特性.早期,Taylor^[21]和Garton等^[22]对这一问题进行了经典研究,受电场作用的液滴和气泡会沿电场方向发生变形,在电场方向上呈现出一个扩展的椭球体形状.陈凤等^[23-24]在蠕动流近似的基础上分析计算了具有不同介电常数工质的气泡在外电场力作用下的受力状况以及气泡内外的速度场分布.计算结果表明电场力的水平分力使气泡沿场强方向伸长, 使得气泡内部的流体运动加剧.并且液体和气体相对介电常数之比越大的工质所受的电场力越大. Yang等^[25]的研究表明电场形式对气泡运动也有影响,垂直电场加速了气泡的上升运动,而水平电场阻碍了上升运动.Rahmat等^[26]研究了电毛细管数与域限制比对气泡（特别是环形气泡）上升的影响,结果表明,当电力足够强时,气泡就会被穿透.对于较小的域限制比率,需要更大的电毛细管数来穿透气泡. Andalib等^[27]通过实验研究了不同电场形式下单个粗气泡形态及运动路径的变化过程.然而实验条件的限制使得多数情况下的气泡行为特性难以通过实验捕捉,目前的多数研究也是依靠数值计算手段实现.且由于静电学与流体动力相互作用的复杂性,相关物性参数对受电场力作用的气泡形态及上升速度的影响的研究尚不充分,对于电场作用下气泡的运动特性及电场的作用机制仍需开展深入系统的研究.

目前已有多种数值方法被开发应用于气液两相流体动力学研究中,包括边界积分法、前跟踪法、流体体积（volume of fluid,VOF）法、水平集（level set,LS）法和光滑粒子流体力学（smoothed particle hydrodynamics,SPH）法等.边界积分法^[28] 和前跟踪法^[11,29] 是精确计算体积力的显式方法,但实现起来非常复杂,且需要消耗大量的计算时间.VOF法、LS法和相场法是一种隐式方法,可以对界面的破裂、聚并等拓扑变化进行常规处理.与相场法^[25]相比,VOF法和LS法在电流体动力学中得到了广泛的应用. VOF方法^[14,26,30]很好地解释了界面的拓扑变化,满足质量守恒,但不能准确计算界面附近的参数. LS方法^[32-34]在保持良好质量守恒的前提下,能准确捕捉两相流体界面.另外,一种利用具有密度和黏度等物理性质的质点来描述流体的SPH方法也可被用于多相流数值研究中.这种方法非常适用于模拟大变形、移动边界和多相混合的动态现象.然而从已有文献可知,目前SPH法^[35]还没有在电流体动力学中得到广泛的验证与应用.

本文借助数值手段对外加电场作用下单个气泡在流体中上升的动力行为进行了模拟研究. 利用水平集方法,通过求解电场方程与Navier-Stokes方程,捕获气液两相界面.该方法能够精确计算界面法向和曲率,平滑界面周围的物理性质,准确地捕捉气液两相流的拓扑变化.之后从网格无关性及与前人研究对比两方面验证了该模型的有效性,通过雷诺数、邦德数和电邦德数等不同参数研究了液体黏度、表面张力和电场力对气泡运动及变形的影响,并给出了不同情况下的气泡形态图.

1 数学描述

本文中的气液两相流体被认为是具有不同介电常数的介质,在没有外加电场的情况下,在流体本体区或界面处没有自由电荷. 用于本文的二维模型如图1所示,初始半径为$R$的气泡位于高$H= 16R$,宽$W= 8R$的矩形计算域中,从气泡中心距离底边$h_{1}= 3 R$处释放,气泡上升至距离顶边$h_{1}$处停止计算.这里选定$R =1$进行模拟,模拟中所用的单位为mm. 计算域的四周边界上应用无滑移壁面边界条件. 在竖直方向上施加电场强度$E_\infty = U_{\rm t} /H$的均匀电场. 电势在顶部边界上为$U_{\rm t}$,在底部边界处为0.

图1

新窗口打开|下载原图ZIP|生成PPT
图1均匀电场中气泡上升的二维模型示意图

Fig. 1Schematic diagram of a two-dimensional model of a rising bubble under a uniform electric field

1.1 水平集方法

水平集方法通过跟踪水平集函数$\phi $的等值线来确定流体界面. 水平集$\phi = 0.5$决定界面的位置. 控制$\phi $的传递和重新初始化的方程为

(1) $ \dfrac{\partial \phi }{\partial t} + {\pmb u} \cdot \nabla \phi = \gamma \nabla \cdot \left( {\varepsilon \nabla \phi - \phi \left( {1 - \phi } \right)\dfrac{\nabla \phi }{\left| {\nabla \phi } \right|}} \right)$

式中,${\pmb u}$为流体速度, $\gamma $和$\varepsilon$为重新初始化参数. 参数$\varepsilon$确定界面周围的层厚度,其中$\phi $从0变为1. $\gamma$参数决定重新初始化的数量.

研究分为相位初始化和时间相关性研究两个阶段. 相位初始化部分求解初始接口到计算域节点$D_{\rm wi}$的距离. 时间相关步骤根据以下表达式使用水平集函数的初始条件

(2) $ \phi = \dfrac{1}{1 + {\rm e}^{\pm D_{\rm wi} / \varepsilon }}$

式中,正号用于界面内部初始域,负号用于界面外部初始域.

由于水平集函数是一个平滑阶跃函数,因此通过以下公式确定全局密度和动力黏度

(3) $ \rho = \rho _{\rm b} + \left( {\rho _{\rm l} - \rho _{\rm b} } \right)\phi$

(4) $ \mu = \mu _{\rm b} + \left( {\mu _{\rm l} -\mu _{\rm b} } \right)\phi$

式中,$\rho_{\rm b}$,$\rho_{\rm l}$,$\mu_{\rm b}$和$\mu_{\rm l}$分别表示气体和周围流体的恒定密度和黏度.

1.2 流体流动方程

对于电场作用下的流体流动,流体流动受惯性力、黏性力、电场力等控制. 在不可混相不可压缩牛顿流体的假设下,将N-S方程表示为

(5) $ \nabla \cdot {\pmb u} = 0$

(6) $ \rho \left( {\dfrac{\partial {\pmb u}}{\partial t} + {\pmb u} \cdot \nabla {\pmb u}} \right) =-\nabla p + \nabla \cdot \mu \left( {\nabla {\pmb u} + \nabla {\pmb u}^{\rm T}} \right) + {\pmb F}_{\rm st} + \rho {\pmb g} +{\pmb F}_{\rm es}$

式中,$\rho $表示密度, ${\pmb u}$是速度,${\pmb g}$为重力加速度,$t$为${\pmb F}_{\rm st}$表示表面张力,${\pmb F}_{\rm es}$表示外加体积力,这里即代表电场力.

表面张力由下式定义

(7) $ {\pmb F}_{\rm st} = \nabla \cdot \left\{ {\sigma \left[ {{\pmb I} + \left( { - {\pmb n}{\pmb n}^{\rm T}} \right)} \right]\delta } \right\}$

式中, $\sigma $为表面张力系数,${\pmb I}$为单位矩阵,${\pmb n}$为界面单位法向,$\delta $为狄拉克$\delta $函数,该函数仅在流体界面为非零值. 界面法向的计算公式如下

(8) $ {\pmb n} = \dfrac{\nabla \phi }{\left| {\nabla \phi } \right|}$

水平集参数$\phi $还用于通过由下式定义的平滑函数来近似$\delta $函数

(9) $ \delta = 6\left| {\phi \left( {1 -\phi } \right)} \right|\left| {\nabla \phi } \right| (9) $

1.3 电场方程

在电流体力学中,由于磁感应效应是可以忽略的,所以电场强度${\pmb E}$被假定为无旋的 ($\nabla \times {\pmb E} = {\bf 0}$). 对于理想的介质材料,电场的控制方程可以写为

(10) $ \nabla \cdot \left( {\varepsilon _0 \varepsilon _{\rm r} {\pmb E}} \right) = {\bf 0}$

式中,$\varepsilon _0 $为真空的介电常数,$\varepsilon _r $为相对介电常数. 在没有任何时变磁场的情况下,电场的卷曲为零.电场可以用电势梯度表示

(11) $ {\pmb E} = -\nabla V $

因此,式(10)可以根据电势$V$重写

(12) $ \nabla \cdot \left( {\varepsilon _0 \varepsilon _{\rm r} \nabla V} \right) = {\bf 0}$

作用在两相流体中的外加体积力${\pmb F}_{\rm es}$即为电场力,由麦克斯韦应力张量的形式给出

(13)$ {\pmb F}_{\rm es} = \nabla \cdot {\pmb T } = \nabla \cdot \Big [{\pmb E}{\pmb D}^{\rm T} - \dfrac{1}{2}({\pmb E} \cdot {\pmb D}) {\pmb I} \Big ]$

式中,${\pmb D}$为电位移场,由下式给出

(14) $ {\pmb D} = \varepsilon _0 \varepsilon _{\rm r} {\pmb E} $

1.4 无量纲参数定义

这里定义了一些特征参数,包括气泡的特征速度$U = u / \sqrt {gR} $,特征时间$T = t / \sqrt {R / g} $以及终端雷诺数$Re^{\ast } = \rho _{\rm l} uR / \mu _{\rm l}$. 在均匀电场中气泡的上升可以用无量纲参数表征,包括雷诺数 (${Re}$)、邦德数 (${Bo}$)、电邦德数 ($Bo_{\rm e}$),以及连续相和分散相的性质比值

(15) $ \left.\begin{array}{l} \rho^ {\ast } = \dfrac{\rho_{\rm b} }{\rho_{\rm l} } , \ \ \mu^{\ast } = \dfrac{\mu_{\rm b} }{\mu_{\rm l} } , \ \ \varepsilon^{\ast } = \dfrac{\varepsilon_{\rm b} }{\varepsilon_{\rm l} } \\ Re = \dfrac{\rho_{\rm l} \sqrt {gR^3} }{\mu_{\rm l} } , \ \ Bo = \dfrac{\rho_{\rm l} gR^2}{\gamma } , \ \ Bo_{\rm e} = \dfrac{\varepsilon_{\rm l} E_\infty ^2 R}{\gamma } \end{array}\right\} $

式中,$R$是气泡半径,$ E_\infty$是未扰动电场强度. 下标b和l分别表示气体和周围液体.

2 模拟结果

2.1 模型验证

为了研究结果的网格依赖性,采用不同的网格数模拟$\rho^{\ast}=0.001$, $\mu ^{\ast }=0.01$, $\varepsilon^{\ast}=0.04$,${Re}=20$,$Bo =0.3$,$Bo_{\rm e}=1$条件下的上升气泡.图2比较了不同网格数下$T =5.74$时的气泡的形状及位置.这里用$\varDelta$表示最大网格尺寸,可以发现,随着$R/\varDelta$的增大,气泡形状及位置的差异逐渐减小,$R/\varDelta =70$之后,气泡形状的差异较小,误差在5%以内,考虑到计算精度和计算机资源成本,这里选用$R/\varDelta=70$的网格尺寸条件进行接下来的模拟.如果区域的宽度远大于气泡的半径,则可以忽略区域的边界效应. 根据研究表明^[36], 当域高和域宽分别小于气泡直径的8倍和3倍时,边界效应对气泡的影响较大,这里设定域高$H= 16 R$,宽$W= 8 R$. 此外,当前2D模拟结果与文献^[19] 的实验结果及对应的3D模拟结果进行了比较. 气泡最终形态与终端雷诺数${Re}^{\ast }$的对比结果分别如表1和图3所示.

图2

新窗口打开|下载原图ZIP|生成PPT
图2不同网格大小下气泡形态及位置

Fig. 2Bubble geometry and position with different grid sizes

Table 1
表 1
表 1当前模拟结果与前人实验结果的对比
Table 1Comparison between current numerical results and previous experimental results

新窗口打开|下载CSV

图3

新窗口打开|下载原图ZIP|生成PPT
图3数值模拟预测的最终气泡形态与前人实验结果^[19]的比较

Fig. 3Comparison of the bubble terminal deformation predicted by numerical simulation with previous experimental results^[19]

可以看出,2D仿真预测的气泡最终形态与前人实验观测结果及3D预测的气泡形态基本一致,且与气泡终端速度相关的终端雷诺数$Re^{\ast}$误差在5%以内,因此2D模型足以用于当前对气泡形态变化及上升速度的数值研究,且比3D模型节省了更多的计算量,验证了该模型的有效性.

2.2 非电场作用下的气泡上升行为

初始气泡形状假设为球形,在本研究中,这一假设将继续成立. 接下来,讨论了流体黏度和表面张力对气泡动力学的孤立影响. 这里采用$\rho^{\ast }=0.001$,$\mu^{\ast }=0.01$条件进行接下来的模拟,通过调整黏度和表面张力系数获取不同${Re}$和${Bo}$,预测的气泡在$h/R =13$处的最终形状如图4所示.

图4

新窗口打开|下载原图ZIP|生成PPT
图4根据不同${Re}$和${Bo}$预测的上升气泡的最终形状

Fig. 4The terminal shape of the rising bubble predicted by different ${Re}$ numbers and ${Bo}$ numbers

一般情况下,当$Re <40$,$Bo <30$时气泡沿直线上升且最终形成沿直线对称的完整的球状或扁椭球状^[12]. 随着${Re}$和${Bo}$的增大,气泡可能会发生破碎与穿透为环状,或者出现路径不稳定的上升的现象,在这种情况下,需要一个三维模型来更好地预测气泡的速度和形状,而不是本文所使用的二维模型. 气泡形状主要受惯性力(即重力)、表面张力和黏性力的控制. 表面张力和黏性力阻止了气泡的变形,并努力使气泡保持球形. 然而,惯性力导致气泡变形. 当惯性力大于黏性力和表面张力时,气泡发生变形. 在初始阶段,气泡的变形和加速上升是由于表面张力,惯性力和黏性力的不平衡. 在达到力平衡后,气泡形状波动微小并稳定上升. 这里通过改变气泡周围流体的黏度及表面张力获得${Re}<40$,${Bo}<30$范围内的不同${Re}$和${Bo}$. 从图4可以发现,在${Re}$和${Bo}$均较小的情况下,气泡更偏近球状,随着${Re}$和${Bo}$增大,气泡最终呈扁长状. 即液体黏度和气液两相间的表面张力越大,气泡在上升过程中克服惯性力,使得气泡维持球形. 随着黏性力及表面张力的减小,气泡不再能维持球状,受惯性力影响发生变形.

2.3 电场强度对气泡形状及速度的影响

在电场作用下,介电液滴在介质流体中的两相界面上受到电场力的作用,会沿电场方向发生拉伸变形^[37]. 定义长径比$Ar = W_{\rm b}/ H_{\rm b} $ 来反映气泡的变形,$W_{\rm b}$和$H_{\rm b}$分别表示气泡的宽度和高度. 通过改变$Bo_{\rm e}$来研究电场强度对气泡运动和变形的影响.可以知道,当${Ar}>1$时气泡发生横向变形,${Ar}$越大,形变程度越大;反之,则发生纵向变形,即沿电场方向拉伸变形,${Ar}$越小,变形程度越大.图 5给出了${Re}=20$,${Bo}=0.3$条件时不同${Bo}_{\rm e}$下气泡长径比${Ar}$随$T$值的变化,气泡初始形状均设为球状. 当${Bo}_{\rme}=0$时,随着$T$值的增大,${Ar}$值逐渐增大,并趋于稳定,${Ar}>1$,气泡沿水平方向变形,最终变形为扁椭球状. 随着${Bo}_{\rm e}$增大,${Ar}$减小,气泡被克服表面张力的电力沿电场力方向拉伸,电场越强,${Ar}$越小,气泡形状拉伸程度越大,且${Ar}$波动变化,较难达到稳定值. 分别选取${Ar}$达到幅值时的不同特征时刻$T$下的气泡形态图,如图6所示.电场力和表面张力之间的相互作用主导了上升的早期阶段,气泡沿电场呈明显拉伸状态,${Ar}$达到一个最小值,之后气泡形状出现竖直和水平方向交替的摇摆不定的振荡形式,且${Bo}_{\rm e}$越大,振荡幅度越大. 随着气泡上升,由于黏滞效应,摇摆运动的振幅随时间而减小,当气泡上升时,惯性力效应进入并与其他机制竞争. 此外可以发现,在施加了电场的情况下,随着气泡垂直上升,不同${Bo}_{\rm e}$下,当前计算工况下的气泡均发生了六次振荡,气泡${Ar}$达到幅值处的特征时间也相差不大,气泡形状发生振荡的时间节点基本一致,不因电场强度的改变而发生过大改变.

图5

新窗口打开|下载原图ZIP|生成PPT
图5气泡长径比${Ar}$随$T$的变化图：${Re}=20$,${Bo}=0.3$

Fig. 5Temporal evolution of the aspect ratio ${Ar}$ of the gas bubble：${Re}=20$,${Bo}=0.3$

图6

新窗口打开|下载原图ZIP|生成PPT
图6不同${Bo}_{\rm e}$下气泡形态演化过程：(a) ${Bo}_{\rm e} =1$,(b) $Bo_{\rm e} =2$,(c) $Bo_{\rm e}=3$,(d) $Bo_{\rm e}=4$

Fig. 6Bubble shape evolution under different ${Bo}_{\rm e}$：(a) ${Bo}_{\rm e} =1$,(b) $Bo_{\rm e} =2$,(c) $Bo_{\rm e}=3$,(d) $Bo_{\rm e}=4$

与此同时,不同$Bo_{\rm e} $下气泡上升特征速度$U$值随$T$值变化的曲线如图7所示. 在没有电场的情况下,上升速度主要受两种力控制：密度差引起的浮力和黏性效应引起的阻力. 然而,流动阻力与气泡的形状密切相关^[11,38]. 其中,当$Ar$较小时,流动阻力较小,上升速度较快,反之亦然. 根据本研究结果可以发现,这种情况也适用于加电情况. 气泡的上升速度是通过气泡质心位置的时间差值得到的. 由图7可见,施加垂直方向的电场导致气泡速度增大. 电场越强时气泡上升得更快,原因主要是由于由电场引起的垂直拉伸导致${Ar}$减小,引起的气泡所经受的阻力减小,气泡速度增大. 由于气泡的不稳定振荡变形带动气泡周围流体的流动,阻力发生变化,上升速度值也随着气泡的不稳定抖动变形出现振荡形式,从图中可见速度$U$发生振荡时间节点与纵横比${Ar}$振荡点几乎一致.

图7

新窗口打开|下载原图ZIP|生成PPT
图7气泡上升速度$U$随$T$的变化图： ${Re}=20$,${Bo}=0.3$

Fig. 7Temporal evolution of the vertical velocity $U$ of the gas bubble：${Re}=20$,${Bo}=0.3$

2.4 表面张力对气泡形状及速度的影响

保持其他参数值不变,改变表面张力,探究不同$Bo$下气泡受电场力作用的上升行为.${Bo}$象征重力和表面张力的比值,较小的$Bo$对应着较大的表面张力,相应的浮力效应较小,不加电场的情况下,$Bo$越小,气泡形状越圆整. 外加电场时,气泡长径比变化如图8所示.对于所有的$Bo$值,气泡上升初始阶段响应电场力,沿电场方向呈明显拉伸状态,$Ar$均小于1.$Bo$为0.1 $\sim$ 0.5时,表面张力效应较大,$Ar$一直小于1,气泡呈长椭球状,$Bo$越小即表面张力越大的情况下,随着气泡的上升,气泡的长径比$Ar$相对越小,即气泡沿电场方向的变形程度越大;此时,$Bo$的减小使得气泡$Ar$值达到最小值的时间节点提前,$Bo$越低时,气泡所受的各种作用力间的相互竞争效应更强,气泡形状振荡越频繁,不稳定性更高.可以看出,当前计算域内,$Bo =0.1$时气泡$Ar$达到振幅的振荡次数几乎是$Bo =0.3$情况下的两倍.此外,如前文所述,表面张力效应增大导致相应的$Ar$值变小,阻力减小,气泡速度$U$越大（如图9）, 且均呈不稳定振荡状态,$Ar$的低峰值点和高峰值点分别对应相同时间节点的$U$的高峰值点和低峰值点.

图8

新窗口打开|下载原图ZIP|生成PPT
图8气泡长径比$Ar$随$T$的变化图：${Re}=20$,$Bo_{\rm e}=2$

Fig. 8Temporal evolution of the aspect ratio $Ar$ of the gas bubble：${Re}=20$,$Bo_{\rm e}=2$

图9

新窗口打开|下载原图ZIP|生成PPT
图9气泡速度$U$随$T$的变化图：${Re}=20$,$Bo_{\rm e}=2$

Fig. 9Temporal evolution of the vertical velocity $U $of the gas bubble：${Re}=20$,$Bo_{\rm e}=2$

$Bo\!=\!1$和$Bo\!=\!10$时,气泡最终长径比$Ar>1$,此时气泡惯性力效应更强,表面张力效应较低,气泡在上升初期响应电场力发生竖直方向拉伸后,受惯性力作用迅速拓展发生水平方向拉伸变形,并且$Ar$逐步趋于稳定值,没有再发生明显的振荡情况,可以看出,$Bo\!=\!1$时气泡后期呈现碗状,$Bo\!=\!10$时,气泡在$T\!=4.504$处呈现为盖状.气泡速度$U$值逐渐增大后也达到相对稳定的值,不再发生明显振荡.并且,$Bo$越小即表面张力越大的情况下,气泡$Ar$值越小,气泡沿电场方向发生拉伸变形程度越大,速度$U$越小.

2.5 液体黏度对气泡形状及速度的影响

黏性力也是影响气泡运动和变形的重要因素.$Re$表征惯性力与黏性力的比值,这里通过改变流体黏度获取不同大小的$Re$,$Re$对气泡形状和速度的影响如图10和图11所示.较大$Re$状态下,$Re$在10$\sim$30范围内时,由于黏性力较小,气泡对电力响应得更快,并因此沿电场方向伸长更多,$Ar$达到一个最小值,气泡上升初期的变形程度最大,随后由于各作用力相互竞争效应,气泡发生不稳定振荡变形,$Re$越大,即黏性力越小的情况下,振荡程度越大,之后变形波动程度逐步减小,并达到相对稳定状态,气泡初期与后期变形程度也相差更大.在$Re =30$时,电场使得气泡沿电场方向发生一定程度的拉伸,气泡在上升后期由于惯性力的主导作用,出现$Ar>1$的状态,即气泡发生横向变形.

图10

新窗口打开|下载原图ZIP|生成PPT
图10气泡长径比$Ar$随$T$的变化图：$Bo =0.3$,$Bo_{\rm e}=2$

Fig. 10Temporal evolution of the aspect ratio $Ar$ of the gas bubble：$Bo =0.3$,$Bo_{\rm e}=2$

图11

新窗口打开|下载原图ZIP|生成PPT
图11气泡速度$U$随$T$的变化图：$Bo =0.3$,$Bo_{\rm e}=2$

Fig. 11Temporal evolution of the vertical velocity $U$ of the gas bubble：$Bo =0.3$,$Bo_{\rm e}=2$

而当黏性力较大时,如$Re=1$和5的情况下,气泡首先响应电场力发生拉伸变形后,黏滞力阻碍了惯性力作用,使得气泡维持稳定的长椭球状,不再发生振荡. 随着$Re$的增大,气泡形态发生不稳定振荡,然而$Ar$达到幅值的时间点接近,相比于表面张力对气泡行为的影响,黏性力的改变似乎不会大幅度影响气泡发生振荡的频率. 气泡上升特征速度$U$值受气泡形状变化的影响也发生改变,相同$Re$下,速度$U$达到幅值的时间点与$Ar$幅值的时间点相对应.然而,不同$Re$下的$U$值大小差异并没有出现前文所述的类似情况,此时$Re$越大,气泡$Ar$值越大,$U$值也越大,这里考虑是黏性力作用影响,$Re$较大的情况下黏性力较小,气泡更容易受惯性力影响快速上升,$Re$较小时,周围液体较大的黏性力阻碍了气泡的上升.

3 结论

本研究对外加电场作用下单个气泡在流体中上升的动力行为进行了数值模拟研究. 利用水平集方法捕捉气液两相界面.通过$Re$,$Bo$和$Bo_{\rm e}$与气泡长径比及上升速度之间的关系研究了$Re<40$,$Bo<30$范围内的液体黏度、表面张力和电场力对气泡运动及变形的影响. 研究得到以下结论.

未加电场的情况下,当液体黏度和表面张力较大时,气泡几乎维持球状;反之,由于表面张力,惯性力和黏性力的不平衡,气泡在初始阶段的变形和加速上升. 在达到力平衡后,气泡形状波动微小并稳定上升.

外加电场对气泡形状和上升速度有显著影响. 气泡沿着电场方向发生拉伸变形,$Bo_{\rm e}$越大,电场强度越大,沿电场方向拉伸更明显,变形程度更大,且由于施加电场力后各作用力之间相互竞争,气泡形状出现竖直和水平方向交替的摇摆不定的振荡形式,且$Bo_{\rm e}$越大,振荡幅度越大. 随着气泡上升,气泡振幅随时间而减小. 此外,发现竖直方向的匀强电场加速了气泡的上升. 电场作用下,保持其他参数不变,随着液体黏度的减小,$Re$增大,气泡上升初期对电场的响应速度更快,气泡初期与后期变形程度相差更大,且气泡上升速度更大,气泡形态出现更为明显的振荡形式,更难达到稳定状态. 液体黏度较大时,气泡响应电场力发生变形后很快维持稳定形态,基本不发生波动,液体较强的黏性力也减缓了气泡的上升速度. 同一$Re$下,随着表面张力系数增大,$Bo$减小,气泡不稳定振荡更频繁. 气泡上升速度受气泡形状变化的影响,发生相同频率的不稳定振荡,且长径比$Ar$越小,气泡上升速度$U$越大.

参考文献原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]

韩宇, 刘志军, 王云峰等. T型微通道反应器内气液两相流动机制及影响因素
力学学报, 2019,51(2):441-449

[本文引用: 1]

( Han

, Liu

Zhijun

, Wang

Yunfeng

, et al. Gas-liquid two-phase flow regimes and impact factors in T-junction microreactor
Chinese Journal of Theoretical and Applied Mechanics, 2019,51(2):441-449 (in Chinese))

[本文引用: 1]

[2]

Yoon

, Koshizuka

, Oka

. Direct calculation of bubble growth, departure, and rise in nucleate pool boiling
International Journal of Multiphase Flow, 2001,27(2):277-298

[3]

Takahashi

, Miyahara

, Nishizaki

. Separation of oily water by bubble column
Journal of Chemical Engineering of Japan, 1979,12(5):394-399

[4]

Al-Shamrani

, James

, Xiao

. Separation of oil from water by dissolved air flotation
Colloids and Surfaces A (Physicochemical and Engineering Aspects), 2002,209(1):15-26

[5]

Pigeonneau

. Mass transfer of a rising bubble in molten glass with instantaneous oxidation--reduction reaction
Chemical Engineering Science, 2009,64(13):3120-3129

[6]

Galvin

, Pratten

, Shankar

, et al. Production of high internal phase emulsions using rising air bubbles
Chemical Engineering Science, 2001,56(21):6285-6293

[7]

Murray

, Dickinson

, Wang

. Bubble stability in the presence of oil-in-water emulsion droplets: Influence of surface shear versus dilatational rheology
Food Hydrocolloids, 2009,23(4):1198-1208

[本文引用: 1]

[8]

Clift

, Grace

, Weber

. Bubbles, Drops, Particles
Courier Corporation, 2005

[本文引用: 2]

[9]

Bonometti

, Magnaudet

. Transition from spherical cap to toroidal bubbles
Physics of Fluids, 2006,18(5):052102

[本文引用: 1]

[10]

Wang

, Zhang

, Yang

, et al. Viscosity effects on the behavior of a rising bubble
Journal of Hydrodynamics, 2010,22(1):81-89

[11]

Hua

, Lou

. Numerical simulation of bubble rising in viscous liquid
Journal of Computational Physics, 2007,222(2):769-795

[本文引用: 3]

[12]

Tripathi

, Sahu

, Govindarajan

. Dynamics of an initially spherical bubble rising in quiescent liquid
Nature Communications, 2015,6:6268

[本文引用: 2]

[13]

王焕然, 李彦鹏, 杨栋 . 黏性液体中单个气泡上升的形状特性
工程热物理学报, 2009,30(9):1492-1494

[本文引用: 1]

( Wang

Huanran

, Li

Yanpeng

, Yang

Dong

. On the shape feature of a single bubble rising in viscous liquids
Journal of Engineering Thermophysics, 2009,30(9):1492-1494 (in Chinese))

[本文引用: 1]

[14]

Chen

, Garimella

, Reizes

, et al. The development of a bubble rising in a viscous liquid
Journal of Fluid Mechanics, 1999,387:61-96

[本文引用: 2]

[15]

Amaya-Bower

, Lee

. Single bubble rising dynamics for moderate Reynolds number using Lattice Boltzmann method
Computers and Fluids, 2010,39(7):1191-1207

[本文引用: 1]

[16]

Alizadeh

, Seyyedi

, Rahni

, et al. Three-dimensional numerical simulation of rising bubbles in the presence of cylindrical obstacles, using Lattice Boltzmann method
Journal of Molecular Liquids, 2017,236:151-161

[本文引用: 1]

[17]

张洋, 陈科, 尤云祥等. 壁面约束对裙带气泡动力学的影响
力学学报, 2017,49(5):1050-1058

[本文引用: 1]

( Zhang

Yang

, Chen

, You

Yunxiang

, et al. Confinement effect on the rising dynamics of a skirted bubble
Chinese Journal of Theoretical and Applied Mechanics, 2017,49(5):1050-1058 (in Chinese))

[本文引用: 1]

[18]

Moore

. The rise of a gas bubble in a viscous liquid
Journal of Fluid Mechanics, 1959,6(1):113-130

[19]

Bhaga

, Weber

. Bubbles in viscous liquids: Shapes, wakes and velocities
Journal of Fluid Mechanics, 2006,105(105):61-85

[本文引用: 3]

[20]

Marco

, Grassi

, Memoli

. Experimental study on rising velocity of nitrogen bubbles in FC-72
International Journal of Thermal Sciences, 2003,42(5):435-446

[本文引用: 1]

[21]

Taylor

. Disintegration of water drops in an electric field
Proceedings A, 1964,280(1382):383-397

[本文引用: 2]

[22]

Garton

, Krasucki

. Bubbles in insulating liquids: stability in an electric field.
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1964,280(1381):211-226

[本文引用: 1]

[23]

陈凤, 宋耀祖, 陈民 . 电场作用下的气泡受力分析
工程热物理学报, 2005,26(S1):146-148

[本文引用: 1]

( Chen

Feng

, Song

Zuyao

, Chen

Min

. Analysis of the electric stresses acting on a bubble in electric field
Journal of Engineering Thermophysics, 2005,26(S1):146-148 (in Chinese))

[本文引用: 1]

[24]

陈凤, 宋耀祖, 陈民 . 电场作用下气泡内外的速度场分析
热科学与技术, 2006,5(2):139-143

[本文引用: 1]

( Chen

Feng

, Song

Zuyao

, Chen

Min

. Analysis of velocity field inside and outside bubble in DC electric field
Journal of Thermal Science and Technology, 2006,5(2):139-143 (in Chinese))

[本文引用: 1]

[25]

Yang

, Li

, Shao

, et al. A phase field numerical study of 3D bubble rising in viscous fluids under an electric field
International Journal of Heat & Mass Transfer, 2014,78(WCE2013):820-829

[本文引用: 2]

[26]

Rahmat

, Tofighi

, Yildiz

. The combined effect of electric forces and confinement ratio on the bubble rising
International Journal of Heat and Fluid Flow, 2017: S0142727X16305203

[本文引用: 2]

[27]

Andalib

, Hokmabad

, Esmaeilzadeh

. Study of a single coarse bubble behavior in the presence of D.C. electric field
Colloids and Surfaces A Physicochemical and Engineering Aspects, 2013,436:604-617

[本文引用: 2]

[28]

Baker

, Meiron

, Orszag

. Boundary integral methods for axisymmetric and three-dimensional Rayleigh-Taylor instability problems
Physica D Nonlinear Phenomena, 1984,12(1):19-31

[本文引用: 1]

[29]

Hua

, Stene

, Lin

. Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method
Journal of Computational Physics, 2008,227(6):3358-3382

[本文引用: 1]

[30]

Zhang

, Yan

, Zu

. Numerical modelling of EHD effects on heat transfer and bubble shapes of nucleate boiling
Applied Mathematical Modelling, 2010,34(3):626-638

[本文引用: 1]

[31]

张洋, 陈科, 尤云祥等. 浮力气泡对水平壁面的回弹动力学特性
力学学报, 2019,51(5):1285-1295

( Zhang

Yang

, Chen

, You

Yunxiang

, et al. Bouncing behaviors of a buoyancy-driven bubble on a horizontal solid wall
Chinese Journal of Theoretical and Applied Mechanics, 2019,51(5):1285-1295 (in Chinese))

[32]

Sussman

, Smereka

. Axisymmetric free boundary problems
Journal of Fluid Mechanics, 1997,341:269-294

[本文引用: 1]

[33]

Norman

, Miksis

. Dynamics of a gas bubble rising in an inclined channel at finite Reynolds number
Physics of Fluids, 2005,17(2):821-833

[34]

M?Hlmann

, Papageorgiou

. Buoyancy-driven motion of a two-dimensional bubble or drop through a viscous liquid in the presence of a vertical electric field
Theoretical and Computational Fluid Dynamics, 2009,23(5):375-399

[本文引用: 1]

[35]

Rahmat

, Tofighi

, Yildiz

. Numerical simulation of the electrohydrodynamic effects on bubble rising using the SPH method
International Journal of Heat and Fluid Flow, 2016,62:313-323

[本文引用: 1]

[36]

Mukundakrishnan

, Quan

, Eckmann

, et al. Numerical study of wall effects on buoyant gas-bubble rise in a liquid-filled finite cylinder
Physical Review E Statistical Nonlinear & Soft Matter Physics, 2007,76(3):036308

[本文引用: 1]

[37]

Sherwood

. Breakup of fluid droplets in electric and magnetic fields
Journal of Fluid Mechanics, 2006,188(188):133-146

[本文引用: 1]

[38]

Ohta

, Imura

, Yoshida

, et al. A computational study of the effect of initial bubble conditions on the motion of a gas bubble rising in viscous liquids
International Journal of Multiphase Flow, 2005,31(2):223-237

[本文引用: 1]