改进的CSF流体表面张力模型

传统连续表面力(continuum surface force,CSF)模型模拟流体表面张力时,流体表面的粒子受到的表面张力都是法向力,很难保证流体表面的光滑性,从而使得表面张力模拟失真;此外,流体表面的粒子不足还导致流体表面密度计算精度降低,模拟稳定性差。文章在CSF模型的基础上,对流体表面施加切向力,使得流体表面更加光滑;对边界粒子进行密度修正,提高了密度计算精度和模拟的稳定性。仿真实验结果表明,该方法模拟的流体表面张力效果更好。

一种模拟激光刻花的二维CSF模型

采用连续表面力的方法 ,建立了模拟激光刻花的二维瞬态模型 .模型考虑了热传导、对流传热和熔池表面形貌变化等因素 ,设计了适合连续表面力方法的热源和虚拟体积力 ,建立了物理模型 ,并利用有限元和多层网格方法求解 ,得出了激光刻花过程中材料内部的等温线。

适用于粒子法的精准连续表面力模型

针对粒子法中表面张力计算的准确性问题,该文对连续表面力（continuum surface force,CSF）模型进行了改进。在采用一种几何法精准识别界面粒子的基础上,曲率由仅需考虑作用域内界面粒子的单位法线面散度计算得到,且表面张力仅作用于界面粒子。圆和椭圆的曲率计算结果表明,改进后的模型在合适的分辨率和光滑长度下可实现较高的曲率计算精度。采用移动粒子半隐式（moving particle semi-implicit,MPS）方法对表面张力作用下的方形液滴振荡和液滴碰撞过程进行了二维单相流动模拟,模拟结果的振荡周期与理论值接近,液滴形状合理且表面光滑,表明改进后的CSF模型可准确地模拟动态算例中的表面张力作用。

二维溃坝的数值模拟研究

为了研究溃坝问题中流体的运动和自由液面的演变,针对简化的块状水坍塌进行了数值模拟。其中,流体控制方程通过有限体积法进行离散,方程的对流项和扩散项分别采用了一阶迎风格式和中心差分格式;对压力与速度耦合的离散方程以IDEAL算法进行求解;以流体体积分数法和PLIC方法进行自由液面捕捉和重构;对两相流问题中的表面张力通过连续表面力(CSF)模型进行转化。基于上,分别对有限空间内的二维溃坝和带障碍物溃坝问题进行了数值模拟,对不同时刻的自由液面进行了捕捉并与实验结果进行了对比,证明了所采用方法的精确性。