粒子法中自由表面粒子识别
|
摘要
针对移动粒子半隐式法MPS(Moving Particle Semi-implicit Method)基于粒子数密度来判断自由表面会出现将内部粒子误判为自由表面粒子的问题,提出了一种结合几何法和体积法的自由表面粒子判定方法。通过对溃坝问题进行数值模拟,结果表明,全新的自由表面粒子判定方法对流体平稳运动以及剧烈运动两种工况,都能准确地判断出自由表面粒子,解决了基于粒子数密度判断方法因粒子分布稀疏产生误判的问题。这种全新的自由表面粒子判定方法对今后采用MPS方法计算两相流问题时,两种介质在界面处的传热传质计算有重要意义。
The conventional surface particle detection method based on the particle number density would lead to the misjudgement of surface particles.To address this issue,this paper proposes a new surface particle tracking method combining geometrical and volumetric principles.Through numerical simulation of dam break,the new surface particles tracking method achieves very good results in both smooth movement fluid and strenuous movement fluid.The simulation results show that this new method is effective in control of surface particle misjudgement,and the issue of misjudgement of surface particles based on the particle number density is solved.This new surface particle detection method has important implications for interfacial heat and mass transfer calculation in two-phase flow problems.
The conventional surface particle detection method based on the particle number density would lead to the misjudgement of surface particles.To address this issue,this paper proposes a new surface particle tracking method combining geometrical and volumetric principles.Through numerical simulation of dam break,the new surface particles tracking method achieves very good results in both smooth movement fluid and strenuous movement fluid.The simulation results show that this new method is effective in control of surface particle misjudgement,and the issue of misjudgement of surface particles based on the particle number density is solved.This new surface particle detection method has important implications for interfacial heat and mass transfer calculation in two-phase flow problems.
引文
[1]Koshizuka S.A particle method for incompressible viscous flow with fluid fragmentation[J].Computational Fluid Dynamics Journal,1995,4:29-46.
[2]Koshizuka S,Oka Y.Moving-particle semi-implicit method for fragmentation of incompressible fluid[J].Nuclear Science and Engineering,1996,123(3):421-434.
[3]Koshizuka S,Ikeda H,Oka Y.Numerical analysis of fragmentation mechanisms in vapor explosions[J].Nuclear Engineering and Design,1999,189(1):423-433.
[4]Lee B H,Park J C,Kim M H,et al.Step-by-step improvement of MPS method in simulating violent freesurface motions and impact-loads[J].Computer Methods in Applied Mechanics and Engineering,2011,200(9-12):1113-1125.
[5]Randles P W,Libersky L D.Smoothed particle hydrodynamics:some recent improvements and applications[J].Computer Methods in Applied Mechanics and Engineering,1996,139(1-4):375-408.
[6]Lee E S,Moulinec C,Xu R,et al.Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method[J].Journal of Computational Physics,2008,227(18):8417-8436.
[7]Tanaka M,Masunage T.Stabilization and smoothing of pressure in MPS method by quasi-compressibility[J].Journal of Computational Physics,2010,229(11):4279-4290.
[8]Dilts G A.Moving least-squares particle hydrodynamics II:Conservation and boundaries[J].International Journal for Numerical Methods in Engineering,2000,48(10):1503-1524.
[9]Haque A,Dilts G A.Three-dimensional boundary detection for particle methods[J].Journal of Computational Physics,2007,226(2):1710-1730.
[10]Shibata K,Masaie I,Kondo M,et al.Improved pressure calculation for the moving particle semi-implicit method[J].Computational Particle Mechanics,2015,2(1):91-108.
[11]Marrone S,Colagrossi A,Le TouzéD,et al.Fast free-surface detection and level-set function definition in SPH solvers[J].Journal of Computational Physics,2010,229(10):3652-3663.
[12]李绍武,熊赞,武霄.MPS水流数学模型在冲击压力计算中的应用[J].天津大学学报,2006,39(9):1031-1036.(LI Shao-wu,XIONG Zan,WU Xiao.Application of a numerical model based on MPS method in the calculation of impact pressure[J].Jouranal of Tianjin University,2006,39(9):1031-1036.(in Chinese))
[13]李绍武,尹振军.无网格数值方法在结构物水动力研究中的应用[J].水科学进展,2004,15(6):739-744.(LI Shao-wu,YIN Zhen-jun.Application of non-grid numerical model to hydrodynamics of structures[J].Advances in Water Science,2004,15(6):739-744.(in Chinese))
[14]席光,项利峰.自由表面流动的移动粒子半隐式模拟方法[J].西安交通大学学报,2006,40(3):249-252,288.(XI Guang,XIANG Li-feng.Numerical simulation of fluid free surface flow using moving particle semi-implicit method[J].Journal of Xi’an Jiaotong University,2006,40(3):249-252,288.(in Chinese))
[15]孙中国,席光,项利峰.不可压缩流动的移动粒子半隐式法数值模拟[J].应用力学学报,2007,24(3):440-442.(SUN Zhong-guo,XI Guang,XIANG Li-feng.Simulation for incompressible flow with moving particle semi-implicit method.[J].Chinses Journal of Applied Mechanics,2007,24(3):440-442.(in Chinese))
[2]Koshizuka S,Oka Y.Moving-particle semi-implicit method for fragmentation of incompressible fluid[J].Nuclear Science and Engineering,1996,123(3):421-434.
[3]Koshizuka S,Ikeda H,Oka Y.Numerical analysis of fragmentation mechanisms in vapor explosions[J].Nuclear Engineering and Design,1999,189(1):423-433.
[4]Lee B H,Park J C,Kim M H,et al.Step-by-step improvement of MPS method in simulating violent freesurface motions and impact-loads[J].Computer Methods in Applied Mechanics and Engineering,2011,200(9-12):1113-1125.
[5]Randles P W,Libersky L D.Smoothed particle hydrodynamics:some recent improvements and applications[J].Computer Methods in Applied Mechanics and Engineering,1996,139(1-4):375-408.
[6]Lee E S,Moulinec C,Xu R,et al.Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method[J].Journal of Computational Physics,2008,227(18):8417-8436.
[7]Tanaka M,Masunage T.Stabilization and smoothing of pressure in MPS method by quasi-compressibility[J].Journal of Computational Physics,2010,229(11):4279-4290.
[8]Dilts G A.Moving least-squares particle hydrodynamics II:Conservation and boundaries[J].International Journal for Numerical Methods in Engineering,2000,48(10):1503-1524.
[9]Haque A,Dilts G A.Three-dimensional boundary detection for particle methods[J].Journal of Computational Physics,2007,226(2):1710-1730.
[10]Shibata K,Masaie I,Kondo M,et al.Improved pressure calculation for the moving particle semi-implicit method[J].Computational Particle Mechanics,2015,2(1):91-108.
[11]Marrone S,Colagrossi A,Le TouzéD,et al.Fast free-surface detection and level-set function definition in SPH solvers[J].Journal of Computational Physics,2010,229(10):3652-3663.
[12]李绍武,熊赞,武霄.MPS水流数学模型在冲击压力计算中的应用[J].天津大学学报,2006,39(9):1031-1036.(LI Shao-wu,XIONG Zan,WU Xiao.Application of a numerical model based on MPS method in the calculation of impact pressure[J].Jouranal of Tianjin University,2006,39(9):1031-1036.(in Chinese))
[13]李绍武,尹振军.无网格数值方法在结构物水动力研究中的应用[J].水科学进展,2004,15(6):739-744.(LI Shao-wu,YIN Zhen-jun.Application of non-grid numerical model to hydrodynamics of structures[J].Advances in Water Science,2004,15(6):739-744.(in Chinese))
[14]席光,项利峰.自由表面流动的移动粒子半隐式模拟方法[J].西安交通大学学报,2006,40(3):249-252,288.(XI Guang,XIANG Li-feng.Numerical simulation of fluid free surface flow using moving particle semi-implicit method[J].Journal of Xi’an Jiaotong University,2006,40(3):249-252,288.(in Chinese))
[15]孙中国,席光,项利峰.不可压缩流动的移动粒子半隐式法数值模拟[J].应用力学学报,2007,24(3):440-442.(SUN Zhong-guo,XI Guang,XIANG Li-feng.Simulation for incompressible flow with moving particle semi-implicit method.[J].Chinses Journal of Applied Mechanics,2007,24(3):440-442.(in Chinese))