基于电流体动力学的气泡浮升行为分析

明德智, 李 强, 刘佳琳, 雷 毛, 郭 旭, 许伟伟, 王振波

基于电流体动力学的气泡浮升行为分析

明德智, 李 强, 刘佳琳, 雷 毛, 郭 旭, 许伟伟, 王振波

(中国石油大学(华东) 新能源学院, 山东 青岛 266580)

为了研究电流体动力学(EHD)在多相流体间传热传质的促进机理，基于开源计算流体力学软件OpenFOAM对电场作用下单个气泡上浮过程进行数值模拟。实现基于interFoam求解器的电场与流场的双向耦合方法，通过与已有试验结果的对比，验证所采用计算方法的准确性。结果表明，加入电场后，微气泡在短时间内发生剧烈形变与振荡，电场力越大，形变与振荡越明显；此外，在电场作用下，气液相表面张力系数过小将会导致气泡进入失稳状态。

电流体动力学(EHD)；气泡；OpenFOAM；两相流；数值分析；界面张力；数值模拟

1 前 言

气液两相流不仅存在于自然界中，还广泛存在于化工[1]、核工[2]、环保[3]等工业生产中，利用电场强化气液相间传热传质有着重要意义。然而该技术并未受到深入研究，直到60年代随着能源紧缺、生态恶化等原因，人们开始重视气液两相强化传热的研究。

目前，电场强化气液两相间传热的研究主要集中在控制气泡的运动与变形，使气泡获得较好的流动状态和上升特性。Taylor[4]和Garton等[5]研究表明，液滴和气泡在电场作用下会沿电场方向变形，在电场方向上气泡呈扩展的椭球形；Cho等[6]对单气泡进行理论和实验研究，结果表明，电场在气泡表面的非均匀分布使得气泡表面张力布局不均匀，导致气泡沿电场平行方向伸长；陈凤等[7-8]在小雷诺数蠕动流假设基础上，分析计算了不同介电常数气泡在外加电场作用下的受力状态以及气泡内外的速度场分布，结果表明，气泡伸长或压缩变形的主要因素取决于气液相的介电常数比值、电导率比值和透水黏性系数比值的相对大小。此外，电场的存在使得气泡内流体运动速度加快；Yang等[9]研究表明，直流电场的施加方向对于气泡的变形影响较大，而水平和垂直电场分别对气泡的上升起抑制和加速作用；王悦柔等[10]通过对比电邦德数，研究了电场力对气泡运动的影响，发现随着电邦德数的增大，气泡达到相对稳定的状态越难。

由于流场和电场耦合的作用，气泡形态的物性影响参数较多，单纯依靠实验难以进行全面比较与验证，数值模拟成为一种重要的研究方法。目前，应用于气液两相流研究的数值方法主要包括有限元法[11]、边界元法[12]、格子Boltzmann法[13]、流体体积(volume-of fluid，VOF)法[14]、水平集(LS)法[15]等。本研究基于OpenFOAM软件提出interEHDFoam求解器，对电场作用下单个气泡在流体中的上浮过程进行了数值模拟，将电场力作为显式体积力添加到流场控制方程中，实现了电场与流场的双向耦合作用，并利用VOF方法准确捕捉气液两相界面；然后，通过与参考文献的实验结果对比验证了该模型的有效性；最后，通过对比电邦德数、厄特沃什数和气液黏度比研究了电场力、液体表面张力系数和气液相黏度比对微气泡形变和运动特性的影响，为两相流强化传热等复杂的电流体动力学研究提供理论支持。

2 数值模型

本研究采用的气液两相模型如图1所示，图中，初始半径为，矩形流体计算域宽为8、高为16，气泡初始形状设定为圆形，气泡中心距离底部边界为2，设定气泡半径=1 mm。假设液滴与周围液体介质互不相溶，计算域边界应用无滑移壁面边界条件，计算域上下边界施加电势差，形成电场方向向下的均匀电场。

图1 均匀电场中气泡与外部流体计算模型示意图

2.1 VOF方法

2.2 电场方程

在电场和流场耦合过程中，产生的动态电流较小，因此磁感应效应可以忽略不计[16]，电场强度可以认为是无旋场，即

式中：为电场强度，V×m-1。定义e为体积电荷密度，C×m-3；e为上下极板电势差，V。电势差和体积电荷密度遵循以下泊松方程：

电荷守恒方程可以表示为

式中：为表面张力系数。

在外加电场作用下，麦克斯韦张量和电场体积力分别表示为

