粒子法中压力振荡的机理研究

朱跃; 姜胜耀; 杨星团; 段日强

doi:10.6052/0459-1879-17-294

粒子法中压力振荡的机理研究

doi: 10.6052/0459-1879-17-294

清华大学核能与新能源技术研究院, 先进核能技术协同创新中心, 教育部先进反应堆工程与安全重点实验室, 北京 100084

基金项目: 国家自然科学基金资助项目(11472155).

详细信息

作者简介:
通讯作者: 朱跃, 博士研究生, 主要研究方向: 反应堆热工水力, 粒子法.E-mail: flyzhuyue@163.com; 段日强, 副研究员, 主要研究方向: 反应堆热工水力, 粒子法. E-mail: duanrq@tsinghua.edu.cn

通讯作者:
朱跃,

朱跃,段日强

计量
- 文章访问数: 1032
- HTML全文浏览量: 151
- PDF下载量: 192
- 被引次数: 0
出版历程
- 收稿日期: 2017-08-29
- 刊出日期: 2018-05-18

MECHANISM ANALYSIS OF PRESSURE OSCILLATION IN PARTICLE METHOD

Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education, Collaborative Innovation Center of Advanced Nuclear Energy Technologyrm, Institute of Nuclear and New Energy Technology, Tsinghua University, Beijing 100084 China

摘要

摘要: 移动粒子半隐式法(moving particle semi-implicit method, MPS)是一种适用于不可压缩流体的无网格方法, MPS方法常应用于自由表面大变形问题.MPS 方法提出至今一直存在着严重的压力振荡问题. 本研究针对MPS 方法中存在的压力振荡现象, 首先将实际的物理问题简化为一维模型, 并从粒子之间相互位置关系的角度说明了MPS 方法中压力波动产生的原因.在采用MPS方法进行模拟时, 加入了粒子碰撞模型, 通过对碰撞系数的选择从而控制粒子之间的相互位置关系.并且对经典的溃坝问题进行了模拟, 结果表明随着碰撞系数的增加, 粒子数密度偏差的波动幅度都会减小, 从而压力振荡的幅度得到了有效的抑制.并且对比了两种不同核函数对压力振荡的影响, 结果表明: 采用高斯核函数时, 压力振荡的幅度更小, 这是因为采用高斯核函数时, 相同的粒子位置波动幅度将会得到较小的粒子数密度偏差的波动.由于在模拟过程中粒子运动的随机性, 这将导致粒子数密度偏差产生随机的波动, 从而产生压力振荡, 因此粒子法中的压力振荡很难彻底消除.
- 移动粒子半隐式法 /
- 压力振荡机理 /
- 碰撞模型 /
- 核函数
Abstract: The moving particle semi-implicit method was developed to simulate incompressible fluid using a meshless method. There was a big problem that the space distribution and the time variation of pressure oscillate drastically in the MPS method. In order to investigate the pressure oscillation in the MPS method, the simplified one-dimensional model was developed. The mechanism of pressure oscillation in the MPS method was illustrated by the movement and the relative position between the center particle and its neighbor particles. The collision model was employed in the simulation and the relative position between particles was controlled by choosing different collision parameters. The classical dam break problem was simulated. With the increase of collision parameter, the fluctuation of the deviation of particle number density decreased. According, the amplitude of pressure oscillation was suppressed. Two different kernel functions were also employed to investigate the pressure oscillations. The results showed that the gauss kernel function improved the stabilization of pressure calculation. It was the reason that the same movement of particles lead to less deviation of particle number density when the gauss kernel function was used. And the randomness of particles motion led to random fluctuation in the deviation of particle number density. As a result, the pressure fluctuation in MPS method occurred and it was difficult to be eliminated.
- moving particle semi-implicit method /
- mechanism of pressure oscillation /
- collision model /
- kernel function

HTML全文

参考文献(1)

