气泡%2f液滴Marangoni迁移数值模拟的研究.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：096030 气泡/液滴Marangoni 迁移数值模拟研究 1 1 2 李震东 ，赵建福 ，秦文韬 1 （中国科学院微重力重点实验室，中国科学院力学研究所，北京 100190) 2 河南濮阳供电公司，濮阳 457000 ） （ （Tel: 010E-mail: jfzhao@imech.ac.cn ） 摘要: 浸没在液体内的气泡或液滴的迁移运动现象广泛存在于诸多工程应用中，是流体力学研究的基本 问题之一。在空间微重力环境中浮力效应的影响被极大削弱，因此由界面温度梯度导致的表面张力分 布不均匀所引发的 Marangoni 效应对液滴、气泡的运动具有重要影响。本文利用 Level Set 方法，模拟 了微重力情况下气泡/液滴的 Marangoni 迁移现象，数值模拟结果表明，由于非线性热对流效应的影响， 无量纲迁移速度随着 Ma 数的增大呈现逐渐减小的趋势。 关键词: Marangoni 迁移，气泡，液滴，Level Set 方法 前言 液滴和气泡在流体介质中的运动对许多物理过程和工业生产活动有着非常重要的 影响。液滴和气泡的运动有多种驱动力，如浮力、界面温度梯度和/或界面浓度梯度引起 的非均匀表面张力等。在空间微重力环境中由于浮力效应的影响被极大削弱，后两种类 型的驱动力在液滴、气泡的运动中具有重要影响。本文研究的是由温度梯度引起的流体 颗粒（液滴和气泡）的运动，即 Marangoni 迁移或热毛细迁移。 由于界面张力会随着温度的不同而变化，若流体颗粒所浸没的母液中存在温度梯 度，作用在流体颗粒表面的界面张力将不再均匀，这种界面张力的不均匀就会在界面上 形成切向作用力，从而驱动周围流体发生流动。流体颗粒内部形成对流胞元，外部则形 成绕流，产生流体颗粒的 Marangoni 迁移现象。 迄今为止，许多学者对这一现象进行了深入研究，对于 Re 数和 Ma 数都非常小的情 况（Re1 ，Ma1 ），Yang 等[1]忽略非线性对流项的影响，得到了描述流体颗粒的 Marangoni 迁移速度的 YGB 公式： 2U V (1) YGB (2 + 3μ 2 μ )(2 + k2 k ) 1 1 对于大 Re 数和 Ma 数的情况，由于复杂的非线性效应，目前尚没有成熟的理论，多依赖 [2] 空间实验 。不过，由于空间实验耗时长、费用大，对该现象的数值模拟研究越来越受 到各国学者的重视 [3] 。但是，对该问题的数值模拟不可避免的会涉及到两相流动，其重 点是要对相界面的运动进行精确追踪。同时，由于相界面会产生迁移、变形甚至破碎拓 扑结构的变化，而对气泡的模拟还要考虑两相间悬殊的物性参数比，使得相应的数值模 * 国家自然科学基金重点项目 资助 拟相当困难。因此，目前绝大部分的模拟研究将液滴、气泡的形状假设为球形且在运动 过程中保持不变，对气泡的模拟也限制物性参数比在相对较小的范围内，与实际情况差 距较大。 在运动界面追踪方法方面，目前已有多种数值方法，如 MAC