Les modèles de séparation d'une émulsion huile-eau sont utilisés pour la simulation des séparateurs gravitationnels huile-eau. Cependant, la plupart des modèles existants ne considèrent pas la coalescence des particules d'huile. Or, afin d'améliorer la précision d'un modèle de séparation, il est préférable de prendre compte le phénomène de coalescence. Nous présentons un modèle eulérien-lagrangien de séparation qui tient compte de la coalescence des particules. Celui-ci est constitué d'un modèle eulérien de séparation couplé à un modèle lagrangien tenant compte de la coalescence. Le développement du modèle lagrangien constitue le point central du travail car l'important est de compléter un modèle eulérien. a priori quelconque, à l'aide d'un modèle lagrangien de coalescence. Pour nous assurer de la validité du modèle eulérien-lagrangien, on simule, à l'aide de ce dernier, la séparation d'une émulsion huile-eau pour un séparateur gravitationnel simple. La comparaison de nos résultats avec ceux d'une expérience physique de ce genre déjà réalisée montre que le modèle reproduit bien la réalité de cette expérience. Antonin Dallaire, étudiant Sylvain Boivin, directeur de recherche Rung Tien Bui, co-directeur de recherche

L'accumulation de givre sur les diverses composantes d'un avion en attente au sol peut causer des pertes de performances aérodynamiques au décollage qui dans nombre de cas se sont avérées néfastes. L'avènement des fluides antigivre, afin de protéger les surfaces portantes de la précipitation a permis de réduire de beaucoup les pertes de portance. Toutefois, des tests en soufflerie réfrigérée dans des conditions de dilution homogène ( eau diluée préalablement au fluide protecteur) ont montré que les résidus de fluide protecteurs eux-mêmes induisent une perte de portance au décollage. Ces pertes sont toutefois acceptables pour un taux inférieur à 5% en vertu de la marge de manoeuvre du pilote et de leur effet temporaire. Cependant dans le cas du vol réel, il s'agit plutôt d'une dilution dynamique, c'est-à-dire que la précipitation vient se diluer dans un fluide concentré. Au basses températures (-20°C) il y a formation d'une mince strate glaciale au dessus du fluide protecteur ce qui est susceptible d'engendrer une perte accrue de portance par rapport à la dilution homogène. Afin de caractériser ces pertes de portance, il a fallu mettre en place une infrastructure expérimentale. Il y a la conception et la réalisation d'une balance aérodynamique basée sur la théorie des déformations élastiques mesurées par jauges de contraintes. Elle possède une ossature rigide tout en étant dotée d'une bonne sensibilité. Il y a la minimisation des effets de parois latérales par l'implantation de murs de séparation qui permettent la bidimensionnalisation du modèle d'aile. Les effets des parois supérieure et inférieure furent compensés par un terme correctif appliqué sur les coefficients de portance. Il y a également l'instauration d'un système de gicleurs simulant la bruine verglaçante de 200 fim de diamètre pouvant opérer jusqu'à - 20 °C. Le tout est géré par un logiciel d'application réalisé par certains membres du LIMA. La simulation en soufflerie fût réalisée sur un modèle symétrique NACA 0018 bidimensionnel. Les tests furent effectués pour une densité de précipitation relativement forte (28,5 g/dm2h) à deux températures (-10°C, -20°C), deux taux de dilution (15%, 30%) et ce pour trois fluides commerciaux de rhéologie non-newtonienne. Le décollage est simulé à l'aide d'une rampe d'accélération de 2,6 m/s2 dans un l'intervalle de vitesse 2 < U < 50 m/s. La rotation de l'aile débute à 40 m/s et elle tourne au taux de 2.9 °/s. Dans chaque condition, trois tests ont été effectués: un test à sec (Courbe de référence), un test en dilution dynamique et un dernier en dilution homogène. À - 20°C et à 30% de dilution, les résultats obtenus pour des angles utilisés lors du décollage (1 à 8°) démontrent une perte de portance accrue pour la dilution dynamique en comparaison à la dilution homogène et ce pour chaque fluide. Les pertes de portance à CLmax dépassent largement la norme de sécurité de 5%. À -10°C, il n'y a pas de perte de portance significative par rapport au test à sec. On conclue que les fluides anti-givres sont efficaces dans bon nombre de situations. Toutefois, le mode de dilution dynamique qui s'apparente au vol réel démontre que l'utilisation de chacun de ces fluides est inacceptable au point de vue des normes de sécurité dans les conditions extrêmes de forte précipitation à basse température.

