A Radiative Transfer Module for Relativistic Magnetohydrodynamics in the PLUTO Code
A Radiative Transfer Module for Relativistic Magnetohydrodynamics in the PLUTO Code
复制标题
PLUTO 代码中相对论磁流体动力学的辐射传输模块
DOI：
10.3847/1538-4365/ab18ff
复制
发表时间：
2019
期刊：
影响因子：
--
通讯作者：
A. Mignone
中科院分区：
文献类型：
--
作者：
J. D. Melon Fuksman;A. Mignone
我们提出了一个求解相对论辐射流体动力学和磁流体动力学方程的数值实现，它被设计为一个独立的模块，位于可免费获得的代码PLUTO中。辐射传递方程在灰色近似下求解，并施加M1闭包，这允许在自由流和扩散极限下处理辐射输运。方程按照隐式-显式格式进行积分，其中辐射-物质相互作用项隐式积分，而输运和所有剩余的源项则通过PLUTO中包含的相同godunov型求解器显式求解。其中，我们介绍了一种新的用于光薄辐射输运的Harten-Lax-van Leer-contact （HLLC）求解器。该代码适用于使用单个处理器或并行架构进行笛卡尔坐标、球面坐标和柱坐标中的多维计算。自适应网格计算也可以通过CHOMBO库实现。该算法的性能通过一系列数值基准来证明，通过研究各种不同的配置，特别强调解在自由流和扩散极限下的行为。
We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer equations are solved under the gray approximation and imposing the M1 closure, which allows the radiation transport to be handled in both the free-streaming and diffusion limits. Equations are integrated following an implicit–explicit scheme, where radiation–matter interaction terms are integrated implicitly, whereas transport and all of the remaining source terms are solved explicitly by means of the same Godunov-type solvers included in PLUTO. Among these, we introduce a new Harten–Lax–van Leer–contact (HLLC) solver for optically thin radiation transport. The code is suitable for multidimensional computations in Cartesian, spherical, and cylindrical coordinates using either a single processor or parallel architectures. Adaptive grid computations are also made possible by means of the CHOMBO library. The algorithm performance is demonstrated through a series of numerical benchmarks by investigating various different configurations with a particular emphasis on the behavior of the solutions in the free-streaming and diffusion limits.