基于无网格法的水中气泡上升运动数值模拟.pdfVIP

中国工程热物理学会 iiic体机械 学术会议论文 编号：067006 基于无网格法的水中气泡． 上升运动数值模拟 孙中国席光王焕然项利锋 西安交通大学能源与动力工程学院710049 电话：029—mai 1．xjtu．edu．cn l：xiguang@IIlai 摘要：本文对无网格法的一种——移动粒子半隐式法(MPS)进行了研究，建立了表面张力和气液两 相流等若干粒子作用模型，编程实现了二维情况下水中气泡自由上升运动的数值模拟；对计算结果进 行分析并与相关实验结果进行了定性比较，模拟结果显示了移动粒子半隐式法在模拟自由表面和气液 ‘ 两相流等涉及到大变形问题的有效性． ， 关键词：移动粒子半隐式法；表面张力；气液两相流 1．引言 气液两相流动在工业生产和日常生活中是极其常见的流动现象．但是由于两相流动 常常伴随着自由表面较剧烈的位移和变形，因此使得传统基于网格的数值模拟方法在模 拟上述流动时面临着网格的动态调整和重组等困难，可能会造成计算误差增大和计算量 的急剧增加，对于较复杂的流动现象甚至无法收敛。 ． 移动粒子半隐式法作为一种新的无网格方法消除了网格对计算的影响，连续介质的 离散化处理使得该方法可以很容易的模拟各种复杂几何大变形，并且不会对计算量产生 任何影响，本文采用该方法对简单的气液两相流问题——水中静止气泡自由上升运动进 行建模和数值模拟。 2．数值方法 压缩流动的无网格数值解法。其基本思想是在求解区域内部及其边界上布置一系列离散 的粒子来代表宏观流体，每个粒子都携带其相应的位置、速度等流动信息，在Lagrange 框架下，通过粒子间相应函数关系的求解来追踪各流体粒子在各个时间层的位置，速度 等信息，以时间层推进的方式，来获得整个流场的动态流动信息。 2．1流动控制方程： ．控制方程满足质量守恒和动量守恒： 坐+刃．u：o (1)(1) —D二t+刀’u。O p{兰=一VP+∥2u+pg+丫 (2) 其中u为流体速度矢量，P为流体密度，P为压力，／z为流体动力粘性系数。g为作用 在流体单位质量上的外力矢量，丫为表面张力矢量。 2．2粒子作用模型 每个粒子通过核函数与其邻点粒子发生相互作用，核函数为： 7 fr “‘ 帕)={等-1 (3) 【0 ，≥‘ 式(3)中r是粒子f和／之间的距离，，=I‘一‘I，re称为粒子作用半径。粒子仅与位 于其粒子作用半径兀范围内的邻点粒子发生相互作用。定义f粒子的粒子数密度为： (以)，=∑w(Il一‘1) (4) 流体密度与粒子数密度成正比。对于不可压缩流体，流体密度保持定常相当于其粒 子数密度维持一个常数矿。 设谚=矽(I)，表示位于I处f粒子的物理参数标量值，则f粒子的梯度矢量和拉普 拉氏算子的粒子作用形式表示为：