Alors que des recherches ont déjà été effectuées sur la résistance en traction et en compression de la glace atmosphérique (givre et verglas), il n'existe à peu près pas de valeurs connues pour la résistance au cisaillement de ce type de glace. La connaissance de ces valeurs et des modèles de comportement rhéologique en cisaillement est nécessaire pour l'analyse du comportement des structures sous l'effet du givre et du verglas, pour la détermination des critères de rupture applicables à la glace atmosphérique et pour l'évaluation des contraintes admissibles dans ce matériau pour les domaines fragile et ductile. L'objectif de cette recherche est de développer une méthode d'évaluation de la résistance au cisaillement de la glace atmosphérique en fonction des conditions météorologiques de formation et du taux de déformation, à partir d'échantillons formés en laboratoire. L'objectif premier est de déterminer les propriétés mécaniques du matériau; cependant, la diversité des essais de cisaillement effectués au cours de la présente étude a été également planifiée pour évaluer les influences des différents paramètres dimensionnels de l'appareil et de l'éprouvette sur les mesures obtenues. C'est pourquoi un seul type de glace a été utilisé pour la fabrication des éprouvettes. Les échantillons de glace atmosphérique sont formés en laboratoire dans des conditions préétablies et contrôlées. Les échantillons obtenus sont plus adéquats que ceux formés dans des conditions atmosphériques naturelles. La glace, du givre dur dont la masse volumique est d'environ 850 kg/m3, est formée à partir de gouttelettes surfondues projetées sur un obstacle (le collecteur) qui est placé dans une soufflerie en circuit ouvert. Les valeurs des paramètres atmosphériques de formation utilisés sont de -14°C pour la température de l'air dans la chambre froide, 16 m/s pour la vitesse de l'air dans la soufflerie, 1,2 g/m3 pour la teneur en eau de l'air et de 40 um pour le diamètre volumique médian des gouttelettes. Les éprouvettes sont usinées à 68 ou 70 mm de diamètre et à 8,10 ou 12 mm d'épaisseur sur un tour parallèle. Un groupe de 6 éprouvettes peut être produit sur le même manchon (glace et collecteur). Finalement, le collecteur est démonté pour permettre la manipulation des éprouvettes. L'évaluation de la résistance au cisaillement est effectuée par une méthode de poinçonnage en utilisant une presse mécanique à vitesse de déplacement contrôlée. L'appareil de cisaillement utilisé a été conçu de façon à permettre le confinement total ou partiel de l'échantillon de glace. Dans ces deux modes de chargement, les 213 ruptures des éprouvettes de glace ont été effectuées à des taux de déformation (y) variant de 0,0001 à 40 sec-1. Des vitesses de déplacement du poinçon de 0,002 à 400 mm/min ont été utilisées pour obtenir ces taux de déformation. En confinement partiel, les contraintes maximales de rupture passent par un maximum dans la zone de transition ductile-fragile du matériau. Elles atteignent 1,5 MPa au taux de déformation de 0,03 sec-1. En confinement total, la contrainte de rupture croît avec le taux de déformation et atteint 8 MPa pour le taux de 40 sec1. Quoiqu'il soit utilisé pour l'identification des essais, il a été possible de démontrer que la valeur du taux de déformation, déterminé par la relation entre le déplacement du poinçon et le jeu radial entre ce dernier et la matrice, ne représente pas le comportement réel de l'éprouvette. Il a été constaté que la valeur du module de cisaillement variait de façon presque proportionnelle avec celle du jeu, ce qui conduit à l'inexactitude du modèle de déformation utilisé. Cette relation existe entre le jeu et les valeurs du module puisque les dimensions utilisées n'ont pas eu d'influence sur les taux de montée en charge de l'éprouvette. Cependant, la dimension du jeu occasionne une légère baisse de la valeur des contraintes maximales lorsque sa dimension est supérieure à 2 mm. En utilisant le modèle de déformation énoncé dans cette étude, la valeur moyenne du module de cisaillement est 2,2 MPa. En se servant de la même méthode que Lavrov (1969) pour déterminer la déformation appliquée à l'éprouvette, les valeurs de ce module se situent entre 12 et 73 MPa pour les essais en confinement total et entre 7 et 40 MPa pour ceux en confinement partiel. Finalement, ces travaux ont permis de démontrer que la valeur de la rigidité de cisaillement brute de l'éprouvette de glace variait proportionnellement avec celle de la contrainte maximale. Cette rigidité représente le taux de montée en charge de la contrainte par mètre de déplacement du poinçon. Elle atteint 7,5 GPa/m au taux de déformation de 40 sec-1 pour les essais effectués en confinement total et 5 GPa/m dans la zone de transition pour ceux effectués en confinement partiel.

