Monte-Carlo method for two component plasmas

The direct simulation method of Monte-Carlo type for Coulomb collisions in the case of two component plasma is considered. The illustrative numerical simulation of the initial distribution relaxation for two sorts of particles in 3D 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 is given for the wide range of numerical parameters.