Laboratoire d'accueil : Laboratoire de dosimétrie biologique (LDB)

Date de début de thèse : 01/10/2009

Nom du doctorant : Ingrid Nosel

Contexte et problématique

Dans les centrales êtectronucléaires actuelles, le combustible se présente sous la forme de pastilles cylindriques d'environ 1 cm de diamètre, empilées dans des tubes métalliques minces. Avant leur mise en réacteur (état dit "vierge"), ces pastilles sont assimilables à un milieu continu - par opposition à milieu discret - constitué par une céramique poreuse et leur empilement peut se concevoir comme une colonne de gros fragments cylindriques. Dès les premiers temps de séjour en réacteur, les pastilles de combustibles se fragmentent sous l'action des différentes sollicitations qu'elles subissent: gradients thermiques et mécaniques, réactions chimiques, irradiation, etc...
L'empilement de pastilles s'approche alors progressivement d'un milieu discret dont la morphologie est complexe.
Que ce soit en fonctionnement normal ou en situations accidentelles, la connaissance du comportement de ce milieu discret est essentielle aux analyses de sûreté. Il conditionne en effet, en cas de rupture du gainage, les potentielles interactions du combustible avec le milieu environnant (interaction thermomécanique avec le réfrigérant, dispersion des produits radioactifs portés par les fragments). Ce comportement va fortement dépendre de la "mésostructure" de l'empilement, c'est à dire de géométrie des fragments, des rapports de taille entre fragments, de leur comportement (lié à sa microstructure), de leurs interactions et du réseau d'effort entre ces fragments (la texture).
Une thèse est actuellement en cours pour comprendre et modéliser le comportement de ce milieu discret en situation accidentelle (Thèse A. MARTIN, conjointe DPAM-LMGC, cofinancée IRSN-CNRS, débutée en 2006). Cette thèse se concentre sur l'écoulement et la stabilité de cette colonne combustible déjà fragmentée et pour ce faire des hypothèses sont formulées a priori sur la mésostructure du milieu discret étudié.
Le travail proposé ici se situe en amont de l'analyse de l'écoulement: il concerne la génération en elle-même de ce milieu discret. La problématique réside donc dans la recherche d'une statistique descriptive pour une mésostructure obtenue par fragmentation d'Un milieu semi·continu sous l'action de gradients thermo-mécaniques. Par semi-continu on entend ici un milieu constitué de gros fragments (les pastilles) qui donneront lieu à des fragments de taille bien inférieure.
A l'heure actuelle, seules des simulations numériques complexes peuvent prétendre apporter une réponse satisfaisante et quantitative à cette problématique. Deux classes de méthodes récentes permettent t'analyse de la multifissuration dynamique sans régularisation et avec une prise en compte de la physique locale: la méthode aux éléments discrets IDEM - Discrete Element Method et la méthode aux éléments de Réseau (LEM - Lattice Element Method).
La méthode DEM consiste à considérer chaque maille ou groupe de mailles d'une discrétisation comme un corps indépendant relié aux autres par des conditions aux limites mixtes. Cette méthode a montré ses potentialités dans la multi·fissuration des matériaux à microstructure complexe tels que les matériaux à gradient de propriétés (Thèse F. Perales, conjointe DPAM-LMGC, cofinancée IRSN-EdF, soutenue en 2005). Elle s'avère cependant coûteuse sur le plan de la simulation numérique.
La méthode LEM consiste à considérer ta discrétisation du milieu sous la forme d'Un réseau de liens 1D. Cette méthode plus simple et plus rapide a été appliquée avec succès dans le cadre de la rupture de l'atbumen de blé (Thèse V. Topin, conjointe CNRS·INRA, cofinancée INRA-Région Languedoc-Roussiton, soutenue en 2008).
Dans cette étude les rôles des interfaces cohésives entre granules d'amidon et matrice protéique et celles avec le réseau des pores ont été mis en évidence par rapport au comportement à la mouture. Un avantage majeur de cette méthode est de permettre la prise en compte simultanée d'une phase discrète et d'une phase continue dans un cadre unique. Le cadre relativement ouvert de cette méthode permet de prendre en compte naturellement les défauts initiaux (pores, fissures). Cette situation n'est qualitativement pas très éloignée de celte de la fragmentation du combustible nucléaire UO2.