Spritz: general relativistic magnetohydrodynamics with neutrinos

Cipolletta, Federico; Kalinani, Jay Vijay; Giangrandi, Edoardo; Giacomazzo, B.; Ciolfi, R.; Sala, Lorenzo; Giudici, Beatrice

doi:10.1088/1361-6382/abebb7

Cited by 16 publications

(16 citation statements)

References 63 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Recently, we introduced a new publicly available GRMHD code called Spritz 1 [25], built to work within the Einstein Toolkit infrastructure [26][27][28]. We further improved and tested the code as described in [29]. It now includes support for high-order numerical schemes, use of finite temperature tabulated EOS, as well as a ray-by-ray neutrino leakage scheme that is based on the ZelmaniLeak public code and which accounts for both neutrino cooling and heating [30,31].…”

Section: And Refs Therein)mentioning

confidence: 99%

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

Kalinani,

Ciolfi,

Kastaun

et al. 2021

Preprint

View full text Add to dashboard Cite

General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) -black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on 'conserved' variables which then need to be converted back into the fundamental ('primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the evolution of a BH-accrection disk system, show that RePrimAnd is able to support highly magnetized, low density environments, even for magnetizations as high as 10 4 , for which the previously used recovery scheme fails.

show abstract

Section: And Refs Therein)mentioning

confidence: 99%

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

Kalinani,

Ciolfi,

Kastaun

et al. 2021

Preprint

View full text Add to dashboard Cite

show abstract

“…At present, the Riemann solver and constrained transport scheme implemented in existing numerical relativity magnetohydrodynamics codes, e.g. [55,[60][61][62][63][64][65][66][67], are based on the HLLE solver [4,59,68]. (An exception is the SpECTRE [69] code, which is based on the discontinuous Galerkin method.)…”

Section: Introductionmentioning

confidence: 99%

Implementation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its application to a binary neutron star merger

Kiuchi¹,

Held²,

Sekiguchi³

et al. 2022

Preprint

View full text Add to dashboard Cite

We implement advanced Riemann solvers HLLC and HLLD [1,2] together with an advanced constrained transport scheme [3] in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We validate our implementation by performing a series of one-and multi-dimensional test problems for relativistic hydrodynamics and magnetohydrodynamics in both Minkowski spacetime and a static black hole spacetime. We find that the numerical solutions with the advanced Riemann solvers are more accurate than those with the HLLE solver [4], which was originally implemented in our code. As an application to numerical relativity, we simulate an asymmetric binary neutron star merger leading to a short-lived massive neutron star both with and without magnetic fields. We find that the lifetime of the rotating massive neutron star formed after the merger and also the amount of the tidally-driven dynamical ejecta are overestimated when we employ the diffusive HLLE solver. We also find that the magnetorotational instability is less resolved when we employ the HLLE solver because of the solver's large numerical diffusivity. This causes a spurious enhancement both of magnetic winding resulting from large scale poloidal magnetic fields, and also of the energy of the outflow induced by magnetic pressure.

show abstract

“…Of these codes, only a few can simulate the merger of magnetized neutron stars with neutrinos and a realistic EoS. Those that can generally use a simplified approximate neutrino scheme called leakage [16][17][18][19] (but more recently also with the M1 formalism [20], although with a simplified temperaturedependent EoS). Here, we extend our code, Mhduet , to allow for tabulated EoS with a leakage scheme to model the neutrinos.…”

Section: Introductionmentioning

confidence: 99%

Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos

Palenzuela,

Liebling,

Miñano

2022

Preprint

View full text Add to dashboard Cite

Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is increasingly important in light of advancements in multi-messenger astronomy generally, and in gravitational wave astronomy in particular (such as the development of thirdgeneration detectors). Several general relativistic codes have been incorporating some of these elements with different levels of realism. Here, we extend our code Mhduet , which can perform large eddy simulations of magnetohydrodynamics to help capture the magnetic field amplification during the merger, and to allow for realistic equations of state and neutrino cooling via a leakage scheme. We perform several tests involving isolated and binary neutron stars demonstrating the accuracy of the code.

show abstract

Spritz: general relativistic magnetohydrodynamics with neutrinos

Cited by 16 publications

References 63 publications

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

Implementation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its application to a binary neutron star merger

Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos

Contact Info

Product

Resources

About