Thèse soutenue

Modélisation et simulation du remplissage de moules verriers : "Prise en compte du transfert radiatif"

FR  |  
EN
Auteur / Autrice : Hoang Quan Nguyen
Direction : Alain DegiovanniBenjamin Rémy
Type : Thèse de doctorat
Discipline(s) : Mécanique et énergétique
Date : Soutenance le 02/10/2009
Etablissement(s) : Vandoeuvre-les-Nancy, INPL
Ecole(s) doctorale(s) : EMMA - Ecole Doctorale Energie - Mécanique - Matériaux
Partenaire(s) de recherche : Laboratoire : Laboratoire d'énergétique et de mécanique théorique et appliquée (Nancy)
Jury : Président / Présidente : Gérard Jeandel
Examinateurs / Examinatrices : Alain Degiovanni, Benjamin Rémy, Gérard Jeandel, Jean-Paul Caltagirone, Didier Delaunay, Franck Humbert
Rapporteurs / Rapporteuses : Jean-Paul Caltagirone, Didier Delaunay

Résumé

FR  |  
EN

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