(3+1)D hydrodynamic simulation of relativistic heavy-ion collisions

We present music, an implementation of the Kurganov-Tadmor algorithm for relativistic 3+1 dimensional fluid dynamics in heavy-ion collision scenarios. This Riemann-solver-free, second-order, high-resolution scheme is characterized by a very small numerical viscosity and its ability to treat shocks and discontinuities very well. We also incorporate a sophisticated algorithm for the determination of the freeze-out surface using a three dimensional triangulation of the hypersurface. Implementing a recent lattice based equation of state, we compute pₓ-spectra and pseudorapidity distributions for Au+Au collisions at s=200 GeV and present results for the anisotropic flow coefficients v₂ and v₄ as a function of both pₓ and pseudorapidity. We were able to determine v₄ with high numerical precision, finding that it does not strongly depend on the choice of initial condition or equation of state.