Related papers: Time integration for neutrino radiation transport …
The ability to accurately approximate trajectories of dynamical systems enables their analysis, prediction, and control. Neural network (NN)-based approximations have attracted significant interest due to fast evaluation with good accuracy…
We present a C++ implementation of a fifth order semi-implicit Runge-Kutta algorithm for solving Ordinary Differential Equations. This algorithm can be used for studying many different problems and in particular it can be applied for…
The development of a neutrino moment based radiative-transfer code to simulate binary neutron-star mergers can easily become an obstacle path because of the numerous ways in which the solution of the equations may fail. We describe the…
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…
We present the implementation of an implicit-explicit (IMEX) Runge-Kutta numerical scheme for general relativistic hydrodynamics coupled to an optically thick radiation field in two existing GR-hydrodynamics codes. We argue that the…
Space discretization of some time-dependent partial differential equations gives rise to systems of ordinary differential equations in additive form whose terms have different stiffness properties. In these cases, implicit methods should be…
We introduce \texttt{Weakhub}, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino-matter, neutrino-photon, and neutrino-neutrino…
To achieve efficient and accurate long-time integration, we propose a fast, accurate, and stable high-order numerical method for solving fractional-in-space reaction-diffusion equations. The proposed method is explicit in nature and…
Simulations in stellar astrophysics involve the coupling of hydrodynamics and nuclear reactions under a wide variety of conditions, from simmering convective flows to explosive nucleosynthesis. Numerical techniques such as operator…
Time integration of Fourier pseudo-spectral DNS is usually performed using the classical fourth-order accurate Runge--Kutta method, or other methods of second or third order, with a fixed step size. We investigate the use of higher-order…
We combine the recent relaxation approach with multiderivative Runge-Kutta methods to preserve conservation or dissipation of entropy functionals for ordinary and partial differential equations. Relaxation methods are minor modifications of…
We present results of explicit asymptotic approximations applied to neutrino--electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary…
The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In…
Some features of the physics of radiation-dominated shock waves are discussed with emphasis on the peculiarities which are important for correct numerical modeling of shock breakouts in supernova. With account of those peculiarities, a…
Many natural processes, such as chemical reactions and wave dynamics, are modeled as production-destruction (PD) systems that obey positivity and linear conservation laws. Classical time integrators do not guarantee positivity and can…
Finite element discretization of time dependent problems also require effective time-stepping schemes. While implicit Runge-Kutta methods provide favorable accuracy and stability problems, they give rise to large and complicated systems of…
A standard approach to solve ordinary differential equations, when they describe dynamical systems, is to adopt a Runge-Kutta or related scheme. Such schemes, however, are not applicable to the large class of equations which do not…
The unified gas-kinetic wave-particle method (UGKWP) has been developed for the multiscale gas, plasma, and multiphase flow transport processes for the past years. In this work, we propose an implicit unified gas-kinetic wave-particle…
We propose a family of integrators, Flow-Composed Implicit Runge-Kutta (FCIRK) methods, for perturbations of nonlinear ordinary differential equations, consisting of the composition of flows of the unperturbed part alternated with one step…
Dynamic systems have a fundamental relevance in the description of physical phenomena. The search for more accurate and faster numerical integration methods for the resolution of such systems is, therefore, an important topic of research.…