Сформулирован (и методически разработан) подход к оптимизации процедуры решения уравнений движения частиц дискретной самосогласованной модели разреженной плазмы на основе сопряжения разных по характеру и вычислительным свойствам динамических схем. Концепция алгоритмически реализована и опробована в математическом (по методу макрочастиц) моделировании низкочастотной кинетической неустойчивости Вайбеля. Показано существенное сокращение времени симуляции при использовании указанной методики оптимизации по сравнению с традиционным счетом на основе только одной (явной или неявной) динамической схемы. В ходе компьютерных экспериментов выявлены особенности поведения вайбелевской раскачки на стадии ее насыщения: присутствие затухающих низкочастотных колебаний плотности энергии магнитного поля - ключевого параметра неустойчивости.
$^1$Кафедра математики Физического факультета МГУ