油滴撞击油膜层内气泡的变形与破裂过程的数值模拟

周剑宏; 童宝宏; 王伟; 苏家磊

doi:10.7498/aps.67.20180133

摘要

旋转工作的机械零部件和机械设备的润滑系统工作过程中普遍存在着油滴和油膜的碰撞行为，这一行为易引起气泡夹带现象.气泡将对油滴撞击油膜时的运动过程和附壁油膜层的形成质量造成不可忽视的影响.基于耦合的水平集-体积分数方法，对油滴撞击含气泡油膜的行为进行数值模拟研究，考察油膜层内气泡的变形运动过程，分析气泡大小和位置等因素对撞击过程中气泡变形特征参数的影响规律，并探讨气泡破裂的动力学机制.研究表明，随着气泡直径的增大，油滴撞击含气泡油膜后气泡会依次出现自由表面破裂、稳定变形以及油膜内部破裂等现象；直径d=20 m的气泡能较稳定地存在于油膜层内，同时该值也是气泡发生自由表面破裂和油膜内部破裂的临界值.此外，气泡所在位置同样对气泡变形历程有一定影响，气泡越接近油膜表面，其变形量越大；位于油膜底层的气泡会附着在壁面上.在自由表面破裂和油膜内部破裂过程中，气泡破裂是由气-液界面不稳定引起的，表面张力对这两种现象起重要作用；而黏性剪切力对油膜内部破裂现象也有着不可忽视的影响.

关键词:

油滴 /
油膜层 /
气泡破裂 /
数值模拟

Abstract

Impact of oil droplet on oil film usually takes place in the lubrication process of rotating mechanical parts and machinery which can easily lead to bubble entrainment. Bubbles have important influences on the motion process of the oil droplet impacting on the oil film and also on the formation quality of the oil film layer. An oil droplet impacting on the oil film which contains a bubble is simulated numerically based on the coupled level set and the method of determining volume fraction. The bubble deformation process in the oil film during an oil droplet impacting on the oil film is investigated by the simulation method. The influences of the bubble size and the bubble position on the bubble deformation characteristic are also analyzed. The dynamic mechanism of the bubble rupture is discussed. The numerical results show that as the oil droplet impacts on the oil film, the bubble may rupture on the free surface, presenting stable deformation, or rupture in the oil film, which is greatly influenced by the bubble size. When the bubble diameter is in a range between 10 m and 20 m, the bubble deformation becomes more serious with the increase of bubble diameter, and the rupture of bubble on the free surface may occur over time. When the bubble diameters are in a range between 20 m and 30 m, the bubble rupture occurs in a short time after the bubble has reached a maximum deformation, and there is no obvious relationship between the maximum bubble deformation and the bubble diameter. The diameter of 20 m is a critical value for a bubble to rupture on a free surface or inside an oil film, with which a bubble can keep stable in an oil film layer. As the bubble position changes, the bubble deformation process changes correspondingly. Under the same impact conditions, bubbles at the top of the oil film are more likely to deform than those in the center of the oil film. Bubbles at the bottom of the oil film have the smallest total deformation and finally attach to the wall. The bubble rupture is caused by the instability of the gas-liquid interface and the surface tension. The viscous shear force also plays an important role when the bubble rupture takes place in the oil film.

Keywords:

作者及机构信息

1.
安徽工业大学机械工程学院, 马鞍山 243032;

2.
合肥工业大学摩擦学研究所, 合肥 230009

通信作者: 童宝宏, bhtong@ahut.edu.cn

基金项目: 国家自然科学基金（批准号：51475135）、清华大学摩擦学国家重点实验室开放基金（批准号：SKLTKF17B01）和安徽工业大学研究生创新基金（批准号：2015040）资助的课题.

Authors and contacts

1.
School of Mechanical Engineering, Anhui University of Technology, Maanshan 243032, China;

2.
Institute of Tribology, Hefei University of Technology, Heifei 230009, China

