A flux-split algorithm applied to relativistic flows

The equations of RFD can be written as a hyperbolic system of conservation laws by choosing an appropriate vector of unknowns, We give an explicit formulation of the full spectral decomposition of the Jacobian matrices associated with the fluxes in each spatial direction, which is the essential ingredient of the techniques we propose in this paper. These techniques are based on the recently derived flux formula of Marquina, a new way to compute the numerical flux at a cell interface which leads to a conservative, upwind numerical scheme. Using the spectral decompositions in a fundamental way, we construct high order versions of the basic First-order scheme described by R. Donat and A. Marquina in (J. Comput. Phys. 125, 42 (1996)) and test their performance in several standard simulations in one dimension. Two-dimensional simulations include a wind tunnel with a flat faced step and a supersonic jet stream, both of them in strongly ultrarelativistic regimes. (C) 1998 Academic Press.