DeepAI AI Chat
Log In Sign Up

Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations

by   Manuel R. Izquierdo, et al.

Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of 3+1 equations closely resembling the hydrodynamic ones. Solving the M1 scheme is still mathematically challenging, since it is necessary to model the radiation-matter interaction in regimes where the evolution equations become stiff and behave as an advection-diffusion problem. Here, we present different global, high-order time integration schemes based on Implicit-Explicit Runge-Kutta (IMEX) methods designed to overcome the time-step restriction caused by such behavior while allowing us to use the explicit RK commonly employed for the MHD and Einstein equations. Finally, we analyze their performance in several numerical tests.


A Kernel Based Unconditionally Stable Scheme for Nonlinear Parabolic Partial Differential Equations

In this paper, a class of high order numerical schemes is proposed to so...

A skeletal high-order structure preserving scheme for advection-diffusion equations

We introduce a nonlinear structure preserving high-order scheme for anis...

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics se...

An implicit symplectic solver for high-precision long term integrations of the Solar System

Compared to other symplectic integrators (the Wisdom and Holman map and ...

Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models

The tracer equations are part of the primitive equations used in ocean m...

An Explicit Mapped Tent Pitching Scheme for Maxwell Equations

We present a new numerical method for solving time dependent Maxwell equ...