Related papers: Formal solutions for polarized radiative transfer.…
Since the early 1970s, inversion techniques have become the most useful tool for inferring the magnetic, dynamic, and thermodynamic properties of the solar atmosphere. The intrinsic model dependence makes it necessary to formulate specific…
Context. Numerical solutions to transfer problems of polarized radiation in solar and stellar atmospheres commonly rely on stationary iterative methods, which often perform poorly when applied to large problems. In recent times, stationary…
In this work, we present approaches to rigorously certify $A$- and $A(\alpha)$-stability in Runge-Kutta methods through the solution of convex feasibility problems defined by linear matrix inequalities. We adopt two approaches. The first is…
In this paper, we study the inverse problem for a class of abstract ultraparabolic equations which is well-known to be ill-posed. We employ some elementary results of semi-group theory to present the formula of solution, then show the…
This paper discusses stochastic numerical methods of Runge-Kutta type with weak and strong convergences for systems of stochastic differential equations in It\^o form. At the beginning we give a brief overview of the stochastic numerical…
Fractional-step methods are a popular and powerful divide-and-conquer approach for the numerical solution of differential equations. When the integrators of the fractional steps are Runge--Kutta methods, such methods can be written as…
Multiphysics systems are driven by multiple processes acting simultaneously, and their simulation leads to partitioned systems of differential equations. This paper studies the solution of partitioned systems of differential equations using…
The evolution of many astrophysical systems is dominated by the interaction between matter and radiation such as photons or neutrinos. The dynamics can be described by the evolution equations of radiation hydrodynamics in which reactions…
We present two new accurate and efficient method to compute the formal solution of the polarized radiative transfer equation. In this work, the source function and the absorption matrix are approximated using quadratic and cubic Bezier…
We consider the inverse problem of reconstructing the optical parameters of the radiative transfer equation (RTE) from boundary measurements in the diffusion limit. In the diffusive regime (the Knudsen number $\mathsf{Kn}\ll 1$), the…
Many important initial value problems have the property that energy is non-increasing in time. Energy stable methods, also referred to as strongly stable methods, guarantee the same property discretely. We investigate requirements for…
We discuss the problem of polarized radiative transfer in general relativity. We present a set of equations suitable for solving the problem numerically for the case of an arbitrary space-time metric, and show numerical solutions to example…
Applied to the master equation, the usual numerical integration methods, such as Runge-Kutta, become inefficient when the rates associated with various transitions differ by several orders of magnitude. We introduce an integration scheme…
Stability of a set of travelling wave solutions to the hyperbolic generalization of the convection-reaction-diffusion equation is studied by means of the qualitative methods and numerical simulation.
We analyze the spectral stability of small-amplitude, periodic, traveling-wave solutions of the Kawahara equation. These solutions exhibit high-frequency instabilities when subject to bounded perturbations on the whole real line. We…
Stability is a key property of both forward models and inverse problems, and depends on the norms considered in the relevant function spaces. For instance, stability estimates for hyperbolic partial differential equations are often based on…
We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined…
The parametric instability arising when ordinary differential equations (ODEs) are numerically integrated with Runge-Kutta-Nystr\"om (RKN) methods with varying step sizes is investigated. It is shown that when linear constant coefficient…
Stabilized methods (also called Chebyshev methods) are explicit methods with extended stability domains along the negative real axis. These methods are intended for large mildly stiff problems, originating mainly from parabolic PDEs. In…
Explicit Runge-Kutta methods are classical and widespread techniques in the numerical solution of ordinary differential equations (ODEs). Considering partial differential equations, spatial semidiscretisations can be used to obtain systems…