A New Fast Monte Carlo Code for Solving Radiative Transfer Equations Based on the Neumann Solution
文献类型:期刊论文
作者 | Yang XL(杨晓林)1,2,3,4; Wang JC(王建成)1,2,3,4; Yang CY(杨初源)1,2,3; Yuan ZL(袁尊理)1,2,3 |
刊名 | ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES |
出版日期 | 2021-06-01 |
卷号 | 254期号:2页码:38 |
ISSN号 | 0067-0049 |
DOI | 10.3847/1538-4365/abec73 |
产权排序 | 第1完成单位 |
文献子类 | Article |
英文摘要 | In this paper, we propose a new Monte Carlo radiative transport (MCRT) scheme, which is based completely on the Neumann series solution of the Fredholm integral equation. This scheme indicates that the essence of MCRT is the calculation of infinite terms of multiple integrals in the Neumann solution simultaneously. Under this perspective, we redescribe the MCRT procedure systematically, in which the main work amounts to choosing an associated probability distribution function for a set of random variables and the corresponding unbiased estimation functions. We select a relatively optimal estimation procedure that has a lower variance from an infinite number of possible choices, such as term-by-term estimation. In this scheme, MCRT can be regarded as a pure problem of integral evaluation, rather than as the tracing of random-walking photons. Keeping this in mind, one can avert some subtle intuitive mistakes. In addition, the d functions in these integrals can be eliminated in advance by integrating them out directly. This fact, together with the optimal chosen random variables, can remarkably improve the Monte Carlo (MC) computational efficiency and accuracy, especially in systems with axial or spherical symmetry. An MCRT code, Lemon (Linear integral Equations' Monte carlo solver based On the Neumann solution; the code is available on the GitHub codebase at https://github.com/yangxiaolinyn/Lemon, and version 2.0 is archived on Zenodo at https://doi.org/10.5281/zenodo.4686355), has been developed completely based on this scheme. Finally, we intend to verify the validation of Lemon; a suite of test problems mainly restricted to a flat spacetime has been reproduced, and the corresponding results are illustrated in detail. |
学科主题 | 天文学 ; 天体物理学 ; 高能天体物理学 |
URL标识 | 查看原文 |
出版地 | TEMPLE CIRCUS, TEMPLE WAY, BRISTOL BS1 6BE, ENGLAND |
WOS关键词 | COVARIANT MAGNETOIONIC THEORY ; INVERSE COMPTON-SCATTERING ; X-RAY-POLARIZATION ; KERR BLACK-HOLE ; SGR A-ASTERISK ; ACCRETION DISKS ; FORMAL SOLUTIONS ; PUBLIC CODE ; LIGHT ; TRANSPORT |
资助项目 | National Natural Science Foundation of China[U2031111] ; National Natural Science Foundation of China[11573060] ; National Natural Science Foundation of China[12073069] ; National Natural Science Foundation of China[11661161010] ; National Natural Science Foundation of China[U1838116] ; National Natural Science Foundation of China[11673060] ; National Natural Science Foundation of China[U1931204] ; Yunnan Natural Science Foundation[2019FB008] |
WOS研究方向 | Astronomy & Astrophysics |
语种 | 英语 |
出版者 | IOP PUBLISHING LTD |
WOS记录号 | WOS:000655269700001 |
资助机构 | National Natural Science Foundation of China[U2031111, 11573060, 12073069, 11661161010, U1838116, 11673060, U1931204] ; Yunnan Natural Science Foundation[2019FB008] |
源URL | [http://ir.ynao.ac.cn/handle/114a53/24395] |
专题 | 云南天文台_高能天体物理研究组 云南天文台_中国科学院天体结构与演化重点实验室 |
通讯作者 | Yang XL(杨晓林) |
作者单位 | 1.Yunnan Observatories, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming, 650216, People's Republic of China 2.Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming, 650216, People's Republic of China 3.Center for Astronomical Mega-Science, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing, 100012, People's Republic of China 4.University of the Chinese Academy of Sciences, Beijing, 100049, People's Republic of China |
推荐引用方式 GB/T 7714 | Yang XL,Wang JC,Yang CY,et al. A New Fast Monte Carlo Code for Solving Radiative Transfer Equations Based on the Neumann Solution[J]. ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES,2021,254(2):38. |
APA | Yang XL,Wang JC,Yang CY,&Yuan ZL.(2021).A New Fast Monte Carlo Code for Solving Radiative Transfer Equations Based on the Neumann Solution.ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES,254(2),38. |
MLA | Yang XL,et al."A New Fast Monte Carlo Code for Solving Radiative Transfer Equations Based on the Neumann Solution".ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES 254.2(2021):38. |
入库方式: OAI收割
来源:云南天文台
其他版本
除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。