La formation de glace sur les avions au sol, qui est nuisible au décollage, est empêchée par l'application de fluides antigivre et dégivrant. La présence de ces fluides sur les ailes, tout comme la glace, au moment de l'accélération diminue la force verticale (portance) appliquée par l'écoulement d'air sur l'aile. Cette diminution est directement reliée à l'épaisseur de la couche limite de l'air qui se développe au-dessus du fluide. Des procédures expérimentales permettent de définir le comportement de l'écoulement de l'air, sur l'aile d'avion, en présence du film fluide à partir des mesures de l'épaisseur de déplacement de la couche limite (EDCL) sur une plaque plane horizontale. Le coût de ces essais force la recherche à s'orienter vers la prédiction analytique basée essentiellement sur les caractéristiques du fluide appliqué. Ce projet est une contribution à la recherche dans ce domaine qui consiste à la détermination de méthode de prédiction des pertes de portance en fonction des caractéristiques du fluide. Le projet consiste en une étude approfondie du comportement pseudo-plastique (viscosité) des fluides antigivre commerciaux associée à celle du développement de la couche limite sur une plaque plane recouverte d'un film fluide antigivre durant une accélération similaire à celle d'un décollage d'un avion de ligne. Un volet analytique s'ajoute aux résultats expérimentaux, d'une part pour obtenir un algorithme de calculs qui permet de prédire la valeur de l'EDCL en turbulence, à partir du profil de vitesses, en conjonction avec l'équation de Von Karman et d'autre part pour prédire les conditions qui amènent l'instabilité hydrodynamique de l'interface air/fluide selon la méthode développée par Yih. Cette étude a permis de définir que les vagues à l'interface air/fluide sont générées à la transition laminaire-turbulent à une vitesse de l'ordre de 15 m/s. De plus, on constate qu'une partie importante du fluide est éliminée lorsque la présence des vagues, qui augmentent le frottement, force le fluide pseudo-plastique à se liquéfier. C'est à ce moment, que l'EDCL est le plus élevé et par conséquent, la perte de portance maximum. Les résultats montrent que le taux d'élimination du fluide au moment de cette instabilité est maximum. De plus, cette étude permet de comparer les hauteurs des vagues à celles des rugosités dynamiques équivalentes. En effet, les vagues ont un impact semblable à celui des rugosités sur une plaque plane rugueuse. Ainsi, lorsque le régime d'écoulement de l'air atteint 65 m/s après 30 s, les effets de hauteur et de densité des vagues sur l'EDCL sont du même ordre. Durant la période d'élimination massique du fluide, au début de l'instabilité et de la formation des vagues, le phénomène rugueux est plus complexe car les vagues forment des amas de fluide qui bloquent l'écoulement de l'air. Une application pratique de cette étude permet de définir que les vitesses de rotation sur le mandrin du viscosimètre nécessaires pour créer la même sollicitation que l'air sur le fluide, appliquées sur l'aile, sont de l'ordre de 0,3 RPM lorsque l'avion est au repos et 200 RPM lorsque l'avion de ligne quitte le sol.

