Parallel Realization of the Hybrid Model Code for Numerical Simulation of Plasma Dynamics

Abstract

The work is devoted to a parallel realization of a hybrid model for study of plasma dynamics in axially symmetric open magnetic traps. The model is based on MHD approximation for electron component of the plasma and on the kinetic approach for ion component. In the model the particle-in-cells method (PIC) with explicit numerical schemes on staggered grids for the cylindrical coordinate system is applied. The permanent injection of the particles, the long characteristic times, the necessary grid resolution and the conditional stability of the method required development of parallel version of the algorithm. In the parallelization we use the mixed decomposition with static load balancing. We present the parallel algorithm as well as results of the computational experiments on plasma dynamics in regime of diamagnetic confinement in an open magnetic trap.