Модификация параллельного алгоритма для численного моделирования открытой магнитной ловушки

Работа посвящена численному моделированию динамики плазмы в диамагнитном режиме в открытой магнит-ной ловушке. Используемая гибридная модель основана на кинетическом описании ионной компоненты и магнитно-гидродинамическом описании электронной компоненты плазмы. Недостатком модели является условная устойчивость и соответствующие требования к временному шагу. На практике измельчение сетки в два раза требует уменьшения временного шага в 6 раз. Для характерных времен плазменных процессов по-рядка 102 обратных ионных циклотронных частот вычислений на сетке 100 × 500 требует несколько дней. При этом более 85 % времени работы программы занимают процедуры обработки частиц, поэтому их эффек-тивная реализация играет ключевую роль в сокращении времени расчетов. В предлагаемом вниманию алго-ритме мы сочетаем динамическую загрузку балансировки и векторизацию вычислений расчета плотности электронов и плотности токов. Представлены результаты численных экспериментов с учетом существенно неравномерного распределения частиц в области и увеличивающегося их числа за счет постоянной инжекции.