摘要
在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.
出处
《计算物理》
EI
CSCD
北大核心
2011年第3期375-384,共10页
Chinese Journal of Computational Physics
基金
国家教育部NCET
国家973(973-61310)
第二炮兵工程学院创新性探索研究(EPXY0806)资助项目