Abstract FR:

On s'interesse dans cette etude a l'analyse numerique des ecoulements instationnaires bidimentionnels de fluides visqueux incompressibles autour d'un obstacle pour les grands nombres de reynolds. On est amene a resoudre numeriquement les equations de navier-stokes incompressibles en formulation vorticite-fonction de courant. Il faut donc determiner une methode numerique qui soit robuste, rapide et parallele. L'etape de convection est resolue par la methode des caracteristiques. On aboutit ensuite a un probleme de stokes generalise (probleme elliptique d'ordre 4) par pas de temps. La determination de la vorticite a la paroi de l'obstacle par une formule asymptotique issue de la methode integrale permet de transformer le probleme de stokes generalise en 4 problemes elliptiques d'ordre 2. Ces operateurs elliptiques sont ensuite discretises par la methode des elements avec joints : c'est une methode de decomposition de domaines sans recouvrement. Et a l'interieur de chaque sous-domaine, la methode des elements finis q1 est appliquee. Le systeme lineaire ainsi obtenu est resolu par la methode du complement de schur. On obtient alors un probleme pose aux interfaces dont la resolution est faite par la methode du gradient conjugue preconditionne, et un probleme interne aux sous-domaines qui est resolu par cholesky. On a ensuite confronte nos resultats avec d'autres resultats numeriques et experimentaux pour des ecoulements autour d'un cylindre et d'un profil a des nombres de reynolds superieurs a 1000. Les mesures de temps sont aussi donnees afin d'estimer les performances de notre solveur parallele.