人工黏性增强三维移动粒子半隐算法稳定性

冯美艳; 黄生洪

doi:10.6052/0459-1879-12-247

人工黏性增强三维移动粒子半隐算法稳定性

冯美艳,
黄生洪^,

中国科学技术大学工程科学学院近代力学系, 合肥 230026

基金项目: 国家自然科学基金(10772174)和国家专项(2013GB113000)资助项目.

详细信息

通讯作者:
黄生洪,副教授,主要研究方向:多相与多场耦合流动.E-mail:hshnpu@ustc.edu.cn

中图分类号: O351
计量
- 文章访问数: 1919
- HTML全文浏览量: 111
- PDF下载量: 600
出版历程
- 收稿日期: 2012-09-10
- 修回日期: 2013-03-14
- 刊出日期: 2013-07-17

ENHANCEMENT OF STABILITY OF 3D MOVING PARTICLE SEMI-IMPLICIT METHOD BY ARTIFICIAL VISCOSITY

Feng Meiyan,
Huang Shenghong^,

Department of Modern Mechanics, University of Science and Technology of China, Hefei 230026, China

Funds: The project was supported by the National Natural Science Foundation of China (10772174) and the National Special Found (2013GB113000).

摘要

摘要: 通过综合比较移动粒子半隐方法(moving-particle semi-implicitmethod, MPS)各种稳定性计算方法, 考虑了多种可能导致计算不稳定的因素, 并首次提出加入人工黏性来抑制非物理压力振荡, 得到了较为稳定的三维MPS算法.采用各种稳定性方案对比模拟了三维立方流体旋转状态下的变形及三维静水压问题, 并进一步与商业软件的流体体积函数方法(volumeof fluid, VOF)模型计算结果对比, 验证了其正确性.结果表明:对三维情形, 单纯应用已有的稳定性算法仍难以满足MPS计算稳定性的要求, 而进一步辅以该文提出的人工黏性方法则可以在准确性及稳定性方面均取得较好的效果, 且计算结果显示, 通过该文提出的开关模式增加人工黏性并不影响压强泊松方程的求解精度, 因而不会对流场的求解产生负面作用.
- 移动粒子半隐方法 /
- 三维 /
- 稳定性 /
- 人工黏性 /
- 无网格方法
Abstract: By the virtue of Lagrangian gridless particle method, the moving particle semi-implicit (MPS) method has been applied in a wide range of engineering applications. However, the built-in instability in computational algorithm has seriously limited its further development and application. In this paper, aiming at developing 3D stabilized MPS method, a comprehensive comparison of stabilizing methods proposed by different researchers has been conducted, in which many elements leading to instability are considered. Then the equations proposed in 2D are extended to 3D and an artificial viscosity model based on smooth particle hydrodynamics (SPH) is proposed to reduce instability in MPS for the first time. With the different combinations of proposed methods and models, the deforming process of a swirling cubic fluid patch and hydrostatic pressure problem are simulated comparatively. The results are also compared with those from volume of fluid (VOF) model conducted by a commercial software. The accuracy and stability of proposed methods and models are validated. It is concluded that: for 3D case, only extending and applying previous stabilizing method is still dissatisfactory, while good effects in accuracy and stability can be obtained by incorporating the artificial viscosity model proposed in present paper. And it is also revealed that no negative effects on exactness of pressure solution will be produced, and thus on the flow field.
- moving particle semi-implicit method /
- three-dimension /
- stability /
- artificial viscosity /
- gridless particle method

HTML全文

参考文献(21)

Koshizuka S, Tamako H, Oka Y. A particle method for incompressible viscous flow with fluid fragmentation. Journal of Computational Fluid Dynamics, 1995, 4: 29-46

Sueyoshi M, Naito S. A study of nonlinear fluid phenomena with particle method (part 2). Journal of Kansai Society Naval Architects, 2002, 237: 181-186 (in Japanese))

Ataie-Ashtiani B, Farhadi L. A stable moving-particle semi-implicit method for free surface flows. Fluid Dynamics Research, 2006,38(4): 241-256

Shao S, Lo EYM. Incompressible SPH method for simulating newtonian and non-newtonian flows with a free surface. Adv Water Resour, 2003, 26: 787-800

潘徐杰,张怀新,孙学尧.移动粒子半隐式法中核函数适用性研究.海洋工程, 2010, 28(4): 70-76 (Pan Xujie,Zhang Huaixin. Study on kernel fimctions of moving particle semi-implicit method. The Ocean Engeering, 2010, 28(4):70-76 (in Chinese))

Tanaka M,Masunaga T. Stabilization and smoothing of pressure in MPS method by quasi-compressibility. Journal of Computational Physics, 2010. 4279-4290

Kondo M, Koshizuka S. Improvement of stability in moving particle semi-implicit method. Int J Numer Methods Fluids, 2011, 65: 638-654

Khayyer A, Gotoh H. Modified moving particle semi-implicit methods for the prediction of 2D wave impact pressure. Coastal Engineering, 2009, 56: 419-440

Khayyer A, Gotoh H. A higher order Laplacian model for enhancement and stabilization of pressure. Applied Ocean Research, 2010, 32: 124-131

Koshizuka S, Oka Y. Moving-particle semi-implicit method for fragmentation of incompressible fluid. Nuclear Science and Engineering, 1996, 123: 421-434

Koshizuka S, Nobe A, Oka Y. Numerical analyses of breaking waves using the moving particle semi-implicit method. Int J Numer Meth Fluids, 1998, 26: 751-769 3.0.CO;2-C">

Khayyer A, Gotoh H. Enhancement of stability and accuracy of the moving particle semi-implicit method. Journal of Computational Physics, 2011, 230: 3093-3118

Price DJ. Smoothed particle hydrodynamics and magnetohydrodynamics. Journal of Computational Physics, 2012, 231: 759-794

黄生洪,冯美艳. 超声速气流中雾化全过程模拟的一体化方法. 见:第十四届全国激波与激波管学术会议文集,黄山,2010 (Huang Shenghong,Feng Meiyan. An integrated method for process simulation of liquid atomization in supersonic flow and its preliminary results. In: Proc. of the 14th Chinese National Symposium on Shock Waves, Huangshan, 2010 (in Chinese))

黄生洪,冯美艳.雾化燃烧模拟一体化方法及应用研究.工程热物理学报,2011, 32(8):1399-1402 (Huang Shenghong, Feng Meiyan. An integrated simulation method of spray combustion and its application. Journal of Engineering Thermophysics, 2011, 32(8):1399-1402 (in Chinese))

Monaghan JJ. Smoothed particle hydrodynamics. Ann Rev Astron Astrophys, 1992, 30: 543-574

Monaghan JJ. SPH and riemann solvers. J Comput Phys, 1997, 136: 298-307

Chow E, Monaghan JJ. Ultrarelativistic SPH. J Comput Phys, 1997, 134: 296-305

Price DJ. Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH. J Comput Phys, 2008, 227: 10040-10057

Fang J, Parriaux A, Rentschler M, et al. Improved SPH methods for simulating free surface flows of viscous fluids. Appl Numer Math, 2009, 59 (2): 251-271

Oger G, Doring M, Alessandrini B, et al. An improved SPH method: Towards higher order convergence. J Comput Phys, 2007, 225 (2): 1472-1492