Моделирование Неравновесной Кинетики И Спектров Сверхновых: Метод Розенброка С Контролем Положительности Решений

Расчет спектров сверхновых на фотосферной стадии требует детального моделирования неравновесной населенности атомных уровней в быстро расширяющейся оболочке. В работе рассматривается фундаментальная задача интерпретации наблюдаемых спектров через решение сверхжесткой системы нестационарных кинетических уравнений в расширяющейся оболочке, где характерные времена атомных процессов на много порядков отличаются от гидродинамического времени. На основе анализа сингулярно возмущенной системы мы вводим эвристическое правило для селективной деактивации малонаселенных уровней, что снижает размерность системы и вычислительные затраты без потери точности. Чрезвычайно большой разброс характерных времен процессов приводит к численной неустойчивости, проявляющейся в нефизичных отрицательных населенностях и “срыве” на ложное решение. Мы представляем модифицированный метод Розенброка четвертого порядка, дополненный процедурой коррекции, гарантирующей положительность решений на основе работ Гольдина и Калиткина. Эта коррекция действует как специализированный контроллер шага, автоматически предотвращая выбор слишком больших шагов интегрирования, которые ведут к неустойчивости. Метод реализован в коде LEVELS и апробирован на модели сверхновой типа II-P SN 2018aoq. Показано, что гибридный подход устойчиво интегрирует систему из более чем 1500 уравнений и позволяет удовлетворительно воспроизводить особенности наблюдаемых оптических спектров, что открывает возможности для уточнения гидродинамических параметров взрыва по спектроскопическим данным.