La thèse est centrée sur la modélisation de la physique du comportement d’un alliage binaire et l’implémentation du meilleur modèle mathématique pour simuler le changement de phase liquide solide en tenant compte de la vitesse de refroidissement, la vitesse de solidification, la ségrégation, la convection naturelle et la surfusion afin d’optimiser la capacité de stockage de chaleur d'un tel matériau. Dans le présent travail, les températures pour lesquelles le changement de phase s'opère sont estimées grâce aux diagrammes des phases et la méthodologie CALPHAD qui retraduisent les différentes phases d'un alliage binaire, y compris la transformation isotherme. Pour cela, la minimisation de l'énergie de Gibbs est résolue dans un code de calcul développé à cette occasion et aboutit à l'identification des phases stables du matériau. Pour un intervalle de température souhaite le code permet d'estimer rapidement la décharge de chaleur pour la composition de l'alliage sélectionné en équilibre ou hors équilibre. Dans la méthode proposée, la vitesse de refroidissement du système permet de calculer la vitesse de solidification. Puis,celle-ci établit la relation entre la cinétique globale et la macrostructure. Basé sur le modèle de non-équilibre local, qui dépend de la variation du coefficient de partition, le degré de surfusion est prédit à partir de la vitesse de refroidissement appliquée. Une étude bibliographique a été réalisée pour amener une comparaison numérique et assurer la capacité de notre méthode à reproduire le changement de phase, en incluant des phénomènes spécifiques tels que la surfusion et la recalescence. / Latent Heat Thermal Energy Storage (LHTES) shows high storage density compared to sensible thermal systems. For high temperature applications, the use of alloys as phase change materials presents many advantages. Principally, varying alloy composition allows controlling the storage\discharge of thermal energy through an expected temperature range (defined by the heat source), and the high thermal conductivity givessuitable heat transfer properties to the system that receives/supplies the energy. However, some systems need a specific temperature range to correctly operate. In such conditions, subcooling (also known as undercooling) and segregation are undesirable phenomena in alloys when they are used as PCM. In thepresent work, we propose a method to predict the latent heat release during phase transformation of a binaryalloy submitted to a controlled cooling rate, including subcooling, segregation and variation of composition.This thesis describes the physical models that apply when heat is released from such a material. We takeinto consideration the cooling rate applied to the PCM, the solidification velocity, convective phenomena,melting temperature and subcooling. In the present work, phase diagrams and the CALPHAD methodologyare used to determine the temperature range for phase change (or constant temperature value for isothermal transformation) by minimizing the Gibbs equilibrium energy. The Gibbs free energy minimization has been implemented in a homemade numerical code. The material can be screened with different compositions for equilibrium or off-equilibrium solidification allowing quick selection of the optimal material for the specific heatsource. In the proposed method, the solidification velocity is obtained from the cooling rate. Then, variationin microstructure is driven by the solidification velocity using the local non-equilibrium diffusion model. Based on the local nonequilibrium model that depends on the partition coefficient variation, the subcooling degree, wich is derived from the applied cooling rate is predicted. A bibliographic study has been carried out and anumerical comparison has been undertaken to ensure the capacity of our code to reproduce the phase change of various materials that include phenomena such as subcooling and recalescence. The results highlight that the cooling rate is one of the most important parameters in the performance of the thermal storage, having a large effect on segregation and subcooling degree. Moreover, we show the influence ofpartition coefficient on the time evolution of solid fraction, considering a constant or a composition-dependent value. We can conclude that the latent heat release can be correctly predicted provided that the method correctly predicts the phase diagram and the variable partition coefficient. This work helps to accelerate the design and development of thermal storage systems and lays the foundation to continue exploring other kinds of materials (e.g. paraffins).

