Immersed Boundary Methods for the Three-Dimensional Modeling of In-Flight Ice Accretion

Le givrage en vol est un danger bien connu pour la sécurité et les performances des aéronefs. Il se produit généralement dans des conditions météorologiques froides, lorsqu’un avion rencontre un nuage surfondu. L’un des principaux effets de ce phénomène est la dégradation des performances aérodynamiques. Malgré le développement de systèmes de protection contre le givrage, les avions commerciaux doivent pouvoir démontrer leur capacité à voler dans de telles conditions. Par conséquent, l’intérêt pour le développement de simulations de givrage en vol s’est accru ces dernières années, avec un focus particulier sur les simulations 3D. Les simulations numériques du givrage reposent généralement sur l’hypothèse que le temps d’exposition de la surface aux conditions de givrage est significativement plus long que le temps caractéristique de l’écoulement aérodynamique. Cette différence d’échelle temporelle permet une stratégie de modélisation modulaire, dans laquelle le processus global est décomposé en une succession d’appels de modules. Tout d’abord, un maillage volumique est généré à partir de la surface considérée afin de résoudre l’écoulement aérodynamique stationnaire et les trajectoires des gouttes surfondues. Ensuite, un bilan thermodynamique est fait à la surface afin de déterminer l’épaisseur locale de glace. Enfin, la géométrie est mise à jour pour intégrer la nouvelle couche de glace formée. Dans une simulation multi-step, le temps d’exposition est subdivisé en plusieurs incréments. Cette approche répète de manière itérative la succession des modules, permettant ainsi à la simulation de prendre en compte l’évolution de la géométrie et son influence sur l’écoulement et l’impact des gouttes, améliorant ainsi la précision de la solution. Cependant, des défis importants apparaissent lors de l’automatisation de ce processus, principalement en raison de la complexité associée à la génération de maillages de haute qualité à chaque itération. L’objectif principal de ce travail est d’intégrer les méthodes de frontières immergées (Immersed Boundary Method, IBM) dans le cadre de simulation du givrage afin d’éliminer la nécessité de générer un maillage volumique à chaque étape, tout en conservant le niveau de précision obtenu avec une approche conforme (Body-Fitted, BF). Ceci permet d’utiliser un seul maillage volumique tout au long des étapes d’une simulation multi-steps, facilitant ainsi l’automatisation du processus. Tout d’abord, une méthode IBM est développée pour les simulations 3D d’écoulement d’air, où l’écoulement non-visqueux est modélisé en résolvant les équations d’Euler. Cependant, comme le bilan thermodynamique à la surface dépend des échanges diffusifs, la solution d’Euler doit être couplée à un solveur de couche limite. Une méthode IBM basée sur un forçage discret à partir d’un maillage conforme à la géométrie initiale lisse a été développée afin de prendre en compte la modification de la géométrie à chaque pas d’accrétion. Ensuite, une IBM pour les équations eulériennes 3D des gouttes est nécessaire. En raison de la nature des équations et des conditions aux limites, une méthode de pénalisation a été choisie. Cette méthode a été validée dans des simulations 3D, démontrant la préservation de la précision de la solution.