Un nouvel point de vue est proposé pour la simulation des plasmas utilisant le modèle cinétique qui couple les équations de Vlasov pour la distribution des particules et les équations de Maxwell pour la contribution des champs électromagnétique. On part du principe que les équations de la Physique sont des objets mathématiques qui mettent en relation des objets géométriques. Afin de conserver les propriètés géométriques des différents objets intervenant dans une équation, on utilise, pour l'étude théorique et numérique, la géométrie différentielle. Il s'avère que toutes les équations de la Physique peuvent s'écrire à l'aide des formes différentielles et que sous ce point de vue celles-ci sont indépendantes du choix des coordonnées. On propose alors une discrétisation des formes différentielles en utilisant les B-splines comme fonctions d'interpolation. Afin d'être cohérent avec la théorie, on proposera également une discrétisation des différentes opérations de la géométrie différentielle agissant sur les formes différentielles. On teste notre schéma tout d'abord sur les équations de Maxwell avec plusieurs conditions aux bords et puisque ce schéma numérique obtenu est indépendant du système de coordonnées, on le teste également lorsque l'on effectue un changement de coordonnées. Enfin, on applique la même méthode sur les équations de Vlasov-Poisson 1D et on propose plusieurs schémas numériques.