Constellation, le dépôt institutionnel de l'Université du Québec à Chicoutimi

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

Fortin Guy. (2003). 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. Thèse de doctorat, Université du Québec à Chicoutimi.

[thumbnail of 17657455.pdf] PDF
20MB

Résumé

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é.

Type de document:Thèse ou mémoire de l'UQAC (Thèse de doctorat)
Date:2003
Lieu de publication:Chicoutimi
Programme d'étude:Doctorat en ingénierie
Nombre de pages:216
ISBN:1412310709
Identifiant unique:10.1522/17657455
Sujets:Sciences naturelles et génie > Génie > Génie des matériaux et génie métallurgique
Département, module, service et unité de recherche:Départements et modules > Département des sciences appliquées > Programmes d'études de cycles supérieurs en ingénierie
Directeur(s), Co-directeur(s) et responsable(s):Laforte, Jean-Louis
Ilinca, Adrian
Mots-clés:Glace--Formation--Simulation par ordinateur, Gouttes--Modèles mathématiques, Verglas, Ice formation--Computer simulation, Drops--Mathematical models, Glaze (Meteorology), THESE, ACCRETION, GLACE, SIMULATION, MODELISATION, MATHEMATIQUE, METHODE, BISSECTRICE, GOUTTE, RUISSELET, OBJET, BIDIMENTIONNEL, GIVRE, BILAN, EAU, DYNAMIQUE, NUMERIQUE, GOUTELETTE, AILE, AVION, VERGLAS
Déposé le:01 janv. 2003 12:34
Dernière modification:03 juin 2011 13:16
Afficher les statistiques de telechargements

Éditer le document (administrateurs uniquement)

Creative Commons LicenseSauf indication contraire, les documents archivés dans Constellation sont rendus disponibles selon les termes de la licence Creative Commons "Paternité, pas d'utilisation commerciale, pas de modification" 2.5 Canada.

Bibliothèque Paul-Émile-Boulet, UQAC
555, boulevard de l'Université
Chicoutimi (Québec)  CANADA G7H 2B1
418 545-5011, poste 5630