高级检索

  • ISSN 1006-3080
  • CN 31-1691/TQ
引用本文:
Citation:

鼓泡塔反应器中两相流动CFD-PBM耦合数值模拟

    作者简介: 贾翔飞(1995-),男,山东济宁人,硕士生,主要研究方向为计算流体力学;
    通讯作者: 陈剑佩, chenjianpei@ecust.edu.cn
  • 中图分类号: TQ519

CFD-PBM Coupled Simulation of Two Phase Flow in Bubble Column Reactor

    Corresponding author: Jianpei CHEN, chenjianpei@ecust.edu.cn
  • CLC number: TQ519

  • 摘要: 针对气液体系,将一个新型气泡破碎模型植入群体平衡模型(PBM)并与欧拉双流体模型进行耦合,对内径为0.44 m、高2.43 m的鼓泡塔中气-液两相流动进行模拟。在对破碎模型的气泡破碎速率和子气泡分布预测值进行验证之后,将计算流体动力学-群体平衡耦合模型(CFD-PBM)在不同表观气速下的气含率、轴向液速模拟值与实验值进行了对比并探究了升力的加入对模拟的影响。结果显示模拟过程中升力的加入使气含率和轴向液速的径向分布更接近实验值。低表观气速和高表观气速下,模型的计算数据与实验数据均吻合良好。同时,随表观气速增大,反应器内的湍动耗散率增加,气泡之间相互作用作用增强,使得反应器内气泡尺寸分布范围变宽。
  • 图 1  鼓泡塔反应器几何模型(a)及网格截面(b)

    Figure 1.  Geometry (a) and mesh (b) of bubble column

    图 2  量纲为一破碎速率预测值与实验值对比

    Figure 2.  Predicted and experimental values of dimensionless breakup rate

    图 3  子气泡分布预测值与实验值

    Figure 3.  Predicted and experimental values of daughter bubble size distribution

    图 5  不同表观气速下气含率径向分布

    Figure 5.  Radial profiles of gas holdup in different superficial gas velocity

    图 6  不同表观气速下液相轴向速度径向分布

    Figure 6.  Radial profiles of liquid axial velocity in different superficial gas velocity

    图 4  不同表观气速下鼓泡塔反应器内气含率及轴向液速分布云图

    Figure 4.  Gas holdup and liquid axial velocity contours of different superficial gas velocity in bubble column

    图 7  不同表观气速下气泡轴向速度径向分布

    Figure 7.  Radial profiles of bubble axial velocity in different superficial gas velocity

    图 8  不同表观气速下湍动耗散率径向分布

    Figure 8.  Radial profiles of turbulent dissipation rate in different superficial gas velocity

    图 9  H/D=3处不同气速下气泡尺寸分布

    Figure 9.  Bubble size distribution at H/D=3 in different superficial gas velocity

    图 10  表观气速下(a)r/R=0和(b)r/R=0.5的气泡尺寸分布

    Figure 10.  Bubble size distribution at (a) r/R=0, and (b) r/R=0.5 in different superficial gas velocity

    图 11  经向位置处气泡尺寸分布

    Figure 11.  Bubble size distribution at different radial positions

  • [1] ROLLBUSCH P, BOTHE M, BECKER M, et al. Bubble columns operated under industrially relevant conditions: current understanding of design parameters[J]. Chemical Engineering Science, 2015, 126: 660-678. doi: 10.1016/j.ces.2014.11.061
    [2] LIANG X F, PAN H, SU Y H, et al. CFD-PBM approach with modified drag model for the gas–liquid flow in a bubble column[J]. Chemical Engineering Research and Design, 2016, 112: 88-102. doi: 10.1016/j.cherd.2016.06.014
    [3] HUANG Q S, ZHANG W P, YANG C. Modeling transport phenomena and reactions in a pilot slurry airlift loop reactor for direct coal liquefaction[J]. Chemical Engineering Science, 2015, 135: 441-451. doi: 10.1016/j.ces.2015.01.003
    [4] VIAL C, PONCIN S, WILD G, et al. A simple method for regime identification and flow characterisation in bubble columns and airlift reactors[J]. Chemical Engineering and Processing: Process Intensification, 2001, 40(2): 135-151. doi: 10.1016/S0255-2701(00)00133-1
    [5] DUDUKOVIĆ M P, LARACHI F, MILLS P L. Multiphase catalytic reactors: A perspective on current knowledge and future trends[J]. Catalysis Reviews, 2002, 44(1): 123-246. doi: 10.1081/CR-120001460
    [6] BHOLE M, JOSHI J, RAMKRISHNA D. CFD simulation of bubble columns incorporating population balance modeling[J]. Chemical Engineering Science, 2008, 63(8): 2267-2282. doi: 10.1016/j.ces.2008.01.013
    [7] 叶庆元, 黄建科, 陈剑佩. 釜式光生物反应器内光分布及光能吸收的CFD研究[J]. 华东理工大学学报(自然科学版), 2016, 42(02): 149-156.
    [8] 刘希祥, 陈剑佩. 螺旋管强化传热的CFD模拟与优化[J]. 华东理工大学学报(自然科学版), 2018, 44(3): 296-302.
    [9] LABORDE-BOUTET C, LARACHI F, DROMARD N, et al. CFD simulation of bubble column flows: Investigations on turbulence models in RANS approach[J]. Chemical Engineering Science, 2009, 64(21): 4399-4413. doi: 10.1016/j.ces.2009.07.009
    [10] SANYAL J, VÁSQUEZ S, ROY S, et al. Numerical simulation of gas–liquid dynamics in cylindrical bubble column reactors[J]. Chemical Engineering Science, 1999, 54(21): 5071-5083. doi: 10.1016/S0009-2509(99)00235-3
    [11] OEY R, MUDDE R, VAN DEN AKKER H. Sensitivity study on interfacial closure laws in two‐fluid bubbly flow simulations[J]. AIChE Journal, 2003, 49(7): 1621-1636. doi: 10.1002/aic.690490703
    [12] TABIB M V, ROY S A, JOSHI J B. CFD simulation of bubble column: An analysis of interphase forces and turbulence models[J]. Chemical Engineering Journal, 2008, 139(3): 589-614. doi: 10.1016/j.cej.2007.09.015
    [13] LIAO Y, MA T, KREPPER E, et al. Application of a novel model for bubble-induced turbulence to bubbly flows in containers and vertical pipes[J]. Chemical Engineering Science, 2019, 202: 55-69. doi: 10.1016/j.ces.2019.03.007
    [14] YANG G, ZHANG H, LUO J, et al. Drag force of bubble swarms and numerical simulations of a bubble column with a CFD-PBM coupled model[J]. Chemical Engineering Science, 2018, 192: 714-724. doi: 10.1016/j.ces.2018.07.012
    [15] GUO X, ZHOU Q, LI J, et al. Implementation of an improved bubble breakup model for TFM-PBM simulations of gas-liquid flows in bubble columns[J]. Chemical Engineering Science, 2016, 152: 255-266. doi: 10.1016/j.ces.2016.06.032
    [16] WANG T, WANG J, JIN Y. A CFD-PBM coupled model for gas-liquid flows[J]. AIChE Journal, 2006, 52(1): 125-140. doi: 10.1002/aic.10611
    [17] XING C, WANG T, GUO K, et al. A unified theoretical model for breakup of bubbles and droplets in turbulent flows[J]. AIChE Journal, 2015, 61(4): 1391-1403. doi: 10.1002/aic.14709
    [18] LUO H, SVENDSEN H F. Theoretical model for drop and bubble breakup in turbulent dispersions[J]. AIChE Journal, 1996, 42(5): 1225-1233. doi: 10.1002/aic.690420505
    [19] 韩路长. 多相反应器内流体颗粒传质、破裂分散和聚并行为的研究 [D]. 湖南湘潭: 湘潭大学; 2010.
    [20] TOMIYAMA A. Struggle with computational bubble dynamics[J]. Multiphase Science and Technology, 1998, 10(4): 369-405. doi: 10.1615/MultScienTechn.v10.i4.40
    [21] LOPEZ DE BERTODANO M. Turbulent bubbly flow in a triangular duct[J]. Nuclear Engineering and Design, 1994, 146(1/3): 43-52
    [22] HAGESAETHER L, JAKOBSEN H A, Svendsen H F. A model for turbulent binary breakup of dispersed fluid particles[J]. Chemical Engineering Science, 2002, 57(16): 3251-3267. doi: 10.1016/S0009-2509(02)00197-5
    [23] PRINCE M J, BLANCH H W. Bubble coalescence and break‐up in air‐sparged bubble columns[J]. AIChE Journal, 1990, 36(10): 1485-1499. doi: 10.1002/aic.690361004
    [24] WANG T, WANG J. Numerical simulations of gas–liquid mass transfer in bubble columns with a CFD–PBM coupled model[J]. Chemical Engineering Science, 2007, 62(24): 7107-7118. doi: 10.1016/j.ces.2007.08.033
    [25] CHESTERS A K. Modelling of coalescence processes in fluid-liquid dispersions: A review of current understanding[J]. Chemical Engineering Research and Design, 1991, 69(A4): 259-270.
    [26] LUO H, SVENDSEN H. Modeling and simulation of binary approach by energy conservation analysis[J]. Chemical Engineering Communications, 1996, 145(1): 145-153. doi: 10.1080/00986449608936473
    [27] HIBIKI T, ISHII M. Two-group interfacial area transport equations at bubbly-to-slug flow transition[J]. Nuclear Engineering and Design, 2000, 202(1): 39-76. doi: 10.1016/S0029-5493(00)00286-7
    [28] HAN L, LUO H, LIU Y. A theoretical model for droplet breakup in turbulent dispersions[J]. Chemical Engineering Science, 2011, 66(4): 766-776. doi: 10.1016/j.ces.2010.11.041
    [29] ANDERSSON R, ANDERSSON B. Modeling the breakup of fluid particles in turbulent flows[J]. AIChE Journal, 2006, 52(6): 2031-2038. doi: 10.1002/aic.10832
    [30] RAZZAGHI K, SHAHRAKI F. Theoretical model for multiple breakup of fluid particles in turbulent flow field[J]. AIChE Journal, 2016, 62(12): 4508-4525. doi: 10.1002/aic.15314
    [31] ZHAO H, GE W. A theoretical bubble breakup model for slurry beds or three-phase fluidized beds under high pressure[J]. Chemical Engineering Science, 2007, 62(1/2): 109-115. doi: 10.1016/j.ces.2006.08.008
    [32] ANDERSSON R, ANDERSSON B. On the breakup of fluid particles in turbulent flows[J]. AIChE Journal, 2006, 52(6): 2020-2030. doi: 10.1002/aic.10831
    [33] CHEN J, GUPTA P, DEGALEESAN S, et al. Gas holdup distributions in large-diameter bubble columns measured by computed tomography[J]. Flow Measurement and instrumentation, 1998, 9(2): 91-101. doi: 10.1016/S0955-5986(98)00010-7
    [34] LEHR F, MILLIES M, MEWES D. Bubble‐size distributions and flow fields in bubble columns[J]. AIChE Journal, 2002, 48(11): 2426-2443. doi: 10.1002/aic.690481103
    [35] VEJRAŽKA J, ZEDNÍKOVÁ M, STANOVSKÝ P. Experiments on breakup of bubbles in a turbulent flow[J]. AIChE Journal, 2018, 64(2): 740-757. doi: 10.1002/aic.15935
    [36] MARTÍNEZ-BAZÁN C, MONTANES J, LASHERAS J C. On the breakup of an air bubble injected into a fully developed turbulent flow. Part 1. Breakup frequency[J]. Journal of Fluid Mechanics, 1999, 401: 157-182. doi: 10.1017/S0022112099006680
    [37] KRISHNA R, WILKINSON P, VAN DIERENDONCK L. A model for gas holdup in bubble columns incorporating the influence of gas density on flow regime transitions[J]. Chemical Engineering Science, 1991, 46(10): 2491-2496. doi: 10.1016/0009-2509(91)80042-W
  • [1] 张圆飞闫宇强宗原曹发海 . 以鲍尔环为内构件的鼓泡塔反应器的流体力学性能. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180427001
    [2] 陈阳刘锦泽于兆臣王智峰欧阳福生 . 重油催化裂化12集总动力学模型的工业模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200316002
    [3] 叶贞成王鑫梅华 . 基于改进差分进化算法的裂解反应动力学系数辨识. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190910002
    [4] 王鹏李涛张海涛房鼎业 . 基于钯炭催化剂的苯酚加氢反应及其本征动力学. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki,1006-3080.20200308002
    [5] 胡志梁于新海赖焕新 . 拟塑性流体在六直叶搅拌器中的模拟与比较. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180517006
    [6] 王宁曹萃文 . 基于XGBoost模型的炼油厂氢气网络动态多输出预测模型. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20181104001
    [7] 潘傲谢明辉夏建业庄英萍 . 基于均龄理论模拟搅拌式反应器的混合时间. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180716001
    [8] 蔡冬莹王跃林段先建胡丹李玉冰胡彦杰 . 预混合高速射流燃烧反应器内温度场的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190221004
    [9] 吴兆伟赵辉吴欣洁刘海峰周鹜蔡小舒 . 液体射流内部气泡对黏性流体气流式雾化的影响. 华东理工大学学报(自然科学版),
    [10] 于苗钱锋胡贵华隆建李天越 . 工业级MIP提升管反应器气固两相流动特性的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913003
    [11] 陆炳祥张德祥 . 低次烟叶热解特性及动力学机理函数. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180716007
    [12] 熊文莉杨家鹏安琦 . 轴承滚子加工精度对高速电主轴动力学性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180613002
    [13] 高炳舒刘士荣 . 基于BoW模型的RGB-D SLAM算法的运动轨迹估计. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180419001
    [14] 张剑超杜文莉覃水 . 基于新型自适应采样算法的催化重整过程代理模型. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180915001
    [15] 马芳芳熊达孙铁栋欧阳福生 . 乙烯装置裂解气压缩机性能预测模型研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191226002
    [16] 李东刘乙奇黄道平 . 基于Tri-training MPLS的半监督软测量模型. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191202008
    [17] 许奕杰王嵘万永菁孙静 . 基于AE-LSTM网络模型的机场周界入侵报警及分类算法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200122001
    [18] 饶毓和凌志浩 . 一种结合主题模型与段落向量的短文本聚类方法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190430001
    [19] 张幸子王晓惠王泽建陈必钦李丹郭美锦储炬庄英萍 . 等离子体作用结合氧限制模型选育辅酶Q10高产菌株. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200227002
    [20] 陈颖呈陈宁 . 基于音频内容和歌词文本相似度融合的翻唱歌曲识别模型. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191029002
  • 加载中
