气液鼓泡塔内多气泡上升运动的直接数值模拟

doi:10.16085/j.issn.1000-6613.2025-0613

化工进展 ›› 2025, Vol. 44 ›› Issue (S1): 51-57.DOI: 10.16085/j.issn.1000-6613.2025-0613

气液鼓泡塔内多气泡上升运动的直接数值模拟

谷嘉进(), 陈彩霞(), 夏梓洪

华东理工大学资源与环境工程学院，上海 200237

收稿日期:2025-04-23 修回日期:2025-05-27 出版日期:2025-10-25 发布日期:2025-11-24
通讯作者: 陈彩霞
作者简介:谷嘉进（2000—），男，硕士研究生，研究方向为多相流模拟。E-mail：1934613557@qq.com。
基金资助:
国家自然科学基金(21276085)

Direct numerical simulation of the rising motion of multiple bubbles in a gas-liquid bubbling tower

GU Jiajin(), CHEN Caixia(), XIA Zihong

School of Resource and Environmental Engineering, East China University of Science and Technology, Shanghai 200237, China

Received:2025-04-23 Revised:2025-05-27 Online:2025-10-25 Published:2025-11-24
Contact: CHEN Caixia

摘要/Abstract

摘要：

回顾了直接数值模拟方法模拟鼓泡塔内气泡上升运动的研究进展，介绍了开源程序Basilisk在网格密度上的优化与传统孤立气泡假设不足以准确描述多气泡行为的问题，阐述了Basilisk在气液界面重建与网格加密上的具体设定，构建了网格密度30+/气泡直径的实时网格。在已有孤立气泡模型基础上构建了多气泡模型，采用该模型进行了单孔鼓泡、多孔鼓泡多气泡上升动力学行为的三维直接数值模拟。结果表明，孔口气速增加后，气泡上升方向的稳定性上升，气泡运动轨迹的振荡幅度减小，气泡间合并频率上升且有部分尚未充分加速的气泡被新生成的气泡追上并发生了纵向合并。多孔鼓泡相互影响的临界孔间距约为气泡直径的3倍。孔口布置方式对气泡相互作用的影响不明显，多孔均匀布置与单孔鼓泡相比，前者气泡上升轨迹更加稳定、振荡幅度更小。气泡上升速度主要受气泡直径大小控制，孔间距和孔口布置等因素的影响并不显著。受液相湍流影响，多气泡上升阻力系数大于孤立气泡上升阻力系数。

关键词: 鼓泡塔多气泡上升, 流体体积法, 直接数值模拟, Basilisk, 自适应网格细化技术

Abstract:

The research progress of simulating the upward movement of bubbles in bubble towers by direct numerical simulation methods was reviewed. The optimization of the grid density in the open source software Basilisk and the limitations of the traditional isolated bubble assumption in accurately describing the behavior of multiple bubbles were introduced. The specific settings of Basilisk in gas-liquid interface reconstruction and grid encryption were detailed, achieving a real-time grid with grid densities above 30 grids per bubble diameter. Based on the existing isolated bubble model, a multi-bubble model was constructed. Three-dimensional direct numerical simulation of the dynamic interactions of multiple bubbles from a single nozzle and multiple nozzles was studied by using this model. The results showed that with increases in inlet gas velocities, the zigzag swings of the bubble became less dramatic, the oscillation amplitude of decreases, while the frequency of bubble coalescence increased. Some bubbles that had not been accelerated sufficiently were caught up with the newly generated bubbles. A critical nozzle spacing for the interaction of multiple bubbles was about 3 times the nozzle diameter. The arrangement effects of multiple nozzle on the bubble interactions were insignificant, with reduced bubble swings and decreased oscillation amplitude compared to a single nozzle. Bubble rising velocities were primarily influenced by bubble sizes, with little effect from nozzle spacing or arrangement. Affected by liquid-phase turbulence, the resistance coefficient of multiple bubbles rising was greater than that of isolated bubbles rising.

Key words: multi-bubble rise in a bubbling tower, volume of fluid (VOF), direct numerical simulation, Basilisk, adaptive mesh refinement technique

中图分类号:

O359⁺.1

谷嘉进, 陈彩霞, 夏梓洪. 气液鼓泡塔内多气泡上升运动的直接数值模拟[J]. 化工进展, 2025, 44(S1): 51-57.

GU Jiajin, CHEN Caixia, XIA Zihong. Direct numerical simulation of the rising motion of multiple bubbles in a gas-liquid bubbling tower[J]. Chemical Industry and Engineering Progress, 2025, 44(S1): 51-57.

图/表 11

表1 气液相物性参数

相态	物质	密度/kg·m^-3	黏度/Pa·s	速度/m·s^-1	表面张力/N·m^-1
气相	空气	1	1.81×10^-5	0.1	0.0728
液相	水	1000	1.0×10^-3	—	0.0728

图1 不同密度下的网格细化示意图

图2 气泡上升速度沿床高变化曲线

图3 孔口气速0.15m/s时两孔鼓泡气泡横向合并

图4 孔口气速0.2m/s时两孔鼓泡气泡纵向合并

图5 不同孔口气速下气泡上升速度随时间的变化

图6 不同孔间距下0.43s的气泡运动

图7 不同孔间距下气泡上升速度的变化

图8 不同孔口布置下0.32s的气泡运动

图9 孔口布置对气泡上升速度的影响

图10 阻力系数与雷诺数的关系

参考文献 22

