Prof. Dr. Thorsten Pöschel
Institute for Multiscale Simulation
Friedrich-Alexander Universität Erlangen-Nürnberg
The moving sand dunes, avalanches and the dense planetary rings are granular systems which consists of countless number of particles. The study of such huge systems is impossible with current particle based methods such as Discrete Element Method (DEM) and Event-Driven Molecular Dynamics (EDMD). However, a hydrodynamical (HD) method is a suitable alternative to study these phenomena. The hydrodynamic equations of the granular flow, describes the system by average field variables like particle number density, flow velocity and temperature. The time evolution of these quantities is governed by hydrodynamic equations. The solution to these equations provides valuable information about the granular systems.
We developed a granular flow solver that is based on C++ and parallelized with MPI standard to solve the hydrodynamic equations of the granular flow. We did several optimizations to increase the performance and parallel efficiency of the Granular Flow Solver. The optimizations include modifying the MPI communication pattern, the inlining of the functions, memory management and eliminating the expensive operations. These optimizations improved the single-core and parallel efficiency of the solver significantly.
The current solver is a powerful tool to simulate a vast range of systems with various boundary conditions including the shaking boundaries which is widely used in the study of granular materials. Figure 1 shows a comparison between the density field of the HD solver and the well-known DEM simulation. The accuracy and performance of the current solver is benchmarked with different test cases. Despite all the approximations in HD model, all tests show a very good agreement with the particle based simulations.