In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and "Rankine-Hugoniot" jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges.

We present a ballistic description of the formation and propagation of the working surface of a relativistic jet. Using simple laws of conservation of mass and linear momentum at the working surface, we obtain a full description of the working surface flow parametrised by the initial velocity and mass injection rate. This allows us to compute analytically the energy release at any time in the working surface. We compare this model with the results obtained numerically through a new hydrodynamical code applied to the propagation of a relativistic fluid in one dimension in order to test the limits of our study. Finally, we compare our analytical results with observed light curves of five long gamma ray bursts and show that our model is in very good agreement with observations using simple periodic variations of the injected velocity profiles. This simple method allows us to recover initial mass discharge and energy output ejected during the burst.