Parmi les scenarios accidentels hypothétiques étudiés dans le cadre des études de sûreté des Réacteurs à Eau Pressurisée (REP), figure l'Accident par Perte de Réfrigérant Primaire (APRP). Dans ce scénario, les gaines (en alliage de zirconium) qui contiennent le combustible nucléaire, sont soumises à une oxydation importante à haute température (T>1200°C) dans de la vapeur d'eau. Les gaines étant la première barrière de confinement des radioéléments, il est primordial qu'elles conservent une certaine ductilité résiduelle après la trempe pour assurer leur intégrité. Ces propriétés sont directement liées aux cinétiques de croissance de la zircone et de la phase alphaZr(O), ainsi qu'au profil d'oxygène dans le métal au cours du transitoire. Dans ce cadre, la compréhension et la modélisation du phénomène d'oxydation et de diffusion de l'oxygène dans les alliages de zirconium à haute température ont fait l'objet de cette thèse. Le modèle cinétique (EKINOX-Zr), développé au cours de cette thèse, est basé sur la résolution numérique d'un problème de diffusion/réaction avec des conditions aux limites sur trois interfaces mobiles : gaz/oxyde, oxyde/alpha-Zr, alpha-Zr/bêta-Zr. Le couplage du code cinétique avec le logiciel ThermoCalc et la Zircobase permet de prendre en compte l'influence des éléments d'alliages (Sn, Fe, Cr, Nb) et de l'hydrogène.
Cette étude s'est plus particulièrement intéressée à deux aspects de l'APRP : l'influence d'une couche de pré-oxyde (formée aux températures nominales de fonctionnement du réacteur) et les effets de l'hydrogène. Grâce au couplage avec la base thermodynamique Zircobase, l'effet de l'hydrogène sur les limites de solubilité de l'oxygène dans les différentes phases peut être pris en compte dans le modèle cinétique. Les simulations permettent ainsi de reproduire les profils d'oxygène mesurés sur différents échantillons pré-hydrurés. L'existence de couche de pré-oxyde de forte épaisseur peut amener à une réduction transitoire de la couche de pré-oxyde dans les premiers instants du palier à haute température sous vapeur d'eau, avant formation de l'oxyde haute température. Une première série de simulations à l'aide du modèle cinétique EKINOX-Zr a permis de reproduire qualitativement ce chemin cinétique transitoire et a montré que cette couche formée à basse température possède des propriétés de diffusion particulières. Des expériences de traceurs 18O ont été réalisées sur des gaines pré-oxydées en autoclave pour étudier la diffusion de l'oxygène à haute température dans ces couches formées à basse température. D'autre part le modèle cinétique EKINOX-Zr a aussi été modifié pour prendre en compte la diffusion des traceurs. La confrontation des expériences avec les calculs a permis d'étudier les propriétés de diffusion particulières de la couche d'oxyde basse-température, ainsi qu'une évolution des propriétés de diffusion de la phase alphaZr(O) formée par dissolution de l'oxyde à haute température. |
One of the hypothetical accident studied in the field of the safety studies of Pressurized light Water Reactor (PWR) is the Loss-Of-Coolant-Accident (LOCA). In this scenario, zirconium alloy fuel claddings could undergo an important oxidation in a high temperature (T>1200°C) steam environment. Cladding tubes constitute in the first confinement barrier of radioelements and then it is essential that they keep a certain amount of ductility after quenching ensuring their integrity. These properties are directly related to the growth kinetics of both oxide and alphaZr(O) phase, and also to oxygen diffusion profile in the cladding tube after the transient. In this context, this work was dedicated to the understanding and the modeling of the both oxidation phenomenon and oxygen diffusion in zirconium based alloys at high temperature. The numerical tool (EKINOX-Zr) used in this thesis is based on a numerical resolution of a diffusion/reaction problem with equilibrium-conditions on three moving boundaries: gas/oxide, oxide/alpha-Zr, alpha-Zr/bêta-Zr. EKINOX-Zr kinetics model is coupled with ThermoCalc software and the Zircobase database to take into account the influence of the alloying elements (Sn, Fe, Cr, Nb) but also the influence of hydrogen.
This study focused on two parts of the LOCA scenario: the influence of a pre-oxide layer (formed in-service) and the effects of hydrogen. Thanks to the link between EKINOX-Zr and the thermodynamic database Zircobase, the hydrogen effects on oxygen solubility limit can be considered in numerical simulations. Thus, simulations can reproduce the oxygen diffusion profiles measured in pre-hydrided samples. The existence of a thick pre-oxide layer on cladding tubes can induce a reduction of the pre-oxide layer before the growth of a high-temperature one during the high temperature dwell under steam. The first simulations performed using the numerical tool EKINOX-Zr demonstrated that this particular kinetics pathway is qualitatively well reproduced bay the modeling. It also showed that the pre-oxide layer (formed at low-temperature) has particular diffusion properties. Consequently, some 18O tracer experiments were done on pre-oxidized cladding tubes (pre-oxidizing done in autoclave) in order to study the diffusion properties of these oxide layers formed at low-temperature. Furthermore, EKINOX-Zr was modified to take into account the diffusion of tracer. Both particular diffusion properties of the low-temperature oxide layer and the evolution of the diffusion properties of the alphaZr(O) phase (formed by reduction of the oxide layer) were demonstrated by comparing experiments and numerical simulations. |