Газодинамическая модель аккреции на нейтронную звезду с учетом вязкости и влияние крупномасштабных вихрей на передачу момента импульса

Работа посвящена построению газодинамической модели аккреционного диска вокруг нейтронной звезды (НЗ). Разработанный многомерный код используется для исследования устойчивости стационарных аксиально симметричных моделей путем проведения эволюционных расчетов в 3D геометрии с учетом вязкости, а также с учетом диффузии излучения в 2D. Показано, что произвольный стационарный аксиально симметричный диск с монотонным уменьшением плотности с цилиндрическим радиусом переходит, благодаря вязкости, торможению и растеканию вещества по НЗ в новую квазистационарную тороидальную конфигурацию. Исследование устойчивости стационарной тороидальной конфигурации подтвердило формирование крупномасштабных вихревых структур при первоначальном периодическом возмущении угловой скорости по азимуту теперь уже с учетом «турбулентной» вязкости, а не схемной. Оказалось, что наличие крупномасштабных структур приводит к ускорению торможения, т. е. к эффективному увеличению вязкости.