La modélisation numérique des décharges plasma joue un rôle important dans la compréhension des mécanismes physiques ou chimiques ayant lieu dans les dispositifs assistés par plasma. Une grande partie de ces mécanismes est déjà prise en compte dans les codes actuels. En revanche, beaucoup d'entre eux ne permettent pas de travailler avec des géométries complexes. Cette limitation provient essentiellement de l'utilisation de maillages structurés, cartésiens. Ceux-ci ne sont pas bien adaptés aux géométries courbes. Les calculs en maillages structurés deviennent rapidement compliqués et spécifiques à une géométrie donnée. Notre travail concerne la modélisation de décharge pour un réacteur de traitement à la pression atmosphérique développé par Dow Corning. Sa configuration complexe ainsi que ses grandes dimensions nous ont incités à faire un nouveau code fonctionnant en maillages non structurés. Celui-ci doit être capable de s'adapter à la présence d'une pointe, d'arrondis et de multiples diélectriques mais aussi permettre le passage rapide à de nouvelles géométries. De plus ses grandes dimensions nécessitent l'utilisation de maillages raffinés uniquement aux endroits nécessaires (pointe, surfaces des diélectriques...). Le modèle mathématique utilisé est basé sur l'équation de Poisson couplée aux équations de transport de type dérive-diffusion. Plusieurs discrétisations numériques ont été testées dans des configurations physiques différentes. Nous présentons et validons les méthodes numériques choisies. Les résultats obtenus pour le réacteur Dow Corning sont alors exposés et commentés.