偏振辐射线性转移问题的数值解 I. 代数表述与定常迭代法

背景。偏振辐射产生与转移的数值建模是太阳与恒星物理研究中的关键任务，并引出了一类可重构为线性系统的离散问题。为求解此类问题，通常依赖高效的定常迭代法。然而，这些方法的收敛性质依赖于具体问题，且当应用于偏振辐射转移问题时，对其收敛条件的严格研究仍然缺乏。目的。在总结偏振辐射数值转移中最广泛使用的迭代方法之后，本文旨在阐明这些方法的收敛如何依赖于不同设计要素，如形式求解器的选择、问题的离散化或阻尼因子的使用。主要目标是从稳定性与收敛速率方面凸显不同迭代方法的优缺点。方法。我们首先引入辐射转移问题的代数表述。该表述使我们能显式组装由不同定常迭代法产生的迭代矩阵，计算其谱半径并推导收敛速率，并检验不同离散化设置、问题参数与阻尼因子的影响。结论。本文所采用的基于偏振辐射线性转移问题完全代数表述的通用方法，为各类迭代格式的收敛速率提供了有用估计。此外，它可催生新颖的求解途径以及对更广泛设置（包括非偏振情形）的分析。