Corresponding author: Tong Bao-Hong, bhtong@ahut.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant No. 51475135), the Tribology Science Fund of State Key Laboratory of Tribology, Tsinghua University, China (Grant No. SKLTKF17B01), and the Graduate Innovation Foundation of Anhui University of Technology, China (Grant No. 2015040).

参考文献

[1]	Peduto D 2015 Ph. D. Dissertation (Karlsruhe: Karlsruhe Institute of Technology)
[2]	Rioboo R, Bauthier C, Conti J, Vou M, de Coninck J 2003 Exp. Fluids 35 648
[3]	Okawa T, Shiraishi T, Mori T 2006 Exp. Fluids 41 965
[4]	Guo J H, Dai S Q, Dai Q 2010 Acta Phys. Sin. 59 2601(in Chinese) [郭加宏, 戴世强, 代钦 2010 物理学报 59 2601]
[5]	Song Y C, Ning Z, Sun C H, Lyu M, Yan K, Fu J 2013 Chin. J. Theor. Appl. Mech. 45 833(in Chinese) [宋云超, 宁智, 孙春华, 吕明, 阎凯, 付娟 2013 力学学报 45 833]
[6]	Liang G T, Guo Y L, Shen S Q 2013 Acta Phys. Sin. 62 024705(in Chinese) [梁刚涛, 郭亚丽, 沈胜强 2013 物理学报 62 024705]
[7]	Bisighini A 2010 Ph. D. Dissertation (Bergamo: University of Bergamo)
[8]	Rein M 1996 J. Fluid Mech. 306 145
[9]	Blanchette F, Bigioni T P 2009 J. Fluid Mech. 620 333
[10]	Chen X, Mandre S, Feng J J 2006 Phys. Fluids 18 051705
[11]	Ray B, Biswas G, Sharma A 2010 J. Fluid Mech. 655 72
[12]	Thoraval M J, Li Y, Thoroddsen S T 2016 Phys. Rev. E 93 033128
[13]	Pumphrey H C, Elmore P A 1990 J. Fluid Mech. 220 539
[14]	Pumphrey H C, Crum L A, Bj rn L 1989 J. Acoust. Soc. Am. 85 1518
[15]	Oguz H N, Prosperetti A 1990 J. Fluid Mech. 219 143
[16]	Zou J, Ji C, Yuan B G, Ren Y L, Ruan X D, Fu X 2012 Phys. Fluids 24 057101
[17]	Wang A B, Kuan C C, Tsai P H 2013 Phys. Fluids 25 1518
[18]	Deng Q, Anilkumar A V, Wang T G 2007 J. Fluid Mech. 578 119
[19]	Thoroddsen S T, Thoaval M T, Takehara K, Etoh T G 2012 J. Fluid Mech. 708 469
[20]	Sigler J, Mesler R 1990 J. Colloid Interface Sci. 134 459
[21]	Saylor J R, Bounds G D 2012 AIChE J. 58 3841
[22]	Sussman M, Puckett E G 2000 J. Comput. Phys. 162 301
[23]	Guo Y L, Wei L, Liang G T, Shen S Q 2014 Int. Commun. Heat Mass. 53 26
[24]	Wang Z, Li Y, Huang B, Gao D M 2016 J. Mech. Sci. Technol. 30 2547
[25]	Ohta M, Kikuchi D, Yoshida Y, Sussman M 2011 Int. J. Multiphase Flow 37 1059
[26]	Fan W, Qi T, Sun Y W, Zhu P, Chen H 2016 Chem. Eng. Technol. 39 1895
[27]	Brackbill J U, Kothe D B, Zemach C 1992 J. Comput. Phys. 100 335
[28]	Cossali G E, Marengo M, Coghe A, Zhdanov S 2004 Exp. Fluids 36 888
[29]	Feonychev A I 2007 J. Eng. Phys. Thermophys. 80 961

施引文献