Ce mémoire présente les résultats d'une étude dont le but est d'évaluer l'influence de quatre paramètres sur les échanges thermiques à l'intérieur d'un four d'homogénéisation de lingots d'aluminium. Cette étude a été réalisée par l'intermédiaire d'un modèle numérique, tri-dimensionnel et en régime transitoire, élaboré à l'aide des logiciels commerciaux suivants: PATRAN (PDA Engineering) comme pré- et post-processeur, FLOTRAN (Compuflo) pour la simulation de l'écoulement des gaz chauds et P/THERMAL, pour la simulation du rayonnement et de la conduction de la chaleur à l'intérieur du lingot, chaleur reçue par convection forcée et/ou par rayonnement. Le modèle permet de quantifier l'apport de chaleur dû au rayonnement et celui dû à la convection forcée. Il permet également d'étudier, sans s'y limiter, l'influence des paramètres suivants sur les échanges thermiques à l'intérieur du four: la distance entre les lingots, le débit massique des gaz, la température des gaz et la longueur des lingots. Les résultats de l'étude montrent que le rayonnement est responsable d'environ 24% de l'apport total de chaleur. Quant au transfert de chaleur par convection forcée, on trouve d'abord qu'il est optimum pour une distance bien précise entre les lingots, ensuite, qu'il varie proportionnellement au débit massique, à la température des gaz et à la longueur des lingots. Il est possible d'obtenir, par calcul, une très bonne évaluation du facteur de proportionnalité en se référant à la méthode de calcul du coefficient de film convectif pour un écoulement libre sur une surface plane. Cette relation entre les résultats numériques et les prédictions empiriques est intéressante car elle permet d'obtenir, sans faire de longues simulations, le coefficient de film moyen à la surface du lingot pour des cas qui n'ont pas été considérés dans cette étude. En fonction des hypothèses de départ et des résultats obtenus, le modèle numérique, élaboré au cours de ces travaux de recherche, est donc apte à simuler les échanges thermiques à l'intérieur d'un four d'homogénéisation.