图(11)
计量
  • 文章访问数:  1967
  • HTML全文浏览量:  606
  • PDF下载量:  26
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-21
  • 网络出版日期:  2019-12-28

鼓泡塔反应器中两相流动CFD-PBM耦合数值模拟

    作者简介:贾翔飞(1995-),男,山东济宁人,硕士生,主要研究方向为计算流体力学
    通讯作者: 陈剑佩, chenjianpei@ecust.edu.cn
  • 1. 华东理工大学 化学工程联合国家重点实验室
  • 2. 煤气化及能源化工教育部重点实验室,上海 200237

摘要: 针对气液体系,将一个新型气泡破碎模型植入群体平衡模型(PBM)并与欧拉双流体模型进行耦合,对内径为0.44 m、高2.43 m的鼓泡塔中气-液两相流动进行模拟。在对破碎模型的气泡破碎速率和子气泡分布预测值进行验证之后,将计算流体动力学-群体平衡耦合模型(CFD-PBM)在不同表观气速下的气含率、轴向液速模拟值与实验值进行了对比并探究了升力的加入对模拟的影响。结果显示模拟过程中升力的加入使气含率和轴向液速的径向分布更接近实验值。低表观气速和高表观气速下,模型的计算数据与实验数据均吻合良好。同时,随表观气速增大,反应器内的湍动耗散率增加,气泡之间相互作用作用增强,使得反应器内气泡尺寸分布范围变宽。

