Simulation de l'accrétion de glace sur un obstacle bidimensionnel par la méthode des bissectrices et par la modélisation des ruisselets et des gouttes de surface

• Main Author: Fortin, Guy
• Format: Others
• Language: fr
• Published: 2003
• Subjects: Génie des matériaux et génie métallurgique
• Online Access: http://constellation.uqac.ca/776/1/17657455.pdf

Le LIMA (Laboratoire International des Matériaux Antigivre) en collaboration avec le CIRA (Italian Aerospace Research Centre) a développé un logiciel simulant l'accrétion de la glace en régimes sec et humide sur un objet bidimensionnel fixe. L'approche utilisée s'appuie sur les travaux de Lozowski pour les bilans énergétiques, sur une étude du comportement du film d'eau, des ruisselets et des gouttes de surface pour le calcul des rugosités et des masses d'eau résiduelle, ainsi que sur une méthode de bissectrice pour l'évolution de la surface de glace. La contribution du CIRA a été de fournir le logiciel pour le calcul des écoulements et de la captation. Le bilan énergétique basé sur la conservation de l'énergie est la sommation de la chaleur latente de fusion, d'évaporation et de sublimation, du réchauffement adiabatique et cinétique, et des pertes de chaleur par convection et conduction, ainsi que de l'évolution thermodynamique de l'eau de son état initial à son état final. La densité de la glace, qui a un impact important sur la simulation, est calculée à partir d'une corrélation empirique développée avec les cylindres tournants. En se basant sur les travaux de Al-Khalil et Hansman, le comportement des gouttes en régimes sec et humide a été décrit analytiquement, ce qui a mené à déterminer la hauteur maximale que peuvent atteindre les gouttes avant mouvement. Cette hauteur, appelée hauteur de mouvement, permet de déterminer l'état de l'eau sur la surface (film, ruisselets ou gouttes), ainsi que la hauteur des rugosités lorsque l'eau existe sous forme de gouttes ou de ruisselets. La hauteur de mouvement est déterminée par l'équilibre entre les forces de cisaillement, induites par les effets aérodynamiques et gravitationnels évalués pour une goutte non déformée, et la force de cisaillement, induite par la tension de surface et la déformation de la goutte. Elle a été validée en laboratoire et la précision obtenue pour la partie aérodynamique et gravitationnelle est de 80%. L'étude de la vague qui se forme sur le film a permis de déterminer la hauteur des rugosités lorsque l'eau existe sous forme de film. La masse d'eau résiduelle est calculée à partir des modèles analytiques élaborés selon l'état de surface et la hauteur de mouvement. Ces modèles, basés sur la physique du processus de croissance et de solidification pour le film, les ruisselets et les gouttes, interprètent la solidification à l'échelle de l'élément de surface. Le modèle pour la masse d'eau arrachée a été construit à partir d'observations numériques, il considère que toute la masse d'eau ruisselante à la surface de l'intrados est arrachée sous l'effet des forces aérodynamiques et/ou gravitationnelles. Finalement, la masse de glace accumulée est additionnée à la surface en appliquant une méthode mathématique basée sur l'aire délimitée par les bissectrices entre les panneaux. Elle permet de simuler l'évolution de la surface du dépôt en additionnant la masse de glace de façon continue, afin d'obtenir les formes complexes observées expérimentalement. Le modèle d'accrétion est validé avec les profils de glace obtenus expérimentalement en soufflerie par Shin et Bond pour un profil d'aile NACA0012 de 0,5334 m de corde, des gouttelettes d'eau surfondues de 20 um, une teneur en eau liquide de 1 g/m3 et une vitesse de 65 m/s. Ces résultats couvrent les deux régimes d'accrétion, sec et humide, dans l'intervalle de température s'échelonnant de -4,4°C à -28,3°C. La rugosité obtenue par simulation est du même ordre de grandeur que celle calculée avec la corrélation empirique développée par Ruff. Des simulations ont démontré que la variation de l'incrément de temps d'accrétion (en le diminuant de moitié ou en le doublant) et de la longueur maximale des panneaux (1%, 2% et 3% de la corde) ont peu d'impact sur la méthode des bissectrices utilisée pour le calcul de la géométrie et sur la hauteur et la distribution des rugosités. L'ajout des modèles analytiques pour le calcul des hauteurs locales des rugosités, des masses d'eau résiduelles et arrachées, ainsi que du modèle de bissectrice au modèle thermique couramment utilisé pour l'accrétion de glace sur les ailes d'avion a amélioré les résultats. Les profils de glace simulés concordent bien avec ceux mesurés en laboratoire, mais, dans la majorité des cas, le volume de glace est légèrement supérieur à celui mesuré.