我们基于蒙特卡罗算法开发了一个新的相对论辐射流体动力学程序。在这个代码中,我们实现了一个新的方案,以在大数据包数的限制下实现时间上的二阶精度,从而解决物质与辐射之间的相互作用。这种高阶时间积分方案的实现方式是确保能量动量守恒达到测地线积分器的精度。辐射过程的空间依赖性,如包传播、发射、吸收和散射,也被考虑到二阶精度。我们根据之前的研究,通过在固定背景度量上解决各种测试问题来验证我们的代码;单区热化、动力学扩散、辐射拖曳、辐射介导激波管、光学厚度极限中的激波管以及Eddington极限问题。我们证明了我们的代码以合理的精度再现了物理上适当的结果,并证明了我们对单区和一维问题的实现确实实现了时间和空间上的二阶精度。
3更多- 收到日期:2022年9月26日
- 2023年1月9日接受
内政部:https://doi.org/10.103/PhysRevD.107.023026
引力、宇宙学和天体物理学