Явный многошаговый алгоритм для моделирования динамики самогравитирующего газа

Для численного моделирования динамики самогравитирующего газа развит многошаговый явный метод и построен пространственно трехмерный код на декартовой сетке в эйлеровых переменных. В основе численного алгоритма лежит метод Белоцерковского–Давыдова (FLIC) с использованием разностной схемы первого порядка аппроксимации с сеточной вязкостью, которая обеспечивает устойчивость численных решений. Алгоритм обладает свойством сходимости на последовательности сгущающихся сеток. Приведены результаты расчетов коллапса изотермического газа, вращения изоэнтропического самогравитирующего газа и его автомодельного расширения. Работа выполнена при поддержке Программы Президиума РАН “Происхождение, строение и эволюция объектов Вселенной”, Программы Президиума РАН “Происхождение биосферы и геобиологическая эволюция”, Интеграционного проекта СО РАН No. 26 “Математические модели, численные методы и параллельные алгоритмы для решения больших задач СО РАН и их реализация на многопроцессорных суперЭВМ”. Расчеты проведены в Сибирском суперкомпьютерном центре РАН на машине с общей памятью SMP16x256.