Le sujet de l'étude concerne la caractérisation des propriétés thermiques de liquides semi-transparents à haute température : les liquides silicatés. La caractérisation de ces matériaux est particulièrement délicate car il faut séparer les contributions des différents modes de transfert thermique (conduction, convection et rayonnement) si l'on veut mesurer des propriétés intrinsèques. Le dispositif expérimental est basé sur une méthode pulsée de type flash. La cellule de mesure, placée au centre d'un four tubulaire, est soumise à une excitation délivrée par un laser continu et l'élévation de température est mesurée à l'aide d'un détecteur infrarouge sur la face opposée à l'excitation. La modélisation du problème direct prend en compte le couplage conducto-radiatif en résolvant l'équation de la chaleur en régime transitoire et l'équation de transfert radiatif (ETR) pour un milieu gris émettant, absorbant et non diffusant à l'aide de la méthode des harmoniques sphériques (approximation P1). La méthode développée permet d'estimer simultanément la diffusivité thermique et un coefficient d'absorption moyen (gris) par méthode inverse. Les résultats expérimentaux obtenus sur des liquides silicatés présentant des propriétés radiatives différentes (i.e., des coefficients d'absorption différents) valident la méthode de caractérisation développée et mise en œuvre dans nos travaux / The study deals with the characterization of the thermal properties of semi-transparent liquids at high temperature: molten silicates. The characterization of such materials is particularly challenging because the contributions of the different heat transfer modes (conduction, convection and radiation) should be accounted for to allow the measurement of intrinsic properties. The experimental setup is based on a transient pulse method derived from the flash method. A heat flux stimulation is generated on the front face of an experimental cell, placed at the center of a tube furnace, with a continuous laser beam and the temperature rise is measured by an infrared detector on the opposite side. The modeling of the direct problem takes into account the conducto-radiative coupling by solving both the heat equation and the radiative transfer equation (RTE) for a gray emitting, absorbing but non-scattering medium with the spherical harmonics method (P1 approximation). The developed method allows to simultaneously estimate the thermal diffusivity and a mean (gray) absorption coefficient by an inverse method. The experimental results obtained on molten silicates with different radiative properties (i.e., different absorption coefficients) validate the characterization method developed and implemented in our work

L’objet de ce travail est de proposer un modèle adapté pour la simulation du remplissage de moules qui réponde au meilleur compromis entre temps de calcul et précision des résultats. La difficulté est double. Il faut prendre en compte le phénomène de remplissage qui est un problème complexe à frontières libres et les spécificités liées au Verre : viscosité fortement thermodépendante et température de fusion élevée qui nécessite de prendre en compte le rayonnement. Le Chapitre I est consacrée à la partie écoulement du Verre liquide. La bibliothèque numérique Aquilon/Thétis, adaptée pour traiter ce type de problèmes et les couplages thermique air/verre/parois, a été utilisée (Méthode V.O.F pour le suivi de l’interface, méthodes de type Lagrangien augmenté/Projection vectorielle pour le couplage Vitesse-Pression). Pour l’aspect radiatif, différentes approches sont proposées : conductivité radiative équivalente (Chapitre II), méthode explicite directe pour la validation (Chapitre III) et méthode d’harmoniques sphériques ou méthode PN (Chapitre IV). Dans le Chapitre V, la méthode PN retenue est validée dans des cas simples et est appliquée ensuite à des cas avec couplage convectif en géométries complexes et obstacles semi-transparents (1D, 2D et 3D, 2D axi-symétrique et milieu non gris). Une version P1 modifiée est présentée. Les résultats sont assez proches de ceux donnés par la méthode P3 avec des temps de calcul modestes. L’intérêt de ce modèle est qu’il est facilement intégrable dans des codes numériques existants : une seule équation différentielle du second ordre stationnaire à résoudre en 3D / The aim of this study is to propose an adapted model for the simulation of mould filling that must be a compromise solution between computational time and results accuracy. The double difficulty is to take into account the filling phenomenon that is a complex problem due to the presence of free boundaries and to the Glass specificities: viscosity that is highly thermal dependant and high melting temperature that requires taking into account radiation effects. Chapter I is devoted to the melting Glass flow. The numerical libraries Aquilon/Thétis, adapted for solving such type of problems and the thermal coupling between Air/Glass/Walls, has been used. (V.O.F method for front tracking, Augmented Lagrangian/Vector Projection methods for solving Pressure/Velocity coupling). For radiative aspect, different approaches are proposed: equivalent radiative conductivity (Chapter II), direct explicit method for validation (Chapter III) and spherical harmonics method or PN method (Chapter IV). In the Chapter V, the selected PN method is validated through simple cases and is then applied in other cases with convective coupling in complex geometries including semi-transparent inclusions (1D, 2D and 3D, 2D axi-symmetric and non grey medium). A P1 modified version is presented. The results are close to those given by P3 method but with reduced computational time. The main interest of this model is that it can be easily implemented in existing numerical codes: a single stationary second order partial differential equation to solve in 3D