Au-dessus d'un taux de combustion seuil ≥ 7 at %, les produits de fission volatils Cs, I, et Te ou métalliques (Mo) sont partiellement relâchés hors du combustible et finissent par constituer une couche de composés de PF qui remplit progressivement le jeu existant entre la périphérie de la pastille et la surface interne de la gaine en acier inoxydable. Nous appelons cette couche JOG pour Joint Oxyde-Gaine. Mon sujet de thèse est axé sur l'étude thermodynamique du système (Cs, I, Te, Mo, O) + (U, Pu) ainsi que sur l'étude de la diffusion de ces produits de fission à travers le combustible vers le jeu combustible-gaine pour former le JOG.L'étude thermodynamique constitue la première étape de mon travail. Sur la base d'une analyse critique des données expérimentales issues de la littérature, les systèmes Cs-Te, Cs-I, Cs-Mo-O ont été modélisés par la méthode CALPHAD. En parallèle, une étude expérimentale a été entreprise pour valider la modélisation CALPHAD du système binaire Cs-Te. Dans une deuxième étape, les données thermodynamiques résultant de la modélisation CALPHAD ont été introduites dans la base de données du code de calcul thermodynamique ANGE (code interne au CEA dérivé du logiciel SOLGASMIX) dont la finalité est le calcul de la composition chimique du combustible irradié. Dans une troisième étape, le code de calcul thermodynamique ANGE (Advanced Numeric Gibbs Energy minimiser) a été couplé avec le code de simulation du comportement thermomécanique du combustible des RNR-Na GERMINAL V2. / For a burn-up higher than 7 at%, the volatile FP like Cs, I and Te or metallic (Mo) are partially released from the fuel pellet in order to form a layer of compounds between the outer surface of the fuel and the inner surface of the stainless cladding. This layer is called the JOG, french acronym for Joint-Oxyde-Gaine.My subject is focused on two topics: the thermodynamic study of the (Cs-I-Te-Mo-O) system and the migration of those FP towards the gap to form the JOG.The thermodynamic study was the first step of my work. On the basis of critical literature survey, the following systems have been optimized by the CALPHAD method: Cs-Te, Cs-I and Cs-Mo-O. In parallel, an experimental study is undertaken in order to validate our CALPHAD modelling of the Cs-Te system. In a second step, the thermodynamic data coming from the CALPHAD modelling have been introduced into the database that we use with the thermochemical computation code ANGE (CEA code derived from the SOLGASMIX software) in order to calculate the chemical composition of the irradiated fuel versus burn-up and temperature. In a third and last step, the thermochemical computation code ANGE (Advanced Numeric Gibbs Energy minimizer) has been coupled with the fuel performance code GERMINAL V2, which simulates the thermo-mechanical behavior of SFR fuel.

The aim of this work was to improve the heat capacity estimation of a material for usage within a CALPHAD-type assessment. An algorithm is derived that estimates the trend of heat capacity with temperature based on zero Kelvin properties and the thermal expansion coefficient at the Debye temperature. The algorithm predicts not only the trend of heat capacity but also the temperature trend of the volume and the bulk modulus, which can be also included in new thermodynamic databases. The algorithm is used to assess thermophysical properties of the intermetallic phases eta (Fe2Al5), epsilon~(Fe5Al8) and tau4 (FeAl3Si2). The heat capacity of the intermetallic phases zeta, eta, theta and epsilon of the Al-Fe system and of tau4 of the Al-Fe-Si system was measured using DSC. For the phases zeta, eta, and theta, a non-linearly increasing heat capacity approaching the melting temperature was observed. In addition, the heat capacity of three bcc-based Al-Fe samples including the B2-->A2 transition were determined. The Al-rich section of the Al-Fe phase diagram was studied using DTA and quenching experiments. The homogeneity ranges of the intermetallic phases were determined using SEM/WDS measurements. Based on own and literature values, a thermodynamic description of the Al-Fe system was assessed including the modelling of A2/B2 ordering and the homogeneity range of all intermetallic phases. In addition, thermodynamic parameters of the Al-Fe-Si, Al-Fe-Mg, and the Fe-Mg-Si system were assessed to obtain a thermodynamic description of the Al-rich side of the Al-Fe-Si-Mg system, which can be used to study phase transitions of typical A356-aluminium alloys.

