Modification of Parallel Algorithm for Numerical Modeling of Open Magnetic Trap

Abstract

The work is devoted to the numerical simulation of plasma dynamics in the diamagnetic regime in an open magnetic trap. The hybrid model we use is based on the kinetic description of the ion and magnetohydrodynamic description of the electron plasma component. The disadvantage of the model is the conditional stability and the corresponding requirements for the time step. In practice, the double mesh nodes increase requires reducing the time step by factor of 6. For the typical times of the plasma processes of the order of 102 reciprocal ion cyclotron frequencies, the calculations on 100 × 500 grid require several days. At the same time, more than 85 % of the program operation time is occupied by the particle processing procedures, so their efficient implementation plays a key role in the decrease of the calculation times. In the proposed algorithm, we combine dynamic load balancing and vectorization of calculations for the electron density and current density computations. The results of numerical experiments are presented taking into account the substantially nonuniform distribution of the particles in the region and their increasing number due to the constant injection.