We present an explicit incompressible smoothed particle hydrodynamics formulation with stabilized pressure distribution and its implementation in a multiple graphics processing unit environment. The pressure Poisson equation is stabilized via both pressure invariance and divergence-free conditions, and its explicit formulation is derived using the first step of the Jacobi iterative solver. Also, we show how to adapt the fixed wall ghost particle for the boundary condition into our explicit approach. Verification and validation of the method include hydrostatic and dam break numerical tests. The computational performance in the multi-GPU environment was notably high with reasonable speedup values compared to our single-GPU implementation. In particular, our code allows simulations with very large number of particles reaching up to 200 million per GPU card. Finally, to illustrate the potential of our formulation in simulating natural disasters, we present a simulation of the famous Fukushima Dai-ichi Power Plant inundation by the tsunami from The Great East Japan Earthquake in 2011, in Japan.

我们提出了一种具有稳定压力分布的显式不可压缩平滑粒子流体动力学公式，并在多图形处理单元环境中实现了该算法。压力泊松方程可通过压力不变性和无散度条件进行稳定，并使用Jacobi迭代求解器的第一步导出其明确的公式。同样，我们展示了如何将边界条件的固定壁重影粒子应用于我们的显式方法。该方法的验证和确认包括静水压力试验和溃坝数值试验。与我们的单GPU实现相比，多GPU环境中的计算性能以合理的加速值显着提高。尤其是，我们的代码允许对非常大量的粒子进行仿真，每个GPU卡最多可以达到2亿个粒子。最后，为说明我们的公式在模拟自然灾害中的潜力，我们对2011年日本东日本大地震海啸淹没的著名福岛第一核电站进行了模拟。