новости космоса
1459 6
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.


(Добавил: Hot Temp)

комментарии
1
Stan65 · 09-06-2021

Эта картинка похожа на параграф инопланетного многомерного языка, который я пытался был постигнуть (безрезультатно, ибо он слишком сложный, кроме того мой земной образовательный бекграунд очень мешает). Многомерные языки надо давать учить ребёнку. До его знакомства с языком людей)

2
Bechmet69 · 09-06-2021

Мы знаем то что мы ничего не знаем, попытку знать никто не отменял

3
Teddy40 · 09-06-2021

Линейные интегральные уравнения решаются и так вполне себе замечательно. Вся беда в том, что реальные уравнения нелинейны.

4
Teddy40 · 09-06-2021

Как они из Linear Integral Equations" Monte Carlo Solver Based on Neumann Solution
получили Lemon?

5
koftochka102 · 10-06-2021

Непонятнамана

6
ostman14 · 10-06-2021

Понятно только Монте-Карло. Там вроде казино есть.

написать комментарий наверх
Для добавления комментария необходимо зарегистрироваться, а затем войти на сайт используя свой логин и пароль.

Если Вы уже зарегистрировались, но забыли пароль - воспользуйтесь нашим разделом восстановления пароля.

© 2002-2021. Все права защищены. AstroNews.ru | Перепечатка любых материалов сайта без разрешения редакции запрещена!