Ecole Doctorale
Physique et Sciences de la Matière
Spécialité
PHYSIQUE & SCIENCES DE LA MATIERE - Spécialité : ENERGIE, RAYONNEMENT ET PLASMA
Etablissement
Aix-Marseille Université
Mots Clés
neutronique,réacteur rapide,schéma de calcul,ASTRID,APOLLO3,méthode 2D/1D
Keywords
neutron transport,fast reactor,calculation scheme,ASTRID,APOLLO3,2D/1D method
Titre de thèse
Développement de schémas de calcul neutronique pour des coeurs nucléaires hétérogènes refroidis au sodium dans le code APOLLO3
Development of Neutronic Calculation Schemes for Heterogeneous Sodium-Cooled Nuclear Cores in the APOLLO3 Code
Date
Vendredi 27 Septembre 2019 à 10:30
Adresse
Faculté de droit d'Aix-Marseille
3 avenue Robert Schuman
13100 Aix-en-Provence Salle du conseil no 1
Jury
Directeur de these |
M. Laurent BUIRON |
Aix Marseille Université - CEA Cadarache |
Rapporteur |
M. Hugues DELORME |
INSTN |
Rapporteur |
M. Alain HéBERT |
Ecole Polytechnique de Montréal |
Examinateur |
Mme Raphaèle HERBIN |
Aix Marseille Université |
Examinateur |
M. Pascal ARCHIER |
CEA Cadarache |
Examinateur |
M. Enrico GIRARDI |
EDF Lab |
Résumé de la thèse
Les réacteurs nucléaires refroidis au sodium offrent des perspectives intéressantes pour la filière nucléaire car ils permettent dexploiter tout le potentiel énergétique de luranium naturel, tout en contribuant à la réduction de la radiotoxicité des déchets nucléaires.
Cependant, la nécessité délever le niveau de sûreté de ces réacteurs aux standards du XXIe siècle tend à augmenter la complexité des curs. En particulier, la recherche dun comportement pardonnant en situation de transitoire non protégé conduit à minimiser la réactivité de vidange sodium : cest lidée du cur CFV dASTRID.
Du point de vue de la modélisation neutronique, cet accroissement du niveau dhétérogénéité géométrique constitue une difficulté supplémentaire. Ainsi, les objectifs de la thèse sont lidentification des principaux phénomènes physiques devant être pris en compte lors du calcul neutronique de curs hétérogènes en spectre rapide, ainsi que le développement de schémas de calcul adaptés dans le code APOLLO3®, développé au CEA.
Après quelques rappels théoriques et méthodologiques, ce document présente une analyse critique des schémas de calcul disponibles dans APOLLO3® pour les réacteurs refroidis au sodium. Cette analyse permet de mettre en évidence la nécessité de simuler, dès létape de préparation des sections efficaces, des modes angulaires du flux qui soient représentatifs de la configuration géométrique du cur.
Pour répondre à ce besoin dans le cadre de géométries présentant une forte hétérogénéité axiale, une approximation 2D/1D à léquation du transport des neutrons 3D est développée.Cette dernière permet de représenter de manière cohérente, et à moindre coût, des effets danisotropie axiale dans des calculs 2D. Une nouvelle modélisation de type traverse de linterface cur / réflecteur est également proposée, ainsi quune méthode de calcul innovante
des barres de contrôle.
Ces méthodes permettent, in fine, de définir un schéma de calcul de référence unique et validé numériquement, permettant de modéliser un large spectre de curs de réacteurs refroidis au sodium, dans APOLLO3® . De façon à élargir la portée de ce travail, le dernier chapitre présente une réflexion sur la notion dadaptabilité du schéma.
Thesis resume
Sodium-cooled nuclear reactors offer interesting perspectives in terms of uranium resources economy and radioactive waste management. In France, the research on the SFR technology took a new start at the beginning of the XXIst century with the ASTRID project.
Simultaneously, the rise of safety standards in the nuclear industry resulted in increasingly complex core concepts, with enhanced natural behavior in accidental situations. In particular, the minimization of the sodium-void reactivity worth led to the CFV concept for ASTRID.
Due to increased levels of spatial heterogeneity, these innovative cores challenge classical neutronic calculation strategies. Hence, the first objective of this thesis is the identification of the main physical phenomena that need to be taken into account when modeling the neutronic behavior of a heterogeneous nuclear core in fast neutron spectrum. The second objective is the development of appropriate calculation schemes in the APOLLO3® code, developed at CEA.
After a brief reminder of neutronic calculation theory and methods, this document presents a critical analysis of the neutronic calculation schemes available in APOLLO3® for sodium-cooled applications. This analysis highlights the necessity to model, during the cross section preparation phase, angular modes of the neutron flux that are representative of the core geometrical configuration.
To meet this need in axially heterogeneous geometries, a 2D/1D approximation to the 3D neutron transport equation is derived and implemented in APOLLO3® . In particular, it is shown that this approximation allows to consistently represent axial angular modes of the flux in 2D calculation domains. Besides, a new traverse model is proposed for the core / reflector radial interface, as well as an innovative control rod calculation method.
The combination of these methods allows to define a unique, and numerically validated, reference calculation scheme in APOLLO3®, suitable for the calculation of a wide range of complex sodium-cooled nuclear cores. The discussion is finally enlarged with a reflection on the adaptability concept for neutronic calculation scheme.