[1]	Peduto D 2015 Ph. D. Dissertation (Karlsruhe: Karlsruhe Institute of Technology)
[2]	Rioboo R, Bauthier C, Conti J, Vou M, de Coninck J 2003 Exp. Fluids 35 648
[3]	Okawa T, Shiraishi T, Mori T 2006 Exp. Fluids 41 965
[4]	Guo J H, Dai S Q, Dai Q 2010 Acta Phys. Sin. 59 2601(in Chinese) [郭加宏, 戴世强, 代钦 2010 物理学报 59 2601]
[5]	Song Y C, Ning Z, Sun C H, Lyu M, Yan K, Fu J 2013 Chin. J. Theor. Appl. Mech. 45 833(in Chinese) [宋云超, 宁智, 孙春华, 吕明, 阎凯, 付娟 2013 力学学报 45 833]
[6]	Liang G T, Guo Y L, Shen S Q 2013 Acta Phys. Sin. 62 024705(in Chinese) [梁刚涛, 郭亚丽, 沈胜强 2013 物理学报 62 024705]
[7]	Bisighini A 2010 Ph. D. Dissertation (Bergamo: University of Bergamo)
[8]	Rein M 1996 J. Fluid Mech. 306 145
[9]	Blanchette F, Bigioni T P 2009 J. Fluid Mech. 620 333
[10]	Chen X, Mandre S, Feng J J 2006 Phys. Fluids 18 051705
[11]	Ray B, Biswas G, Sharma A 2010 J. Fluid Mech. 655 72
[12]	Thoraval M J, Li Y, Thoroddsen S T 2016 Phys. Rev. E 93 033128
[13]	Pumphrey H C, Elmore P A 1990 J. Fluid Mech. 220 539
[14]	Pumphrey H C, Crum L A, Bj rn L 1989 J. Acoust. Soc. Am. 85 1518
[15]	Oguz H N, Prosperetti A 1990 J. Fluid Mech. 219 143
[16]	Zou J, Ji C, Yuan B G, Ren Y L, Ruan X D, Fu X 2012 Phys. Fluids 24 057101
[17]	Wang A B, Kuan C C, Tsai P H 2013 Phys. Fluids 25 1518
[18]	Deng Q, Anilkumar A V, Wang T G 2007 J. Fluid Mech. 578 119
[19]	Thoroddsen S T, Thoaval M T, Takehara K, Etoh T G 2012 J. Fluid Mech. 708 469
[20]	Sigler J, Mesler R 1990 J. Colloid Interface Sci. 134 459
[21]	Saylor J R, Bounds G D 2012 AIChE J. 58 3841
[22]	Sussman M, Puckett E G 2000 J. Comput. Phys. 162 301
[23]	Guo Y L, Wei L, Liang G T, Shen S Q 2014 Int. Commun. Heat Mass. 53 26
[24]	Wang Z, Li Y, Huang B, Gao D M 2016 J. Mech. Sci. Technol. 30 2547
[25]	Ohta M, Kikuchi D, Yoshida Y, Sussman M 2011 Int. J. Multiphase Flow 37 1059
[26]	Fan W, Qi T, Sun Y W, Zhu P, Chen H 2016 Chem. Eng. Technol. 39 1895
[27]	Brackbill J U, Kothe D B, Zemach C 1992 J. Comput. Phys. 100 335
[28]	Cossali G E, Marengo M, Coghe A, Zhdanov S 2004 Exp. Fluids 36 888
[29]	Feonychev A I 2007 J. Eng. Phys. Thermophys. 80 961

