Numerical simulation of \(N\)-body gravitational dynamics by PIC method: a parallel implementation. (Russian. English summary) Zbl 1027.85003

Summary: The evolution of self-gravitating systems such as the accretion discs is of great interest to astrophysics. The aim of this work is to create a parallel program for the accretion disc simulation on high-performance multiprocessor computers. The disc structure formation is an \(N\)-body problem in a self-consistent gravity field. A good approximation to the problem is the Vlasov-Liouville kinetic equation. In the present work, the equation is solved by the PIC method. The main difficulty here is the evaluation of the gravitational potential which is given by the 3D Poisson equation. The parallel scheme of the algorithm was designed for the MIMD computers in an assembly technology. This means that the program is assembled of minimal fragments, each being a ready-made program containing potential values and the particles from one or more grid layers. The values of a grid potential are uniformly distributed among the processor elements uniformly in the radial direction. As the potential evaluation takes the main time, the distribution of particles is of minor importance here. Test computations conducted on the ICT cluster of Pentium-III workstations showed the linear acceleration as compared to the sequential version.


85-08 Computational methods for problems pertaining to astronomy and astrophysics
85A05 Galactic and stellar dynamics
68W10 Parallel algorithms in computer science
70F10 \(n\)-body problems