基于粒子分解的SPH并行算法研究与应用

计算机工程与科学 ›› 2022, Vol. 44 ›› Issue (06): 964-970.

基于粒子分解的SPH并行算法研究与应用

许晓阳，王斯棋

（西安科技大学计算机科学与技术学院，陕西西安 710054）

收稿日期:2021-04-20 修回日期:2021-10-11 接受日期:2022-06-25 出版日期:2022-06-25 发布日期:2022-06-17
基金资助:
国家自然科学基金(12071367);陕西省青年科技新星项目(2019KJXX-012)

Research and application of a SPH parallel algorithm based on particle decomposition

XU Xiao-yang,WANG Si-qi

(School of Computer Science and Technology,Xi’an University of Science and Technology,Xi’an 710054,China)

Received:2021-04-20 Revised:2021-10-11 Accepted:2022-06-25 Online:2022-06-25 Published:2022-06-17

摘要/Abstract

摘要： 作为一种典型的拉格朗日型无网格数值方法，光滑粒子流体动力学（SPH）方法在模拟自由表面流问题时具有天然优势。但是，该方法计算量大、耗时长，为此提出了一种基于粒子分解的SPH并行算法。该算法将所有粒子平均分配到各个进程进行计算，每个时间步通信仅调用一次发送、接收和广播函数，因此易于实现且可扩展性较好。应用该并行算法对二维溃坝流和三维液滴冲击液膜问题进行数值模拟,结果表明：该并行算法能显著减少模拟所消耗的计算时间，有利于进行三维大规模计算问题的数值模拟；当粒子数大于百万时，最大加速比可达30以上。

关键词: 光滑粒子流体动力学, 自由表面流, 并行算法, 粒子分解, 数值模拟

Abstract: As a typical meshless numerical method, Smoothed Particle Hydrodynamics (SPH) method has natural advantages in modeling free surface flows. However, this method is computationally heavy and time-consuming. Therefore, this paper proposes a SPH parallel algorithm based on particle decomposition. The algorithm distributes all particles to each process evenly for calculation, and the sending, receiving as well as broadcasting functions are called only once at each time step of communication, which is easy to implement and has good scalability. This parallel algorithm is used to carry out the numerical simulations of 2D dam-break flow and 3D droplet impact onto a liquid film. The results show that this SPH parallel algorithm can significantly reduce the simulation time and is conducive to 3D large-scale computation. The maximum speedup can reach more than 30 when the number of particles is greater than one million.

Key words: smoothed particle hydrodynamics;free surface flow;parallel algorithm, particle decomposition;numerical simulation

许晓阳, 王斯棋. 基于粒子分解的SPH并行算法研究与应用[J]. 计算机工程与科学, 2022, 44(06): 964-970.

XU Xiao-yang, WANG Si-qi. Research and application of a SPH parallel algorithm based on particle decomposition[J]. Computer Engineering & Science, 2022, 44(06): 964-970.

参考文献［22］

［1］	Lucy L B.A numerical approach to the testing of the fissionhypothesis［J］.Astronomical Journal,1977,82:1013-1024.
［2］	Gingold R A,Monaghan J J.Smoothed particle hydrodynamics theory and application to non-spherical stars［J］.Monthly Notices of the Royal Astronomical Society,1977,181:375-389.
［3］	Monaghan J J.Simulating free surface flows with SPH［J］.Journal of Computational Physics,1994,110(2):399-406.
［4］	Shao S,Lo E Y M.Incompressible SPH method for simulat- ing Newtonian and non-Newtonian flows with a free surface［J］.Advances in Water Resources,2003,26:787-800.
［5］	Khayyer A, Shimizu Y,Gotoh H,et al.A coupled incompressible SPH-Hamiltonian SPH solver for hydroelastic FSI corresponding to composite structures［J］.Applied Mathematical Modelling,2021,94:242-271.
［6］	Zhang M Y, Deng X L,Shen Z J.A robust sharp interface method for SPH［J］.Engineering Analysis with Boundary Elements,2019,106:275-285.
［7］	Meng Z F,Wang P P,Zhang A M,et al.A multiphase SPH model based on Roe’s approximate Riemann solver for hydraulic flows with complex interface［J］.Computer Methods in Applied Mechanics and Engineering,2020,365:112999.
［8］	Long T, Yang P Y, Liu M B.A novel coupling approach of smoothed finite element method with SPH for thermal fluid structure interaction problems［J］.International Journal of Mechanical Sciences,2020,174:105558.
［9］	Li G Y,Ma X J,Zhang B W,et al.An integrated smoothed particle hydrodynamics method for numerical simulation of the droplet impacting with heat transfer［J］.Engineering Analysis with Boundary Elements,2021,124:1-13.
［10］	Xu X Y,Ouyang J,Liu W M,et al.SPH simulations of 2D transient viscoelastic flows using Brownian configuration fields［J］.Journal of Non-Newtonian Fluid Mechanics,2014,208-209:59-71.
［11］	Xu X Y,Yu P.A multiscale SPH method for simulating transient viscoelastic flows using bead-spring chain model［J］.Journal of Non-Newtonian Fluid Mechanics,2016,229:27-42.
［12］	Wang Pei, Hong Tao. Parallel computation of three- dimensional smoothed particle hydrodynamics［J］.Chinese Journal of Computational Physics,2006,23(4):431-435.（in Chinese）
［13］	Cherfils J M,Pinon G,Rivoalen E.JOSEPHINE:A parallel SPH code for free-surface flows［J］.Computer Physics Communications,2012,183(7):1468-1480.
［14］	Wu Jian-song,Li Na,Liu Wen-yu,et al.Numerical modeling of flooding current flow over the underground complex staircases with a GPU-based SPH method［J］.Journal of Safety and Environment,2019,19(3):774-779.（in Chinese）
［15］	Fan Xiao-kang,Xia Ze-yu,Long Si-fan,et al. SIMD optimization of smoothed particle hydrodynamics based on ARM SVE［J］.Computer Engineering & Science,2021,43(6):989-996.（in Chinese）
［16］	Liang Lan-bo,Jin A-fang,Wen Teng-teng.Application of SPH-GPU parallel computing in wind and sand flow［J］.Computer Engineering and Applications,2022,58(1):248-254.（in Chinese）
［17］	Liu M B, Liu G R. Smoothed particle hydrodynamics (SPH):An overview and recent developments［J］.Archives of Computational Methods in Engineering,2010,17(1):25-76.
［18］	Xu X Y,Yu P.Modeling and simulation of injection molding process of polymer melt by a robust SPH method［J］.Applied Mathematical Modelling,2017,48:384-409.
［19］	Xu X Y,Yu P.Extension of SPH to simulate non-isothermal free surface flows during the injection molding process［J］.Applied Mathematical Modelling,2019,73:715-731.
［20］	Morris J P,Fox P J,Zhu Y.Modeling low Reynolds number incompressible flows using SPH［J］.Journal of Computational Physics,1997,136:214-226.
［21］	Fang J N,Owens R G,Tacher L,et al. A numerical study of the SPH method for simulating transient viscoelastic free surface flows［J］.Journal of Non-Newtonian Fluid Mechanics,2006,139(1/2):68-84.
［22］	Chen Guo-liang,An Hong,Chen Ling,et al.Practice of parallel algorithms［M］.Beijing:Higher Education Press,2004.（in Chinese）
	附中文参考文献：
