修正表面张力算法的SPH方法及其实现  

Smoothed Particle Hydrodynamics Method with Modified Surface Tension and Its Implementation

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

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

出　　处：《计算物理》2011年第3期375-384,共10页Chinese Journal of Computational Physics

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

摘　　要：在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到表面张力修正方程组;通过半圆形算例测试方法和Morris方法在边界定位、法向计算和曲率计算等影响表面张力关键因素的求解精度,研究曲率计算中应采用的光滑长度值.模拟初始方形液滴在表面张力作用下的自然变化过程,并与Morris方法及VOF有限体积法进行对比,表明方法精度较高,稳定性较好.最后,模拟水溶液中两个油滴的互溶过程.Modified equations for surface tension are derived by modifying normal and curvature with corrected smoothed particle method（CSPM）.It is based on smoothed particle hydrodynamics（SPH） method with surface tension proposed by Morris.Both Morris and our method are tested via a semicircular problem.Factors that affect accuracy are investigated including surface definition,normal and curvature calculation.Smoothed length in curvature calculation is also confirmed reasonable.Furthermore,formation of a liquid drop with initial square shape under surface tension is simulated.Compared with Morris method and grid-based volume of fluid method,it is proved that the accuracy of our method is higher and particle distribution is more homogeneous.Finally,coalescence process of two oil drops in water under surface tension is simulated.

关 键 词：光滑粒子流体动力学 表面张力算法 CSPM方法 CSF方法 曲率 

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