Développements et validation de calculs à énergie continue pondérés par l'importance

par Guillaume Truchet

Thèse de doctorat en Mécanique des fluides, procédés, énergétique

Sous la direction de Alain Santamarina et de Pierre Leconte.

Le président du jury était Daniel Heuer.

Le jury était composé de Alain Santamarina, Pierre Leconte, Joachim Miss, Laurence Villatte, Fausto Malvagi.

Les rapporteurs étaient Paul Reuss, Wilfred Van Rooijen.


  • Résumé

    L'un des enjeux actuel de la neutronique concerne la propagation rigoureuse des incertitudes d'entrée (e.g. données nucléaires, tolérances de fabrications, etc.) aux résultats finaux calculés par les codes (e.g. keff, taux de réaction, etc.). Pour propager les incertitudes, il est de coutume de faire l'hypothèse de petites variations autour d'une référence et de calculer, dans un premier temps, des profils de sensibilités. Or, les codes Monte-Carlo, qui se sont imposés comme des références de calcul, ne possèdent pas -- ou n'ont intégré que très récemment -- un moyen direct de calculer des sensibilités et donc de réaliser un calcul précis d'incertitudes. Les approches déterministes, elles, permettent le calcul de ces sensibilités mais introduisent parfois de très fortes hypothèses, notamment sur la géométrie.Le premier objectif de se travail de thèse est d'introduire dans le code Monte Carlo du CEA de transport des neutrons, TRIPOLI-4, des méthodes à même de calculer des profils de sensibilités du keff aux données nucléaires ou à toute autre perturbation. Pour cela, il a d'abord été nécessaire de mettre en place le calcul du flux adjoint d'un milieu critique. Pour la première fois, et grâce aux développements informatiques de ce travail, il a été possible de calculer dans un cas réel, concret, et applicatif, des spectres de flux adjoints en un point quelconque d'un réacteur. Ceci a été réalisé à l'aide de la probabilité itérée de fission (Iterated Fission Probability ou IFP) qui assimile le flux adjoint à l'importance d'un neutron dans un réacteur exactement critique. Ce calcul de flux adjoint a, par la suite, ouvert la porte au premier développement d'une méthode de calcul de « perturbations exacte » en Monte Carlo, théorie qui permet de s'affranchir des hypothèses de petites variations, et qui ouvre la porte à certaines applications jusqu'alors difficiles à analyser.Au delà de l'analyse poussée de la méthode IFP et de son application au calcul de flux adjoint, cette thèse propose également, d'obtenir dans le code TRIPOLI-4, les paramètres cinétiques d'un réacteur pondérés par le flux adjoint ou bien des aires de migration. A cette fin, l'implémentation reprend et améliore un algorithme déjà développé par la communauté scientifique pour estimer des perturbations au premier ordre.

  • Titre traduit

    Development and validation of continuous energy adjoint-weighted calculations


  • Résumé

    A key issue in nowadays Reactor Physics is to propagate input data uncertainties (e.g. nuclear data, manufacturing tolerances, etc.) to nuclear codes final results (e.g. keff, reaction rate, etc.). In order to propagate uncertainties, one typically assumes small variations around a reference and evaluates at first sensitivity profiles. Problem is that nuclear Monte Carlo codes are not -- or were not until very recently -- able to straightforwardly process such sensitivity profiles, even thought they are considered as reference codes.First goal of this PhD thesis is to implement a method to calculate keff-sensitivity profiles to nuclear data or any perturbations in TRIPOLI-4, the CEA Monte Carlo neutrons transport code. To achieve such a goal, a method has first been developed to calculate the adjoint flux using the Iterated Fission Probability (IFP) principle that states that the adjoint flux at a given phase space point is proportional to the neutron importance in a just critical core after several power iterations. Thanks to our developments, it has been made possible, for the fist time, to calculate the continuous adjoint flux for an actual and complete reactor core configuration. From that new feature, we have elaborated a new method able to forwardly apply the exact perturbation theory in Monte Carlo codes. Exact perturbation theory does not rely on small variations which makes possible to calculate very complex experiments.Finally and after a deep analysis of the IFP method, this PhD thesis also reproduces and improves an already used method to calculate adjoint weighted kinetic parameters as well as reference migrations areas.


Il est disponible au sein de la bibliothèque de l'établissement de soutenance.

Consulter en bibliothèque

La version de soutenance existe

Où se trouve cette thèse\u00a0?

  • Bibliothèque : Service Interétablissement de Documentation. LLSH Collections numériques.
  • Bibliothèque : Université Savoie Mont Blanc (Chambéry-Annecy). Service commun de la documentation et des bibliothèques universitaires. Bibliothèque électronique.
  • Bibliothèque : Service interétablissements de Documentation. STM. Collections numériques.
Voir dans le Sudoc, catalogue collectif des bibliothèques de l'enseignement supérieur et de la recherche.