Les instabilités fluides peuvent dégrader le confinement du plasma au sein des tokamaks. Étant données les échelles spatio-temporelles, on choisit les modèles fluides obtenus à partir de la dérivation des modèles cinétiques. On dérive plusieurs modèles hiérarchiques de la MagnétoHydroDynamique (MHD) et en particulier les modèles de la MHD réduite du Current Hole et de l'équilibre de Grad-Shafranov. Une des difficulté de l'ensemble de ces modèles est de respecter l'équation modélisant l'absence de monopôles magnétiques. Pour assurer cette condition en tout point du domaine, le champ magnétique est réécrit avec un potentiel vecteur. L'utilisation de potentiels fait apparaître des équations faisant intervenir des dérivées d'ordre supérieurs. La stratégie numérique développée est l'utilisation de la méthode des éléments finis avec des éléments C1 de Bell. Sur un maillage non structuré, ces éléments ont l'intérêt de présenter une base réduite définir exclusivement avec des variables aux noeuds du maillage. Les modèles de MHD réduite du Current Hole et de Grad-Shafranov ont été résolus avec ces éléments. La résolution du cas test de Grad-Shafranov avec les conditions de bords exactes a permis d'obtenir l'ordre optimale de 5. La résolution du système du Current Hole avec ces éléments, validée par l'obtention du paramètre η1/3, a permis l'observation de développement d'instabilités en dents de scies. / Fluid instabilities can degrade plasma confinement in tokamaks. Given the spatial and temporal scales, we choose the fluid models obtained from the derivation of kinetic models. We derived several hierarchical models of MagnetoHydroDynamic (MHD) and in particular models of reduced MHD like the Current Hole and the Grad-Shafranov equilibrium. One of the difficulty of all these models is to respect the absence of magnetic monopoles equation. To ensure this condition at any point, the magnetic field is rewritten with a vector potential. The use of vector portential implies that higher order derivatives appear in the equation. The numerical strategy is developed using the finite element method with C1 Bell's elements. On a unstructured mesh, these have the advantage to present a reduced basis with degrees of freedom defined exclusively on the nodes of the mesh. The reduced MHD models of the Current Hole and Grad-Shafranov have thus been resolved with these elements. The resolution of a Grad-Shafranov test case with exact boundary conditions yields the optimal order of 5. The resolution of the Current Hole system with thesse elements has been validated by obtaining physical parameter η1/3 and allowed the observation of the development of sawtooth instabilities.