Fast Iterative Techniques for Polarized Radiative Transfer in Spherically Symmetric Moving Media

• Published in: The Astrophysical journal Vol. 903; no. 1; pp. 6 - 15
• Main Authors: Megha, A., Sampoorna, M., Nagendra, K. N., Anusha, L. S., Sankarasubramanian, K.
• Format: Journal Article
• Language: English
• Published: Philadelphia The American Astronomical Society 01.11.2020; IOP Publishing
• Subjects: Astronomical models; Astrophysics; Atmospheric models; Computational methods; Convergence; Iterative methods; Line spectra; Perturbation; Radiative transfer; Radiative transfer equation; Radiative transfer simulations; Solar atmosphere; Spectropolarimetry; Stellar atmospheres; Velocity distribution
• ISSN: 0004-637X; 1538-4357
• DOI: 10.3847/1538-4357/abb6f4

For a more precise modeling of polarized spectral lines formed in extended and expanding stellar atmospheres, the solution of the radiative transfer equation for the Stokes vectors must be obtained in a spherical geometry rather than in a planar geometry. In this paper, we present the modern iterative techniques based on operator perturbation to solve the spherically symmetric polarized radiative transfer equation with velocity fields. We consider scattering on a two-level atom and account for partial frequency redistribution. An accurate numerical solution to such problems requires the use of spatial grids with higher resolution. Consequently, Jacobi-based methods lead to slower convergence rate. The convergence rate can be improved by a factor of 2 or more when fast iterative schemes based on Gauss-Seidel (GS) and successive overrelaxation (SOR) methods are used over the Jacobi-based method. Here we present the Jacobi, GS, and SOR iterative techniques for solving the abovementioned problem, and discuss their convergence behavior.