Analyse des erreurs induites par une modélisation simplifiée sur l’évolution des combustibles REP Impact des fuites neutroniques dans les calculs cellules

Les études de scénarios d'un parc électronucléaire, ainsi que les études de sûreté, sont essentielles pour explorer les différentes stratégies du nucléaire du futur. Pour mener à bien ces études, il est nécessaire d'estimer le temps d'irradiation d'un combustible donné, ainsi que sa composition isotopique pendant la campagne de production d'électricité. Ces estimations reposent sur des simulations de réacteurs nucléaires, dont les calculs d'évolution doivent être les plus représentatifs possible. Les schémas de calcul classiquement utilisés s'effectuent en deux étapes : un calcul cellule pour résoudre l'équation du transport des neutrons (de type déterministe ou Monte Carlo) suivi d'un calcul cœur (déterministe). Le calcul cellule est une simulation d'évolution d'un assemblage dans des conditions infinies. À partir de ce calcul, des sections efficaces homogénéisées et condensées, ainsi que des grandeurs de diffusion, sont calculées comme données d'entrée pour l'étape suivante, celle du calcul cœur. Le calcul cellule est donc une étape fondamentale et celui-ci doit être le plus représentatif possible d'un assemblage du cœur. Or, les approximations à la base de ce modèle sont nombreuses, plus particulièrement les fuites neutroniques sont négligées. L'objectif de ce travail est d'étudier les effets physiques de fuites neutroniques et de quantifier les biais associés par rapport à une simulation infinie. Dans une première partie, la problématique des fuites neutroniques axiales est étudiée. Dans ce cas, les fuites de neutrons provoquent une variation forte du spectre neutronique localisée dans les derniers centimètres de l'assemblage ainsi qu'une variation plus faible mais globale sur l'ensemble de l'assemblage. Une deuxième partie est dédiée aux fuites radiales de neutrons. L'effet des assemblages voisins, ainsi que le comportement particulier des assemblages en position périphérique sont étudiés et les biais de composition en fin d'irradiation sont quantifiés. Un calcul d'évolution d'un réacteur très simplifié permet de visualiser, dans une dernière partie, l'ensemble des effets physiques observés et qui impactent l'évolution de l'irradiation. De nombreuses approximations du calcul cellule restent à explorer, comme le suivi de réactivité, par l'intermédiaire de la concentration du poison de neutrons thermiques solubilisé dans le modérateur ou présent dans le combustible. Cependant, la détermination des phénomènes physiques à prendre en compte pour le calcul cellule représente une première étape indispensable vers une amélioration de la représentativité du calcul cellule, voire conduire à des nouvelles méthodes de simulation d'un cœur du réacteur. À terme, les quantifications des biais liés aux fuites neutroniques serviront à estimer l'incertitude sur les compositions isotopiques du combustible en fin d'irradiation. Ces incertitudes, propagées dans les études de scénarios, permettront de quantifier le degré de validité des résultats obtenus.