L'industrie de la production primaire de l'aluminium a besoin de nouveaux développements technologiques afin d'augmenter la rentabilité des usines. L'amélioration de l'efficacité énergétique des usines s'avère une voie prometteuse pour atteindre ce but. Ce mémoire aborde l'étude de la problématique des pertes ohmiques reliées à la connexion acier-fonte-carbone de l'anode. La modélisation de ce problème requiert le couplage des champs thermo-électro-mécanique (TEM) puisque l'augmentation de la température dans l'anode influence la dilatation thermique des matériaux qui influence la pression de contact et, a fortiori, la transmission du courant électrique entre celles-ci. À l'aide de la méthode des éléments finis, il a été possible de créer un modèle numérique générique d'une anode complète utilisée dans les cuves Hall-Héroult qui tient compte de la résistance électrique de contact, du déplacement et du contact mécanique, de la conductance thermique et de la génération de chaleur par effet Joule. Le modèle générique complété, une étude de cas a été réalisée pour étudier les effets de la réduction du diamètre des rondins d'acier du tripode insérés à l'intérieur des tourillons. Avec une méthode statistique, il a été possible d'établir que la variation du diamètre des rondins extérieurs était le paramètre principal expliquant la variation des pertes ohmiques dans l'anode. Les résultats ont démontré qu'une diminution du diamètre des rondins de 180 mm (modèle de base) lorsque neufs à 160 mm suite à leur corrosion augmentait du potentiel électrique dans l'anode de l'ordre de 10 mV en comparaison avec le modèle de base. Une fois validé, ce modèle pourra être utile pour la réalisation de plusieurs travaux, comme par exemple étudier les effets de la modification de la forme et du nombre d'ailettes dans le tourillon, de la composition chimique de la fonte ou de la configuration de l'anode et d'en comparer les effets sur le potentiel électrique par rapport à un modèle base.

Le revêtement d'une cuve d'électrolyse est constitué en grande partie de matériaux à base de carbone. La caractérisation thermomécanique de ces matériaux est essentielle afin de modéliser leur comportement lors du préchauffage et en phase de production d'aluminium. Une vaste étude expérimentale a ainsi été élaborée afin de caractériser les cathodes de carbone et se veut une collaboration entre l'Université Laval et Carbone Savoie. Cette étude porte sur la détermination du module d'Young et de l'en-dommagemciit et ce, pour des températures allant de 25°C à 1000'C. Les échantillons cylindriques utilisés ont un diamètre de 50 mm, soit des échantillons de très grandes tailles considérant les contraintes techniques dans un contexte de caractérisation à haute température. Ce mémoire présente donc les résultats obtenus sur des échantillons de cathode de type graphitique pour des températures allant de la température ambiante jusqu'à la température d'opération.

Une étude expérimentale du scellement d'un ensemble anodique utilisé dans une cuve d'electrolyse permet de mieux connaître les enjeux thermiques et mécaniques rencontrés lors de la solidification du collet de fonte entre le rondin d'acier et le tourillon en carbone. Trois tourillons d'une anode cuite sont mesurés avec une machine à mesurer tridimensionnelle pour obtenir une représentation haute précision en trois dimensions de l'interface de carbone. Ces mesures sont ensuite comparées à la surface extérieure des collets de fonte solidifiés et refroidis pour connaître la dimension de l'espace d'air à l'interface fonte/carbone. Trente-neuf thermocouples, placés dans une configuration stratégique permettent la reconstruction du champ de température dans les trois matériaux depuis le moment de la coulée jusqu'au refroidissement à température ambiante. Ceci rend possible l'évaluation des coefficients de transfert de chaleur aux interfaces fonte/carbone, fonte/acier et carbone/acier à l'aide d'un modèle thermique. De plus, une analyse métallographique est associée aux courbes de refroidissement enregistrées. En parallèle, un modèle thermomécanique tridimensionnel est développé pour prédire la formation de l'espace d'air à l'interface fonte/carbone. La solidification de la fonte implique un retrait qui dépend de la composition chimique du métal, des changements de phase, des déformations d'origine thermique et des déformations viscoplastiques dans la phase semi-solide et dans l'état solide à haute température. Le logiciel de thermodynamique Thermo-Cale est utilisé pour calculer le parcours de solidification de la fonte et obtenir les propriétés thermiques nécessaires dans le modèle. Le solveur d'éléments finis commercial ABAQUS calcule la solution thermomécanique en tenant compte de la géométrie complexe d'un tourillon et des propriétés des matériaux dépendantes de la température. Pour valider n le modèle, la prédiction de l'espace d'air est comparée aux mesures expérimentales et aux données enregistrées par les thermocouples. Il a été trouvé que l'espace d'air ne peut être décrit par une simple relation linéaire en fonction de l'épaisseur de fonte. L'analyse métallographique révèle une variation significative de la microstructure à travers le collet de fonte et de la macroségrégation du carbone. Il est conclu que la meilleure évaluation des propriétés thermiques est faite en combinant le calcul à l'équilibre de Thermo-Cale, le calcul du module de Scheil et les observations expérimentales des courbes de refroidissement. Le modèle thermomécanique reproduit bien les mesures de températures, mais ne tient pas compte de la viscoplasticité à cause de problèmes de convergence, ce qui donne une prédiction incohérente de l'espace d'air. Toutefois, il est possible d'observer des tendances susceptibles de se produire dans une anode. Il est recommandé d'utiliser une nouvelle approche de modélisation en traitant le problème en plusieurs étapes indépendantes accompagnées de validations expérimentales avant l'intégration d'un modèle complet. La programmation en UMAT dans ABAQUS devrait être priviligiée. Les travaux expérimentaux futurs devraient inclure des essais mécaniques sur la fonte sur toute la plage de température concernée ainsi que des mesures d'espace d'air en temps réel avec des LVDT pendant la solidification d'un anneau de fonte. D'après les mesures d'espace d'air, il est aussi recommandé que les résultats des modélisations thermo-électro-mécaniques antérieures à ce projet soient réévalués selon les nouvelles données disponibles.