2.3 流体运动控制方程

电场作用下的气液两相流动受到多种力的耦合作用，主要包括惯性力、黏性力、电场力、表面张力等。对于不可压缩流体，考虑重力及电场力源项情况，流体的连续性方程和动量方程可分别表示为

2.4 求解方法

本研究提出的interEHDFoam求解器基于OpenFOAM 6.0版本，在Navier-Stokes方程中添加电场力源项并求解。首先，利用VOF方法计算网格内的相场值，得到气液界面的位置与形状，以及网格内流体的物性参数；其次，求解电场方程式(7)、(8)，得到计算域内的电场分布与电荷分布，利用方程式(11)求得电场作用于流体单元网格的体积力；最后，将电场体积力代入N-S方程，计算连续性方程式(12)和动量方程式(13)，进行速度与压力的耦合求解。interEHDFoam求解器基于有限体积法，速度和压力的耦合采用PIMPLE算法，计算过程中设定库朗特数0.5，计算流程见图2。图中，total为计算时间，s；D为时间步长，s。

图2 InterEHDFoam求解过程

3 计算模型验证

3.1 无量纲参数定义

为方便描述气泡的运动特性，用电邦德数e表示气泡所受电场力与界面表面张力关系，e越大，电场力影响越大；定义厄特沃什数表示气泡所受浮力和界面表面张力关系，厄特沃什数越小，表面张力系数越大，如式(15)、(16)所示。

图3 气泡长宽比D定义

3.2 interEHDFoam求解器准确性验证

为了验证求解器的准确性，本研究建立了如图4所示的计算域。计算域内充满乙醇，设置边界条件Air-in为速度入口，Out为压力出口。Upwall与Downwall距离200 mm，施加电场后计算域内形成竖直向下的均匀电场，工况条件及物性参数参考文献[17]。当e=0.34时，气泡形状与长宽比随时间变化如图5、表1所示。由图5、表1可知，数值模拟与参考文献[17]实验测得的气泡形状相似度较高，气泡长宽比接近，表明interEHDFoam能够较准确地捕捉到电场作用下的气液相界面形状与位置。

图4 求解器精确性验证模型

Fig.4 Model for accuracy verification of solver

图5 模拟与参考文献[17]中实验图像对比

表1 当前模拟结果与参考文献[17]中实验结果对比

表2 网格无关性验证结果

3.3 网格无关性验证

为避免因网格密度所导致的差异性结果，分别采用20万、25万、27万、30万、40万网格数目对图1中矩形计算域进行网格无关性验证，当e=1.9时，得到稳定后的气泡长宽比见表2。计算域网格数量达到30万时，气泡的长宽比相对误差已经降至1%，说明继续加密网格对计算结果的影响已经很小，因此本研究模拟计算中，采用网格数量30万较为合适。

4 模拟结果与分析

4.1 电场力对气泡运动特性的影响

本研究在其他物性参数不变的条件下，通过改变电压大小模拟了6种不同e情况下气泡的变形情况与速度场分布规律，如图6所示。当e=0时计算域内没有电场存在，由于计算域高度较小，气泡在上浮过程中压力变化可以忽略不计，因此气泡长宽比保持不变。当流体域中加入电场后，随着e的增加，气泡在稳定前均出现振荡情况。并且施加的电压越大，气泡趋于稳定的时间越短，气泡稳定后的长宽比越大。

如图7、8所示为当e=3.9时电压分布图与电场强度分布图，图中为沿竖直中心线的纵坐标，由于气液界面处的介电常数梯度变化较大，流场中电压和电势会在该位置发生突变。气液界面处的电荷主要聚集于气泡界面的顶端和底端，并相对于电势方向产生正向偶极子，气泡界面处感应电荷受到竖直方向的电场力作用，使得气泡被拉成瘦长形。从图7、8中还可以看出，在=0.004～0.010 s时，气泡顶端移动距离较小，气泡的底端回弹明显，说明电场在气泡表面的非均匀分布使得气泡表面张力布局趋于不均匀，不均匀的表面张力布局对气泡底端影响更大，使气泡底端形状相较于顶端更加“尖锐”。

图6 不同电邦德数气泡长宽比随时间变化图

图7 计算域竖直中心线电压随坐标变化图

图8 计算域竖直中心线电场强度随坐标变化图

