Context. The numerical modeling of the generation and transfer of polarized radiation is a key task in solar and stellar physics research and has led to a relevant class of discrete problems that can be reframed as linear systems. In order to solve such problems, it is common to rely on efficient stationary iterative methods. However, the convergence properties of these methods are problemdependent, and a rigorous investigation of their convergence conditions, when applied to transfer problems of polarized radiation, is still lacking. Aims. After summarizing the most widely employed iterative methods used in the numerical transfer of polarized radiation, this article aims to clarify how the convergence of these methods depends on different design elements, such as the choice of the formal solver, the discretization of the problem, or the use of damping factors. The main goal is to highlight advantages and disadvantages of the different iterative methods in terms of stability and rate of convergence. Methods. We first introduce an algebraic formulation of the radiative transfer problem. This formulation allows us to explicitly assemble the iteration matrices arising from different stationary iterative methods, compute their spectral radii and derive their convergence rates, and test the impact of different discretization settings, problem parameters, and damping factors. Results. Numerical analysis shows that the choice of the formal solver significantly affects, and can even prevent, the convergence of an iterative method. Moreover, the use of a suitable damping factor can both enforce stability and increase the convergence rate. Conclusions. The general methodology used in this article, based on a fully algebraic formulation of linear transfer problems of polarized radiation, provides useful estimates of the convergence rates of various iterative schemes. Additionally, it can lead to novel solution approaches as well as analyses for a wider range of settings, including the unpolarized case.
Context. The polarization signals produced by the scattering of anistropic radiation in strong resonance lines encode important information about the elusive magnetic fields in the outer layers of the solar atmosphere. An accurate modeling of these signals is a very challenging problem from the computational point of view, in particular when partial frequency redistribution (PRD) effects in scattering processes are accounted for with a general angle-dependent treatment. Aims. We aim at solving the radiative transfer problem for polarized radiation in nonlocal thermodynamic equilibrium conditions, taking angle-dependent PRD effects into account. The problem is formulated for a two-level atomic model in the presence of arbitrary magnetic and bulk velocity fields. The polarization produced by scattering processes and the Zeeman effect is considered. Methods. The proposed solution strategy is based on an algebraic formulation of the problem and relies on a convenient physical assumption, which allows its linearization. We applied a nested matrix-free GMRES iterative method. Effective preconditioning is obtained in a multifidelity framework by considering the light-weight description of scattering processes in the limit of complete frequency redistribution (CRD). Results. Numerical experiments for a one-dimensional (1D) atmospheric model show near optimal strong and weak scaling of the proposed CRD-preconditioned GMRES method, which converges in few iterations, independently of the discretization parameters. A suitable parallelization strategy and high-performance computing tools lead to competitive run times, providing accurate solutions in a few minutes. Conclusions. The proposed solution strategy allows the fast systematic modeling of the scattering polarization signals of strong resonance lines, taking angle-dependent PRD effects into account together with the impact of arbitrary magnetic and bulk velocity fields. Almost optimal strong and weak scaling results suggest that this strategy is applicable to realistic 3D models. Moreover, the proposed strategy is general, and applications to more complex atomic models are possible.
We consider the space-time discretization of the diffusion equation, using an isogeometric analysis (IgA) approximation in space and a discontinuous Galerkin (DG) approximation in time. Drawing inspiration from a former spectral analysis, we propose for the resulting space-time linear system a multigrid preconditioned GMRES method, which combines a preconditioned GMRES with a standard multigrid acting only in space. The performance of the proposed solver is illustrated through numerical experiments, which show its competitiveness in terms of iteration count, run-time and parallel scaling.
Context. The numerical modeling of the generation and transfer of polarized radiation is a key task in solar and stellar physics research and has led to a relevant class of discrete problems that can be reframed as linear systems. In order to solve such problems, it is common to rely on efficient stationary iterative methods. However, the convergence properties of these methods are problemdependent, and a rigorous investigation of their convergence conditions, when applied to transfer problems of polarized radiation, is still lacking. Aims. After summarizing the most widely employed iterative methods used in the numerical transfer of polarized radiation, this article aims to clarify how the convergence of these methods depends on different design elements, such as the choice of the formal solver, the discretization of the problem, or the use of damping factors. The main goal is to highlight advantages and disadvantages of the different iterative methods in terms of stability and rate of convergence. Methods. We first introduce an algebraic formulation of the radiative transfer problem. This formulation allows us to explicitly assemble the iteration matrices arising from different stationary iterative methods, compute their spectral radii and derive their convergence rates, and test the impact of different discretization settings, problem parameters, and damping factors. Results. Numerical analysis shows that the choice of the formal solver significantly affects, and can even prevent, the convergence of an iterative method. Moreover, the use of a suitable damping factor can both enforce stability and increase the convergence rate. Conclusions. The general methodology used in this article, based on a fully algebraic formulation of linear transfer problems of polarized radiation, provides useful estimates of the convergence rates of various iterative schemes. Additionally, it can lead to novel solution approaches as well as analyses for a wider range of settings, including the unpolarized case.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.