A multi-species solver for plasma at thermodynamical equilibrium is developed. A numerical scheme, based on Roe’s approximate Riemann solver, is implemented with some modification regarding averaged quantities. Treatments for perfect and real gas have been carried out. The real gas approach considers changes in the composition of the gas in local thermodynamics equilibrium assumptions. These models have been applied to the classical shock tube problem involving two different gases and air as a multi-species gas.