水中上浮气泡动态特性研究-物理学报.PDFVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 18 (2014) 184701 水中上浮气泡动态特性研究 李帅 孙龙泉 张阿漫 (哈尔滨工程大学船舶工程学院, 哈尔滨 150001) ( 2014 年3 月12 日收到; 2014 年4 月10 日收到修改稿) 针对水中的上浮气泡, 计入表面张力的作用, 然后将气泡边界层内黏性效应分为法向和切向两部分, 其中 附加法向黏性应力通过Young-Laplace 关系考虑; 附加切向黏性力是基于黏性耗散能量等效原理, 引入黏性 修正压力代替. 首先建立了轴对称和三维上浮气泡边界元模型, 将数值结果与理论值和实验值进行对比分析, 有良好的符合度, 验证了数值模型的有效性; 然后针对毫米量级上浮气泡的平衡速度与形态, 讨论了气泡初始 条件、表面张力和黏性对气泡上浮过程中动力学行为的影响; 最后, 提出了一种处理三维上浮气泡融合的数值 方法, 计算结果与实验现象符合良好, 并且能够反映气泡融合后的复杂现象细节. 关键词: 上浮气泡, 黏性, 平衡速度, 融合 PACS: 47.55.dd, 47.15.km, 47.11.Df DOI: 10.7498/aps.63.184701 广到一般情况. 1 引 言 实验方法能够直接获得气泡上浮速度与形态 特征, 是一种有效的研究手段. Duineveld3 研究了 上浮气泡在船舶与海洋工程中发挥着重要的 直径为0.33—1 mm 的气泡在纯净水中的上浮特性, 作用, 比如可燃冰开采中需要收集甲烷气体(气 研究了气泡获得的平衡速度和形状系数以及 泡)、船舶减阻中的微气泡、水下抗爆抗冲击的气泡 之间的关系, 发现能够以稳定形态和速度上浮的 帷幕; 在现代医学上, 采用气泡在血管中携带药物; 最大韦伯数不超过3.2. Raymond 和Rosant4 以及 另外, 气泡还在化工、生物、核能利用等领域都有相 Zenit 和Magnaudet5 都是采用不同液体或者向水 关应用. 所以, 对气泡进行基础性研究对工程应用 中加入不同比例的化学物质来改变流体的密度、黏 能够起到参考作用, 具有重要的学术和工程意义. 性以及表面张力系数等参数, 探究了不同 和 目前对气泡上浮过程的研究方法有理论分析、实验 情况下的气泡阻力系数、形态和上浮速度等物理量 测量和数值模拟, 研究的内容包括气泡的上浮速 的变化规律与内在联系. Wu 和Gharib6 发现直径 度、形态、破裂以及融合特性等. 在0.1—0.2 cm 之间的气泡上浮过程中基本呈圆球 Moore1 计入边界层内黏性作用, 推导了球状 形或椭球形, 当气泡体积不变时, 具有较大上浮速 气泡的上浮阻力为12 (其中 为气泡半径, 度的气泡一般为椭球形; 另外, 直径小于1.5 mm 的 为气泡上浮速度, 为流体的动力黏性系数); 后 气泡上浮路径为直线, 当气泡尺寸更大时, 气泡的 来, Moore2 又在气泡为椭球体的假设下, 推导了 上浮路径将呈“Z” 形或者螺旋形. Amirnia 等7 采 气泡形状系数 与韦伯数 的函数关系, 并预测 用黄原胶和羟甲基化纤维素作为液体介质, 发现气 大于3.745 时, 气泡形态和上浮速度将变得不 泡上浮的平衡速度随着气泡尺寸变化也是连续变 稳定. 不过理论研究的适用范围往往较窄, 不能推