An Immersed Boundary Method Approach Using Hierarchical and Overlapping Grids for Unsteady Aerodynamics

La simulation numérique d'écoulements compressibles autours d'objets mobiles comprend de nombreuses applications, comme la séparation de magasins à grande vitesse, le suivi de débris, les jets à un nombre de Mach élevé, et l'écoulement à l'intérieur des disjoncteurs haute-tension. La modélisation numérique de l'ensemble de ces applications se heurte principalement à certains défis : les complexités géométriques, l'écoulement compressible à haute vitesse, et de la prise en compte du mouvement des objets. A ce sujet, cette recherche vise à relever ces défis en développant un algorithme pour simuler les écoulements compressibles sur des géométries complexes en mouvement relatif. Ce travail présente le développement d'un code numérique basé sur l'intégration de la méthode des frontières immergées (Immersed Boundary Method, IBM) avec des maillages superposés (Overset Grid) hiérarchiques pour simuler des problèmes d'écoulement compressible non visqueux. L'approche proposée réduit la complexité de la simulation numérique par rapport aux approches conventionnelles avec des maillages ajustés aux frontières. La méthode des frontières immergées choisie est la méthodologie (Sharp Interface, SI), où la géométrie est représentée avec précision en fonction de la résolution de la grille Eulérienne. Cette méthode est adaptée pour représenter des géométries rigides. La résolution de la grille est contrôlée par un raffinement local basé sur l'approche hiérarchique linéaire. Les quadtrees linéaires sont une structure de données efficace pour représenter les grilles hiérarchiques pour une représentation précise de géométries complexes avec un faible coût de traitement. Des contrôles de grille supplémentaires sont présentés, en tant que préadaptation du maillage, et pour adapter localement le domaine de calcul aux variations physiques auparavant étudiées. Une approche de marquage robuste est présentée, basée sur le parcours de la géométrie, comme l'une des contributions de cette thèse. Cette approche est basée sur un bouclage sur la courbe géométrique plutôt que sur l'ensemble de la grille Eulérienne pour identifier la frontière immergée. Cette approche prend en charge l'identification de géométries complexes sur des maillages cartésiens et hiérarchiques. Les informations de la géométrie immergée générées par l'approche de marquage sont transférées vers d'autres parties de l'algorithme pour intégrer le module des maillages superposés et le solveur fluide. La prise en compte du mouvement des objets est réalisée par l'utilisation de maillages superposés où la topologie immergée et la grille se déplacent dans un repère solidaire. Cette approche permet d'effectuer un marquage unique et constant du maillage même avec le déplacement de l'objet. En effet, on élimine le besoin de refaire le marquage de la grille à chaque pas de temps. Les applications avec des objets mobiles sont jusqu'à présent limitées aux mouvements de translation sans collisions entre les parois. Une approche de reconstruction implicite par la méthode des moindres carrés a été mise en oeuvre permettant de reconstruire la solution d'écoulement. Ce schéma prend en compte précisément la reconstruction des conditions aux limites pour les équations d'Euler au niveau des cellules d'interface avec la frontière immergée, ainsi que la reconstruction entre les maillages superposés pour transférer la solution d'écoulement entre celles via une couche de cellules interpolées autour de la frontière de la superposition. Le schéma de reconstruction a été vérifié analytiquement et numériquement pour les deux cas (pour l'application des conditions aux limites et la reconstruction pour le transfert de la solution entre les maillages superposés). Un résoluteur d'écoulement à schéma de volumes finis a été implémenté sur la base d'un schéma de Roe, où les calculs des flux sont généralisés à la discrétisation de l'espace Cartésien et hiérarchique. La méthode mise en oeuvre est du premier ordre dans l'espace et dans le temps, avec une intégration temporelle explicite. De plus, ce résoluteur est mis à jour avec la formulation Arbitrary Lagrangian-Eulerian (ALE) pour prendre en compte les vitesses de grille introduites par le module de maillages superposés. Un ensemble de cas de test de vérification sont menés et classés par type de configuration spatial (maillage unique ou maillages superposés), type d'écoulement (régime permanent ou instationnaire) et le mouvement de la géométrie (stationnaire, en mouvement à vitesse uniforme ou en mouvement accéléré) . Les cas tests réalisés vérifient l'intégration de chaque module du code développé et évaluent l'ordre de précision des schémas implémentés. Les résultats de l'algorithme développé sont limités à l'ordre de convergence du résoluteur (premier ordre) dans le cas d'écoulements sans discontinuités, et un ordre de convergence autour de (0.5) est atteint pour les cas d'écoulements discontinus ce qui concorde avec des configurations similaires disponibles dans la littérature. Cette recherche constitue, comme l'une des principales contributions de ce travail, la première version d'un algorithme entièrement basé sur l'approche IBM intégrée avec les maillages superposés hiérarchiques. De plus, ce jalon est extensible pour prendre en charge les mouvements complexes qui incluent la collision des géométries, les implémentations de schémas d'ordre élevé et les extensions 3D. La réalisation de ces extensions est en cours par l'équipe de recherche permettant de repousser les limitations de la méthodologie et du code développés et d'ouvrir de nouvelles perspectives pour des contributions scientifiques supplémentaires.