欢迎访问过程工程学报, 今天是

过程工程学报 ›› 2018, Vol. 18 ›› Issue (6): 1178-1186.DOI: 10.12034/j.issn.1009-606X.217435

• 流动与传递 • 上一篇    下一篇

流经突然扩张的微通道的单气泡界面形态演变的数值模拟

鲍允州,吴忱韩,张莹*   

  1. 南昌大学机电工程学院,江西 南昌 330031
  • 收稿日期:2017-12-12 修回日期:2018-03-16 出版日期:2018-12-22 发布日期:2018-12-19
  • 通讯作者: 张莹 yzhan2033@163.com
  • 基金资助:
    毛细管内脉动两相流动的演化规律和机理研究;南昌大学研究生创新专项

Numerical simulation of interface evolution of single bubble rising through sudden expansion channel

Yunzhou BAO, Chenhan WU, Ying ZHANG*   

  1. School of Mechanical and Electrical Engineering, Nanchang University, Nanchang, Jiangxi 330031, China
  • Received:2017-12-12 Revised:2018-03-16 Online:2018-12-22 Published:2018-12-19
  • Contact: Ying ZHANG yzhan2033@163.com

摘要: 采用界面追踪方法(FTM)研究了狭窄通道内受浮力上升的单气泡经过突然扩张通道时的运动形态,考察了狭窄通道和宽阔通道的通道宽度比及Eotvos数(Eo)和Morton数(Mo)对气泡形态和上升速度的影响. 结果表明,在狭窄通道、扩张部位和宽阔通道内气泡形态呈不同特性,气泡在离开扩张处时高宽比最大. 随通道宽度比增大,气泡最大高宽比减小. Eo数较小时,减小Mo数,气泡高宽比减小,气泡上升速度增大;增大Eo数,气泡上升时形变更严重. 改变初始时刻气泡相对于狭窄通道的位置对气泡形态有影响,初始位置靠近单侧壁面时,气泡横向变形较严重,在两侧壁面之间摆动上升.

关键词: 浮升气泡, Front Tracking Method方法, 界面形态, 数值模拟, 流体动力学

Abstract: A single bubble driven by buoyance rising through a sudden expansive channel was simulated by Front Tracking Method (FTM). Firstly, the influence of the width ratio of the narrow channel and the wide channel on shape and rising velocity of bubble was investigated. It was found that the morphology of bubble interface presented different characteristics in the three parts, i.e., the narrow channel, the expansion part and the wide channel, respectively. When the bubble left the expansion, the surrounding fluid flowed from both sides of the bubble into the bottom area, causing both sides of the bubble to be squeezed, the bubble reached the highest aspect ratio. With the increase of the channel width ratio, the flow field of the bubble was less inhibited by both sides of the wall, and the maximum aspect ratio of the bubble decreased. When Wn/We reached over 4/6, the influence of expansion segment can be ignored. Then the influence of Eotvos number Eo and Mo number were analyzed. At low Eo, when the Mo number reduced, the bubble was more likely to deform, the aspect ratio of the bubble decreased and the bubble?s rising speed increased. When the Eo number increased, the deformation became more serious when the bubble rised. Changing the position of the bubble relative to the narrow channel at the initial moment also affected the shape of the bubble. When the initial position of the bubble was close to the single side wall, due to the low flow velocity of the surrounding fluid near the wall, the bubbles would be stretched horizontally while moving, and the bubble deformed badly in the horizontal direction. Asymmetrical pressure distributions occurred on the left and right sides of the bubble level. The pressure on the side of the bubble near the wall surface was lower, making the movement path of the bubble in the narrow channel appear to swing left and right. After leaving the expansion part, the movement of bubbles gradually stabilized due to the absence of wall constraint on the bubbles. When the bubble was far away from the expansion part, the influence of the wall surface decreased. After 2.3 times from the expansion part, it can be considered that the bubble enters a stable rising state.

Key words: buoyancy-driven bubble, Front Tracking Method, bubble shape, direct numerical simulation, fluid dynamics