A new general relativistic magnetohydrodynamics code for dynamical spacetimes

Abstract

We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code has been developed with the main objective of studying astrophysical scenarios in which both, high magnetic fields and strong gravitational fields appear, such as the magneto-rotational collapse of stellar cores, the collapsar model of GRBs, and the evolution of neutron stars. The code is based on an existing and thoroughly tested purely hydrodynamics code and on its extension to accommodate weakly magnetized fluids (passive magnetic field approximation). The numerical code we present here is based on high-resolution shock-capturing schemes to solve the GRMHD equations together with the flux constraint transport method to ensure the solenoidal condition of the magnetic field. Since the astrophysical applications envisaged do not deviate much from spherical symmetry, the conformal flatness condition approximation is used for the formulation of the Einstein equations. In addition, the code can handle several equations of state, from simple analytical expressions to microphysical tabulated ones. In this paper we present stringent tests of our new GRMHD numerical code, which show its ability to handle all aspects appearing in the astrophysical scenarios for which the code is intended, namely relativistic shocks, highly magnetized fluids, and equilibrium configurations of magnetized neutron stars. As an application, magneto-rotational core collapse simulations of a realistic progenitor are presented, comparing the results with our previous finding in the passive magnetic field approximation.

18 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

References in corpus (17)

• Producing ultra-strong magnetic fields in neutron star mergers
• Merger of binary neutron stars to a black hole: Disk mass, short gamma-ray bursts, and quasinormal mode ringing
• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• Relativistic neutron star merger simulations with non-zero temperature equations of state I. Variation of binary parameters and equation of state
• WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics
• 3D Collapse of Rotating Stellar Iron Cores in General Relativity including Deleptonization and a Nuclear Equation of State
• Magnetorotational collapse of massive stellar cores to neutron stars: Simulations in full general relativity
• WHAM: A WENO-based general relativistic numerical scheme I: Hydrodynamics
• Generic Gravitational Wave Signals from the Collapse of Rotating Stellar Cores
• Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
• Waveless Approximation Theories of Gravity
• Axisymmetric simulations of magnetorotational core collapse: Approximate inclusion of general relativistic effects
• Rotating Collapse of Stellar Iron Cores in General Relativity
• Magnetohydrodynamics of Neutrino-Cooled Accretion Tori around a Rotating Black Hole in General Relativity
• Dynamics of magnetized relativistic tori oscillating around black holes
• General relativistic simulations of pasive-magneto-rotational core collapse with microphysics
• Towards relativistic simulations of magneto-rotational core collapse