[1]	姚斡维, 郑学波. 基于界面追踪法的受限通道内气泡运移特性研究[J]. 热能动力工程, 2024, 39(6): 113-122.
	YAO Weiwei, ZHENG Xuebo. Research on bubble migration characteristics in confined channel based on interface tracking method[J]. Journal of Engineering for Thermal and Power, 2024, 39(6): 113-122.
[2]	郝文杰, 韩晋玉, 刘瑶, 等. 基于VOF方法的同轴气泡聚并数值模拟研究[J]. 原子能科学技术, 2025: 1-06. .
	HAO Wenjie, HAN Jinyu, LI Yao, et al. Numerical simulation of coaxial bubble merger based on VOF method[J]. Atomic Energy Science and Technology, 2025: 1-06. .
[3]	胡波, 庞明军. 特征时间对剪切稀化流体气泡上浮特性的影响[J]. 化工进展, 2021, 40(5): 2440-2451.
	HU Bo, PANG Mingjun. Effect of characteristic time on bubble floating characteristics of shear-thinning fluid[J]. Chemical Industry and Engineering Progress, 2021, 40(5): 2440-2451.
[4]	CIONCOLINI Andrea, MAGNINI Mirco. Solitary bubbles rising in quiescent liquids: A critical assessment of experimental data and high-fidelity numerical simulations, and performance evaluation of selected prediction methods[J]. Physics of Fluids, 2025, 37(2): 021304.
[5]	CIONCOLINI Andrea, MAGNINI Mirco. Shapes and rise velocities of single bubbles in a confined annular channel: Experiments and numerical simulations[J]. Fluids, 2021, 6(12): 437.
[6]	LIU Zehua, FARSOIYA Palas Kumar, PERRARD Stéphane, et al. Direct numerical simulation of bubble rising in turbulence[J]. Journal of Fluid Mechanics, 2024, 999: A11.
[7]	BHAGA D, WEBER M E. Bubbles in viscous liquids: Shapes, wakes and velocities[J]. Journal of Fluid Mechanics, 1981, 105: 61-85.
[8]	TRIPATHI Manoj Kumar, SAHU Kirti Chandra, GOVINDARAJAN Rama. Dynamics of an initially spherical bubble rising in quiescent liquid[J]. Nature Communications, 2015, 6: 6268.
[9]	GERLACH D, ALLEBORN N, BUWA V, et al. Numerical simulation of periodic bubble formation at a submerged orifice with constant gas flow rate[J]. Chemical Engineering Science, 2007, 62(7): 2109-2125.
[10]	NAHRA Henry K, KAMOTANI Y. Bubble formation from wall orifice in liquid cross-flow under low gravity[J]. Chemical Engineering Science, 2000, 55(20): 4653-4665.
[11]	HUA Jinsong, LOU Jing. Numerical simulation of bubble rising in viscous liquid[J]. Journal of Computational Physics, 2007, 222(2): 769-795.
[12]	GUMULYA M, JOSHI J B, UTIKAR R P, et al. Characteristics of energy production and dissipation around a bubble rising in water[J]. Chemical Engineering Science， 2019, 193: 38-52.
[13]	田航. 基于数值模拟的水下气泡运动特性及其流场压力分析[D]. 沈阳: 沈阳工业大学, 2022.
	TIAN Hang. Motion characteristics of underwater bubbles and analysis of flow field pressure based on numerical simulation[D]. Shenyang: Shenyang University of Technology, 2022.
[14]	MURADOGLU Metin, TRYGGVASON Gretar. Simulations of soluble surfactants in 3D multiphase flow[J]. Journal of Computational Physics, 2014, 274: 737-757.
[15]	LOISY Aurore, NASO Aurore, SPELT Peter D M. Buoyancy-driven bubbly flows: Ordered and free rise at small and intermediate volume fraction[J]. Journal of Fluid Mechanics, 2017, 816: 94-141.
[16]	GUMULYA M, UTIKAR R P, PAREEK V K, et al. Dynamics of bubbles rising in pseudo-2D bubble column: Effect of confinement and inertia[J]. Chemical Engineering Journal, 2021, 405: 126615.
[17]	POPINET Stéphane. Gerris: A tree-based adaptive solver for the incompressible Euler equations in complex geometries[J]. Journal of Computational Physics, 2003, 190(2): 572-600.
[18]	YAN Hongjie, ZHANG Heyang, ZHANG Huimin, et al. Three-dimensional dynamics of a single bubble rising near a vertical wall: Paths and wakes[J]. Petroleum Science, 2023, 20(3): 1874-1884.
[19]	CHEN Caixia, FAN Liang-Shih. Discrete simulation of gas-liquid bubble columns and gas-liquid-solid fluidized beds[J]. AIChE Journal, 2004, 50(2): 288-301.
[20]	LIU Liu, YAN Hongjie, ZHAO Guojian. Experimental studies on the shape and motion of air bubbles in viscous liquids[J]. Experimental Thermal and Fluid Science, 2015, 62: 109-121.
[21]	BROWN Phillip P, LAWLER Desmond F. Sphere drag and settling velocity revisited[J]. Journal of Environmental Engineering, 2003, 129(3): 222-231.
[22]	LOTH E. Quasi-steady shape and drag of deformable bubbles and drops[J]. International Journal of Multiphase Flow, 2008, 34(6): 523-546.

气液鼓泡塔内多气泡上升运动的直接数值模拟

Direct numerical simulation of the rising motion of multiple bubbles in a gas-liquid bubbling tower

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 11

参考文献 22

相关文章 1

编辑推荐

Metrics

本文评价