We present a second-order upwind numerical scheme for equations of relativistic hydrodynamics with a source term. A new non-linear Riemann solver is constructed. Solution of a Riemann problem on a cells boundary is based on exact relations in case of zero tangential velocities. In this sense our solver is "exact". In case of non-zero tangential velocities a reasonable approximation is made that allows to avoid solution of very complicated exact equations. The scheme is tested on several one- and three-dimensional solutions demonstrating a good performance for both strong shocks and strong rarefactions.

Simple, self-similar, elliptic solutions of non-relativistic fireball hydrodynamics are presented, generalizing earlier results for spherically symmetric fireballs with Hubble flows and homogeneous temperature profiles. The transition from one dimensional to three dimensional expansions is investigated in an efficient manner.