Méthodes d'éléments finis pour le problème de changement de phase en milieux composites

Dans ces travaux de thèse on s’intéresse au développement d’un outil numérique pour résoudre le problème de conduction instationnaire avec changement de phase dans un milieu composite constitué d’une mousse de graphite infiltrée par un matériau à changement de phase tel que le sel, dans le contexte du stockage de l’énergie thermique solaire.Au chapitre 1, on commence par présenter le modèle sur lequel on va travailler. Il estséparé en trois sous-parties : un problème de conduction de chaleur dans la mousse, un problème de changement de phase dans les pores remplis de sel et une condition de résistance thermique de contact entre les deux matériaux qui est traduite par une discontinuité du champ de température.Au chapitre 2, on étudie le problème stationnaire de conduction thermique dans un milieu composite avec résistance de contact. Ceci permet de se focaliser sur la plus grande difficulté présente dans le problème qui est le traitement de la condition de saut à l’interface.Deux méthodes d’éléments finis sont proposées pour résoudre ce problème : une méthode basée sur les éléments finis Lagrange P1 et une méthode hybride-duale utilisant les éléments finis Raviart-Thomas d’ordre 0 et P0. L’analyse numérique des deux méthodes est effectuée et les résultats de tests numériques attestent des efficacités des deux méthodes [10]. Les matériaux à changement de phase qu’on étudie dans le cadre de cette thèse sont des matériaux pures, par conséquent le changement de phase s’effectue en une valeur de température fixe qui est la température de fusion. Ceci est modélisé par un saut dans la fonction fraction liquide et par conséquent dans la fonction enthalpie du matériau. Cette discontinuité représente une difficulté numérique supplémentaire qu’on propose de surmonter en introduisant un intervalle de régularisation autour de la température de fusion.Cette procédure est présentée dans le chapitre 3 où une étude analytique et numérique montre que l’erreur sur la température se comporte comme " en dehors de la zone de mélange, où " est la largeur de l’intervalle de régularisation. Cependant, à l’intérieur l’erreur se comporte comme p " et on montre que cette estimation est optimale. Cette diminution de vitesse de convergence est due à l’énergie qui reste bloquée dans la zone de mélange [58].Dans le chapitre 4 on présente quatre des schémas les plus utilisés pour le traitement de la non-linearité due au changement de phase: mise à jour du terme source, linéarisation de l’e