<p>Combination of the Calphad method and theoretical calculations provides new possibilities for the study of materials science. This work is a part of the efforts within the CCT project (Centre of Computational Thermodynamics) to combine these methods to facilitate modeling and to extend the thermodynamic databases with critically assessed volume data. In this work, the theoretical calculations refer to first-principles and Debye-Grüneisen calculations. The first-principles (i.e. ab initio) electronic structure calculations, based on the Density- Functional Theory (DFT), are capable of predicting various physical properties at 0 K, such as formation energy, volume and bulk modulus. The ab initio simulation software, VASP, was used to calculate the binding curves (i.e. equation of state at 0 K) of metallic elements, cubic carbides and nitrides. From the binding curves, the equilibrium volumes at 0 K were calculated for several metastable structures as well as stable structures. The vibrational contribution to the free energy was calculated using the Debye-Grüneisen model combined with first-principles calculations. Two different approximations for the Grüneisen parameter, γ, were used in the Debye-Grüneisen model, i.e. Slater’s and Dugdale-MacDonald’s expressions. The thermal electronic contribution was evaluated from the calculated electronic density of states. The calculated thermal expansivities for metallic elements, cubic carbides and nitrides were compared with Calphad assessments. It was found that the experimental data are within the limits of the calculations using the two approximations for γ. By fitting experimental heat capacity and thermal expansivity around Debye temperatures, we obtained optimal Poisson’s ratio values and used them to evaluate Young’s and Shear moduli. In order to reach a reasonable agreement with the experiments, it is necessary to use the logarithmic averaged mass of the constitutional atoms. The agreements between the calculations and experiments are generally better for bulk modulus and Young’s modulus than that for shear modulus. A new model describing thermodynamic properties at high pressures was implemented in Thermo-Calc. The model is based on an empirical relation between volume and isothermal bulk modulus. Pure Fe and solid MgO were assessed using this model. Solution phases will be considered in a future work to check the model for compositional dependence.</p>

Combination of the Calphad method and theoretical calculations provides new possibilities for the study of materials science. This work is a part of the efforts within the CCT project (Centre of Computational Thermodynamics) to combine these methods to facilitate modeling and to extend the thermodynamic databases with critically assessed volume data. In this work, the theoretical calculations refer to first-principles and Debye-Grüneisen calculations. The first-principles (i.e. ab initio) electronic structure calculations, based on the Density- Functional Theory (DFT), are capable of predicting various physical properties at 0 K, such as formation energy, volume and bulk modulus. The ab initio simulation software, VASP, was used to calculate the binding curves (i.e. equation of state at 0 K) of metallic elements, cubic carbides and nitrides. From the binding curves, the equilibrium volumes at 0 K were calculated for several metastable structures as well as stable structures. The vibrational contribution to the free energy was calculated using the Debye-Grüneisen model combined with first-principles calculations. Two different approximations for the Grüneisen parameter, γ, were used in the Debye-Grüneisen model, i.e. Slater’s and Dugdale-MacDonald’s expressions. The thermal electronic contribution was evaluated from the calculated electronic density of states. The calculated thermal expansivities for metallic elements, cubic carbides and nitrides were compared with Calphad assessments. It was found that the experimental data are within the limits of the calculations using the two approximations for γ. By fitting experimental heat capacity and thermal expansivity around Debye temperatures, we obtained optimal Poisson’s ratio values and used them to evaluate Young’s and Shear moduli. In order to reach a reasonable agreement with the experiments, it is necessary to use the logarithmic averaged mass of the constitutional atoms. The agreements between the calculations and experiments are generally better for bulk modulus and Young’s modulus than that for shear modulus. A new model describing thermodynamic properties at high pressures was implemented in Thermo-Calc. The model is based on an empirical relation between volume and isothermal bulk modulus. Pure Fe and solid MgO were assessed using this model. Solution phases will be considered in a future work to check the model for compositional dependence.

