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

计算物理 ›› 2011, Vol. 28 ›› Issue (3): 375-384.

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

强洪夫, 陈福振, 高巍然

第二炮兵工程学院201 室, 陕西西安 710025

收稿日期:2010-05-10 修回日期:2010-09-28 出版日期:2011-05-25 发布日期:2011-05-25
作者简介:强洪夫(1965-),男,江苏武进,教授,博士,博导,国家973项目技术首席专家,从事材料与结构的失效行为和高性能数值模拟研究,E-mail:Qiang@263.net
基金资助:
国家教育部NCET;国家973(973-61310);第二炮兵工程学院创新性探索研究(EPXY0806)资助项目

A modified PUFF equation of state is proposed with considenation of both non-linear volume change and anisotropic strength of material under compression and dilation.Taking 2-D impact of orthotropic material and thermal shock wave under X-ray radiation as examples,a finite element program is developed.The mean stresses obtained by traditional and modified PUFF equation of state are compared.They are obviously different at low compression and tension states.With increase of stress,they tends to be consistent gradually.

QIANG Hongfu, CHEN Fuzhen, GAO Weiran

201 Staff Room, Xi'an Hi-Tech Institute, Xi'an 710025, China

Received:2010-05-10 Revised:2010-09-28 Online:2011-05-25 Published:2011-05-25

摘要/Abstract

摘要： 在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到表面张力修正方程组;通过半圆形算例测试方法和Morris方法在边界定位、法向计算和曲率计算等影响表面张力关键因素的求解精度,研究曲率计算中应采用的光滑长度值.模拟初始方形液滴在表面张力作用下的自然变化过程,并与Morris方法及VOF有限体积法进行对比,表明方法精度较高,稳定性较好.最后,模拟水溶液中两个油滴的互溶过程.

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

Abstract: 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.

Key words: smoothed particle hydrodynamics, surface tension method, CSPM, CSF, curvature

中图分类号:

强洪夫, 陈福振, 高巍然. 修正表面张力算法的SPH方法及其实现[J]. 计算物理, 2011, 28(3): 375-384.

QIANG Hongfu, CHEN Fuzhen, GAO Weiran. A modified PUFF equation of state is proposed with considenation of both non-linear volume change and anisotropic strength of material under compression and dilation.Taking 2-D impact of orthotropic material and thermal shock wave under X-ray radiation as examples,a finite element program is developed.The mean stresses obtained by traditional and modified PUFF equation of state are compared.They are obviously different at low compression and tension states.With increase of stress,they tends to be consistent gradually.[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(3): 375-384.

[1]	宫翔飞, 杨基明, 张树道. 使用AMR型背景网格的并行SPH方法[J]. 计算物理, 2016, 33(2): 183-189.
[2]	韩亚伟, 强洪夫. 改进的物理粘性SPH方法及其在溃坝问题中的应用[J]. 计算物理, 2012, 29(5): 693-699.
[3]	王安文, 徐绯, 张岳青. SPH方法在液固撞击数值模拟中的应用[J]. 计算物理, 2012, 29(4): 525-533.
[4]	郑兴, 马庆位, 段文洋. K2_SPH方法及二维破碎波的模拟[J]. 计算物理, 2012, 29(3): 317-325.
[5]	程俊霞. 高能炸药爆速曲率关系的-种计算方法[J]. 计算物理, 2011, 28(6): 817-824.
[6]	郑兴, 段文洋. K2_SPH方法及其对二维非线性水波的模拟[J]. 计算物理, 2011, 28(5): 659-666.
[7]	解立强. 高亏格膜泡形状的数值计算[J]. 计算物理, 2009, 26(5): 763-766.
[8]	徐金中, 汤文辉. 高速碰撞SPH方法模拟中的初始光滑长度和粒子间距[J]. 计算物理, 2009, 26(4): 548-552.
[9]	刘学哲, 沈智军, 岳晶岩. 曲率与运动界面发展[J]. 计算物理, 2008, 25(6): 668-672.
[10]	强洪夫, 高巍然. 完全变光滑长度SPH法及其实现[J]. 计算物理, 2008, 25(5): 569-575.
[11]	经荥清, 单文态, 罗经宇. 谱计算中的曲率反馈[J]. 计算物理, 1992, 9(S1): 597-601.
[12]	密柯, 章本照. 含尘离心风机叶轮磨损研究的一种二维边界层求解方法[J]. 计算物理, 1992, 9(4): 473-476.

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

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 12

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