广义相对论守恒流体动力学程序中基于M1闭合处理辐射的半隐式格式
高能天体物理现象
2013-01-16 v2
摘要
本文描述了一种在多维广义相对论守恒流体动力学程序中包含辐射的数值格式。该方法采用协变形式的M1闭合方案来闭合辐射矩，并对辐射源项进行半隐式处理，以同时应对光学薄和光学厚两种状态。该格式已在一个守恒型广义相对论辐射流体动力学程序KORAL中实现。通过一系列测试问题展示了程序的稳健性，包括辐射相对论激波管、静态辐射压支撑大气、阴影、弯曲时空中的光束以及辐射邦迪吸积。讨论了M1闭合相对于爱丁顿闭合和通量限制扩散等其他方法的优势，并指出了其局限性。
引用
@article{arxiv.1212.5050,
title = {Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes},
author = {Aleksander Sadowski and Ramesh Narayan and Alexander Tchekhovskoy and Yucong Zhu},
journal= {arXiv preprint arXiv:1212.5050},
year = {2013}
}
备注
19 pages, 15 figures, as published in MNRAS