Новый решатель Пуассона повысил точность моделирования коллапса массивных звёзд
Группа исследователей представила новый многоклеточный решатель уравнения Пуассона, предназначенный для задач численной релятивистики. Работа опубликована на arXiv и описывает реализацию метода, который позволяет повысить точность моделирования гравитационного коллапса массивных звёзд.
Новый алгоритм основан на многоклеточной сетке (multi-grid), что ускоряет решение эллиптических уравнений. Разработчики протестировали его на нескольких начальных задачах: двухдырочная чёрная дыра, статическая нейтронная звезда, равномерно вращающаяся нейтронная звезда и гравитационно коллапсирующая массивная звезда.
Основная демонстрация проводилась для звезды массой 9 солнечных. Моделирование учитывало перенос нейтринного излучения и релятивистскую гидродинамику. Симуляция велась до момента отскока ядра, с использованием нового решателя для сохранения ограничений при перестроении сетки.
Результаты показали, что барионная масса, масса ADM (Арновитта-Дезера-Мизнера) и момент импульса ADM сохраняются с точностью порядка 0,001% для массы и 0,01–0,1% для момента. Это значительно повышает надёжность долгосрочных симуляций.
Разработка важна для астрофизики: точное моделирование коллапса звёзд помогает лучше понять механизмы взрывов сверхновых и образование нейтронных звёзд или чёрных дыр. Новый решатель может стать стандартным инструментом в численной релятивистике.
Полный текст статьи доступен на arXiv под идентификатором 2606.12542.
