No abstract

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 between matter particles and radiation quanta couples the hydrodynamic equations to those of radiative transfer, see Munier & Weaver (1986a), Munier & Weaver (1986b). The numerical treatment has to account for their potential stiffness (e.g., in optically thick environments). In this article, we will present a new method to numerically integrate these equations in a stable way by using minimally implicit Runge-Kutta methods. With these methods, the inversion of the implicit operator can be done analytically. We also take into account the physical behavior of the evolved variables in the limit of the stiff regime. We will show the results of applying this method to the reactions between neutrinos and matter in core-collapse supernovae simulations.