论文部分内容阅读
辐射转移问题一直以来都是天体物理学中重要的内容之一,而求解辐射转移方程是辐射转移理论的核心问题。求解转移方程有助于我们得到恒星大气的物理结构和化学组成,以及热动力学量(如温度,压强,密度等)的信息,这些结果对研究恒星演化具有非常重要的意义。本文主要对恒星大气层中辐射转移方程进行研究,并对其进行计算所涉及的物理量予以详细介绍,其中包括辐射强度,光学深度,相函数等等。近些年,由于计算机的快速发展,其运算速度及内存有了极大的提高,早期计算时间的限制已不再成为问题,为此辐射转移方程的直接法有必要再次被讨论,直接法的优势在于推导简单方便,程序也容易运行,所以它是很有效的一种数值计算方法。我们首先给出了一个具有4阶局部截断误差的高精度算法,应用此算法进行求解,并将数值结果与精确解及文献[23]的结果进行了对比,肯定了该算法的优势所在。其次,通过利用一般化的数值微分和数值积分格式来近似转移方程的导数和积分部分,然后,我们离散所得到的方程,将辐射转移方程转化成了具有一般性的一个线性方程组。随后我们给出了直接法的一般误差估计式和计算量分析,充分说明直接法目前仍然是求解辐射转移方程的很有效的数值方法。最后,我们还进一步给出辐射转移方程直接法的一个收敛准则,这个结论是很重要的,对今后人们判断新的直接法是否能够正确求解辐射转移方程具有指导作用。同时我们在前人求解的例子下,采用相同的参数设定,对不同数值格式加以计算和验证,将具体的数值计算结果与精确解进行了比较,结果充分说明该法则的可行性。