Modeling of the solid-liquid phase change. Application to the latent heat thermal energy storage Suitable for concentrated solar power plant

L'un des principaux leviers technologiques permettant le développement industriel de process de production énergétique renouvelable et à haute efficacité, consiste en l'élaboration d'une solution innovante de stockage de l'énergie. Ce système de stockage doit permettre de lisser la période de production et ainsi de suivre au plus près les besoins des consommateurs. Parmi les solutions existantes, le stockage thermique par chaleur latente présente de nombreux avantages qui font qu'aujourd'hui il fait l'objet de plusieurs travaux de recherche et de développement. Cette technologie est basée sur le principe que certaines classes de matériaux, appelés matériaux à changement de phase (MCP), libèrent (transition liquide/solide) ou accumulent (transition solide/liquide) de l'énergie lorsqu'ils sont soumis à un changement de phase. En amont du développement d'un design de stockage, il est essentiel de comprendre et de maitriser les processus thermiques entrant en jeu lors des phases de fusion et de solidification du matériau et cette compréhension passe par le développement de modèles numériques adaptés aux problématiques rencontrées. Dans le cadre de ce manuscrit, la filière technologique qui nous intéresse est celle des centrales solaires à concentration. Porté par l'ADEME dans le cadre du projet STARS (Stockage Thermique appliqué à l'extension de pRoduction d'énergie Solaire thermodynamique), le travail réalisé au sein du LaTEP consiste à analyser les performances d'une solution de stockage via la modélisation de cette dernière en considérant les phénomènes thermiques et hydrauliques. Le travail de modélisation est effectué à l'aide du logiciel de CFD libre de droit OpenFOAM dans lequel est développé et implémenté, par le laboratoire, un module dédié au problème qui nous concerne, à savoir la résolution eulérienne (maillage fixe) des équations de conservation pour un fluide incompressible, en présence d'un changement de phase solide-liquide dominé par des mouvements convectifs (convection-dominated phase change). Concernant les problèmes de transition de phase, diverses méthodes mathématiques et numériques ont été développées pour rendre compte finement de la physique de ces phénomènes. Après avoir effectué une revue de ces dernières dans la première partie du manuscrit, nous avons sélectionné deux formulations que nous avons implémenté dans OpenFOAM. Une fois ce travail réalisé nous avons taché de comparer les résultats renvoyés par ces différentes formulations en les confrontant aux résultats expérimentaux disponibles dans la littérature. Cela nous a permis d'une part de nous conforter dans l'utilisation de nos solveurs et sur la pertinence des résultats obtenus avec ces derniers et d'autre part de mettre en évidence les écarts entre les solutions renvoyées par chaque formulation. Fort de ce constat, nous avons souhaité évaluer l'impact de l'équation d'état utilisée pour relier l'enthalpie et la température, indispensable à la fermeture thermodynamique du système d'équations. Cette comparaison s'est faite par la simulation d'un échangeur type stockage thermique (simulations en 2D) et par l'analyse des performances de ce dernier lors des phases de stockage, de déstockage et au cours de plusieurs séries de cycles. Les résultats obtenus nous ont permis de conclure sur l'importance d'une bonne caractérisation des MCP afin de pouvoir modéliser leur comportement au plus juste via la formulation mathématique et la loi d'état la plus adaptéeA major technological lever to the industrial development of renewable energy production processes with high efficiency, is the development of an innovative solution to store the energy. This storage device should help to smooth the production period and to follow closely the demand. Among the existing solutions, latent heat thermal storage has many advantages that make today it is the subject of several research and development projects. This technology is based on the principle that certain classes of material, called phase change materials (PCMs), release (during liquid to solid transition) or accumulate (during solid to liquid transition) energy when subjected to a phase change. Upstream of the development of a new storage design, it is essential to understand and master the thermal processes involved in the melting and solidification phase of the material and this knowledge comes through the development of numerical models adapted to the problems encountered. In the particular context of this Phdthesis, the technological process that interests us is that of CSP (Concentrated Solar Power). Funded by ADEME under the STARS Project (Thermal STorage Applied to the expansion of pRoduction of thermodynamic Solar energy), the work done by the LaTEP is to analyze the performance of a storage solution by modeling the latter, considering the thermal and hydraulic phenomena. The modeling work is done with the free source OpenFOAM CFD software in which is developed and implemented by the laboratory, a new module dedicated to the problem we are concerned, namely the resolution of Eulerian (fixed grid) conservation equations for an incompressible fluid in the presence of a solid-liquid phase change dominated by convective motions. Regarding the phase transition problems, various mathematical and numerical methods have been developed to finely consider the physics of these phenomena. After conducting a review of methods in the first part of the Phd thesis, we selected two formulations that we have implemented in OpenFOAM. Once this work done, we have managed to compare the results returned from these formulations by comparing them with experimental results available in the literature and also with analytical cases. This allowed us firstly to strengthen us in the use of our solvers and the accuracy of the obtained results and secondly to highlight the differences between the solutions returned by each formulation. After that, we wanted to assess the impact of the equation of state used to connect the enthalpy to the temperature, essential for closing the thermodynamic equations. This comparison was made by the simulation of a thermal storage exchanger (2D simulation) and by analyzing the performance of this latter during the charge phase, the discharge one and during several series of cycles. The obtained results allowed us to conclude about the importance of a good characterization of PCM in order to model their behavior as accurately via the mathematical formulation and the most suitable state la