BHAC-QGP: three-dimensional MHD simulations of relativistic heavy-ion collisions, I. Methods and tests

We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions in the presence of electromagnetic fields. It is derived from the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes in a general-relativistic magnetohydrodynamical description. As the original Black Hole Accretion Code, BHAC-QGP benefits from the use of Adaptive Mesh Refinement (AMR), which allows us to dynamically adjust the resolution where necessary, and makes use of time-dependent Milne coordinates and the ultrarelativistic equation of state, $P = e/3$. We demonstrate that BHAC-QGP accurately passes a number of systematic and rigorous tests.