作　　者： (杨秀峰）; (刘谋斌）;

机构地区： Department of Mechanical Engineering Iowa State University Ames IA 50011 USA

出　　处： 《物理学报》 2017年第16期155-163,共9页

摘　　要： 提出了一种适用于模拟多相流的光滑粒子法,该方法对密度方程在交界面处的离散格式进行了修正以适应多相流所涉及的大密度比问题,在不同相粒子之间施加了很小的排斥力以防止粒子穿透交界面,并采用了最新发展的双曲型光滑函数以消除应力不稳定问题.应用该多相流光滑粒子法模拟研究了单模态和多模态瑞利-泰勒不稳定问题.通过与文献中结果的对比研究表明:在模拟瑞利-泰勒不稳定问题时,本文方法的结果明显优于文献中的大部分光滑粒子法模拟结果,与Grenier等(2009 J.Comput.Phys.228 8380)的结果相当,但本文方法比Grenier等的方法简单方便.对于单模态瑞利-泰勒不稳定问题,研究了交界面的形态,涡结构的演化过程以及贯穿深度随时间的变化关系.对于多模态瑞利-泰勒不稳定问题,研究了交界面演化过程中小尺度结构合并成大尺度结构的过程,水平方向的平均密度随高度的变化关系,以及贯穿深度随时间的变化关系. In this paper, we present a smoothed particle hydrodynamics （SPH） method for modeling multiphase flows. The multiphase SPH method includes a corrective discretization scheme for density approximation around the fluid interface to treat large density ratio, a small repulsive force between particles from different phases to prevent particles from unphysically penetrating fluid interface, and a newly-developed hyperbolic-shaped kernel function to remove possible stress instability. This multiphase SPH method is then used to study the single- and multi-mode Rayleigh-Taylor instability problems. A comparison between our results with the results from existing literature shows that our results are obviously better than most available results from other SPH simulations. The present results are close to those by Grenier et al. while the present multiphase SPH method is simpler and easier to implement than that in the work by Grenier et al. （Grenier, et al. 2009 J. Comput. Phys. 2288380）. For the single-mode Rayleigh-Taylor instability, the evolutions of the interface pattern and vortex structures, and the penetration depth each as a function of time are investigated. For the multi-mode Rayleigh-Taylor instability, the merging of small structures into a large structure during the evolution of the interface is studied. The horizontal average density and the penetration each as a function of height are also studied.

关 键 词： 瑞利 泰勒不稳定性 光滑粒子法 多相流