Étude de conditions limites de sortie pour la méthode des éléments finis

La résolution numérique des équations de Navier-Stokes nécessite l'utilisation d'un domaine de calcul discrétisé. Quelque soit la méthode de détermination des champs de vitesse et de pression du fluide, l'augmentation du nombre de Reynolds implique un agrandissement indubitable de la taille du domaine. Ce qui revient à accroître la durée de la simulation ainsi que la demande en ressources informatiques. Au contraire, dans le cas où la taille du domaine de calcul est petite, les frontières influenceront l'écoulement interne et pollueront le signal et les données de la simulation. Nous présentons une famille de conditions limites diffusives de sortie utile lors de simulations numériques d'écoulements fluides incompressibles sur des domaines sévèrement tronqués. La forme générale de ces conditions limites assure la stabilité des calculs et la précision des résultats, même lorsque l'écoulement provoque des tourbillons importants et des reflux de fluide à la frontière de sortie. Le code numérique employé pour traiter ce genre de cas est basé sur un schéma implicite fortement couplé et sur la méthode des éléments spectraux. La discrétisation temporelle est bâtie sur une heuristique adaptative en pas de temps et en ordre. Ce logiciel est interne à l'École Polytechnique de Montréal et en continuel développement. Le projet de ce mémoire est de mettre en oeuvre la famille de conditions limites et d'analyser les résultats de simulations à haut nombre de Reynolds d'un écoulement autour d'un cylindre carré en repère absolu et en repère relatif en translation. Dans un premier temps, nous vérifions et validons le code complet pour ensuite optimiser le choix des valeurs des paramètres de configuration. Puis, nous introduisons la théorie énergétique derrière les conditions limites diffusives et les implémentons numériquement. Nous effectuons alors de nombreux tests en repère absolu à divers nombres de Reynolds et sur différentes longueurs de domaine afin de comparer les résultats avec des valeurs provenant de la littérature. Finalement, nous simulons le phénomène des vibrations induites par vortex (VIV) autour d'un cylindre carré à Re = 500 sur un domaine sévèrement tronqué. L'analyse poussée de l'influence des conditions limites de sortie sur un tel phénomène amène à la conclusion suivante. Quelque soit la condition limite choisie, si la longueur du domaine est trop faible, alors les mouvements longitudinaux du solide seront forcément limités. En revanche, une légère perte en précision est acceptable face aux gains en temps et en ressources obtenus par la réduction de la taille du domaine de calcul et par l'imposition de conditions limites de sortie.