ISPH中的固壁边界的处理.pdfVIP

第l0卷 第6期 v01．10NO．6 2011年 6月 Jun．20l1 ISPH中的固壁边界的处理 梁朝湘 (广西广播 电视大学 梧州市分校，广西 梧州 543002) 摘 要 ：介绍了ISPH(incompressibleSmoothedParticleHydrodynamics，ISPH)的原理，以及SPH的一些基本原理和 基本公式后 ，侧重介绍 ISPH 的 固壁边界 的处理 方法 。通过模拟 Couette流，并与 WCSPH(WeaklyCompressible SPH)的结果及级数解作对比，验证 了这种处理方法的有效性 。 关键词 ：ISPH；不可压缩流体 ；边界处理；数值模拟 中图分类号 ：TP301 文献标识码 ：A 文章编号 ：1672—7800(2011)06一o021—02 式可描述如下： 0 引言 + ： 0 (1) l0 L／t SPH(SmoothedParticleHydrodynamics)是一种基于 Du = 一 P + 2-V (2) 粒子 的无 网格数值模拟方法 。它最初是 由Lucy、Gingold U ‘ p p 这里 ，p是流体粒子密度；t为时间 为粒子速度；P 和Monaghan用于求解可压缩 的天体物理 问题 ，后来 由 为压力； 为重力加速度 。本文忽略重力加速度的影响。 Monaghan和 Takeda扩展到求解不 可压缩流 问题 。在 SPH 中，粒子的任一函数值可通过对其紧支域 内所有粒 ISPH 的计算方法包括两步：第一步是预估步，是一个 子相对应 的函数值进行加权平均来近似 ，而边界外没有粒 不执行不可压缩特性的显式积分 。此时不考虑压力 的影 子 ，所 以边界上或邻近边界处的粒子的函数值在进行加权 响。这一步，不可压缩特性未能得到满足。流体密度将偏 平均时会被边界截断，导致求解结果错误，这个 问题称之 离初始密度 。因此，第二步修正步用于调节粒子密度 回到 为粒子缺陷问题 。因此，当支持域不完全在研究 区域中 初始密度 。在修正过程 中，压力的影响将被考虑 。由此可 时，就需要进行边界处理 ，如在边界上设置带有特殊性质 获得如下的压力泊松方程： 的粒子 。边界处理对于 ISPH (incompressibleSmoothed cJ0*一 去口-／5r ㈥ ParticleHydrodynamics)来说更加重要 ，因为它是靠边界 粒子 的压力梯度来 阻止粒子穿越边界的。 P 为时刻 +1时的粒子压力 ，P 为第一步所得 的 对 固壁边界的处理主要有 3种方法 ，一种是用边界虚 粒子的流体密度 。△t是时间步长。(3)式的左边使用下 粒子的方法 ，这种边界虚粒子对靠近边界的内部粒子施加 面相应的SPH公式，可 以离散化为线性方程组。求解这 短距离辐射排斥力 ，阻止粒子穿越边界 ；另一种是用假想 个方程组 ，就可以得到压力 。 的边界上和边界外 的虚粒子模拟边界 ，靠这些虚粒子产生 的压力梯度来阻止粒子