偏振辐射转移的形式解 III：刚性与不稳定性
太阳与恒星天体物理
2018-05-02 v1
摘要
偏振辐射转移方程的高效数值逼近具有挑战性，因为该常微分方程组表现出刚性，可能导致数值不稳定性。这会对形式解算子的准确性产生负面影响，通常需要很小的步长才能获得物理解。本文对偏振光辐射转移方程的形式解算子进行了稳定性分析，识别了不稳定性问题，并提出了实用的补救措施。特别是，Runge-Kutta 方法稳定性分析的假设与限制发挥了关键作用。在此基础上，概述并测试了一种合适且实用的形式解算子。此外，还展示了与标量辐射转移方程的深刻比较。
引用
@article{arxiv.1803.10945,
title = {Formal solutions for polarized radiative transfer. III. Stiffness and instability},
author = {Gioele Janett and Alberto Paganini},
journal= {arXiv preprint arXiv:1803.10945},
year = {2018}
}