Modeling of Thin Water Films on Swept Wings in Icing Condition

Les outils de simulation numérique du givrage sont d'une importance significative pour une conception optimale de l'avion et sa sécurité en vol. Des logiciels de givrage de nouvelle génération sont en cours d'implantation au sein de l'industrie aéronautique, amenant des outils de prédiction précis et applicable à tout type de géométrie. Néanmoins, une modélisation adéquate des films d'eau caractéristique de la glace de type verglas est requise pour améliorer les prédictions actuelles et pour répondre aux récentes normes de certification concernant les gouttelettes surfondues de grande taille (communément nommées SLD). Ce projet de recherche vise à modéliser correctement la glace de type verglas sur des profils d'ailes en condition de givrage. Pour atteindre cet objectif, trois aspects principaux sont considérés : l'évaluation adéquate du transfert de chaleur par convection, une bonne représentation des nuages de gouttelettes et une modélisation juste de la physique des films d'eau. De plus, les méthodes utilisées doivent être compatible pour la simulation sur des ailes en flèche infinies. En premier lieu, le calcul du coefficient de transfert de chaleur convectif est étudié afin de remplacer les méthodes semi-empiriques actuellement utilisées. Une méthode basée sur la solution des équations de Navier-Stokes est proposée afin de rendre le calcul applicable à tout type géométrie et à tout type d'écoulement. De plus, afin d'obtenir une meilleure représentation du nuage de gouttelettes et du même coup une meilleure prédiction des limites de glace, une distribution de tailles de gouttelettes est ajouté. Il est démontré que cette modification est nécessaire autant pour les tailles standards que pour les grandes tailles de gouttelettes. Enfin, l'évaluation de la force de traînée appliquée sur les gouttelettes est modifiée pour considérer la déformation de la gouttelette puisque celle-ci n'est plus strictement sphérique en régime SLD, elle tend plutôt vers la forme d'un disque. En second lieu, un modèle thermodynamique basé sur des équations aux dérivées partielles pour l'évaluation de l'épaisseur de glace, l'épaisseur du film d'eau et de la température de surface est implémenté. Ce modèle suppose un film d'eau se déplaçant selon l'effet de la friction de l'air et donne des résultats très similaires à un modèle algébrique. Le modèle est résolu de manière instationnaire permettant l'ajout de phénomènes transitoires dans la modélisation du film d'eau. Un cas pathologique est démontré sur une géométrie manufacturée comportant des zones de recirculations où une accumulation d'eau peut être observée. Cette situation est rencontrée autant pour le modèle algébrique que pour le modèle aux dérivées partielles.