Тестирование магнитогидродинамического кода на задачах естественной конвекции и геодинамо

В рамках метода контрольного объема разработан программный код для численного решения задач неидеальной магнитной гидродинамики вязкой несжимаемой жидкости на структурированных разнесенных сетках в сферических координатах. При дискретизации уравнения индукции магнитного поля использован алгоритм ограниченного переноса (Сonstrained Transport Algorithm) и схема QUICK с методом отложенной коррекции для аппроксимации конвективных членов. Для решения уравнений гидродинамики использован алгоритм SIMPLER. Программный код разработан для моделирования естественной конвекции и гидромагнитного динамо во вращающемся шаре или сферическом слое. Представлены результаты решения тестовых задач естественной конвекции и геодинамо с вакуумными граничными условиями, демонстрирующие достаточно точное соответствие результатам эталонных расчетов. Программное обеспечение разработано для ускорителей вычислений, поддерживающих технологию CUDA, с использованием набора расширений к языку программирования Фортран.