基于SPH方法的表面张力修正算法及其应用

针对传统SPH方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组。应用本文方法模拟了水溶液中初始方形油液滴在表面张力作用下的自然变化过程,与VOF有限体积法进行对比,结果非常吻合;最后,基于本文方法模拟了二维水溶液中两个油滴的碰撞破碎过程。