Hosted by: Chun Shen

Relativistic hydrodynamic simulations play a key role in exploring the QGP bulk property and the QCD phase transition from analyses of high-energy heavy-ion collisions at RHIC and LHC. From the intensive study based on relativistic viscous hydrodynamic models with event-by-event initial fluctuations, we can extract detailed information of the bulk feature of the QGP such as transport coefficients and the QCD equations of states. In the quantitative analyses of the QGP property, high-precision numerical treatment on the hydrodynamic calculation is important. Recently, we developed a new 3+1 dimensional relativistic viscous hydrodynamics code in Cartesian coordinates. In the algorithm, we use a Riemann solver based on the two-shock approximation which is stable under existence of large shock waves. We extend the algorithm in Cartesian coordinates to that in Milne coordinates so that we can efficiently apply it to the analyses of relativistic heavy-ion collisions. We check the correctness of the numerical algorithm by comparing numerical calculations and analytical solutions in various problems for ideal and viscous fluids. The new numerical scheme is stable even with small numerical viscosity, which is very important to discuss the physical viscosities at RHIC and LHC.