An approach to the optimization of the procedure for solving the equations of motion of particles of a discrete self-consistent model of a rarefied plasma based on the conjugation of dynamic schemes of different nature and computational properties is developed. The concept was algorithmically implemented and tested in mathematical modeling (using the PIC method) of the low-frequency kinetic Weibel instability. A significant reduction in the simulation time is shown when using the specified optimization technique in comparison with traditional computation based on only one (explicit or implicit) dynamic scheme. During computer experiments, the features of the behavior of the Weibel process at the stage of its saturation were revealed: the presence of damped low-frequency oscillations of the energy density of the magnetic field, the key parameter of instability.
$^1$Department of Mathematics, Faculty of Physics, Moscow State University