用格子boltzmann方法模拟液滴合并过程 simulation of droplets coalescence process by lattice boltzmann method.pdfVIP

第6卷第3期 热 科 学 与 技 术 V01．6No．3 2007年9月 ofThermaland Journal Science Sep．2007 Technology 文章编号：1671-8097(2007)03—0198—06 用格子Boltzmann 方法模拟液滴合并过程 李维仲， 李 爽 (大连理工大学能源与动力学院，辽宁大连 116024) 摘要：用格 Boltzmann方法中的伪势模型对两个液滴的合并过程进行了数值模拟。详细研究了两液滴能否 合并临界尺寸和液滴合并过程中液桥的形成与演化的关系，并研究了表面张力对合并速度的影响。研究结果 发现当两个液滴之间距离小于2倍界面厚度时，两个液滴在不受外力的作用下能够自动合并；液桥的宽度与 演化时间有一定的指数关系；表面张力越大，合并速度越快，这个结果与前人的理论预测和实验结果一致。 关键词：格子Boltzmann方法；伪势模型；两相流；液滴合并 中图分类号：0359 文献标识码：A 0 引 言 1 格子Boltzmann模型 在冷却塔的设计、喷雾系统的设计以及自然 本文使用二维九速正方形格子模型，即 界中雨滴的形成过程中，对于液滴合并动力学的 D2Q9模型，如图1所示，图中每个格点的速度离 研究有显著的重要性。因此，在过去的几十年，液 散成9个方向，即速度矢量为 滴的合并过程成为大量的理论研究、实验研究和 eo一(O，O) 数值模拟的主要对象。格子Bohzmann[13方法是e，=(cos0。，sin0i)c(护i一(i一1)7c／2，i=1，…， 发展迅速的一种计算流体力学新方法。它与以宏 4) 观连续方程为基础的传统计算流体力学方法不 ei一~／2(cos Oi，sin 同，格子Boltzmann方法是基于流体微观模型和 7r／4，i一5，…，8) 介观动理论方法。由于格子Boltzmann方法的微 观粒子特性，它可以方便地描述相间的相互作用。 格子Boltzmann方法在多相流方面的研究经历了 三个主要阶段：1991年Gunstensen等提出的着色 ，＼ ‘／ 模型[2_3]，1993年Shan和Chen等提出的基于均场 理论的伪势模型[4曲3和Swift等提出的自由能模 ∥＼| 型[10|。本文使用Shan和Chen等提出的伪势模型 模拟了液滴合并过程，寻找液滴合并条件以及液 图1 D2Q9模型示意图 Schematicof model 桥形成过程，并通过对不同液滴合并的模拟，掌握 Fig．1 drawingD2Q9 不同介质或不同表面张力对该过程的影响。为相 粒子的时空演化方程为 关工业设备的设计和运行提供理论依据和参考。 ^(z+e，，t+1)一fi(z，z)一 收稿日期：2007—07—03；修回日期：2007—08—11． 基金项目：国家自然科学基金资助项目． 作者简介：李维仲(1956一)，男．博士，教授．博士生导师 万方数据 第3期 李维仲等：用格子Boltzmann方法模拟液滴合并过程 199 一丛生尘j型，i一。～8 (