[1]	Koshizuka S. A particle method for incompressible viscous flow with fluid fragmentation. Comput Fluid Dynamics J, 1995, 4: 29-46
[2]	Koshizuka S. Moving-particle semi-implicit method for fragmentation of incompressible fluid. Nuclear Science ＆ Engineering the Journal of the American Nuclear Society, 1996, 123(3): 421-434
[3]	Koshizuka S,Ikeda H,Oka Y. Numerical analysis of fragmentation mechanisms in vapor explosions. Nuclear Engineering and Design, 1998, 189(1-3): 423-433
[4]	Koshizuka S, Oka Y. Application of moving particle semi-implicit method to nuclear reactor safety. Computational Fluid Dynamics, 2001, 9: 366-375
[5]	Xie H, Koshizuka S, Oka Y. Simulation of drop deposition process in annular mist flow using three-dimensional particle method. Nuclear Engineering and Design, 2005, 16(235): 1687-1697
[6]	Tian W, Ishiwatari Y, Ikejiri S, et al.Numerical computation of thermally controlled steam bubble condensation using moving particle semi-implicit (MPS) method. Annals of Nuclear Energy, 2010, 1(37): 5-15
[7]	Yoon HY, Koshizuka S, Oka Y. A mesh-free numerical method for direct simulation of gas-liquid phase interface. Nuclear Science and Engineering, 1999, 133: 1-9
[8]	Yoon HY, Koshizuka S, Oka Y, Direct calculation of bubble growth, departure and rise in nucleate pool boiling. International Journal of Multiphase Flow, 2011, 27(2): 277-298
[9]	Chen R, Chen L, Guo K, et al.Numerical analysis of the melt behavior in a fuel support piece of the BWR by MPS. Annals of Nuclear Energy, 2017, 102: 422-439
[10]	Li Y, Chen R, Guo K, et al.Numerical analysis of the dissolution of uranium dioxide by molten zircaloy using MPS method. Progress in Nuclear Energy, 2017, 100: 1-10
[11]	Jena D, Biswal KC. A numerical study of violent sloshing problems with modified MPS method. Journal of Hydrodynamics, Ser.B, 2017, 29(4): 659-667
[12]	Shibata K, Koshizuka S. Numerical analysis of shipping water impact on a deck using a particle method. Ocean Engineering, 2007, 34(3-4): 585-593
[13]	Sueyoshi M, Kashiwagi M, Naito S. Numerical simulation of wave-induced nonlinear motions of a two-dimensional floating body by the moving particle semi-implicit method. Journal of Marine Science and Technology, 2008, 13(2): 85-94
[14]	Shibata K, Koshizuka S, Tanizawa K. Three-dimensional numerical analysis of shipping water onto a moving ship using a particle method. Journal of Marine Science and Technology, 2009, 14(2): 214-227
[15]	Guo K, Chen R, Li Y, et al.Numerical investigation of the fluid-solid mixture flow using the focus code. Progress in Nuclear Energy, 2017, 97: 197-213
[16]	Khayyer A, Gotoh H. Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure. Coastal Engineering, 2009, 56(4): 419-440
[17]	Lee B, Park J, Kim M, et al.Step-by-step improvement of MPS method in simulating violent free-surface motions and impact-loads. Computer Methods in Applied Mechanics and Engineering, 2011, 200(9-12): 1113-1125
[18]	Sueyoshi M Naito S. A study of nonlinear fluid phenomena with particle method (Part 2). Journal of Kansai Society Naval Architects, 2002, 2002(237): 181-186
[19]	Khayyer A, Gotoh H. A higher order Laplacian model for enhancement and stabilization of pressure calculation by the MPS method. Applied Ocean Research, 2010, 32(1): 124-131
[20]	Kondo M, Koshizuka S. Improvement of stability in moving particle semi-implicit method. International Journal for Numerjcal Methpds in Fluids, 2011, 65(6): 638-654
[21]	Tanaka M, Masunaga T. Stabilization and smoothing of pressure in MPS method by quasi-compressibility. Journal of Computational Physics, 2010, 229(11): 4279-4290
[22]	Shibata K, Masaie I, Kondo M, et al.Improved pressure calculation for the moving particle semi-implicit method. Computational Particle Mechanics, 2015, 2(1): 91-108
[23]	潘徐杰, 张怀新.移动粒子半隐式法晃荡模拟中的压力震荡现象研究.水动力学研究与进展, 2008, 23(4): 453-463
[23]	(Pan Xujie, Zhang Huaixin. A study on the oscillations appearing in pressure calculation for sloshing simulation by using moving-particle semi-implicit method. Chinese Journal of Hydrodynamics, 2008, 23(4): 453-463 (in Chinese))
[24]	Duan G, Chen B, Koshizuka S, et al.Stable multiphase moving particle semi-implicit method for incompressible interfacial flow. Computer Methods in Applied Mechanics and Engineering, 2017, 318: 636-666
[25]	郑仙佩, 顾声龙, 任立群.基于SWE-SPH的下游为湿河床的溃坝水流模拟.青海大学学报, 2016, 34(5): 40-44
[25]	(Zheng Xianpei, Gu Shenglong, Ren Liqun. Simulation of dam-break flow with wet bed at downstream based on SWE-SPH. Journal of Qinghai University, 2016, 34(5): 40-44 (in Chinese))
[26]	李婧文, 陈昌平, 孙家文等.基于溃坝模型的SPH方法光滑函数模拟.中国海洋平台, 2017, 32(2): 34-40
[26]	(Li Jingwen, Chen Changping, Sun Jiawen, et al.Smoothing fution simulation of SPH method based on dam-break model. China Offshore Platform, 2017, 32(2): 34-40 (in Chinese))
[27]	刘汉涛, 干湧, 常建忠等.溃坝问题的粒子方法模拟.应用力学学报, 2016, 33(4): 547-552
[27]	(Liu Hantao, Gan Yong, Chang Jianzhong, et al.Dam break simulation employing particle methods. Chinese Journal of Applied Mechanics, 2016, 33(4): 547-552 (in Chinese))
[28]	陈健云, 李静, 徐强等.溃坝涌浪及其对重力坝影响的数值模拟.浙江大学学报(工学版), 2016, 50(11): 2164-2169
[28]	(Chen Jianyun, Li Jing, Xu Qiang, et al.Numerical simulation of Landslide surge and its effects on gracity dam. Journal of Zhejiang University(Engineering Science), 2016, 50(11): 2164-2169 (in Chinese))
[29]	许晓阳, 任胜章, 邓方安.三维溃坝流的光滑粒子动力学方法模拟.计算力学学报, 2016, 33(5): 676-682
[29]	(Xu Xiaoyang, Ren Shengzhang, Deng Fang’an. Numerical simulation of 3D dam-break flow using smoothed particle hydrodynamics. Chinese Journal of Computational Mechanics, 2016, 33(5): 676-682 (in Chinese))
[30]	Daneshvar FA, Rakhshandehroo GR, Talebbeydokhti N. New modified gradient models for MPS method applied to free-surface flow simulations. Applied Ocean Research, 2017, 66: 95-116
[31]	Hu C, Kashiwagi M. A CIP-based method for numerical simulations of violent free-surface flows. Journal of Marine Science and Technology, 2004, 9: 143-157