(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(T)-spectra and pseudorapidity distributions for Au+Au collisions at root s = 200 GeV and present results for the anisotropic flow coefficients nu(2) and nu(4) as a function of both p(T) and pseudorapidity eta. We were able to determine v4 with high numerical precision, finding that it does not strongly depend on the choice of initial condition or equation of state.