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

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

Research and application of a SPH parallelalgorithm based on particle decomposition

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.