一种新的半隐式法压力缓解方法.docxVIP

一种新的半隐式法压力缓解方法 目前，船舱和海洋应用的数值模拟方法主要包括网格法和无网格法。与网格法相比，无网格法的优点是建模简单，适合复杂边界条件和复杂流动，并采用lagore描述来消除对边界条件的离散。mps（move-semi-iii）方法是一种新型的非网络法。在水流数值的模拟过程中，这种方法将水分解成颗粒状，这对水流的流动和自由表面的模拟具有独特的优势。因此，在提出后，由于它的快速发展，mps和其他无网格方法的实际计算中存在了严重的压力振荡问题，这严重影响了该方法的发展。sueyoshi等人重视pms法本身的压力衰减，并提出了四种不同的方法来缓解压力扰动，对提高压力得到了一定的作用，但收效甚微。因此，由于采用人工压力收缩法，提出了一种新的压力缓解法。潘徐杰等人对振动中的压力扰动进行了研究，比较了核心函数和自由表面的条件，并使用了平均面积时间法进行了研究。 本文根据MPS的计算特性, 将数密度守恒和速度散度为零条件共同作为不可压缩流体的保证条件, 并在原压力计算的泊松方程中加入了速度散度为零约束条件的作用项, 以此缓解压力的振荡性, 并在计算完成后再加以时间平均的优化.比较溃坝和船舱横摇晃荡的计算结果可知, 压力振荡现象缓解较好, 并且计算量也未增加, 计算效率较高, 该方法值得推广与应用. 1 初始数密度模型 MPS法利用Lagrangian形式的Navier-Stokes方程作为控制方程. 质量守恒方程: dρdt=0(1)dρdt=0(1) 动量守恒方程: 式中:ρ为流体密度;t为时间;u为速度矢量;P为压力;ν为运动黏度;F为体积力.控制方程中的一阶微分项利用梯度模型来离散, 如下式所示: 控制方程中的二阶微分 (梯度) 项利用拉普拉斯模型来离散, 如下式所示: 式中:f为粒子物理量;d为空间维数;n0为在模拟初始求得的粒子初始数密度;xi和xj为粒子i和粒子j的坐标;w为核函数;λ为参数, λ=∑j≠iw(|xj?xi|)|xj?xi|2∑j≠iw(|xj?xi|)(5)λ=∑j≠iw(|xj-xi|)|xj-xi|2∑j≠iw(|xj-xi|)(5) 粒子数密度表示在核函数的作用范围内, 某粒子周围邻居核函数数值之和.粒子i的数密度为 ni=∑j≠iw(|xi?xj|)(6)ni=∑j≠iw(|xi-xj|)(6) 利用粒子数密度模型来计算式 (3) 、 (4) 中梯度模型和拉普拉斯模型中左端的数值.在式 (3) ～ (6) 中都有核函数的作用, 核函数为粒子间相互作用关系的表达式, 在本文中选取的核函数为 w(r)=?????????????????23?4(rre)2+4(rre)3,43?4rre+4(rre)2?43(rre)3,0,0≤r≤re2re2r≤rerrew(r)={23-4(rre)2+4(rre)3,0≤r≤re243-4rre+4(rre)2-43(rre)3,re2r≤re0,rre 式中:r为2个粒子间的距离;re为核函数作用距离.MPS法中压力计算是通过压力泊松方程求解, 如下式: 式中:k为循环次数;〈n*〉i为粒子i的不同时间步的数密度. 所有粒子的压力泊松方程可用共轭梯度法 (ICCG) 求解.MPS法的计算流程如图1所示. 2 加标计算的粒子数密度法 MPS法应用到实际计算领域中的缺陷:虽然压力整体的变化趋势是可以接受的, 但振荡性过于明显, 超出了可用的范围.Sueyoshi等也在研究MPS的过程中发现了此问题并提出了一些改进方法, 这些改进方法在一定程度上解决了压力振荡的问题, 但增加了计算量, 因此, 应用并不是十分广泛. MPS法中利用粒子的数密度守恒来保证流体的不可压缩性, 并主要利用压力泊松方程式 (7) 求解压力.但是每次计算的粒子数密度的具体数值并不是理论值1, 因此需要进行修正;而在压力计算过程中的同一时间步中并没用修正的方法, 因此数密度的差异造成了压力的振荡性很大.而在计算流体力学里, 流体的不可压缩性可以通过速度散度为零的条件来保证, 如下式所示. 将此条件代入压力泊松方程, 得到压力与速度的关系: 式中,u*为计算过程中的粒子中间速度. 将式 (9) 和原数密度条件同时代入泊松方程, 可得改进方法的压力泊松方程: 式中,γ为松弛因子, 可根据不同的计算条件来修改, 以达到很好的计算效果.式 (10) 利用速度散度为零的条件和数密度条件共同求解压力, 缓解了由数密度不同造成的压力振荡. 3 比较分析的算术示例 基于本文提出的缓解压力振荡的改进MPS法, 分别对溃坝模型和液舱横摇模型进行改进前、改进后的对比. 3.1 坍塌 3.1.1 计算步长及效率 水箱内长为0.72 m, 内高为0.35 m, 水箱共布置内中外三层粒子, 粒子间距为2.5 mm