不可压流体自由表面流动的SPH数值模拟

根据SPH方法的原理提出了一套模拟流体自由表面流动的方法,场变量及其导数通过核函数插值求取,不需进行差分,也不需要网格,时间积分采用分数步长法,避免了不可压条件带来的压力计算不稳定问题。对流项通过粒子的运动求解,完全消除了数值扩散和自由表面模糊问题。可以模拟飞溅、融合等复杂自由表面现象,并对水坝坍塌这一典型自由表面流动问题进行了模拟,模拟结果与试验吻合良好。

粘性不可压流场数值模拟的无网格方法

采用无网格伽辽金法(elementfreeGalerkinmethod,EFG)对二维不可压粘性流场进行了数值模拟,主要计算了几个经典的流体力学问题:管道流、绕台阶流以及方腔驱动流。计算中采用标准的Galerkin离散,使用规则分割产生的背景网格来形成积分单元,并与有限元、有限差分法计算所得的速度剖面图进行了比较,结果相差很小。更重要的是,与有限元法相比,无网格虽然计算量比较大,但是它能够处理结点分布畸形甚至重合的情况,从而为解决多相流动直接数值模拟中网格重构的问题奠定基础。

用无网格法求解不同Re下圆柱绕流问题

采用无网格伽辽金法对不同Re下二维不可压粘性流场圆柱绕流问题进行了数值模拟。该计算不仅预测了不同条件下的流动特性,而且预测了流体作用在圆柱上的曳力和升力。这些力通过积分圆柱表面的压力和粘性力获得。结果显示,当Re大于80时,圆柱所受曳力和升力开始振荡。计算所得的脱落频率(Sr)与前人的结果吻合良好。