Параллельный алгоритм решения задач гравитационной физики, основанный на декомпозиции области

Разработан параллельный алгоритм для решения нестационарных трехмерных задач бесстолкновительной динамики свыше $10^9$ тел в самосогласованном гравитационном поле в декартовых координатах. При распараллеливании алгоритма применен метод декомпозиции вычислительной области, что привело к созданию параллельной программы для суперкомпьютеров с распределенной памятью. Полученную программу можно использовать для расчетов на сетках, содержащих более $1000^3$ узлов. Представлено решение одной из таких задач динамики развития неустойчивости, которое получено на суперЭВМ с помощью созданной программы. Работа выполнена по проекту программ Президиума РАН “Происхождение и эволюция звезд и галактик” (координатор академик А.А. Боярчук) и интеграционному проекту СО РАН № 26 (координатор академик Б.Г. Михайленко).