Au cours de la production de l'aluminium via le procédé de Hall-Héroult, le bain gelé, obtenu par solidification du bain électrolytique, joue un rôle significatif dans le maintien de la stabilité de la cellule d'électrolyse. L'objectif de ce travail est le développement d'un modèle numérique bidimensionnel afin de prédire le profil du bain gelé dans le système biphasé bain liquide/bain gelé, et ce, en résolvant trois problèmes physiques couplés incluant le problème de changement de phase (problème de Stefan), la variation de la composition chimique du bain et le mouvement de ce dernier. Par souci de simplification, la composition chimique du bain est supposée comme étant un système binaire. La résolution de ces trois problèmes, caractérisés par le mouvement de l'interface entre les deux phases et les discontinuités qui ont lieu à l'interface, constitue un grand défi pour les méthodes de résolution conventionnelles, basées sur le principe de la continuité des variables. En conséquence, la méthode des éléments finis étendus (XFEM) est utilisée comme alternative afin de traiter les discontinuités locales inhérentes à chaque solution tandis que la méthode de la fonction de niveaux (level-set) est exploitée pour capturer, implicitement, l'évolution de l'interface entre les deux phases. Au cours du développement de ce modèle, les problématiques suivantes : 1) l'écoulement monophasique à densité variable 2) le problème de Stefan couplé au transport d'espèces chimiques dans un système binaire sans considération du phénomène de la convection et 3) le problème de Stefan et le mouvement du fluide qui en résulte sont investigués par le biais du couplage entre deux problèmes parmi les problèmes mentionnées ci-dessus. La pertinence et la précision de ces sous-modèles sont testées à travers des comparaisons avec des solutions analytiques ou des résultats obtenus via des méthodes numériques conventionnelles. Finalement, le modèle tenant en compte les trois physiques est appliqué à la simulation de certains scénarios de solidification/fusion du système bain liquide-bain gelé. Dans cette dernière application, le mouvement du bain, induit par la différence de densité entre les deux phases ou par la force de flottabilité due aux gradients de température et/ou de concentration, est décrit par le problème de Stokes. Ce modèle se caractérise par le couplage entre différentes physiques, notamment la variation de la densité du fluide et de la température de fusion en fonction de la concentration des espèces chimiques. En outre, la méthode XFEM démontre sa précision et sa flexibilité pour traiter différents types de discontinuité tout en considérant un maillage fixe. / During the Hall-Héroult process for smelting aluminium, the ledge formed by freezing the molten bath plays a significant role in maintaining the internal working condition of the cell at stable state. The present work aims at building a vertically two-dimensional numerical model to predict the ledge profile in the bath-ledge two-phase system through solving three interactive physical problems including the phase change problem (Stefan problem), the variation of bath composition and the bath motion. For the sake of simplicity, the molten bath is regarded as a binary system in chemical composition. Solving the three involved problems characterized by the free moving internal boundary and the presence of discontinuities at the free boundary is always a challenge to the conventional continuum-based methods. Therefore, as an alternative method, the extended finite element method (XFEM) is used to handle the local discontinuities in each solution space while the interface between phases is captured implicitly by the level set method. In the course of model building, the following subjects: 1) one-phase density driven flow 2) Stefan problem without convection mechanism in the binary system 3) Stefan problem with ensuing melt flow in pure material, are investigated by coupling each two of the problems mentioned above. The accuracy of the corresponding sub-models is verified by the analytical solutions or those obtained by the conventional methods. Finally, the model by coupling three physics is applied to simulate the freezing/melting of the bath-ledge system under certain scenarios. In the final application, the bath flow is described by Stokes equations and induced either by the density jump between different phases or by the buoyancy forces produced by the temperature or/and compositional gradients. The present model is characterized by the coupling of multiple physics, especially the liquid density and the melting point are dependent on the species concentration. XFEM also exhibits its accuracy and flexibility in dealing with different types of discontinuity based on a fixed mesh.

