Geodesic
Numerical simulation of N-body gravitational dynamics by PIC method: a parallel implementation
Sibirskij žurnal vyčislitelʹnoj matematiki, Tome 6 (2003) no. 1, pp. 25-36 
V. A. Vshivkov; V. È. Malyshkin; A. V. Snytnikov; V. N. Snytnikov. Numerical simulation of N-body gravitational dynamics by PIC method: a parallel implementation. Sibirskij žurnal vyčislitelʹnoj matematiki, Tome 6 (2003) no. 1, pp. 25-36. http://geodesic.mathdoc.fr/item/SJVM_2003_6_1_a2/
@article{SJVM_2003_6_1_a2,
     author = {V. A. Vshivkov and V. \`E. Malyshkin and A. V. Snytnikov and V. N. Snytnikov},
     title = {Numerical simulation of {N-body} gravitational dynamics by {PIC} method: a~parallel implementation},
     journal = {Sibirskij \v{z}urnal vy\v{c}islitelʹnoj matematiki},
     pages = {25--36},
     year = {2003},
     volume = {6},
     number = {1},
     language = {ru},
     url = {http://geodesic.mathdoc.fr/item/SJVM_2003_6_1_a2/}
}
TY  - JOUR
AU  - V. A. Vshivkov
AU  - V. È. Malyshkin
AU  - A. V. Snytnikov
AU  - V. N. Snytnikov
TI  - Numerical simulation of N-body gravitational dynamics by PIC method: a parallel implementation
JO  - Sibirskij žurnal vyčislitelʹnoj matematiki
PY  - 2003
SP  - 25
EP  - 36
VL  - 6
IS  - 1
UR  - http://geodesic.mathdoc.fr/item/SJVM_2003_6_1_a2/
LA  - ru
ID  - SJVM_2003_6_1_a2
ER  - 
%0 Journal Article
%A V. A. Vshivkov
%A V. È. Malyshkin
%A A. V. Snytnikov
%A V. N. Snytnikov
%T Numerical simulation of N-body gravitational dynamics by PIC method: a parallel implementation
%J Sibirskij žurnal vyčislitelʹnoj matematiki
%D 2003
%P 25-36
%V 6
%N 1
%U http://geodesic.mathdoc.fr/item/SJVM_2003_6_1_a2/
%G ru
%F SJVM_2003_6_1_a2

Voir la notice de l'article provenant de la source Math-Net.Ru

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 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 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 readymade 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.

[1] Bertin G., Dynamics of Galaxies, University Press, Cambridge, 2000 | MR | Zbl

[2] Snytnikov A. V., “Modelirovanie evolyutsii galakticheskogo diska metodom chastits v yacheikakh”, Trudy konferentsii molodykh uchenykh, posvyaschennoi 10-letiyu IVT SO RAN. Matematicheskoe modelirovanie, T. 2, SO RAN IVT, Novosibirsk, 2000, 145–147

[3] Berezin Yu. B., Vshivkov V. A., Metod chastits v dinamike razrezhennoi plazmy, Nauka, Novosibirsk, 1980 | MR

[4] Khokni R., Istvud Dzh., Chislennoe modelirovanie metodom chastits, Mir, M., 1987

[5] Grigorev Yu. N., Vshivkov V. A., Chislennye metody “chastitsy-v-yacheikakh”, Nauka, Novosibirsk, 2000 | MR

[6] Boris J. P., “Relativistic plasma simulation - optimization of a hybrid code”, Proc. Fourth Conf. Num. Sim. Plasmas, Washington, 1970, 3–67

[7] Ilin V. P., Metody nepolnoi faktorizatsii dlya resheniya algebraicheskikh sistem, Fizmatlit, M., 1995 | MR

[8] Sellwood J. A., “Spiral Structure as a Recurrent Instability”, Astrophysical Dynamics – in commemoration of F. D. Kahn, Kluwer, Dordrecht, 1999

[9] Kraeva M. A., Malyshkin V. E., “Implementation of PIC Method on MIMD Multicomputers with Assembly Technology”, Proc. HPCN Europe 1997, LNCS, 1255, Springer Verlag, 1997, 541–549

[10] Kraeva M. A., Malyshkin V. E., “Assembly Technology for Parallel Realization of Numerical Models on MIMD-Multicomputers”, The special issue of the International Journal on Future Generation Computer Systems devoted to Parallel Computing Technologies, 17 (2001), 755–765 | Zbl

[11] Kraeva M. A., Malyshkin V. E., “Dinamicheskaya balansirovka zagruzki v realizatsii PIC-metoda na MIMD multikompyuterakh”, Programmirovanie, Vyp. 1, Nauka, M., 1999, 47–53