Luis Antón

Numerical 3+1 General Relativistic Magnetohydrodynamics: A Local Characteristic Approach

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e., Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic MHD equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic MHD in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick disks accreting onto a black hole and subject to the magnetorotational instability.

RELATIVISTIC MAGNETOHYDRODYNAMICS: RENORMALIZED EIGENVECTORS AND FULL WAVE DECOMPOSITION RIEMANN SOLVER

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wave front in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However, its relative efficiency increases in multidimensional simulations.

A new general relativistic magnetohydrodynamics code for dynamical spacetimes

We present a new numerical code that solves the general relativistic magneto-hydrodynamical (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within a 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 hydrodynamical code and on its extension to accommodate weakly magnetized fluids (passive magnetic-field approximation). These codes have been applied in the past to simulate the aforementioned scenarios with increasing levels of sophistication in the input physics. The numerical code we present here is based on high-resolution shockcapturing schemes to solve the GRMHD equations, which are cast in first-order, flux-conservative hyperbolic form, together with the flux constraint transport method to ensure the solenoidal condition of the magnetic field. Since the astrophysical applications envisaged do not deviate significantly from spherical symmetry, the conformal flatness condition approximation is used for the formulation of the Einstein equations; this has repeatedly shown to yield very good agreement with full general relativistic simulations of corecollapse supernovae and the evolution of isolated neutron stars. 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, magnetorotational core-collapse simulations of a realistic progenitor are presented and the results compared with our previous findings in the passive magnetic-field approximation.

Cosmic Microwave Background Maps Lensed by Cosmological Structures: Simulations and Statistical Analysis

A method for ray-tracing through n-body simulations has been recently proposed. It is based on a periodic universe covered by simulation boxes. Photons move along appropriate directions to avoid periodicity effects. Here an improved version of this method is applied to simulate lensed CMB maps and maps of lens deformations. Particle mesh n-body simulations with appropriate boxes and resolutions are used to evolve the nonlinear inhomogeneities until the present time. The resulting maps are statistically analyzed to look for deviations from Gaussianity. These deviations are measured, for the first time, using correlations for configurations of n directions (n ≤ 6). A wide range of angular scales are considered. Some interesting prospects are outlined.

Upwind Relativistic MHD Code for Astrophysical Applications

We describe the status of devolpment of a 2.5D numerical code to solve the equations of ideal relativistic magnetohydrodynamics. The numerical code, based on high-resolution shock-capturing techniques, solves the equations written in conservation form and computes the numerical fluxes using a linearized Riemann solver. A special procedure is used to force the conservation of magnetic flux along the evolution.

A Divergence-Free High-Resolution Code for MHD

We describe a 2.5D numerical code to solve the equations of ideal magnetohydrodynamics (MHD). The numerical code, based on high-resolution shock-capturing (HRSC) techniques, solves the equations written in conservation form and computes the numerical fluxes using a linearized Riemann solver. A special procedure is used to force the conservation of magnetic flux along the time.