Formal solutions for polarized radiative transfer. III. Stiffness and instability

Abstract

Efficient numerical approximation of the polarized radiative transfer equation is challenging because this system of ordinary differential equations exhibits stiff behavior, which potentially results in numerical instability. This negatively impacts the accuracy of formal solvers, and small step-sizes are often necessary to retrieve physical solutions. This work presents stability analyses of formal solvers for the radiative transfer equation of polarized light, identifies instability issues, and suggests practical remedies. In particular, the assumptions and the limitations of the stability analysis of Runge-Kutta methods play a crucial role. On this basis, a suitable and pragmatic formal solver is outlined and tested. An insightful comparison to the scalar radiative transfer equation is also presented.

References in corpus (4)

• DELO-Bezier formal solutions of the polarized radiative transfer equation
• PORTA: A three-dimensional multilevel radiative transfer code for modeling the intensity and polarization of spectral lines with massively parallel computers
• Formal Solutions for Polarized Radiative Transfer. I. The DELO Family
• Formal Solutions for Polarized Radiative Transfer. II. High-order Methods

Cited by in corpus (12)

• STiC -- A multi-atom non-LTE PRD inversion code for full-Stokes solar observations
• Modeling the scattering polarization of the solar Ca i 4227 Å line with angle-dependent partial frequency redistribution
• Formal Solutions for Polarized Radiative Transfer. IV. Numerical Performances in Practical Problems
• The transfer of polarized radiation in resonance lines with partial frequency redistribution, J-state interference, and arbitrary magnetic fields. A radiative transfer code and useful approximations
• Covariant polarized radiative transfer on cosmological scales for investigating large-scale magnetic field structures
• Numerical solutions to linear transfer problems of polarized radiation I. Algebraic formulation and stationary iterative methods
• Empirical determination of atomic line parameters of the 1.5 $\mu$m spectral region
• Discontinuities in numerical radiative transfer
• Assessment of the CRD approximation for the observer's frame RIII redistribution matrix
• Diagnostic potential of wavelength-integrated scattering polarisation signals of the solar He II Ly-alpha line
• A numerical approach for modelling the polarisation signals of strong resonance lines with partial frequency redistribution. Numerical applications to two-term atoms and plane-parallel atmospheres
• Fast and accurate approximation of the angle-averaged redistribution function for polarized radiation