Discontinuities in numerical radiative transfer
arXiv:1903.08891 · doi:10.1051/0004-6361/201833984
Abstract
Observations and magnetohydrodynamic simulations of solar and stellar atmospheres reveal an intermittent behavior or steep gradients in physical parameters, such as magnetic field, temperature, and bulk velocities. The numerical solution of the stationary radiative transfer equation is particularly challenging in such situations, because standard numerical methods may perform very inefficiently in the absence of local smoothness. However, a rigorous investigation of the numerical treatment of the radiative transfer equation in discontinuous media is still lacking. The aim of this work is to expose the limitations of standard convergence analyses for this problem and to identify the relevant issues. Moreover, specific numerical tests are performed. These show that discontinuities in the atmospheric physical parameters effectively induce first-order discontinuities in the radiative transfer equation, reducing the accuracy of the solution and thwarting high-order convergence. In addition, a survey of the existing numerical schemes for discontinuous ordinary differential systems and interpolation techniques for discontinuous discrete data is given, evaluating their applicability to the radiative transfer problem.
References in corpus (8)
- 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
- IRIS: A Generic Three-Dimensional Radiative Transfer Code
- Supersonic Downflows in a Sunspot Light Bridge
- Formal Solutions for Polarized Radiative Transfer. I. The DELO Family
- Formal solutions for polarized radiative transfer. III. Stiffness and instability
- Formal Solutions for Polarized Radiative Transfer. IV. Numerical Performances in Practical Problems
- Formal Solutions for Polarized Radiative Transfer. II. High-order Methods