气泡在上浮过程中主要受浮力、曳力和惯性力作用，其中曳力受气泡形状尤其是气泡的长宽比影响较大[18]。无电场时，气泡可以在极短时间内达到稳定状态匀速上升，上浮过程中的形状和速度保持不变。如图9所示为加入电场后，在不同e下，气泡上浮速度up随时间的变化曲线，从图中可以看出，气泡稳定前的上浮速度随着e的增加而增加，主要原因是电场力拉伸气泡，使得气泡的曳力减小，在未稳定前气泡加速上浮，并且气泡的长宽比越大，初始加速度越大，达到稳定时的速度也越大。

图9 不同电邦德数下气泡上浮速度

4.2 电场作用下表面张力和黏性力对气泡运动特性的影响

为了进一步研究电场作用下表面张力对气泡上浮形状和运动状态的影响，在其他物性参数不变的情况下，图10、11给出了e=3.9(e=10 000 V)时不同数下气泡上浮过程中气泡体积率云图，图中l为液相体积分数。由图11可知，在高数下，气泡将不再保持椭圆形而是呈纺锤形运动，长宽比不再适合描述气泡形态。并且在上浮过程中，气泡受到电场力作用被继续拉伸直至进入失稳状态[6,19]。

图10 低Eo数下气泡云图

图11 高Eo数下气泡云图

图12 低Eo数下气泡长宽比随时间变化图

图12给出了=0.056、=0.56时气泡的长宽比变化，从图中可以发现，加入电场后，气泡出现明显的振荡现象，当=0.56时气泡仅振荡一次，0.010 s后保持稳定的长宽比；而当=0.056时气泡出现了两次振荡，在0.006 s后保持稳定的长宽比。由上述可知，电场力作用下气液2相表面张力系数对维持气泡形态有重要作用，表面张力系数越大，气液相界面越具有“弹性”，气泡保持圆形度的能力越高。

通过改变液相的黏度来改变气液相之间的黏度比，研究电场作用下气泡振荡时所造成的流场变化。当e=3.9(e=10 000 V)时，假设其他参数不变，1为默认黏度比，取2=101，3=1001。

图13给出了气泡稳定前(=0.003～0.012 s)气泡内部及周围流体的电场、速度矢量图。由图13(a)可知，电场作用下不同气液黏度比对气泡形状造成的影响变化较小，但气泡内部以及周围的电场出现了不同程度扰动。黏度比越小，由图13(b)可以看出，气泡底部短时间内(0.003～0.006 s)产生的微射流也会更明显。

图13 不同黏度比对气泡的影响

5 结 论

[1] 霍磊, 时骋, 计世俊. 化工装置中气液两相流管道的设计探讨[J]. 化工设计, 2019, 29(5): 17-20.

HUO L, SHI P, JI S J. Discussion on the design of gas-liquid two-phase flow pipeline in chemical plant [J]. Chemical Engineering Design, 2019, 29(5): 17-20.

[2] 毕景良. 微尺度核态池沸腾传热传质特性及机理研究[D]. 北京: 清华大学, 2014.

BI J L. Micro-scale nucleate pool boiling heat and mass transfer characteristics and mechanisms [D]. Beijing: Tsinghua University, 2014.

[3] 张敏, 宋昭峥, 孙珊珊, 等. 微纳米气浮技术用于炼化污水的深度处理[J]. 环境工程学报, 2016, 10(2): 599-603.

ZHANG M, SONG Z Z, SUN S S,. Advanced treatment of refining sewage with a new micro/nano-bubble flotation technology [J]. Chinese Journal of Environmental Engineering, 2016, 10(2): 599-603.

[4] TAYLOR G. Disintegration of water drops in an electric field [J]. Proceedings A, 1964, 280(1382): 383-397.

[5] GARTON G G, KRASUCKI Z. Bubbles in insulating liquids: stability in an electric field [J]. Proceedings of the Royal Society of London, 1964, 280(1381): 211-226.

[6] CHO H J, KANG I S, KWEON Y C,. Study of the behavior of a bubble attached to a wall in a uniform electric field [J]. International Journal of Multiphase Flow, 1996, 22(5): 909-922.

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

CHEN F, SONG Y Z, CHEN M. Analysis of the electric stresses acting on a bubble in electric field [J]. Journal of Engineering Thermophysics, 2005, 26(S1):146-148.

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

CHEN F, SONG Y Z, CHEN M. Analysis of velocity field inside and outside bubble in DC electric field [J]. Journal of Thermal Science and Technology, 2006, 5(2): 139-143.

[9] YANG Q, LI B Q, SHAO J,A phase field numerical study of 3D bubble rising in viscous fluids under an electric field [J]. International Journal of Heat and Mass Transfer, 2014, 78: 820-829.

