Implementation of the Kinetic Plasma Code with Locally Recursive non-Locally Asynchronous Algorithms

25th IUPAP Conference on Computational Physics

物理学;计算机科学

Perepelkina, A.Yu.^1 ; Levchenko, V.D.^1,2 ; Goryachev, I.A.^1

Keldysh Institute of Applied Mathematics, Moscow, Russia^1

National Research Nuclear University MEPhI, Moscow, Russia^2

关键词: Asynchronous algorithms;    Dependency graphs;    Filamentation instabilities;    Laser interaction;    Laser pulse amplitudes;    Multiscale problem;    Numerical experiments;    Particle in cell codes;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/510/1/012042/pdf
DOI  :  10.1088/1742-6596/510/1/012042

学科分类：计算机科学（综合）

来源: IOP

【 摘 要 】

Numerical simulation is presently considered impractical for several relevant plasma kinetics problems due to limitations of computer hardware even with the use of supercomputers. To overcome the existing limitations it is suggested to develop algorithms which would effectively utilize the computer memory subsystem hierarchy to optimize the dependency graph traversal rules. The ideas for general cases of numerical simulation and implementation of such algorithms to particle-in-cell code is discussed in the paper. This approach enables the simulation of previously unaccessible for modeling problems and the execution of series of numerical experiments in reasonable time. The latter is demonstrated on a multiscale problem of the development of filamentation instability in laser interaction with overdense plasma. One variant of the simulation with parameters typical for simulations on supercomputers is performed with the use of one cluster node. The series of such experiments revealed the dependency of energy loss on incoming laser pulse amplitude to be nonmonotonic and reach over 4%, an interesting result for research of fast ignition concept.