An in-house developed, 2D/3D unstructured CFD solver has been extended to deal with a mixture of thermally perfect gases in chemical non-equilibrium. The Euler equations have been coupled with a state-to-state kinetic model for argon plasma. The spatial discretization uses compact stencil Residual Distribution Schemes and shock waves can be modelled using either shock-capturing or shock-fitting. Promising results have been obtained using the shock-fitting approach for a 2D hypersonic flow past the fore-body of a circular cylinder
