[en] In a wall-confined environment, the plasma behaviour can be identified by dividing the domain in two main regions: the bulk, where the quasi-neutrality prevails (the charge density is close to zero) and the sheath region near the walls, where positive space charge is built. This phenomenon is due to a higher mobility of electrons, creating an electric potential structure that accelerates the ions towards the surface and repels the electrons. In this work we present the effort carried on in order to obtain high order simulations of a mixture of argon plasma: in particular we present a novel physical model (the multicomponent equations in contrast with the most common multifluid approach) and we show the first results of the application of the Discontinuous Galrkin framework, detailing challenges and solutions adopted. Simulations show encouraging results for future developments involving more complex geometries, more variegated mixtures and more accurate time schemes.