Exponential discontinuous solution of the multispecies relativistic heavy-ion transport equations

Linear schemes applied to charged-particle transport problems demonstrate high-order accuracy but under certain conditions can also produce negative solutions. However, the recently developed nonlinear exponential discontinuous (ED) method has been shown to produce accurate strictly positive solutions for positive sources in neutral particle transport applications. We have applied this method to the solution, in space and energy, of the multispecies transport equations for relativistic heavy ions. The solution includes projectile fragments resulting from collisions between projectile ions and atoms in the target medium. The primary ion and all fragments are treated using the straight-ahead approximation under which the fragments continue on with the same velocity as the original projectile