Abstract. Computer simulation of dynamic of plasma target for highly efficient neutralization of powerful negative ion beams is considered. The plasma is confined within a magnetic trap with multipole magnetic walls. Mathematical model is based on the Boltzmann equation for the distribution functions for ions and electrons and system of the Maxwell's equations for the selfconsistent electromagnetic fields. The combination of the modified PIC-method in the cylindrical R-Z coordinates and the MonteCarlo methods is used to solve these equations. The complex nature of the processes studied, and also the need of calculation of trajectories of billions of particles required the use scalable parallel algorithm. The use of modern supercomputers has allowed to calculate plasma dynamics, to determine plasma streams both on the walls of the trap and through end holes.