等直径微液滴碰撞过程的改进光滑粒子动力学模拟-物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 22 (2013) 224701 等直径微液滴碰撞过程的改进光滑粒子动力学模拟* 12 1† 1 蒋涛 陆林广 陆伟刚 1) ( 扬州大学数学科学学院, 扬州大学水利科学与工程学院, 扬州 225002 ) 2) ( 西北工业大学应用数学系, 西安 710129 ) ( 2013年5月1 日收到; 2013年8月15 日收到修改稿) 运用一种改进光滑粒子动力学(SPH) 方法模拟了相溶和不相溶两种情况下的等直径微液滴碰撞动力学过程. 为提高传统SPH 方法的数值精度和稳定性, 采用一种不涉及核导数计算的核梯度改进形式; 为处理液滴界面张力采 用修正的van der Waals 表面张力模型. 通过模拟牛顿液滴碰撞聚并变形过程并与相关文献或试验结果进行对比, 验 证了改进SPH 方法模拟微液滴碰撞过程的可靠性. 随后, 研究了基于van der Waals 模型相溶聚合物微液滴碰撞聚并 变形过程及不相溶微液滴碰撞后的反弹、分离过程, 讨论了碰撞过程中碰撞速度、碰撞角度、密度比等参数对碰撞 变形过程的影响, 分析了流体桥、旋转角度等因素的变化情况. 关键词: 光滑粒子动力学, 微液滴, 聚合物液滴, 碰撞 PACS: 47.11.j, 44.10.+i, 47.85. g DOI: 10.7498/aps.62.224701 滴碰撞变形过程时存在诸多困难, 更不易将其应用 1 引言 到聚合物微液滴碰撞变形过程的模拟. 此外, 由文 献 [14, 15] 知, 对聚合物液滴的碰撞动力学行为进 12 微液滴 (主要指较小尺度下的液滴) 碰撞变 行实验研究亦存在诸多困难, 且相关的数值模拟研 形过程是工业上产品制造过程和日常生活中常见 究也鲜有报道. 的现象, 它包含复杂的流动过程, 如液滴的融合、 因此, 本文采用一种无网格粒子方法 1620 旋转和分离等现象38 , 其变形过程受多种因素的 20 —— 光滑粒子动力学(SPH) 方法 来对相溶聚合 影响, 如碰撞速度、密度比、两液滴直径比及碰撞 物微液滴的碰撞聚并过程以及不相溶微液滴的碰 角度等. 微液滴碰撞的流动机理研究对化工工程、 撞分离过程进行模拟研究. SPH 方法是一种纯无网 环境工程以及石油产业等有着重要意义. 格粒子方法, 它最初提出是为求解三维无边界空间 目前, 已有很多学者对两个相溶微液滴 的碰撞动力学行为进行了实验 36910 或数值 的行星运动问题. 在 1994 年开始被应用于自由表 21 模拟 1241113 研究. 大部分实验研究中考虑了 面流问题 , 后来逐渐被应用于黏性流动、黏弹 真空条件下液滴的碰撞动力学行为, 数值模拟中研 性流动、多相流动等问题2225 . 与网格类方法相 究了液滴变形过程中表面张力的处理问题. 然而,