论文部分内容阅读
本文基于高分辨率激波格式(High resolution shock capturing schemes，HRSC)、辐射扩散近似(Flux-limited diffusion，FLD)和局域热动平衡(Local thermody-namic equilibrium，LTE)假设设计了一种多维辐射磁流体程序。HRSC是一种守恒代数方法，可确保总能量守恒。这里我们采用三阶Runge-Kutte格式和HLLE黎曼问题求解法。在FLD近似下，可以忽略辐射能量流的演化方程，辐射能量流正比于辐射能量密度梯度。因此，在FLD近似下辐射磁流体程序容易被实施，尽管FLD近似在处理光薄区域内高度各向异性的辐射场时精度较低。本文还利用经典的磁流体问题和辐射激波对程序做了测试，结果表明，程序具备足够的精度和效率求解辐射扩散近似下的辐射主导的天体物理过程。