Développement d'une méthode de Monte Carlo dépendante du temps et application au réacteur de type CANDU-6

La résolution de l'équation de Boltzmann demeure une étape importante dans la prédiction du comportement d'un réacteur nucléaire. Malheureusement, la résolution de cette équation présente toujours un défi pour une géométrie complexe (réacteur) tout comme pour une géométrie simple (cellule). Ainsi, pour prédire le comportement d'un réacteur nucléaire,un schéma de calcul à deux étapes est nécessaire. La première étape consiste à obtenir les paramètres nucléaires d'une cellule du réacteur après une étape d'homogénéisation et condensation. La deuxième étape consiste en un calcul de diffusion pour tout le réacteur en utilisant les résultats de la première étape tout en simplifiant la géométrie du réacteur à un ensemble de cellules homogènes le tout entouré de réflecteur. Lors des transitoires (accident), ces deux étapes sont insuffisantes pour pouvoir prédire le comportement du réacteur. Comme la résolution de l'équation de Boltzmann dans sa forme dépendante du temps présente toujours un défi de taille pour tous types de géométries,un autre schéma de calcul est nécessaire. Afin de contourner cette difficulté, l'hypothèse adiabatique est utilisée. Elle se concrétise en un schéma de calcul à quatre étapes. La première et deuxième étapes demeurent les mêmes pour des conditions nominales du réacteur. La troisième étape se résume à obtenir les nouvelles propriétés nucléaires de la cellule à la suite de la perturbation pour les utiliser, au niveau de la quatrième étape, dans un nouveau calcul de réacteur et obtenir l'effet de la perturbation sur le réacteur. Ce projet vise à vérifier cette hypothèse. Ainsi, un nouveau schéma de calcul a été défini. La première étape de ce projet a été de créer un nouveau logiciel capable de résoudre l'équation de Boltzmann dépendante du temps par la méthode stochastique Monte Carlo dans le but d'obtenir des sections efficaces qui évoluent dans le temps. Ce code a été utilisé pour simuler un accident LOCA dans un réacteur nucléaire de type CANDU-6. Les sections efficaces dépendantes du temps ont été par la suite utilisées dans un calcul de diffusion espace-temps pour un réacteur CANDU-6 subissant un accident de type LOCA affectant la moitié du coeur afin d'observer son comportement durant toutes les phases de la perturbation. Dans la phase de développement, nous avons choisi de démarrer avec le code OpenMC, développé au MIT,comme plateforme initiale de développement. L'introduction et le traitement des neutrons retardés durant la simulation ont présenté un grand défi à surmonter. Il est important de noter que le code développé utilisant la méthode Monte Carlo peut être utilisé à grande échelle pour la simulation de tous les types des réacteurs nucléaires si les supports informatiques sont disponibles.