HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
نویسندگان
چکیده
We describe a conservative, shock-capturing scheme for evolving the equations of general relativistic magnetohydrodynamics. The fluxes are calculated using the Harten, Lax, and van Leer scheme. A variant of constrained transport, proposed earlier by Tóth, is used to maintain a divergence free magnetic field. Only the covariant form of the metric in a coordinate basis is required to specify the geometry. We describe code performance on a full suite of test problems in both special and general relativity. On smooth flows we show that it converges at second order. We conclude by showing some results from the evolution of a magnetized torus near a rotating black hole. Subject headings: accretion, accretion disks, black hole physics, Magnetohydrodynamics: MHD, Methods: Numerical
منابع مشابه
An efficient shock-capturing central-type scheme for multidimensional relativistic flows II. Magnetohydrodynamics
A third order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented and validated against several numerical tests. The simple and efficient central scheme described in Paper I (Del Zanna & Bucciantini 2002) for relativistic hydrodynamics is here extended to the magnetic case by following the strategies prescribed for classical M...
متن کاملAn HLLC Solver for Relativistic Flows – II. Magnetohydrodynamics
An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The HLLC solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper (Mignone & Bodo 2004) to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks ...
متن کاملPrimitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between “conserved” variables such as momentum and energy density and “primitive” variables such as restmass density, internal energy, and components of the four-velocity. The forward transformation (primitive to conserved) has a closed-form solution, but the inverse transforma...
متن کامل0 70 11 09 v 2 5 J un 2 00 7 WhiskyMHD : a new numerical code for general relativistic magnetohydrodynamics
The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of γ-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativi...
متن کاملar X iv : g r - qc / 0 70 11 09 v 1 1 9 Ja n 20 07 WhiskyMHD : a new numerical code for general relativistic magnetohydrodynamics
The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of γ-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativi...
متن کامل