［12］	王裴,洪滔.三维光滑粒子流体动力学并行计算程序［J］.计算物理,2006,23(4):431-435.
［14］	吴建松,李娜,刘文雨,等.基于SPH方法与GPU并行计算的复杂阶梯流数值模拟［J］.安全与环境学报,2019,19(3):774-779.
［15］	范小康,夏泽宇,龙思凡,等.基于ARM SVE的光滑粒子流体动力学SIMD加速方法［J］.计算机工程与科学,2021,43(6):989-996.
［16］	梁岚博,金阿芳,闻腾腾.SPH-GPU并行计算在风沙流中的应用［J］.计算机工程与应用,2022,58(1):248-254.
［22］	陈国良,安虹,陈崚,等.并行算法实践［M］.北京:高等教育出版社,2004.

[1]	刘杰, 龚春叶, 杨博, 郭晓威, 甘新标, 李胜国, 李超, 陈旭光, 肖调杰, 穆利安, 宋敏, 赵冬勇, 鞠羽中. YH-ACT：热工流体力学并行应用程序[J]. 计算机工程与科学, 2021, 43(01): 58-69.
[2]	李雯, 迟利华, 张会, 张哲, 刘杰, . 潜艇使用自航式声诱饵防御鱼雷模型并行计算方法研究[J]. 计算机工程与科学, 2020, 42(12): 2163-2168.
[3]	王俊吉1,2,朱朝艳1,2,4,陈建军1,2,郑澎3，5,徐权3. 基于无锁原子操作的多线程并行Delaunay三角化算法[J]. 计算机工程与科学, 2018, 40(05): 773-779.
[4]	章子凯，武继刚，姜文超，刘竹松. 容错处理器阵列的多逻辑列并行重构算法[J]. 计算机工程与科学, 2018, 40(01): 24-33.
[5]	李娟1，汤德佑1，傅娟2，3. 一种基于蚁群算法的生物序列并行比对方法[J]. 计算机工程与科学, 2017, 39(09): 1610-1616.
[6]	宋威，吉红蕾，李晋宏. 一种高效用项集并行挖掘算法[J]. J4, 2015, 37(03): 422-428.
[7]	杨菊，吴卓锋，王刚，刘晓光. 多核平台PAML并行算法研究[J]. J4, 2013, 35(9): 15-19.
[8]	周琛，陈振杰，张帅. 基于边界代数法的矢量栅格化并行算法设计与实现[J]. J4, 2013, 35(4): 37-41.
[9]	王多雷，龙汉，陈璇. 组合打击下一类面目标易损性的并行计算研究[J]. J4, 2013, 35(4): 42-46.
[10]	刘立国，莫锦军，袁乃昌. 基于“天河一号”的并行电磁计算求解器[J]. J4, 2013, 35(3): 43-47.
[11]	焦善飞1,2，何晨2，豆育升2，唐红1,2. 短程力分子模拟在Hadoop上的实现及优化[J]. J4, 2013, 35(2): 26-31.
[12]	龚春叶1,2,3,包为民1,闵昌万1,张烨琛1，刘杰3. 分数阶微分方程的一种细粒度数据级并行算法[J]. J4, 2013, 35(11): 76-79.
[13]	徐莹，徐磊，姜恺. 三维NavierStokes方程分步法的并行算法在异构平台上实现初探[J]. J4, 2012, 34(9): 33-39.
[14]	徐莹，徐磊. 三维Navier-Stokes方程的差分-谱方法混合法在GPU上的实现与优化[J]. J4, 2012, 34(8): 53-58.
[15]	寇大治1，梁好均2. 自洽场理论研究嵌段聚合物的并行算法实现[J]. J4, 2012, 34(8): 108-113.