Study of higher harmonics based on (3+1)-d relativistic viscous hydrodynamics

We develop a state-of-the-art numerical scheme for solving the relativistic hydrodynamic equation, based on relativistic Godunov method. We find that our numerical algorithm has less artificial viscosity compared to other numerical methods (SHASTA, KT and NT) in a grid size, which is important for physical viscosity analyses in relativistic heavy ion collisions. We employ our new numerical scheme to investigation of higher harmonics which contains detail information of the quark-gluon plasma (QGP) property.