[10] 王悦柔, 王军锋, 刘海龙. 电场作用下气泡上升行为特性的数值计算研究[J]. 力学学报, 2020, 52(1): 31-39.

WANG Y R, WANG J F, LIU H L. Numerical simulation on bubble rising behaviors under electric field [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(1): 31-39.

[11] FENG J Q, SCOTT T C. A computational analysis of electrohydrodynamics of a leaky dielectric drop in an electric field [J]. Journal of Fluid Mechanics, 1996, 311: 289-326.

[12] LAC E, HOMSY G M. Axisymmetric deformation and stability of a viscous drop in a steady electric field [J]. Journal of Fluid Mechanics, 2007, 590: 239-264.

[13] 黄伟峰, 李勇, 刘秋生. 格子Boltzmann方法在电流体动力学中的应用: 均匀电场中液滴的变形和失稳[J]. 科学通报, 2007, 11: 1232-1236.

HUANG W F, LI Y, LIU Q S. Application of lattice Boltzmann method in electrohydrodynamics: Droplet deformation and instability in uniform electric field [J]. Chinese Science Bulletin, 2007, 11: 1232-1236.

[14] HIRT C W, NICHOLS B D. Volume of fluid (VOF) method for the dynamics of free boundaries [J]. Journal of Computational Physics, 1981, 39(1): 201-225.

[15] HO Y K, YONG W K. Homogeneous nucleation and macroscopic growth of gas bubble in organic solutions [J]. International Journal of Heat and Mass Transfer, 1998, 41(4/5): 757-767.

[16] MELCHER J R, TAYLOR G I. Electrohydrodynamics: A review of the role of interfacial shear stresses [J]. Annual Review of Fluid Mechanics, 1969, 1(1): 111-146.

[17] ZHANG W, WANG J, LI B,EHD effects on periodic bubble formation and coalescence in ethanol under non-uniform electric field [J]. Chemical Engineering Science, 2019, 215: 115451.

[18] 闫盛楠. 鼓泡流化床不规则形状颗粒气固两相流动特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2014.

YAN S N. Investigation on gas-solid two-phaseflow characteristics of irregularlyshaped particle in a bubblingfluidized bed [D]. Harbin: Harbin Institute of Technology, 2014.

[19] MORIYA S, ADACHI K, KOTAKA T. Deformation of droplets suspended in viscous media in an electric field 1. Rate of deformation [J]. Langmuir, 1986, 2(2): 155-160.

Investigation of rising bubble dynamics based on electrohydrodynamics

MING De-zhi, LI Qiang, LIU Jia-lin, LEI Mao, GUO Xu, XU Wei-wei, WANG Zhen-bo

(College of New Energy, China University of Petroleum, Qingdao 266580, China)

A single bubble floating process under electric field was numerically simulated based on OpenFOAM software to study the mechanism of electrohydrodynamics (EHD) promoted multiphase heat and mass transfer. The bidirectional coupling method of electric field and flow field was realized based on interFoam solver and the accuracy of the calculation method was verified by comparing with existing test results. The results show that microbubbles undergo severe deformation and oscillation in a short time after the addition of electric field, and the deformation and oscillation are more obvious with the increase of electric field force. In addition, the bubble is more unstable under electric field if the surface tension coefficient of gas-liquid interphase is too small.

electrohydrodynamics; bubble; OpenFOAM; two-phase flow; surface tension; numerical simulation

1003-9015(2022)06-0793-08

TQ021.1

A

10.3969/j.issn.1003-9015.2022.06.003

2021-09-15；

2022-01-19。

国家自然科学基金(51506225，52176050)；山东省自然科学基金 (ZR2020ME174)；山东省重点研发项目(2018GHY115018)；中央高校基本科研业务费专项资金(18CX02129A)。

明德智(1996-)，男，山东泰安人，中国石油大学(华东)硕士生。

李强，E-mail：liq@upc.edu.cn

明德智, 李强, 刘佳琳, 雷毛, 郭旭, 许伟伟, 王振波. 基于电流体动力学的气泡浮升行为分析[J]. 高校化学工程学报, 2022, 36(6): 793-800.

：MING De-zhi, LI Qiang, LIU Jia-lin, LEI Mao, GUO Xu, XU Wei-wei, WANG Zhen-bo. Investigation of rising bubble dynamics based on electrohydrodynamics [J]. Journal of Chemical Engineering of Chinese Universities, 2022, 36(6): 793-800.