Le but de ce travail est de concevoir une méthode permettant de calculer le coefficient transitoire de transfert de chaleur à l'interface moule/métal lors de la solidification de l'aluminium pur commercial. La méthode de calcul proposée consiste à déterminer la température à la surface de l'aluminium, la température à la paroi interne du moule et le flux de chaleur traversant l'interface moule/métal. L'obtention de ces valeurs se fait par l'utilisation complémentaire d'une expérimentation et d'un modèle numérique. L'expérimentation permet d'obtenir les conditions initiales et aux limites du modèle numérique. Le montage expérimental se compose d'un moule cylindrique (hauteur de 22.5 cm, rayon interne de 5 cm et rayon externe de 9 cm) permettant d'obtenir une solidification unidirectionnelle et radiale, d'un système d'alimentation, et de capteurs de température. Le modèle numérique repose sur les équations de conduction de chaleur, développées en coordonnées cylindriques pour la direction radiale. Le moule et l'aluminium solide sont traités avec les équations de conduction de chaleur, et l'aluminium liquide et diphasique sont traités avec la méthode enthalpique de Voiler. Diverses techniques ont été développées pour combiner l'expérimentation et le modèle numérique: une méthode de lissage permettant d'éliminer le bruit contenu dans les données expérimentales; une méthode spécifique de calcul des conditions aux limites à partir des valeurs expérimentales; et, enfin, une méthode de correction du modèle unidimensionnel qui prend en compte la perte de chaleur induite par la dissymétrie verticale. Les résultats obtenus permettent de valider l'expérimentation et le modèle numérique, et, aussi, d'évaluer les causes d'erreurs dues aux capteurs utilisés. Le coefficient de transfert de chaleur à l'interface moule/métal montre un comportement transitoire jusqu'à 200 s puis quasi stationnaire par la suite. La valeur maximale du coefficient est de l'ordre de 2400 W/m2Kpour un écart de température de l'ordre de 500°C entre la paroi interne du moule et la surface de l'aluminium liquide. La valeur du coefficient apparaît comme relativement indépendante des conditions initiales de température (dans le domaine d'utilisation industrielle). La situation expérimentale obtenue s'avère correctement unidimensionnelle dans le moule. L'approximation unidimensionnelle est moins bonne dans l'aluminium, mais la méthode de calcul reste suffisamment précise, de l'ordre de 4%, dans la partie transitoire. La méthode décrite peut être utilisée pour déterminer les valeurs du coefficient en fonction des conditions de moulage comme les températures initiales du moule et de l'aluminium, le type et l'épaisseur du revêtement recouvrant les parois du moule. Après modification, elle pourra être aussi utilisée pour les alliages.

L'étude des événements de givrage atmosphérique sur les réseaux de transport de l'énergie électrique constitue une préoccupation importante des habitants des régions de climat froid, telles que la province de Québec. Ceci inclut la compréhension de plusieurs phénomènes complexes en vue de réduire les risques de dommages aux réseaux de transmission existants et à ceux qui seront être implantés dans le futur. Dans cette perspective, il est important de comprendre ces phénomènes par l'analyse des événements de givrage atmosphériques antérieurs et de les modéliser. La présente recherche, effectuée dans le cadre des travaux de la Chaire industrielle NSERC/HYDRO-QUÉBEC/UQAC sur le givrage atmosphérique des équipements des réseaux électriques (CIGELE), avait pour but d'analyser les événements passés tels qu'enregistrés par le réseau de mesure SYGIVRE d'Hydro-Québec, et de créer un modèle de prédiction en temps réel de l'évolution des tempêtes de verglas pour différentes régions du Québec. Les informations météorologiques disponibles dans la base de données SYGIVRE, mesurées à l'aide d'un givromètre amélioré (ice rate meter), couvrent 6 ans, à partir de 1992, ceci pour un total de 28 stations. Une technique par réseaux de neurones a été choisie comme base pour le travail de modélisation. L'analyse exploratoire des événements de verglas contenus dans la base de données SYGIVRE constitue la première étape de l'étude. Cette étude a consisté à extraire de la base de données les tempêtes de verglas, à faire ressortir les caractéristiques de chaque paramètre de la base de données et à rendre les données adaptées au traitement par réseaux de neurones. Quatre groupes de stations ayant un rapport géographique et météorologique sont impliqués dans la prédiction, chacun étant traité séparément. Un premier modèle prédictif basé sur la technique des réseaux de neurones a été conçu en vue d'étudier l'évolution des tempêtes en se basant sur la relation spatiale qui existe entre les stations d'un groupe. La variable à prédire à chaque point pour une station est une variable dichotomique qui prend la valeur 1 s'il y a un événement de givre et 0 sinon. Ce modèle constitue une grande amélioration comparativement à des modèles antérieurs décrits dans la littérature, basés sur la méthode de régression logistique. Afin d'obtenir un modèle plus réaliste, un second modèle a été créé pour prédire en temps réel le poids de la glace accumulée sur la structure. Le modèle utilise les informations binaires avec des variables météorologiques, et la prédiction est ajustée pour évaluer différents temps futurs. Les résultats ont montré que la meilleure performance possible du modèle peut être réalisée en ajoutant au modèle de prédiction les variables de température et du poids de la glace. Toutefois, le modèle détecte encore difficilement les événements de verglas extrême. En conclusion, il est apparu qu'inclure un plus grand nombre d'années d'observations aux deux modèles devrait améliorer la performance de la prédiction, puisque ceci augmenterait la quantité des informations à étudier. Il s'agit d'ailleurs d'une condition essentielle aux prédictions par réseaux de neurones.

Ce mémoire décrit sur une nouvelle méthode pour briser la glace en utilisant les ondes de choc. Le travail effectué a consisté à concevoir, fabriquer et caractériser un dispositif pouvant focaliser des ondes de choc dans l'air. Trois réflecteurs elliptiques ont été utilisés de façon à focaliser des ondes de choc sur un disque de glace. Les ondes de choc sont générées par une décharge électrique d'environ 10 kV dans l'air à -4,5°C au foyer primaire du réflecteur. Les échantillons de glace qui ont un diamètre de 8,5 cm et une épaisseur variant de 3 à 10 mm, sont placés au foyer secondaire. L'efficacité des réflecteurs a été caractérisée dans les directions axiales et radiales à l'aide d'un capteur de pression commercial. La pression maximale enregistrée au point focal est de 400 kPa (près de quatre fois la pression atmosphérique) avec un temps de montée de 2 us. La rupture complète a été observée sur des échantillons d'une épaisseur maximale de 9,2 mm après un seul passage de l'onde de choc. Les échantillons dont l'épaisseur est comprise entre 9,2 et 10 mm se fracturent après plusieurs coups et finalement ceux dont l'épaisseur est supérieure à 10 mm ne se fracturent pas. Le processus de rupture semble être reliée à un processus de flexion en domaine fragile. Une hypothèse, reliée à la présence des bulles dans les échantillons de glace, a été proposée pour expliquer la rupture de ceux-ci par plusieurs coups successifs.

Les réseaux de transmissions électriques et de télécommunications des pays froids sont soumis à des conditions météorologiques extrêmes caractérisées par des tempêtes de givre ou de verglas. En janvier 1998, la région de Montréal (Québec, Canada) a été paralysée par une tempête de verglas d'une durée exceptionnelle entraînant la rupture mécanique de nombreux câbles de transport d'électricité et même de lignes. Un site de givrage naturel comportant deux lignes électriques expérimentales sur le Mont Valin, au nord de Chicoutimi (Québec, Canada), fournit une base de données contenant douze saisons de givrage, de la saison 86-87 à la saison 97-98. Ces lignes expérimentales ne sont pas alimentées en courant électrique. La première ligne installée sur le site est un câble témoin d'une portée de 31,4 mètres. La seconde, appelée ligne expérimentale, est composée d'un câble de garde, d'un câble Bersimis et d'un faisceau de quatre conducteurs Bersimis. Sa portée est de 96,515 mètres. Chaque saison de givrage se décompose en événements de glace atmosphérique, un événement pouvant comprendre trois phases successives: une phase d'accrétion de la glace, une phase de maintien de la glace et une phase de délestage de la glace. La phase de maintien de la glace est un concept nouveau qu'il faut définir clairement à l'aide de critères. Cette phase caractérise l'état de stabilité de la glace qui se maintient de façon aléatoire sur les câbles. En général, cette phase est la plus longue des trois phases; elle peut durer de quelques heures à quelques mois. Nous distinguons habituellement trois grands types de glace d'origine atmosphérique susceptibles d'entraîner des surcharges importantes sur les lignes aériennes par dépôt et formation éventuelle de manchons autour des conducteurs : le givre, le verglas et la neige mouillée. Dans l'atmosphère, les nuages, constitués de très petites gouttelettes d'eau et parfois de cristaux de glace, sont à l'origine de la formation de la glace atmosphérique. Sur le site du Mont Valin, les câbles sont principalement recouverts de dépôts de givre et de verglas, parfois de neige mouillée. La neige mouillée doit être différenciée de la neige sèche. Pour le verglas et le givre, la rigidité en torsion d'un conducteur ou d'un câble a un effet sur la quantité de glace qui peut s'y déposer et sur la forme du manchon de glace. Sur un conducteur toronné, la rigidité en torsion du conducteur influence la forme axiale ou cylindrique de l'accumulation de neige alors que la rigidité en torsion n'a aucun effet sur un conducteur lisse recouvert de neige. Des critères ont été définis afin d'identifier le début et la fin de la phase de maintien de chaque événement des douze saisons de givrage et afin d'identifier les différents types de glace (givre, verglas, neige mouillée) et les mélanges de ces types de glace présents sur les câbles pendant la phase de maintien. Dans le cas des mélanges, l'ordre successif de formation de chaque type de glace (givre, verglas, neige mouillée) a été pris en compte. L'utilisation de ces critères a permis de créer une base de données sur Excel pour chaque type de câble, donc pour le câble Bersimis, le câble du faisceau, le câble de garde et le câble témoin. Ces données se rapportent à la phase de maintien et à la charge maximale atteinte pendant la phase de maintien pour 75 événements du câble Bersimis, 75 événements du câble du faisceau, 85 événements du câble de garde et 62 événements du câble témoin. Quatorze mélanges ont été identifiés. Une analyse de régression multiple a été réalisée sur la charge de glace maximale des événements pendant la phase de maintien en fonction des variables suivantes : le type de glace ou de mélange correspondant à la charge maximale, la durée entre le début de la phase de maintien et l'instant où la charge maximale est atteinte, la charge de glace en fin de phase d'accrétion et le taux d'accrétion moyen pendant la phase d'accrétion. Le résultat attendu est de savoir, à partir de l'analyse de régression linéaire, quels paramètres influencent le plus la formation de la charge maximale d'un événement de glace atmosphérique. Une description statistique de la charge maximale a aussi été effectuée par type de câble et par type de glace et de mélange. La phase de maintien a été analysée en considérant les paramètres suivants : les types de glace et les mélanges présents sur les câbles pendant la phase de maintien, la température moyenne de l'air pendant la phase de maintien, la variation de la charge de glace pendant la phase de maintien, la vitesse moyenne et maximale du vent perpendiculaire à la ligne expérimentale pendant la phase de maintien, la durée de la phase de maintien, les délestages de glace pendant la phase de maintien et les types de délestages en fin de phase de maintien. En conclusion, l'analyse de régression a montré que le type de glace ou de mélange correspondant à la charge maximale n'a pas d'influence sur la valeur de la charge maximale. La durée de la phase de maintien diffère d'un câble à l'autre. Ceci est lié au fait que les charges de glace maximales et les délestages de glace différent d'un câble à l'autre. L'analyse statistique est un outil très intéressant qu'il faut savoir interpréter avec précaution en tenant toujours compte des phénomènes physiques qui se rapportent au sujet traité.

Ce projet s'inscrit dans des études à approfondir la connaissance du givrage afin d'adapter les éoliennes pour opérer dans les conditions nordiques. L'objectif principal est de développer des techniques pour évaluer, à long terme, l'impact du givrage sur le fonctionnement d'un projet éolien et d'optimiser des techniques de dégivrage. Plusieurs études complexes doivent être effectuées avant de proposer un modèle précis capable de réaliser une simulation tridimensionnelle du givrage sur des pales des éoliennes en rotation. Pour atteindre les objectifs, nous avons utilisé des logiciels commerciaux de CFD et l'interface conviviale de MS-Excel pour valider un cas test des simulations de l'accumulation de glace sur un cylindre. Ces méthodes sont destinées à être appliquées à la simulation en 3D de l'accumulation de glace sur les éoliennes, étant donné que l'étude du givrage du cylindre, pour lequel des résultats analytiques et expérimentaux nombreux sont disponibles pour validation, est fondamentale et pré-requis pour ce domaine de recherche. La capacité d'un objet à capturer des gouttelettes d'eau dans un écoulement est appelé l'efficacité de collection. Une évaluation de l'efficacité de collection sur un cylindre a été simulée en utilisant une approche Eulérienne basée sur les modèles de turbulence multiphasiques dans ANSYS-CFX. Les résultats ont été validés avec des résultats des approches Lagrangiennes et Eulériennes dans FLUENT ainsi qu'avec des résultats expérimentaux obtenus à partir d'études antérieurs. Les résultats ont démontré l'efficacité des modèles multiphasiques du logiciel CFX à simuler les fractions volumiques d'eau et à définir la zone de collection et les limites d'impact autour du cylindre. Pareillement, un logiciel pour calculer les trajectoires des gouttelettes d'eau dans un écoulement d'air interceptées par un cylindre a été réalisé sous Excel avec un code VBA (Visual Basic for Applications). Cette interface permet de démontrer les différents scénarios pouvant aider à valider les simulations avec CFX. Les deux simulations réalisées avec Excel et CFX ont démontré une cohérence entre le comportement des lignes de courant, des trajectoires et des forces agissant sur les gouttelettes. - This project has been conducted in the context of intensifying the knowledge of icing on wind turbines to be adapted in Nordic conditions. The main objective is to develop techniques to evaluate the long-term impact of icing on the performance of wind turbine projects and to optimize the de-icing techniques. Several upstream complex studies are to be conducted prior to propose a precise model capable of achieving a three-dimensional simulation of icing on wind turbine rotating blades. In order to overcome these difficulties, we made use of the high capacity of commercial CFD software, together with the friendly user interface of MS-Excel. The validation of these tools for the ice accretion past a cylinder is the prerequisite for applying these tools for the 3D simulation of icing over wind turbines, given that the study of cylinder icing is fundamental in several areas of icing research. The ability of an object to capture water droplets presented in an air flow is called collection efficiency. Evaluation of the collection efficiency over a cylinder has been simulated using Eulerian approach based on multi-phase turbulence models in CFX. The results have been validated with those of Lagrangian and Eulerian approaches in FLUENT as well as with experimental results obtained from previous studies. The results showed that the use of ANSYS-CFX multi-phase models was effective in simulating water volume fractions through the domain and in defining the collection zone and the impingement limits around the cylinder. In addition, a code to calculate the trajectories of supercooled water droplets in the air intercepted by a cylinder has been achieved using MS-Excel sheets supported with VBA (Visual Basic for Applications). This interface has been used to demonstrate different scenarios that can help to validate similar simulations using CFX. In both Excel and CFX simulations, the streamlines and trajectories have demonstrated similar behaviour which is consistent with the forces acting on water droplets.