La présente thèse présente le développement d'un premier modèle d'électrode poreuse pour une cathode Hall-Héroult ainsi que le développement de nouvelles méthodes numériques. Les nouvelles méthodes numériques sont l'implantation du critère d'électroneutralité par contrainte de type multiplicateur de Lagrange ainsi que l'application d'un équilibre thermodynamique ionique par la méthode de pénalisation. Une revue des divers processus physico-chimiques associés à la cathode est présentée. Le modèle d'électrode poreuse est développé avec une formulation ionique pour le bain comportant cinq espèces. Le modèle est développé avec une méthode d'homogénéisation pour le volume poreux. Les lois de transport sont formulées en fonction des concentrations ioniques. Les concentrations ioniques d'équilibre ont été calculées par un modèle thermodynamique simple formulé selon des données thermodynamiques empiriques. Le modèle prend en compte la migration, la diffusion et la convection des espèces ioniques. Les divers paramètres de transport ont été calés avec des mesures empiriques de conductivité conjointement aux concentrations ioniques d'équilibre. Le modèle prend en compte la réaction électrochimique de production d'aluminium au niveau des pores du carbone. De plus, un équilibre ionique est imposé entre certaines espèces. La méthode par éléments finis a été utilisée pour effectuer les calculs numériques. De nouvelles méthodes numériques d'implantation du critère d'électroneutralité et d'un équilibre ionique sont présentées. Les méthodes ont été utilisées avec succès et permettent d'éviter d'éliminer une équation de conservation. La méthode de pénalisation s'est avérée comparable aux méthodes d'élimination par le biais d'un terme source mais les résultats sont partiellement différents du modèle principal procédant par élimination d'équations de conservation. Certaines méthodes par élimination négligent en partie la contribution au transport provenant des termes associés aux équations de conservation éliminées. Par la suite, des résultats de simulation du début de l'électrolyse sont montrés. Le modèle donne des résultats conformes aux ordres de grandeur des mesures expérimentales.

