PhasicFlow: A parallel, multi-architecture open-source code for DEM simulations

• Published in: Computer physics communications 2023-10, Vol.291, p.108821, Article 108821
• Main Author: Norouzi, H.R.
• Format: Article
• Language: English
• Subjects: CUDA; Discrete element method; GPU; Granular flow; OpenMP; Parallel computing
• ISSN: 0010-4655
• DOI: 10.1016/j.cpc.2023.108821

PhasicFlow is an open-source code for DEM simulations developed in C++, licensed under GPL V-3.0 [www.github.com/PhasicFlow/phasicFlow]. It is parallelized based on shared-memory model and can be executed on multi-core CPUs or GPUs. The philosophy of its development is to provide an efficient, easy-to-use, and easy-to-extend tool for a wide range of hardware including cheap to high-end desktop computers. Here, the code structure, the most important features, its performance for simulating small to large problems (250,000 to 80,000,000 particles) on various hardware are presented. In all cases, the execution time varied linearly with number of particles that is ideal for DEM. The computation load per particle was almost constant for all problem sizes except for small simulations (higher). For CPU simulations, a good scaling was observed with number of cores, and scaling was closer to ideal state for larger simulations. The performance of PhasicFlow was compared with similar tools, wherever possible. It was shown that the performance of PhasicFlow is much higher than the similar tools in CPU execution mode and it is around 10–15% faster than similar tools in GPU execution mode. This new tool broadens the application of computer simulation to larger problems and allows researchers to perform DEM simulations on a wider range of hardware. Program title: PhasicFlow CPC Library link to program files:https://doi.org/10.17632/fs7w7jb8mr.1 Developer's repository link:www.github.com/PhasicFlow/phasicFlow Licensing provisions: GNU General Public License 3 Programming language: C++ Nature of problem: The code covers the simulation of granular matter in powder, gas-solid and liquid-solid flows. The formulation is based on the soft-sphere discrete element method (DEM) which is a powerful tool for simulating solid powder flows. Solid particles are considered distinct bodies that may interact with solid surfaces, other particles and fluid around them. This method can be purely used for investigating various flow problems or can be coupled with proper fluid flow simulation code to be used in multiphase flow systems. The main goal here is to provide a powerful tool for researchers to perform large-scale powder and multiphase flows on cheaper and more accessible hardware. Solution method: The solution method is based on the high-order explicit integration of equations of motions for discrete particles over time to obtain their linear and rotational position and velocity. Tria