Let be the Jacobian matrix associated with one of the fluxes of the original system, and the vector of unknowns. Then, the locally constant matrix , depending on and (the left and right state defining the local Riemann problem) must have the following four properties:
Once a matrix satisfying Roe's conditions has been obtained for every numerical interface, the numerical fluxes are computed by solving the locally linear system. Roe's numerical flux is then given by
where , , and are the eigenvalues and the right and left eigenvectors of , respectively (p runs from 1 to the number of equations of the system).
Roe's linearization for the relativistic system of equations in a general spacetime can be expressed in terms of the average state [49, 50]
where g is the determinant of the metric tensor . The role played by the density in case of the Cartesian non-relativistic Roe solver as a weight for averaging, is taken over in the relativistic variant by k, which apart from geometrical factors tends to in the non-relativistic limit. A Riemann solver for special relativistic flows and the generalization of Roe's solver to the Euler equations in arbitrary coordinate systems are easily deduced from Eulderink's work. The results obtained in 1D test problems for ultra-relativistic flows (up to Lorentz factors 625) in the presence of strong discontinuities and large gravitational background fields demonstrate the excellent performance of the Eulderink-Roe solver .
Relaxing condition 3 above, Roe's solver is no longer exact for shocks but still produces accurate solutions, and moreover, the remaining conditions are fulfilled by a large number of averages. The 1D general relativistic hydrodynamic code developed by Romero et al.  uses flux formula (26) with an arithmetic average of the primitive variables at both sides of the interface. It has successfully passed a long series of tests including the spherical version of the relativistic shock reflection (see Section 6.1).
Roe's original idea has been exploited in the so-called local characteristic approach (see, e.g., ). This approach relies on a local linearization of the system of equations by defining at each point a set of characteristic variables, which obey a system of uncoupled scalar equations. This approach has proven to be very successful, because it allows for the extension to systems of scalar nonlinear methods. Based on the local characteristic approach are the methods developed by Marquina et al.  and Dolezal & Wong , which both use high-order reconstructions of the numerical characteristic fluxes, namely PHM  and ENO  (see Section 9.4).
|Numerical Hydrodynamics in Special Relativity
Jose Maria Martí and Ewald Müller
© Max-Planck-Gesellschaft. ISSN 1433-8351
Problems/Comments to firstname.lastname@example.org