Related papers: Formal solutions for polarized radiative transfer.…
When integrating the radiative transfer equation for polarized light, the necessity of high-order numerical methods is well known. In fact, well-performing high-order formal solvers enable higher accuracy and the use of coarser spatial…
The numerical computation of reliable and accurate Stokes profiles is of great relevance in solar physics. In the synthesis process, many actors play a relevant role: among them the formal solver, the discrete atmospheric model, and the…
The discussion regarding the numerical integration of the polarized radiative transfer equation is still open and the comparison between the different numerical schemes proposed by different authors in the past is not fully clear. Aiming at…
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…
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…
A time discretization method is called strongly stable, if the norm of its numerical solution is nonincreasing. It is known that, even for linear semi-negative problems, many explicit Runge--Kutta (RK) methods fail to preserve this…
We prove that Runge-Kutta (RK) methods for numerical integration of arbitrarily large systems of Ordinary Differential Equations are linearly stable. Standard stability arguments -- based on spectral analysis, resolvent condition or strong…
Stabilized Runge-Kutta methods are especially efficient for the numerical solution of large systems of stiff nonlinear differential equations because they are fully explicit. For semi-discrete parabolic problems, for instance, stabilized…
In this technical note a general procedure is described to construct internally consistent splitting methods for the numerical solution of differential equations, starting from matching pairs of explicit and diagonally implicit Runge-Kutta…
Stiff and chaotic differential equations are challenging for time-stepping numerical methods. For explicit methods, the required time step resolution significantly exceeds the resolution associated with the smoothness of the exact solution…
Radiative transfer is essential in astronomy, both for interpreting observations and simulating various astrophysical phenomena. However, self-consistent line radiative transfer is computationally expensive, especially in 3D. To reduce the…
We investigate the well-posedness of the radiative transfer equation with polarization and varying refractive index. The well-posedness analysis includes non-homogeneous boundary value problems on bounded spatial domains, which requires the…
Explicit stabilized methods are an efficient alternative to implicit schemes for the time integration of stiff systems of differential equations in large dimension. In this paper, we derive explicit stabilized integrators of orders one and…
This paper deals with stability of classical Runge-Kutta collocation methods. When such methods are embedded in linearly implicit methods as developed in [12] and used in [13] for the time integration of nonlinear evolution PDEs, the…
The reliability and precision of numerically solving stochastic non-Markovian equations by standard numerical codes, more specifically, with the fourth-order Runge-Kutta routine for solving differential equations, is gauged by comparing the…
Highly precise numerical solutions to the radiative transfer equation with polarization present a special challenge. Here, we establish a precise numerical solution to the radiative transfer equation with combined Rayleigh and isotropic…
The linearization principle states that the stability (or instability) of solutions to a suitable linearization of a nonlinear problem implies the stability (or instability) of solutions to the original nonlinear problem. In this work, we…
A new formalism is introduced for the transfer of polarized radiation. Stokes parameters are shown to be four-vectors in a Minkowski-like space and, most strikingly, the radiative transfer equation (RTE) turns out to be an infinitesimal…
The radiative transfer equation (RTE) for polarized light accepts a convenient exponential solution when the absorption matrix commutes with its integral. We characterize some of the matrix depth variations which are compatible with the…
An explicit stabilized additive Runge-Kutta scheme is proposed. The method is based on a splitting of the problem in severely stiff and mildly stiff subproblems, which are then independently solved using a Runge-Kutta-Chebyshev scheme. The…