Alumino-silicate systems are of great interest for materials scientists and geochemists. Thermodynamic knowledge of these systems is useful in steel and ceramic industries, and for understanding geochemical processes. A popular and efficient approach used to obtain a self-consistent thermodynamic dataset is called CALPHAD. It couples phase diagram information and thermochemical data with the assistance of computer models. The CALPHAD approach is applied in this thesis to the thermodynamic modelling and assessments of the CaO-Al2O3-SiO2, MgO-Al2O3-SiO2 and Y2O3-Al2O3-SiO2 systems and their subsystems. The compound energy formalism is used for all the solution phases including mullite, YAM, spinel and halite. In particular, the ionic two sub-lattice model is applied to the liquid solution phase. Based both on recent experimental investigations and theoretical studies, a new species, AlO2-1, is introduced to model liquid Al2O3. Thus, the liquid model corresponding for a ternary Al2O3-SiO2-M2Om system has the formula (Al+3,M+m)P (AlO2-1,O-2, SiO4-4,SiO20)Q, where M+m stands for Ca+2, Mg+2 or Y+3. This model overcomes the long-existing difficulty of suppressing the liquid miscibility gap in the ternary systems originating from the Al2O3-free side during the assessments. All the available and updated experimental information in these systems are critically evaluated and finally a self-consistent thermodynamic dataset is achieved. The database can be used along with software for Gibbs energy minimization to calculate any type of phase diagram and all thermodynamic properties. Various phase diagrams, isothermal and isoplethal sections, and thermochemical properties are presented and compared with the experimental data. Model calculated site fractions of species are also discussed. All optimization processes and calculations are performed using the Thermo-Calc software package. / QC 20100607

The recent emphasis on magnesium alloys can be appreciated by following the research push from several agencies, universities and editorial efforts. With a density equal to two-thirds of Al and one-thirds of steel, Mg provides the best opportunity for lightweighting of metallic components. However, one key bottleneck restricting its insertion into industrial applications is low strength values. In this respect, Mg-Y-Nd alloys have been promising due to their ability to form strengthening precipitates on the prismatic plane. However, if the strength is compared to Al alloys, these alloys are not attractive. The primary reason for low structural performance in Mg is related to low alloying and microstructural efficiency. In this dissertation, these terminologies are discussed in detail. A simple calculation showed that the microstructural efficiency in Mg-4Y-3Nd alloy is 30% of its maximum potential. Guided by the definitions of alloying and microstructural efficiency, the two prime objectives of this thesis were to: (i) to use thermomechanical processing routes to tailor the microstructure and achieve high strength in an Mg-4Y-3Nd alloy, and (ii) optimize the alloy chemistry of the Mg-rare earth alloy and design a novel rare—earth free Mg alloy by Calphad approach to achieve a strength of 500 MPa. Experimental, theoretical and computational approaches have been used to establish the process-structure-property relationships in an Mg-4Y-3Nd alloy. For example, increase in strength was observed after post aging of the friction stir processed/additive manufactured microstructure. This was attributed to the dissolution of Mg2Y particles which increased the alloying and microstructural efficiency. Further quantification by numerical modeling showed that the effective diffusivity during friction stir processing and friction stir welding is 60 times faster than in the absence of concurrent deformation leading to the dissolution of thermally stable particles. In addition, the investigation on the interaction between dislocations and strengthening precipitate revealed that, specific defects like the I1 fault aid in the accelerated precipitation of the strengthening precipitate in an Mg-4Y-3Nd alloy. Also, the effect of external field (ultrasonic waves) was studied in detail and showed accelerated age hardening response in Mg-4Y-3Nd alloy by a factor of 24. As the bottleneck of low strength is addressed, the answers to the following questions are discussed in this dissertation: What are the fundamental micro-mechanisms governing second phase evolution in an Mg-4Y-3Nd alloy? What is the mechanical response of different microstructural states obtained by hot rolling, friction stir processing and friction stir additive manufacturing? Is defect engineering critical to achieve high strength Mg alloys? Can application of an external field influence the age hardening response in an Mg-4Y-3Nd alloy? Can a combination of innovative processing for tailoring microstructures and computational alloy design lead to new and effective paths for application of magnesium alloys?