English Abstract

  • 我国的自然资源呈现贫油少气富煤的现象,每年需要进口大量的石油。为了保障我国的能源安全,我国对煤液化技术进行了技术储备并且成功将该技术工业化。煤液化技术的关键设备是液化反应器,其本质上属于鼓泡塔反应器。该种反应器具有良好的混合性能并且生产能力大、结构简单、易于操作[1-2]

    煤液化反应器中涉及到氢气(气)、供氢溶剂(液)、煤颗粒(固)的流动问题,但反应器中的煤粉颗粒较小并且煤粉浓度较低,煤粉颗粒实际上是伴随着液相一起运动[3],因此通常将煤粉颗粒与液相作为混合均匀的浆相处理,将气-液-固三相流动简化为气-浆两相流动。根据反应器操作的表观气速不同,反应器内的状态可以分为:均匀鼓泡状态、过渡状态、湍动鼓泡态。在均匀鼓泡状态下,鼓泡塔内的气泡尺寸小而均匀,气含率与表观气速呈现近似线性关系[4]。工业生产过程中的鼓泡塔反应器一般处于湍动鼓泡态[5]。在湍动鼓泡态下,流动结构由于气泡大量的聚并和破碎行为而变得非常复杂。强烈的聚并和破碎行为使得鼓泡塔内的气泡尺寸分布范围较宽并对鼓泡塔内对于反应器设计和放大有关的气含率、液速、热量和质量传递速率等参数产生影响[1]。因此,反应器内的气泡破碎聚并以及基础流体力学行为对反应器的设计放大以及安全运行有着重要影响。

    在对反应器进行设计和放大时,计算流体动力学(CFD)逐渐成为一种重要的工具[2, 6-8]。在综合考虑计算量和计算精度的前提下,通常采用Eulerian-Eulerian方法[9]对鼓泡塔进行模拟。使用该方法研究气液鼓泡塔的流体力学行为关键在于合理的描述气-液相间作用力和气泡的大小分布。相间作用力包括曳力、升力、壁面润滑力及湍流扩散力。曳力是由气液之间的摩擦产生的,曳力的大小决定了气泡在液相中的上升速度并最终影响鼓泡塔内的气含率的大小,而升力、壁面润滑力和湍流扩散力会影响气含率在径向上的分布。

    由于气泡与气泡之间存在复杂相互作用,鼓泡塔内的气泡存在局部气泡直径分布。早期的研究者通常将该气泡直径分布简化为一个恒定的当量气泡直径[10-12]。对于均匀鼓泡态下,气泡之间的相互作用较弱,气泡的尺寸分布较窄,气泡直径均一,通过调节曳力模型及气泡直径会获得与实验值较吻合的结果。但在湍动鼓泡态下,气泡之间强烈的相互作用造成较宽的气泡的尺寸分布,上述对气泡直径的简化无法获得与实验值吻合的结果。

    由于气泡之间的相互作用对鼓泡塔内的流体力学参数有重要影响,群体平衡方程逐渐被应用于描述鼓泡塔内的气泡尺寸的变化及分布,并将其与计算流体动力学进行耦合来对鼓泡塔进行模拟[13-16]。多年来,人们做出了大量的努力来建立合理的聚并及破碎模型[16-19],并希望将其应用于CFD-PBM耦合模型中。

    在本文中,一个新颖的气泡破碎模型[19]将被植入群体平衡模型并与欧拉双流体模型进行耦合。首先在ANSYS Fluent平台中将CFD模型所需要的相间作用力模型通过用户自定义函数(UDF)加入,然后将PBM模型中聚并、破碎模型源项通过UDF植入并与CFD模型耦合建立CFD-PBM模型。采用该模型对实验室尺寸的圆柱鼓泡塔进行二维轴对称模拟并与文献中的实验数据进行比较。

    • 在欧拉双流体模型中,离散相被当做拟流体,并与连续相相互渗透。在这里液相被视为连续相,而气相作为离散相。气相和液相的连续性方程及动量方程分别如式(1)和式(2)所示,其中P′为修正压力,μeff,i代表有效黏度,下标i代表气相或者液相(i=l, g)。

      文献[9]中提到,采用重整化群k-εRNG k-ε)湍流模型可以更好地描述鼓泡塔中的流动。我们采用RNG k-ε湍流模型计算液相中的湍动能和湍流耗散率。

    • 曳力的大小会影响气泡在液相中的上升速度,并进一步影响鼓泡塔内气含率的大小,曳力通常采用式(3)计算[14]

      其中Kgl代表两相之间的动量交换系数,Cdi代表曳力系数。

    • 液相中速度梯度的存在会使气泡受到一个横向升力,这会对流体力学参数的径向分布产生影响,升力系数的计算采用Tomiyama模型[20]

      其中,Eod数是根据气泡的变形后的长轴计算出的修正Eötvös数。

    • 气泡受到的壁面润滑力采用Tomiyama模型[20]

    • 湍流耗散力采用Lopez de Bertodano模型[21]

    • 在鼓泡塔中,存在着数量非常庞大的气泡,这些气泡会发生聚并和破碎,最终使鼓泡塔中的气泡存在一定的尺寸分布。群体平衡方程(PBE,式(15))以及合理的聚并、破碎源项的引入使得鼓泡塔中的气泡尺寸分布得以进行描述。

      式(15)本质上为气泡数密度函数$n\left( {{v_i},t} \right)$的传输方程[22],其与气泡体积v和时间t有关。

      Si为第i组气泡中的PBM方程的源项[22],在目前的研究工作中,成核及生长速率被忽略,只考虑聚并生成、聚并消失、破碎生成及破碎消失项(BC,iDC,iBB,iDB,i)。各源项的计算采用文献[22]的方法。

    • 鼓泡塔中引起气泡之间相互碰撞的机理通常有三种,分别是液相湍流引起的碰撞、气泡上升速度差异引起的碰撞以及大气泡尾涡引起的碰撞。但不是每次气泡碰撞都会使两个气泡发生聚并,所以通过引入气泡聚并效率的概念来表示两个气泡聚并的可能性。最终,可以通过式(17)来计算直径为didj的气泡的聚并速率。鼓泡塔中总的聚并速率可以认为是三种机理聚并速率的线性叠加。

      其中$w\left( {{d_i},{d_j}} \right)$$P\left( {{d_i},{d_j}} \right)$分别为碰撞频率和聚并效率。

      (1)湍流涡引起的聚并

      ① 由湍流涡引起的碰撞频率

      液相湍流引起的碰撞是通过类比气体分子动理论来获得的[23],对其进行修正后,碰撞频率可以表示为

      其中$\dfrac{{{\alpha _{{\rm{g}},\max}}}}{{{\alpha _{{\rm{g}},\max}} - {\alpha _{\rm{g}}}}}$是考虑到气泡体积的影响而对碰撞频率的修正,Γij为考虑气泡湍动距离的影响而进行的修正[24]

      ② 湍流涡引起的聚并效率

      根据文献[25]提出的液膜理论,气泡的相互作用时间和气泡聚并所需时间的相对大小决定着气泡聚并的概率。文献[25-26]给出了直径为didj的气泡的聚并概率。

      其中,$W{e_{ij}} = \dfrac{{{\rho _{\rm l}}{d_i}{u_{ij}}}}{\sigma }$${\xi _{ij}} = {d_i}/{d_j}$

      (2)由气泡尾涡引起的聚并

      ① 碰撞频率

      鼓泡塔中气泡在尺寸较大时会由球形变成球帽形,并在气泡的尾部形成涡流。此时若其他气泡进入该大气泡的尾涡影响区域,就会被该气泡的尾涡夹带而加速上升。该过程会使大气泡和夹带气泡发生碰撞并发生聚并,文献[16]提出尾涡引起的气泡碰撞频率为:

      由于小气泡的尾涡相对较小,因此这里只考虑大气泡尾涡引起的聚并,所以在碰撞频率中加入因子Θ

      ② 聚并效率

      尾涡引起的聚并效率采用[27]提出的关联式。

      (3)由于气泡上升速度不同引起的聚并

      ① 碰撞频率

      由气泡上升速度不同引起的聚并模型与液相湍流引起的聚并模型相似[23],其表达式为:

      ② 聚并效率

      对于由气泡上升速度不同造成的聚并,文献[14, 24]将聚并效率设置为为0.5.

    • 对于气泡的破碎过程,我们只考虑湍流涡引起的破碎,并且破碎过程为二元破碎,即每个母气泡破碎成两个子气泡,故气泡破碎速率可采用式(26)进行计算[16, 18]

      积分下限${\xi _{\min}} = {\lambda _{\min}}/d$,其中λmin为破碎气泡的最小湍流尺度[16]${\lambda _{\min}} = 11.4 - 31.4{\left( {\nu _{\rm l}^3/\varepsilon } \right)^{1/4}}$。对于积分上限ξmax,文献[28]认为气泡的破碎主要由λ≤3d的涡造成,所以此处ξmax=3。

      (1) 气泡破碎碰撞频率

      目前,一些研究者[17, 28-30]认为大尺寸的湍流涡会把部分能量传递给气泡从而使气泡发生破碎。因此将λd的涡和λ>d的涡与气泡的碰撞频率分别进行计算[19]

      (2)气泡破碎碰撞效率${P_{\rm b}}\left( {d,{f_{\rm v}},\xi } \right)$

      ${P_{\rm b}}\left( {d,{f_{\rm v}},\xi } \right)$为尺寸为λ的湍流涡与直径为d的气泡发生碰撞时以破碎比fv发生破碎的概率。破碎是否能够发生与参与碰撞的湍流涡所携带的能量有关,若该能量大于气泡破碎所需的临界能量,则破碎发生。

      通过假设气-液体系中湍流涡湍动速度可以用正态分布[18]描述可推导得气泡破碎的概率密度函数${P_{\rm b}}\left( {d,{f_{\rm v}},\xi } \right)$

      其中χ为量纲为一湍流涡能量,χc气泡破碎的量纲为一临界能量。

      气泡破碎时必须同时满足能量约束和压力约束[1-4]。当λ>d[5],湍流涡中的能量只能部分作用于气泡并导致气泡的破碎,因此在量纲为一湍流涡能量中引入修正因子Ψ

      气泡破碎时必须同时满足能量约束和压力约束[13, 15-16, 31]。当λ>d[28],湍流涡中的能量只能部分作用于气泡并导致气泡的破碎,因此在无量纲湍流涡能量中引入修正因子Ψ

      在对气泡和液滴的破碎过程进行了实验研究时[32],高速摄像机拍摄到气泡更倾向于破碎成两个不同大小的子气泡。气泡在破碎之前,两部分内的压力不同,会产生气泡内的流体重分布颈部的气体流动会产生动压pd。文献[19]给出了动压pd的计算方法。

      对气泡颈处的受力进行分析最终得到压力破碎因子为

      能量约束条件的推导可见文献[18]

    • 基于文献[33]设置相关实验并进行二维轴对称模拟,采用空气-水作为试验介质,鼓泡塔的相关几何模型见图1(a)。鼓泡塔网格采用ANSYS ICEM CFD划分并对进口处以及壁面处进行加密处理。最终经过网格无关化检验后确定四边形网格数量为9581,计算过程中的时间步长取0.002 s。气体均匀地从底部r/R<0.9的中心区域网格进入(图1(b))。实验的气速范围为0.02~0.1 m/s,覆盖了均匀鼓泡流和湍动鼓泡流。所有计算均在ANSYS Fluent平台上进行,压力-速度方程使用SIMPLE算法处理,动量方程采用QUICK格式,体积分数方程等均采用一阶迎风格式离散,松弛因子均为默认值。模拟中所用到的相间作用力模型均采用UDF加入Fluent求解平台。

      图  1  鼓泡塔反应器几何模型(a)及网格截面(b)

      Figure 1.  Geometry (a) and mesh (b) of bubble column

      鼓泡塔反应器的入口及出口边界条件为速度进口和压力出口,壁面边界条件为无滑移边界。PBE方程的求解采用离散法,同时破碎和聚并模型均采用UDF实现。本文气泡尺寸为$0.5 \le {d_{\rm b}} \le 40.3mm$并将气泡尺寸分为20组,各组气泡大小的计算如式(35)所示。

    • 采用CFD-PBM耦合模型对鼓泡塔反应器内的气含率及轴向液速分布情况进行模拟,气泡聚并及破碎模型通过UDF接口植入ANSYS Fluent中。其中鼓泡塔反应器入口表观气速分别为0.02 m/s和0.1 m/s,并与相应表观气速下的实验数据[33]进行对比。

    • 图2示出了破碎模型预测的量纲为一破碎速率与实验测定的量纲为一破碎速率,图中横坐标表示气泡量纲为一直径(d*),纵坐标表示气泡量纲为一破碎速率,其计算方法如式(36)、(37)所示

      图  2  量纲为一破碎速率预测值与实验值对比

      Figure 2.  Predicted and experimental values of dimensionless breakup rate

      上式中LT分别为长度尺度和时间尺度[34],可通过式(38)、(39)进行计算

      图2中直线代表模型预测的量纲为一破碎速率,数据点为文献[35-37]中的实验测量值。从图2中可以看出,在气泡量纲为一直径较大时略小于实验值,模型预测值与实验值基本吻合。

      图3比较了子气泡分布的模型预测值与实验测定值,其中实验数据来自文献[32]中的实验测定。实线为模型预测的子气泡分布概率密度曲线,灰色直方图为概率密度曲线在不同区间上的积分值,红色直方图为相应区间上的实验测量值。从图中可以看出由于模型采用了压力约束条件,在破碎比fv趋于0时子气泡分布的概率密度(β)曲线也趋于0。同时,从图3中也可以看出,模型预测的结果为气泡更倾向于不均匀破碎,这也与实验测定的结果相吻合[32]

      图  3  子气泡分布预测值与实验值

      Figure 3.  Predicted and experimental values of daughter bubble size distribution

    • 升力是由于液相速度梯度差使得气泡在径向方向上受到的力,其决定了气泡在径向上的运动方向,升力的加入会使大气泡向反应器中心移动,因此会对气含率和轴向液速在径向上的分布产生影响。图5图6示出了添加升力对气含率和轴向液速在径向方向上的分布的影响,其中实线为未添加升力时的气含率、轴向液速径向分布。从图5中可以看出,当模拟不加入升力时,气含率在径向上的分布更加平坦,在鼓泡塔中心处模拟值小于实验测量值,在靠近壁面处远大于实验测量值。图6中可以看出,升力的加入使得径向方向上速度梯度变大,模拟值与实验值更加吻合。总体来说,升力对于气泡的鼓泡塔内流体力学参数的径向分布十分重要,升力的加入使径向上气含率的分布与实验值更加吻合。

      图  5  不同表观气速下气含率径向分布

      Figure 5.  Radial profiles of gas holdup in different superficial gas velocity

      图  6  不同表观气速下液相轴向速度径向分布

      Figure 6.  Radial profiles of liquid axial velocity in different superficial gas velocity

    • 图4中气含率云图和图5中气含率径向分布曲线图可以看出,鼓泡塔内气含率沿径向方向逐渐降低。由于鼓泡塔内存在液相的循环,液体在中心处上升而在壁面处下降,这种液相的循环使得气含率在鼓泡塔的中心较高,而在壁面处较低。从气含率云图与径向分布图可以看出,低表观气速(0.02 m/s)下,气含率的径向分布比较均匀,但在高表观气速(0.1 m/s)下,鼓泡塔中心与壁面处差异较大。

      图  4  不同表观气速下鼓泡塔反应器内气含率及轴向液速分布云图

      Figure 4.  Gas holdup and liquid axial velocity contours of different superficial gas velocity in bubble column

      图5中显示了CFD模拟值与实验值,可以看出,在高表观气速下,鼓泡塔壁面处的模拟值大于实验值;低表观气速(0.02 m/s)下,模拟值在径向方向上略小于实验值。从整体上看,在低表观气速下,目前模型的计算值低估了气含率。

    • 图4中轴向液速云图及图6中轴向液速径向分布图可以看出,鼓泡塔内的轴向液速沿径向方向逐渐降低。由图4可以看出,在径向方向上,中心处的轴向液速最大,沿径向方向,轴向液速逐渐降低并变为负值。此时的负值表示液体的流动方向变为与原来相反的方向。造成这种现象的原因是鼓泡塔内的气相在浮力的驱动下沿鼓泡塔上升,进而造成鼓泡塔内的液相产生这种轴向循环。

      图6中将现有模型的模拟值与实验数据进行了比较。在高表观气速下,也即处于湍动鼓泡态下[9],模拟值与实验值吻合的良好,但在中心处略微低估了轴向液速。在均匀鼓泡态下,模拟值与实验值大致吻合,但在鼓泡塔中心处,模拟值要高于实验值。

      图7为鼓泡塔内不同表观气速下气泡上升速度的径向分布,可以看出,随表观气速增加,气泡的上升速度增加,并且气泡轴向速度径向分布的不均匀性也在增加。结合图6图7也可以看出,气液之间的相对速度变大。这是因为随着表观气速的增加,反应器状态由均匀鼓泡态变为湍动鼓泡态[4],气泡大小分布变的更宽,反应器内大气泡数量增加[16],大气泡尾涡产生的加速效应使得气液之间的相对速度变大。

      图  7  不同表观气速下气泡轴向速度径向分布

      Figure 7.  Radial profiles of bubble axial velocity in different superficial gas velocity

    • 气泡的聚并与破碎速率均与湍动耗散率有关,并且气泡的聚并破碎速率会对局部气泡直径产生影响并进而影响气液相间作用力的大小,因此合理的描述鼓泡塔内的液相湍流耗散率对于流场的准确模拟非常重要。

      图8为不同表观气速下,液相的湍流耗散率在径向上的计算结果。从此图中可以看出,湍动耗散率与表观气速呈正比关系。同时,由于壁面的存在以及流体粘性的作用使得壁面处液相的速度梯度较大,因此壁面处的湍流耗散率较大。

      图  8  不同表观气速下湍动耗散率径向分布

      Figure 8.  Radial profiles of turbulent dissipation rate in different superficial gas velocity

    • 图9为鼓泡塔在H/D=3处截面平均气泡尺寸分布,从图中可以看出,随着表观气速增加,H=1.32 m处的截面平均气泡尺寸分布变得更加低矮,同时大气泡的份额增加。

      图  9  H/D=3处不同气速下气泡尺寸分布

      Figure 9.  Bubble size distribution at H/D=3 in different superficial gas velocity

      图10H/D=3处,r/R=0和r/R=0.5处不同表观气速下的气泡尺寸分布。从图10(a)可以明显看出,随着表观气速增加,鼓泡塔由均匀鼓泡态转变为湍动鼓泡态,鼓泡塔中心位置的气泡尺寸分布变的更宽,大气泡的份额增加。这主要是由于湍动鼓泡态下,较大的表观气速造成鼓泡塔内的气泡数目增加。单位体积内气泡数量的增加使得气泡之间的相互作用更加频繁,在液相湍流的诱导下气泡聚并产生更多的大气泡,而大气泡数量的增加也使得更多的气泡因为气泡尾涡效应和气泡上升速度差效应而发生聚并。另一方面,鼓泡塔内聚并产生的大气泡拥有更大的碰撞截面积,因此会有更多的湍流涡体与其发生相互作用。与此同时,高表观气速下湍动耗散率增加使得湍流涡体的湍动能增加,最终表现为气泡的破碎速率增加。这些因素的综合作用下使得气泡尺寸的分布范围变宽。同时,在升力的作用下,大气泡趋向于向鼓泡塔中心运动,故我们可以看到高表观气速下,r/R=0处气泡尺寸分布较宽且大气泡数量较多。从图10(b)中可以看出,沿径向方向远离鼓泡塔中心,高表观气速与低表观气速的气泡尺寸分布差异在减小。但高表观气速下气泡尺寸分布范围依然大于低表观气速。

      图  10  表观气速下(a)r/R=0和(b)r/R=0.5的气泡尺寸分布

      Figure 10.  Bubble size distribution at (a) r/R=0, and (b) r/R=0.5 in different superficial gas velocity

      图11为表观气速为0.1 m/s时不同径向位置处的气泡尺寸分布,从图中可以非常清晰的看出,沿径向方向靠近鼓泡塔中心,大气泡的份额在增加,小气泡的份额在减少。

      图  11  经向位置处气泡尺寸分布

      Figure 11.  Bubble size distribution at different radial positions

    • 本文将一个新型的破碎模型通过UDF植入ANSYS Fluent 19.0软件中,并将聚并和破碎模型与双流体模型进行双向耦合建立一个通用的CFD-PBM耦合模型,在考虑了不同机理对气泡破碎和聚并的影响以及气液相间作用力后对处于均匀鼓泡态与湍动鼓泡态下的鼓泡塔流体力学参数进行研究。通过数值模拟获得以下结论:

      (1)在对鼓泡塔进行模拟时,升力是否加入对于气含率和轴向液速等流体力学参数的轴向分布有重要影响,升力的加入使模拟数据在径向上的分布与实验测量值吻合得更好。

      (2)在对破碎模型的气泡破碎速率和子气泡分布预测值进行验证之后,将模型在不同表观气速下的气含率、轴向液速模拟值与实验值进行了对比。结果显示CFD-PBM耦合模型的计算数据与实验数据吻合较好,这也证明了该模型可以用于湍动鼓泡态下鼓泡塔的流体力学模拟。

      (3)鼓泡塔内液相湍动耗散率随表观气速的增加而增加。同时,由于壁面的存在以及流体粘性的作用使得壁面处液相的速度梯度较大,因此壁面处的湍流耗散率较大。

      (4)随表观气速增加,鼓泡塔内的气泡数目增加,单位体积内气泡数量的增加使得气泡之间的相互作用更加频繁,气泡聚并速率增加。另一方面,鼓泡塔内聚并产生的大气泡拥有更大的碰撞截面积,因此会有更多的湍流涡体与其发生相互作用。与此同时,高表观气速下湍动耗散率增加使得湍流涡体的湍动能增加,最终表现为气泡的破碎速率增加。这些因素的综合作用下使得气泡尺寸的分布范围变宽。

      符号说明:

      $ {b\left( {d,{f_{\rm v}}} \right)}$−气泡破碎速率,1/s

      CL−升力系数

      CTD−湍流耗散力系数

      CW−壁面润滑力系数

      $ {c\left( {{d_i},{d_j}} \right)}$−气泡聚并核函数,m3/s

      db−气泡直径,m

      dc−临界气泡直径,m

      Eo−Eötvös数

      Eod−根据气泡的变形后的长轴计算出的修正Eötvös数

      Fij−相间作用力,N/m3

      FD−曳力,N/m3

      FL−升力,N/m3

      Fwl−壁面润滑力,N/m3

      FTD−湍流耗散力,N/m3

      fv−气泡破碎体积分数

      g−重力加速度,m/s2

      hb,ij−直径为didj的气泡之间的平均距离,m

      k−湍动能,m2/s2

      lbt,ij−直径为didj的气泡之间的平均相对湍动距离,m

      nW−单位向量

      Re−雷诺数

      ub−气泡速度,m/s

      uslip−气液相间滑移速度,m/s

      α−相含率

      βfvd)−子气泡大小分布

      ρ−密度,kg/m3

      ε−湍流耗散率,m2/s3

      σ−表面张力,kg/m2

      下角标:

      i−第i组气泡

      g−气相

      l−液相

(11)  参考文献 (37) 相关文章 (20)

目录

    /

    返回文章