Le procédé de l’électrolyse Hall-Héroult demeure la méthode principale pour la production mondiale de l’aluminium primaire depuis son invention en 1887. L'utilisation de cathodes mouillables au lieu de cathodes de carbone usuelles a été proposée afin de réduire de plus de 10% de la consommation d'énergie électrique du procédé, ce qui constitue plus de 35% de coûts de la production de l'aluminium. Cependant, à cause des conditions sévères qui prévalent dans le bain d'électrolyse, la fabrication d'une bonne cathode mouillable a été un défi au cours des 60 dernières années et aucune cathode mouillable commerciale n’est encore disponible sur le marché mondial. Dans ce projet, une nouvelle céramique poreuse a été développée par frittage sans pression de TiB2 avec additif de Ti-Fe pré-allié. Ce matériau possède les propriétés requises pour servir de cathode mouillable. Dans cette étude doctorale, le frittage en phase liquide sans pression a été choisi comme méthode de consolidation permettant la fabrication de grandes pièces à un coût relativement bas. Des essais ont été réalisés afin de comprendre l'effet de différentes conditions de traitement y compris la composition d'additif, la température de frittage, le temps de broyage, et la pré-alliage des additifs sur les propriétés physiques, mécaniques et métallurgiques ainsi que le comportement de mouillage et la stabilité des spécimens dans l'aluminium liquide. Après l'ajustement des paramètres de procédé, le matériau sélectionné a été fabriqué par le mélange de TiB2 en poudre avec 10% en poids d'additif 7Ti-3Fe pré-allié dans le broyeur à billes à haute énergie pendant 30 min, suivi par un pressage à 150 MPa et un frittage sous atmosphère de Ar/H2 à 1650°C pendant 1 h. Une microstructure sans fissures avec une distribution uniforme de pores, une densité maximale relative de 91%, une résistance à la flexion de 300 MPa et une résistivité électrique de 54 μΩ.cm ont donc été obtenues. Une goutte d’aluminium a très bien mouillé la surface de l'échantillon et une solidification isotherme s'est produite lors de sa pénétration due à l'interaction avec les additifs métalliques et la formation des phases TiAl3 et Fe4Al13. Malgré la dissolution des additifs métalliques, le matériau développé a montré une excellente stabilité après exposition dans l'aluminium fondu à 960°C pour une durée maximale de 5 jours tout en maintenant sa forme, et aucun signe d’expansion ou de gonflement n’ont été observés. Les analyses microstructurales ont révélé la formation de ponts de TiB2 entre les particules, en présence de phase liquide de Ti-Fe au cours du frittage, et donc la formation d’un squelette TiB2 qui est la cause de la stabilité du matériau développé dans l'aluminium liquide. Par conséquent, ce matériau est proposé en tant qu’un candidat fiable pour l'application en tant que cathodes mouillables dans la production d'aluminium. / Hall-Héroult electrolysis process has been the major method for world production of primary aluminum since its invention in 1887. The use of wettable cathodes instead of usual carbon cathodes has been proposed to reduce more than 10% of the electrical energy consumption of the process which constitutes more than 35% of the aluminum production costs. However, due to the severe conditions of the electrolysis bath, the fabrication of a proper wettable cathode has been a challenge during the last 60 years and no commercial wettable cathode is available in the world market yet. In this project, a novel porous ceramic by pressureless sintering of TiB2 with pre-alloyed Ti-Fe additives was developed. This material showed to meet the required properties to be used as wettable cathode. In this doctoral study, the pressureless sintering in the presence of liquid phase was selected as the consolidation method allowing the fabrication of large parts at relatively lower temperatures and costs. Experimental efforts were made in order to understand the effect of different processing conditions including additive composition, sintering temperature, milling time and pre-alloying of additives on the physical, mechanical and metallurgical properties as well as wetting behavior and stability in liquid aluminum of specimens. Based on the results of the adjustment of processing parameters, the selected material was fabricated by mixing of TiB2 powder and 10 wt% pre-alloyed 7Ti-3Fe additive in high energy ball mill for 30 min, compacting under the pressure of 150 MPa to prepare the green parts, and sintering under Ar/H2 atmosphere at 1650C for 1 h. Uniform crack-free microstructure with even distribution of pores as well as maximum relative density of 91%, bending strength of 300 MPa and electrical resistivity of 54 µΩ.cm were accordingly obtained. Aluminum drop wetted the surface of the specimen very well and isothermal solidification occurred during its penetration due to the interaction with the metallic additives and the formation of TiAl3 and Fe4Al13 phases. Despite of the dissolution of metallic additives, this material showed excellent stability after being exposed to molten aluminum at 960C for up to 5 days by maintaining its shape and no sign of expansion or swelling was observed. Microstructural investigation revealed the precipitation of inter-particle bridges of TiB2 nature in the presence of Ti-Fe liquid phase during sintering forming a TiB2 skeleton, which is the cause of the stability of the developed material in liquid aluminum. This material is proposed as a reliable candidate for application as wettable cathodes in aluminum smelting.

L’industrie de l’aluminium vit des bouleversements économiques, d’où l’importance de réduire les coûts d’opération. Le présent projet se penche sur les pertes ohmiques du scellement anodique. Les essais ont été complétés entre les conditions ambiantes jusqu’à charge et température maximales permises par le banc d’essais. Les travaux furent complétés sur un banc d’essais novateur - permettant une stabilisation rapide de la température et une faible oxydation de la portion carbonée de l’échantillon - et des échantillons (acier-fonte-carbone) fabriqués pour reproduire la réalité. Ensuite, une régression non linéaire a été utilisée pour modéliser le comportement mesuré en laboratoire. L’analyse statistique démontre que le modèle représente adéquatement le comportement des interfaces à mesure que la charge appliquée et la température s’intensifient. La comparaison avec la littérature démontre une similitude, toutefois les résultats diffèrent de façon significative. Ceci s’explique par l’utilisation de matériaux différents et par le changement de méthode expérimentale. / The aluminum industry faces economical headwinds, thus needing to reduce its operating costs. The project studies specifically energy losses at anode sealing location. Experiments were conducted from room conditions up to the maximum capacity of the bench test. Laboratory work was completed using tri-material samples (made in such a way to replicate accurately reality) and an innovative bench test using magnetic induction as a source of heat that enables minimal heat-up time and sample oxidation during experiments. Non-linear regression was then used to retrieve a model from laboratory results. Analysis showed that a power model represents accurately the interfaces behavior from room to operation conditions. Comparison with literature showed that the order of magnitude is the same but results are not similar. This observation can be explained by the use of slightly different materials and also due to the use of a different experimental procedure.