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

Предложена параллельная реализация метода RKDG (Runge-Kutta Discontinuous Galerkin) для системы уравнений идеальной магнитной гидродинамики в двумерной осесимметричной постановке на неструктурированных треугольных сетках. Разработан алгоритм бездивергентной реконструкции магнитного поля для цилиндрической системы координат, позволяющий получать физически адекватные результаты расчетов с высоким порядком точности. Показана высокая масштабируемость созданного параллельного кода на многопроцессорных установках, проведен анализ эффективности распараллеливания программы на кластере К-100 Института прикладной математики им. М.В. Келдыша РАН. Представлены и обсуждены результаты моделирования магниторотационной неустойчивости в аккрецирующей оболочке протозвезды.