Near optimal angular quadratures for polarised radiative transfer

Abstract

In three-dimensional (3D) radiative transfer (RT) problems, the tensor product quadratures are generally not optimal in terms of the number of discrete ray directions needed for a given accuracy of the angular integration of the radiation field. In this paper, we derive a new set of angular quadrature rules that are more suitable for solving 3D RT problems with the short- and long-characteristics formal solvers. These quadratures are more suitable than the currently used ones for the numerical calculation of the radiation field tensors that are relevant in the problem of the generation and transfer of polarised radiation without assuming local thermodynamical equilibrium (non-LTE). We show that our new quadratures can save up to about 30\,\% of computing time with respect to the Gaussian-trapezoidal product quadratures with the same accuracy.

arXiv comment7 pages, 3 figures, 2 tables; accepted for publication in A&A

References in corpus (4)

• A publicly available simulation of an enhanced network region of the Sun
• Distributing many points on spheres: minimal energy and designs
• PORTA: A three-dimensional multilevel radiative transfer code for modeling the intensity and polarization of spectral lines with massively parallel computers
• The Hanle and Zeeman polarization signals of the solar Ca II 8542 Å line

Cited by in corpus (11)

• Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics
• Flare Kernels May be Smaller than You Think: Modelling the Radiative Response of Chromospheric Plasma Adjacent to a Solar Flare
• Novel framework for the three-dimensional NLTE inverse problem
• Atomic line radiative transfer with MCFOST I. Code description and benchmarking
• Radiance Cascades: A Novel High-Resolution Formal Solution for Multidimensional Non-LTE Radiative Transfer
• Improved near optimal angular quadratures for polarised radiative transfer in 3D MHD models
• 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
• The uncertainty of magnetic fields in 3D non-local thermodynamic equilibrium inversions
• Numerical evaluation of orientation averages and its application to molecular physics
• Irregular grids for 3D NLTE radiative transfer in stellar atmospheres
• A single frequency approach to nonequilibrium modeling of the chromosphere