论文标题
具有高阶方案的基于蒙特卡洛的相对论辐射流体动力学代码
A Monte-Carlo based relativistic radiation hydrodynamics code with a higher-order scheme
论文作者
论文摘要
我们根据蒙特卡洛算法开发了一种新的相对论辐射流体动力学代码。在此代码中,我们实施了一个新方案,以在大量数据包号码的限制中实现二阶精度,以解决物质与辐射之间的相互作用。这种高阶时间集成方案的实施方式是以确保能量弹药保护到地球积分器的精度的方式。辐射过程的空间依赖性(例如数据包传播,发射,吸收和散射)也被考虑到二阶精度。在先前的研究之后,我们通过解决各种测试问题来验证我们的代码;单区热化,动力扩散,辐射拖动,辐射介导的冲击管,光学厚极限的冲击管以及爱丁顿的极限问题。我们表明,我们的代码以合理的准确性重现了身体上适当的结果,并且还证明了我们对一区和一维问题的实施确实可以实现时空的二阶精度。
We develop a new relativistic radiation hydrodynamics code based on the Monte-Carlo algorithm. In this code, we implement a new scheme to achieve the second-order accuracy in time in the limit of a large packet number for solving the interaction between matter and radiation. This higher-order time integration scheme is implemented in the manner to guarantee the energy-momentum conservation to the precision of the geodesic integrator. The spatial dependence of radiative processes, such as the packet propagation, emission, absorption, and scattering, are also taken into account up to the second-order accuracy. We validate our code by solving various test-problems following the previous studies; one-zone thermalization, dynamical diffusion, radiation dragging, radiation mediated shock-tube, shock-tube in the optically thick limit, and Eddington limit problems. We show that our code reproduces physically appropriate results with reasonable accuracy and also demonstrate that the second-order accuracy in time and space is indeed achieved with our implementation for one-zone and one-dimensional problems.