一种改进的SPH方法及在不可压缩流模拟中的应用.pdfVIP

一种改进的SPH方法及其在不可压缩流模拟中的应用牛 周浩，汤文辉，徐志宏，冉宪文 国防科学技术大学理学院工程物理研究所，长沙410073 本文提出一种改进的SPH方法，其基本思想是在已知每个粒子状态的条件下，采用广 义有限差分方法 GFDM 计算每个粒子有关物理量的空间一阶和二阶导数 如速度散度、 压力梯度和速度的拉普拉斯算子等 ，然后根据控制方程更新粒子的位置、速度、压力等物 理量。即使粒子分布不均匀，如果某物理量为线性分布，本方法也能够精确求解其在各个粒 子处的空间一阶和二阶导数。与传统SPH方法和有限粒子法 FPM 相比，本方法中不需 要核函数，且过程中没有数值积分，因而过程更加简单，易于编程实现。二维的Poiseuille 流和Couette流的模拟结果显示，本方法精度与稳定性都优于SPH。本方法也适合于静态问 题的模拟。 在模拟Poiseuille流和Couette流时，当雷诺数为1．25X10-2，SPH方法和改进方法都能 给出较好结果，但SPH方法的相对误差约为改进方法的6倍。当雷诺数为1．25X103时，SPH 方法由于不稳定而无法求解，但改进方法仍然能给出精确结果。图1和图2为雷诺数为1250 内。 ～ ． 一 一 ” Ⅶ ¨ ： ^l、 “ ∞ ∞ ，，，．．，：：＼、、一 ∞ 一0 ?∥|l| 一 了、。一、、。、弋一 鬲、≥一 一～么一 ∥∥≯、|一…乏＼一 between between Fig．1．comparison velocityprofilesFig．2．comparisonvelocityprofiles obtainedournew methodandobtained new method our and using improved using improved solutionsfor2DPoiseuilleflow for2D solutionsCouetteflow analytical analytical 《Re 1250 Re 1250 改进方法也适合于静态问题的模拟。图3为一杯水在重力季 2，-9．7046，0 ms-2作用 下的压力分布图。 fieldofafluidina Fig．3．hydrostaticpressure tankunder gravity