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

Modified algorithm for surface tension in smoothed particle hydrodynamics and its application

作　　者：强洪夫[1] 陈福振[1] 高巍然[1] 

机构地区：[1]第二炮兵工程学院201室,西安710025

出　　处：《计算力学学报》2011年第B04期37-42,共6页Chinese Journal of Computational Mechanics

基　　金：国家教育部NCET资助项目;国家973资助项目(973-613102);第二炮兵工程学院创新性探索研究资助项目(EPXY0806)

摘　　要：针对传统SPH方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组。应用本文方法模拟了水溶液中初始方形油液滴在表面张力作用下的自然变化过程,与VOF有限体积法进行对比,结果非常吻合;最后,基于本文方法模拟了二维水溶液中两个油滴的碰撞破碎过程。The continuum surface force （CSF）model for surface tension, which is employed in the traditional smoothed particle hydrodynamics （SPH） method, is very low accurate when handling the problem with large deformation and surface tension. The reason is that the deficiency of particles near the boundary and sharp-angled areas, and it caused large error of curvature calculation in the traditional SPH method. In this paper modified equations for surface tension were derived by modifying normal and curvature with corrective smoothed particle method （CSPM）. It was based on the SPH method with surface tension proposed by Morris. Formation of an oil drop under surface tension from initial square shape in water was simulated with the new method. Simultaneous, two-dimensional simulation of breaking up process of two colliding oil drops in water under surface tension was also simulated.

关 键 词：光滑粒子流体动力学 表面张力算法 CSF方法 曲率 液滴碰撞破碎 

分 类 号：O35[理学—流体力学]