A Fourth-Order Finite Volume Scheme for Resistive Relativistic Magnetohydrodynamics
Authors:
Andrea Mignone,
Vittoria Berta,
Marco Rossazza,
Matteo Bugli,
Giancarlo Mattia,
Luca Del Zanna,
Lorenzo Pareschi
Abstract:
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advance…
▽ More
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advanced as cell volume-averages. Spatial accuracy employs 5th-order accurate WENO-Z reconstruction from point values (as described in a companion paper) to obtain left and right states at zone interfaces. Explicit flux evaluation is carried out by solving a Riemann problem at cell interfaces, using the Maxwell-Harten-Lax-van Leer with contact wave resolution (MHLLC). Time stepping is based on the implicit-explicit (IMEX) Runge-Kutta (RK) methods, of which we consider both the 3rd-order strong stability preserving SSP3(4,3,3) and a recent 4th-order additive RK scheme, to cope with the stiffness introduced by the source term in Ampere's law. Numerical benchmarks are presented in order to assess the accuracy and robustness of our implementation.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
Authors:
Manuel R. Izquierdo,
Lorenzo Pareschi,
Borja Miñano,
Joan Massó,
Carlos Palenzuela
Abstract:
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…
▽ More
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.
△ Less
Submitted 29 May, 2023; v1 submitted 31 October, 2022;
originally announced November 2022.