Dans cette thèse nous nous intéressons a l'approximation numérique des équations incompressibles de Navier-Stokes évoluant dans un domaine en mouvement par la méthode des éléments spectraux et des intégrateurs en temps d'ordre élève. Dans une première phase, nous présentons la méthode des éléments spectraux et les outils de base pour effectuer des discrétisations spectrales du type Galerkin ou Galerkin avec intégration numérique (G-NI). Nous couvrons un large éventail de possibilités concernant les éléments de référence, fonctions de base, points d'interpolation et points de quadrature. Dans cette approche, l'intégration et la différentiation des fonctions polynomiales est faite numériquement grâce a l'aide d'ensembles de points convenables. En ce qui concerne la différenciation, nous présentons une étude numérique des points qui doivent être utilisés pour atteindre une meilleure stabilité numérique (parmi les choix que nous avons actuellement). Deuxièmement, nous introduisons les équations incompressibles stationnaires et non-stationnaires de Stokes et de Navier-Stokes et son approximation spectrale. Dans le cas non-stationnaire, nous introduisons une combinaison de la méthode Backward Différentiation Formula (BDF) et une formule d'extrapolation du même ordre pour l'intégration par rapport au temps. Une fois les équations discrétisées, un système linéaire doit être résolu pour obtenir la solution approchée. Dans ce contexte, nous resolvons ce système avec un préconditionneur par blocs. Nous montrons que le préconditionneur est optimal par rapport au nombre d'iterations utilisées par la méthode GMRES dans le cas stationnaire, mais pas dans le cas non-stationnaire. Une autre alternative est d'utiliser les méthodes de factorization algébrique de type Yosida et séparer le calcul de la vitesse et de la pression. Un cas test est présente pour déterminer les proprietes de convergence de ce type de méthodes dans notre contexte. Troisièmement, nous 'tendons les algorithmes développés dans le cas ou le domaine est fixe au cadre de la formulation Arbitraire Lagrange-Euler (ALE). La question de la définition d'une carte ALE d'ordre élevé est aborée. Cela permet de construire un domaine de calcul qui est d'ecrit avec des éléments courbes. Un cas test utilisant une méthode directe et les méthodes Yosida-q pour resoudre le système linéaire est présente pour montrer les ordres de convergence de la méthode proposée. Finalement, nous appliquons la méthode développée pour résoudre une un problème d'interaction fluide-structure pour un exemple simple bidimensionnel d'hémodynamique. Nous considérons deux approches: une implicite entièrement couplée et une semi-implicite.