The new direct simulation method of Monte-Carlo type (DSMC) for Coulomb collisions in the case of two component plasma is considered. A brief literature review and preliminary information concerning the problem are given. Further the idea that lies in the basis of the method is discussed and its scheme is provided. The illustrative numerical simulation of the initial distribution relaxation for one and two sorts of particles in 3D case in the velocity space is performed. Simulation results are compared with the numerical results based on the completely conservative finite difference schemes for the Landau-Fokker-Planck equation. Estimation of calculation accuracy obtained from numerical results is given.

Mathematical modelling in actual problems of science and technics