Новый код на основе метода Монте-Карло для решения уравнений лучистого переноса
AstroNews.ru - 09 Июня 2021 11:16:41
В новом исследовании Сяолинь Ян (Yang Xiaolin) и его коллеги из Юньнаньской астрономической обсерватории Академии наук Китая разработали новый быстрый код под названием Lemon (Linear Integral Equations" Monte Carlo Solver Based on Neumann Solution), который предназначен для точного решения уравнений лучистого переноса. Схема построения этого кода основана на линейном интегральном уравнении и серии его решений, называемой серией Ньюмана.
Лучистый перенос представляет собой широко распространенный в астрофизике процесс, который играет большую роль как в теоретических исследованиях, так и в практических наблюдениях. Для решения уравнений лучистого переноса были предложены различные методы, среди которых метод Монте-Карло является наиболее важным и широко используемым численным методом из-за его простоты и высокой производительности.
Обычная реализация метода Монте-Карло (или схема отслеживания фотонов), однако, имеет неотъемлемый недостаток, состоящий в том, что в результате проведения большого количества вычислений достигается не очень устойчивое решение с высокими величинами отклонений, поэтому большая часть вычислительных мощностей расходуется впустую.
Для преодоления этого недостатка Сяолинь Ян и его коллеги предложили новую схему, в которой они предлагают использовать для решения уравнений лучистого переноса интегральное уравнение и его ньюмановское решение вместо отслеживания фотонов.
У этой новой схемы есть значительные преимущества. Она предполагает учет вклада фотонов в итоговый результат на каждом участке рассеяния, что существенно повышает эффективность и точность расчетов. В результате описанный выше недостаток метода удается устранить. Метод может быть использован как для уравнений лучистого переноса, в которых присутствует поляризация излучения, так и для уравнений, описывающих неполяризованное излучение, в рамках единого подхода. Процедура расчета упрощается, если система имеет осевую или сферическую симметрию. Кроме того, метод может быть применен напрямую для решения любых дифференциально-интегральных уравнений с корректно указанными начальными или граничными условиями.
Код Lemon разработан полностью на этих новых принципах и написан на языке FORTRAN 90. Он находится в открытом доступе и может быть скачан по ссылке: github.com/yangxiaolinyn/Lemon
Исследование появилось на сервере предварительных научных публикаций arxiv.org.