[1]	聂腾飞, 徐强, 罗欣怡, 洪奥越, 曹泽水, 郭烈锦. 水分解中氧气泡生长动力学研究. 物理学报, 2025, 74(10): 108201. doi: 10.7498/aps.74.20250014
[2]	张超, 布龙祥, 张智超, 樊朝霞, 凡凤仙. 丁二酸-水纳米气溶胶液滴表面张力的分子动力学研究. 物理学报, 2023, 72(11): 114701. doi: 10.7498/aps.72.20222371
[3]	赵豪, 吴志豪, 胡晓红, 凡凤仙, 苏明旭. 外加液滴条件下固体细颗粒声凝并特性. 物理学报, 2023, 72(6): 064702. doi: 10.7498/aps.72.20221912
[4]	潘文韬, 文琳, 李姗姗, 潘振海. Y型微通道内气泡非对称破裂行为的数值研究. 物理学报, 2022, 71(2): 024701. doi: 10.7498/aps.71.20210832
[5]	林茜, 谢普初, 胡建波, 张凤国, 王裴, 王永刚. 不同晶粒度高纯铜层裂损伤演化的有限元模拟. 物理学报, 2021, 70(20): 204601. doi: 10.7498/aps.70.20210726
[6]	叶欣, 单彦广. 疏水表面振动液滴模态演化与流场结构的数值模拟. 物理学报, 2021, 70(14): 144701. doi: 10.7498/aps.70.20210161
[7]	梁煜, 关奔, 翟志刚, 罗喜胜. 激波汇聚效应对球形气泡演化影响的数值研究. 物理学报, 2017, 66(6): 064701. doi: 10.7498/aps.66.064701
[8]	陈大伟, 王裴, 孙海权, 蔚喜军. 爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究. 物理学报, 2016, 65(2): 024701. doi: 10.7498/aps.65.024701
[9]	周祥曼, 张海鸥, 王桂兰, 柏兴旺. 电弧增材成形中熔积层表面形貌对电弧形态影响的仿真. 物理学报, 2016, 65(3): 038103. doi: 10.7498/aps.65.038103
[10]	马理强, 苏铁熊, 刘汉涛, 孟青. 微液滴振荡过程的光滑粒子动力学方法数值模拟. 物理学报, 2015, 64(13): 134702. doi: 10.7498/aps.64.134702
[11]	王勇, 林书玉, 莫润阳, 张小丽. 含气泡液体中气泡振动的研究. 物理学报, 2013, 62(13): 134304. doi: 10.7498/aps.62.134304
[12]	刘邱祖, 寇子明, 韩振南, 高贵军. 基于格子Boltzmann方法的液滴沿固壁铺展动态过程模拟. 物理学报, 2013, 62(23): 234701. doi: 10.7498/aps.62.234701
[13]	陈石, 王辉, 沈胜强, 梁刚涛. 液滴振荡模型及与数值模拟的对比. 物理学报, 2013, 62(20): 204702. doi: 10.7498/aps.62.204702
[14]	耿少飞, 唐德礼, 邱孝明, 聂军伟, 于毅军. 霍尔漂移对阳极层霍尔等离子体加速器电离效率的影响. 物理学报, 2012, 61(7): 075210. doi: 10.7498/aps.61.075210
[15]	刘小平, 范广涵, 张运炎, 郑树文, 龚长春, 王永力, 张涛. 量子阱垒层掺杂变化对双波长LED调控作用研究. 物理学报, 2012, 61(13): 138503. doi: 10.7498/aps.61.138503
[16]	马理强, 刘谋斌, 常建忠, 苏铁熊, 刘汉涛. 液滴冲击液膜问题的光滑粒子动力学模拟. 物理学报, 2012, 61(24): 244701. doi: 10.7498/aps.61.244701
[17]	马理强, 常建忠, 刘汉涛, 刘谋斌. 液滴溅落问题的光滑粒子动力学模拟. 物理学报, 2012, 61(5): 054701. doi: 10.7498/aps.61.054701
[18]	张运炎, 范广涵. 不同掺杂类型的GaN间隔层和量子阱垒层对双波长LED作用的研究. 物理学报, 2011, 60(1): 018502. doi: 10.7498/aps.60.018502
[19]	耿少飞, 唐德礼, 赵杰, 邱孝明. 圆柱形阳极层霍尔等离子体加速器的质点网格方法模拟. 物理学报, 2009, 58(8): 5520-5525. doi: 10.7498/aps.58.5520
[20]	杨玉娟, 王锦程, 张玉祥, 朱耀产, 杨根仓. 采用多相场法研究三维层片共晶生长的厚度效应. 物理学报, 2008, 57(8): 5290-5295. doi: 10.7498/aps.57.5290

计量

文章访问数: 7904
PDF下载量: 318
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板