Etude numérique et expérimentale de la fissuration des combustibles nucléaires oxydes et de la séparation de l’interface oxyde-gaine
Le CEA mène des études et expertises sur les combustibles nucléaires céramiques à base de dioxyde d’uranium (UO2). Des schémas numériques avancés pour la simulation prédictive du comportement de ces combustibles sont développés et s’appuient sur une démarche d’amélioration continue des modèles et des lois relatives aux propriétés physiques des matériaux. Les combustibles nucléaires sont des céramiques poreuses, dont la microstructure dépend de leur procédé de fabrication. Afin de garantir le confinement des produits de fission et le bon refroidissement du combustible, des gaines en alliage métallique sont disposées autour des combustibles.
Au cours de leur utilisation en réacteur, on observe une interaction mécanique entre le combustible et la gaine et potentiellement un accrochage, ce qui favorise le refroidissement du combustible en réduisant la résistance thermique de l’interface. Des fissures peuvent également apparaitre au sein du combustible créant, selon leur orientation, des barrières thermiques ou des chemins d’écoulement privilégiés pour les produits de fission créés lors de l’irradiation. Ces gaz vont ensuite exercer une pression sur les faces intérieures de la gaine ce qui peut, dans certaines conditions, entrainer une séparation/un décollement de l’interface combustible gaine et potentiellement l’apparition d’une lame de gaz entre le combustible et la gaine ce qui aurait un impact à la fois sur la thermique, et sur le comportement des produits de fission volatils.
Le but de cette thèse est donc de développer une démarche de caractérisation et modélisation de l’interface combustible-gaine fondée sur l’expérience, permettant notamment de prédire les conditions pouvant conduire à son décollement. Le travail consistera tout d’abord en la simulation numérique d’essais de décollement d’interface existants, ainsi que de l’effet d’une pression de gaz à l’intérieur du combustible, en présence ou non de fissures au sein du combustible. Il est envisagé de décrire numériquement l’interface par un modèle de zone cohésive, régi par la contrainte maximale, ainsi que par l’énergie dissipée pour la création de fissure. Les simulations permettront de comprendre le rôle des chacun de ces paramètres dans le décollement de la gaine. En fonction de ces résultats, il est envisagé de concevoir de nouveaux essais pour identifier les paramètres des lois de comportement mises en œuvre.
Ce travail sera basé au Département d’Etude des Combustibles de l’Institut IRESNE (CEA-Cadarache) pour les aspects spécifiques au comportement du combustible, et mené en étroite collaboration avec le Laboratoire de Mécanique Paris-Saclay pour l’étude des interfaces. Ainsi la personne travaillant sur ce sujet bénéficiera d’un environnement scientifique riche et stimulant et aura, en outre, la responsabilité de proposer, développer, réaliser et interpréter des expériences sur combustible nucléaire. Les compétences et connaissances acquises par le candidat seront valorisées à travers la rédaction de publications dans des journaux scientifiques internationaux et la présentation de ses travaux dans des conférences internationales.
Influence du dopage au chrome du combustible UO2 sur la spéciation des produits de fission en conditions accidentelles
Le développement des réacteurs nucléaires s’inscrit dans une démarche d’amélioration de la sûreté, avec par exemple le déploiement de combustibles nucléaires à propriétés améliorées vis-à-vis de leur comportement en conditions accidentelles : les combustibles nucléaires dits E-ATF (Enhanced Accident Tolerant Fuel). Parmi les combustibles E-ATF envisagés, le combustible UO2 dopé avec Cr2O3 est développé par l’opérateur industriel FRAMATOME. En revanche, très peu de données existent sur le comportement des produits de fission d’un combustible dopé Cr en conditions accidentelles.
La thèse propose de mettre au point un procédé de synthèse d’un combustible UO2 dopé Cr simulant le combustible irradié pour étudier le comportement des éléments (Cr et produits de fission) en température et sous différentes pressions partielles d’oxygène. La méthodologie repose sur une approche expérimentale couplant synthèse de matériaux modèles et caractérisation chimique approfondie, complétée par une approche théorique (calculs thermodynamiques) permettant de dimensionner les séquences thermiques et conforter les mécanismes réactionnels proposés.
La thèse sera réalisée au CEA de Cadarache (France), au sein de l’IRESNE (Institut de recherche sur les systèmes nucléaires pour la production d'énergie bas carbone). Le(La) doctorant(e) sera accueilli(e) dans un laboratoire dédié à l’étude des composés à base d’uranium du Département d’étude des combustibles (DEC). Selon les procédés de densification choisis, des expériences de plus ou moins longue durée pourront être menées dans d’autres laboratoires en France ou en Europe.
La thèse sera réalisée au CEA de Cadarache (France), au sein de l’IRESNE (Institut de recherche sur les systèmes nucléaires pour la production d'énergie bas carbone). Le(La) doctorant(e) sera accueilli(e) dans un laboratoire dédié à l’étude des composés à base d’uranium du Département d’étude des combustibles (DEC). Selon les procédés de densification choisis, des expériences de plus ou moins longue durée pourront être menées dans d’autres laboratoires en France ou en Europe.
Le doctorant aura l’opportunité de se former à des techniques pointues de caractérisation des sciences des matériaux céramiques, d’accéder à des expériences sur grands instruments (synchrotron) et de participer à des échanges avec le monde académique (CNRS, Universités, JRC). Il pourra valoriser ses travaux à travers des publications et des participations à congrès.
A l’issue de cette thèse, le doctorant aura acquis des compétences en science des matériaux et en caractérisation du solide qu’il pourra mettre à profit dans différents domaines des matériaux, ainsi qu’une expérience dans le milieu nucléaire d’intérêt pour l’industrie nucléaire.
Implémentation d’algorithmes parallèles sur GPU pour les simulations du combustible nucléaire sur supercalculateurs exaflopiques
Dans un contexte où le calcul haute performance (HPC) est en perpétuelle évolution, le design des nouveaux supercalculateurs tend à intégrer toujours plus d’accélérateurs ou de cartes graphiques (GPUs), qui fournissent l’essentiel de la puissance de calcul de la plupart des supercalculateurs actuels. En raison de leurs différences architecturales par rapport aux unités centrales de calcul (CPUs) et des environnements logiciels en constante évolution, les GPUs posent de profonds défis de programmation. Une utilisation efficace de leur puissance de calcul demande une refonte des algorithmes et logiciels de simulation existants pour atteindre un parallélisme massif.
Le CEA a développé la plateforme de calcul PLEIADES dédiée à la simulation du comportement des combustibles, depuis la fabrication jusqu’au comportement en réacteur, puis lors du stockage. Elle inclut une parallélisation en mémoire distribuée MPI permettant une parallélisation sur plusieurs centaines de cœurs. Cette plateforme répond aux exigences des partenaires du CEA que sont EDF et Framatome, mais il est nécessaire de l’adapter pour les nouvelles infrastructures de calcul. Proposer une solution flexible, portable et performante pour les simulations sur des supercalculateurs équipés de GPUs est d'un intérêt majeur afin de capturer des physiques toujours plus complexes sur des simulations comportant des domaines de calcul toujours plus grands.
Dans ce cadre, la thèse visera d’élaborer puis évaluer différentes stratégies de portage de noyaux de calculs sur GPU ainsi que l’utilisation de méthodes de répartition dynamique de la charge adaptés aux supercalculateurs GPUs actuels et futurs. Le candidat s’appuiera sur des outils développés au CEA comme les solveurs thermo-mécaniques MFEM-MGIS [1,2] ou MANTA [3]. Les solutions logicielles et algorithmes parallèles proposés avec cette thèse permettront à terme la réalisation de grands calculs de modélisation multi-physique en 3D du comportement des crayons combustibles sur des supercalculateurs comportant des milliers de cœurs de calcul et des GPUs.
Le candidat travaillera au sein du Laboratoire de développement des Outils de Calcul Scientifique (OCS) combustibles PLEIADES (LDOP) au département d'études des combustibles (DEC - Institut IRESNE, CEA Cadarache). Il sera amené à évoluer dans une équipe pluridisciplinaire composée de mathématiciens, physiciens, mécaniciens et informaticiens. A terme, les contributions de la thèse visent à enrichir la plate-forme numérique pour la simulation de combustibles nucléaires PLEIADES.
Références : [1] MFEM-MGIS - https://thelfer.github.io/mfem-mgis/; [2] Th. Helfer, G. Latu. « MFEM-MGIS-MFRONT, a HPC mini-application targeting nonlinear thermo-mechanical simulations of nuclear fuels at mesoscale ». IAEA Technical Meeting on the Development and Application of Open-Source Modelling and Simulation Tools for Nuclear Reactors, June 2022, https://conferences.iaea.org/event/247/contributions/20551/attachments/10969/16119/Abstract_Latu.docx, https://conferences.iaea.org/event/247/contributions/20551/attachments/10969/19938/Latu_G_ONCORE.pdf; [3] O. Jamond et al. «MANTA : un code HPC généraliste pour la simulation de problèmes complexes en mécanique », https://hal.science/hal-03688160
Essais thermomécaniques jusqu’à de hautes températures sur une céramique nucléaire irradiée
Ces travaux de thèse s’inscrivent dans le cadre des études sur les interactions pastille-gaine dans les crayons combustibles utilisés dans les réacteurs des centrales nucléaires à eau pressurisé. L’exploitant doit assurer et démontrer l’intégrité de l’ensemble des crayons, en toutes situations. Or, les contraintes mécaniques subies par la gaine, première barrière de sûreté, sont reliés aux propriétés viscoplastiques du combustible. Il faut donc connaître ces propriétés et leurs évolutions en fonctionnement.
Le sujet de thèse proposé s’intéresse particulièrement à la caractérisation en cellule de haute activité de combustibles irradiés en réacteur nucléaire. Une des difficultés majeures réside dans le fait que les combustibles issus d’une irradiation en réacteur sont multi-fissurés, ce qui rend leurs caractérisations mécaniques particulièrement compliquées. Toutefois, une thèse en cours (2022-25) a permis de (i) concevoir une machine spécifique d’essais thermomécaniques, (ii) qualifier en partie ce dispositif, (iii) mettre en place des outils et méthode d’extraction d’échantillon fissuré, (iv) ainsi qu’une modélisation intégrale du système (jumeau numérique).
La thèse proposée ici correspond à la suite de ces travaux et sera construite en quatre étapes sur trois plateformes expérimentales disponibles au CEA
1. Prise en main et amélioration des outils numériques et expérimentaux existant,
2. Implémentation du dispositif en cellule de haute activité sur un four existant,
3. Réalisation d’essai thermomécanique sur combustible irradié ; ce qui constituera en soi une première mondiale dans ces conditions.
Les essais demanderont des post-traitements dédiés basés sur une discussion essai-simulation. En premier lieu, cela permettra de décolérer les phénomènes et de se focaliser sur le comportement viscoplastique du matériau. Une fois la base expérimentale suffisamment étoffée et interprétée, il sera alors possible de conforter ou revoir les lois de comportement du combustible. Un lien avec la microstructure des matériaux pourrait être abordé.
Tout au long de ces étapes, le thésard s’appuiera sur les compétences et expertises de différents laboratoires du Département d'études des Combustible (Institut IRESNE, CEA Cadarache) et sur un environnement collaboratif académique. Cette thèse s'inscrit également dans le cadre du projet européen OPERA HPC et en est un enjeu majeur.
Le doctorant devra présenter un goût marqué pour l’approche expérimentale et quelques facilités pour l’utilisation d’outils numériques. Des connaissances en science des matériaux sont le minimum requis. Au long des trois années de doctorat, le doctorant améliorera ses compétences multiphysiques en conception de dispositifs expérimentaux et comportement des matériaux à haute température, ainsi qu’en simulation numérique, ce qui facilitera son insertion professionnelle. Au regard de l’ampleur de la tâche, il serait préférable, mais non obligatoire, que le candidat ait suivi préalablement le stage de master lié à ce sujet (voir par ailleurs).
Modélisation thermo-chemo-mécanique du frittage : effet de l’atmosphère et de la densification différentielle sur le retrait des pastilles
Les combustibles de dioxyde d’uranium (UO2), utilisés dans les centrales nucléaires sont des céramiques, dont le frittage en phase solide est une étape-clé de la fabrication. L’étape de frittage consiste en un traitement thermique sous pression partielle contrôlée de O2 permettant de consolider, densifier le matériau et faire grossir les grains de UO2. La densification induit un retrait macroscopique de la pastille. Si le compact (poudre comprimée par pressage avant le frittage) admet de fortes hétérogénéités de densité, une différence de densification dans la pastille peut avoir lieu entraînant un retrait différentiel et l’apparition de défauts.
Cette thèse se consacre à la mise en place d’un modèle thermo-chemo-mécanique du frittage pour simuler l’impact de la composition et les propriétés physiques de l’atmosphère sur la densification du combustible à l’échelle de la pastille. Cette échelle nous permettra de considérer les gradients de densité issus du pressage, mais également de prendre en compte la cinétique de diffusion d’oxygène impactant localement la vitesse de densification qui elle-même impactera le processus de transport. Une simulation multiphysique est nécessaire pour simuler le couplage de ces phénomènes.
Ce travail de thèse sera mené au sein du Laboratoire commun MISTRAL (Aix-Marseille Université/CNRS/Centrale Marseille et l'institut IRESNE du CEA-Cadarache). Le doctorant valorisera ses résultats au travers de publications et participations à des congrès et aura acquis de solides compétences qui sont recherchées et valorisables dans un grand nombre de domaines académiques et industriels.
Simulation de l'évolution des microstructures de dislocations dans UO2 : impact de la montée des dislocations à haute température
La neutralité carbone passe par le développement de systèmes de production d’énergie bas carbone incluant le nucléaire. L’analyse de sûreté du fonctionnement des réacteurs nucléaires porte sur le confinement des produits de fission dans toutes les situations de fonctionnement, avec notamment l’intégrité de la première barrière composée des éléments combustibles. Pour les concepts de type crayon, constitués d’un empilement de pastilles combustibles dans une gaine métallique, le comportement mécanique des pastilles en dioxyde d'Uranium (UO2) joue un rôle important dans l'évaluation de l'intégrité de la gaine. Ainsi, en situation de transitoire de puissance, le contact combustible-gaine accroît les sollicitations mécaniques de la gaine et le fluage du combustible peut permettre une accommodation des déformations de gonflement réduisant ainsi les contraintes appliquées à la gaine. Un des enjeux porte sur la compréhension et la prédiction de ce phénomène de fluage de l’UO2 avec les mécanismes qui le pilotent à l’échelle microstructurale polycristalline, notamment impliquant les dislocations.
L’objectif de la thèse sera de construire une méthode de simulation indispensable à l’acquisition de résultats de référence en support à la modélisation multi-échelle du comportement mécanique du combustible à haute température fortement dépendant du phénomène de montée des dislocations. Ce type de démarche de simulation et les résultats qui seront obtenus seront particulièrement novateurs et n’ont encore jamais été mis en œuvre dans le cas des combustibles oxyde pour lesquels l’évolution de la microstructure de dislocation à également un impact fort sur le comportement des produits de fission gazeux en plus des aspects mécaniques étudiés dans la thèse. Pour cela le doctorant développera un schéma de calcul, basée sur le couplage entre un code de dynamique des dislocations (NUMODIS) et un code de résolution des équations aux dérivées partielles non linéaires par FFT (AMITEX-FFTP). Ceci permettra de décrire l’évolution d’une microstructure de dislocations sous l’effet de la montée des dislocations (NUMODIS) induite par la diffusion des lacunes (AMITEX-FFTP). Ensuite, des simulations basées sur cette approche permettront de quantifier les phénomènes de restauration de la densité des dislocations stockées avec l’effet des mécanismes de montée dans différentes configurations (températures, contraintes…). Ce travail permettra in fine d’améliorer et valider la modélisation micromécanique existante et mise en œuvre dans la plateforme de simulation PLEIADES du CEA.
Cette thèse sera réalisée dans le cadre d’un co-encadrement entre le Département d'Etude des Combustibles (Institut IRESNE, CEA Cadarache) et Le Département de recherche sur les Matériaux et la Physico-chimie (Institut ISAS, CEA Saclay), et d’une collaboration avec l’IM2NP d’Aix Marseille Université. Les travaux de thèse seront menés au sein des laboratoires LM2C (Cadarache) et LC2M (Saclay) dans un environnement donnant accès à une grande expertise sur la modélisation multiéchelle des matériaux. Les travaux de recherche seront valorisés par des publications et des participations à des conférences internationales dans le domaine des matériaux.
Solveur Intégrodifférentiel HPC Parallèle pour la Dynamique des Dislocations
Contexte : La compréhension du comportement des métaux à forts taux de déformation [4] (entre 104 et 108 s-1) représente un défi scientifique et technologique considérable. Cette déformation irréversible (plastique) est due à la présence de défauts linéaires d'alignement cristallin : les dislocations, qui interagissent via le champ élastique à longue portée et par interactions de contact.
Actuellement, le comportement des métaux à forts taux de déformation ne sont accessibles expérimentalement que par chocs laser. D’où la nécessité d’un outil de simulation. Deux grands types d’approches sont possibles : la dynamique moléculaire, et les simulations élastodynamiques. Cette thèse s’inscrit dans le second type d’approche, capitalisant sur nos travaux récents [1, 2] qui ont permis les premières simulations numériques de l’équation de Peierls-Nabarro Dynamique (PND) [5]. Celle-ci décrit des phénomènes intervenant à l’échelle de la dislocation.
PND est une équation intégrodifférentielle non-linéaire qui présente une double difficulté : la non-localité en temps et en espace des opérateurs. Nous l’avons simulée pour la première fois grâce à une stratégie numérique efficace [1], issue de [6]. Mais la nature mono-processeur de son implantation actuelle constitue un verrou, limitant fortement la taille du système et l’étude de son comportement en temps long.
Sujet de thèse : Les objectifs de cette thèse sont de deux natures :
- Numérique. Sur la base algorithmique développée dans [1], implémenter un solveur HPC (Calcul Haute Performance) parallélisé en espace et en temps, avec mémoire distribuée.
- Physique. Grâce au code développé, éclaircir des points cruciaux relatifs à la phénoménologie des dislocations en régime dynamique rapide. L’exploitation des résultats numériques requerra des techniques de traitement de données et de statistiques - potentiellement assistées par de l’IA.
En fonction de l’avancement, il sera possible d’appliquer la méthode numérique développée aux phénomènes de fissuration dynamique [3].
Profil du candidat : Le sujet de thèse proposé est pluridisciplinaire, à la croisée des chemins entre simulation numérique, physique des dislocations et de la propagation de fissures, et traitement statistique. Le candidat devra d’abord posséder une solide formation en calcul scientifique appliqué aux équations aux dérivées partielles et un gout prononcé pour les applications physiques. La maîtrise du C++, avec des compétences en OpenMP et MPI seraient fortement appréciées. Des connaissances en mécanique des milieux continus seraient aussi vu comme un plus.
La thèse se déroulera au Département d'Etudes des Combustibles (Institut IRESNE, CEA/DES, centre de Cadarache), avec des déplacements réguliers en région parisienne pour la collaboration avec le CEA/DAM et le CEA/DRF.
[1] Pellegrini, Josien, Shock-driven motion and self-organization of dislocations in the dynamical Peierls model, soumis.
[2] Josien, Etude mathématique et numérique de quelques modèles multi-échelles issus de la mécanique des matériaux. Thèse. (2018).
[3] Geubelle, Rice. J. of the Mech. and Phys. of Sol., 43(11), 1791-1824. (1995).
[4] Remington et coll., Metall. Mat. Trans. A 35, 2587 (2004).
[5] Pellegrini, Phys. Rev. B, 81, 2, 024101, (2010).
[6] Lubich & Schädle. SIAM J. on Sci. Comp. 24(1), 161-182. (2002).
Dynamique multiéchelle d’une structure élancée avec singularités de frottement : application à un assemblage de combustible
La modélisation dynamique de structures complexes peut nécessiter la prise en compte de phénomènes intervenant à des échelles très différentes. Or, une modélisation fine de ce type de structures entraîne généralement des coûts de calculs prohibitifs. La modélisation multiéchelle se présente alors comme une solution alternative à cette problématique en tenant compte de chaque phénomène à l’échelle la plus adéquate.
Nous nous intéressons ici à des structures élancées soumises à des sollicitations mécaniques qui conduisent à des contacts frottants entre la structure et les éléments de maintien. Le comportement des structures élancées est en général représenté par des modèles de type poutre équivalente, mais la prise en compte précise du contact/frottement local nécessite des modèles 3D massifs.
L’originalité du travail proposé ici est de bâtir une approche multiéchelle et multimodèle efficace entre modèles poutres et massifs qui permette de prendre en compte localement le contact frottant de structures élancées. Nous nous orientons ainsi vers l’utilisation de méthodes multigrilles (ou multiniveaux) locales qui permettent naturellement un couplage multiéchelle non intrusif. La précision de ces méthodes repose alors sur le choix des opérateurs de transfert entre échelles, qui devront être définis avec soin. Il faudra également prendre en compte la non compatibilité des maillages soutenant les modèles sur les différentes échelles pertinentes. Ainsi, le modèle final sera un modèle de poutre enrichi permettant de prendre en compte des phénomènes de contact locaux.
Le modèle développé sera comparé à des résultats expérimentaux obtenus lors de campagnes d’essais déjà réalisées, et à des solutions numériques de référence, de complexité croissante, destinées à valider finement la pertinence de l’approche multiéchelle proposée.
Le potentiel fort des approches multiéchelles visées, appliqué dans ce sujet au domaine du nucléaire, pourra être valorisé par le candidat à d’autres problématiques industrielles telles que celles de l’aéronautique ou encore de l’automobile.
La thèse se déroulera dans le cadre du laboratoire commun MISTRAL entre le CEA et le LMA (laboratoire de mécanique et d’acoustique) de Marseille. Le doctorant réalisera la majeure partie de sa thèse au CEA au sein de l'institut IRESNE (Cadarache) dans les équipes spécialisées autour des méthodes numériques et de la modélisation dynamique de structures complexes. Il se rendra régulièrement à Marseille pour échanger avec les encadrants universitaires.
De l’Angström au micron : un modèle d’évolution microstructurale du combustible nucléaire dont les paramètres sont calculés à l’échelle atomique
La maîtrise du comportement des gaz de fission dans le combustible nucléaire (oxyde d’uranium) est un enjeu industriel important puisque leur relâchement ou leur précipitation limite l'utilisation du combustible à forts taux de combustion. Or ces phénomènes sont fortement influencés par l’évolution microstructurale du matériau aux défauts générés par l’irradiation (création de défauts ponctuels, agrégations de ceux-ci en cavités et bulles de gaz ou en boucles ou lignes de dislocation…). La dynamique d’amas (DA) est un modèle de type cinétique chimique permettant de décrire la nucléation/croissance des amas de défauts, leur contenu en gaz et le relâchement de celui-ci. Le modèle utilisé est paramétré à partir de données de base calculées à diverses échelles (ab initio, potentiels empiriques, Monte Carlo). Ce modèle rend déjà compte d’expériences de recuit d’UO2 implanté en atomes de gaz de fission et a mis en évidence le fort impact des défauts d’irradiation sur le relâchement gazeux. L’objectif de la thèse est d’une part d’améliorer le modèle et ses paramètres d’entrée, notamment le taux de création de défauts d’irradiation, et d’autre part d’étendre son domaine de validation en le confrontant à de nombreuses expériences issues de thèses récemment soutenues au département (mesure de relâchement gazeux par recuit d’échantillons implantés via un accélérateur d’ions, observation de cavités, bulles de gaz et boucles de dislocation par microscopie électronique à transmission, caractérisation du dommage par spectrométrie d’annihilation de positons). Le candidat sera donc amené à faire évoluer certains des sous-modèles constitutifs de la DA, interpréter et simuler l’ensemble des expériences disponibles. En parallèle cela permettra d’affiner la paramétrisation du modèle.
Ce sujet de modélisation présente l’intérêt pour le candidat d’articuler à une dimension “théorique” (amélioration du modèle), ainsi que de physique numérique (simulation en Dynamique Moléculaire de cascades de déplacements) une dimension “expérimentale” (interprétation d’expériences déjà réalisées, voire conception et suivi de nouvelles expériences). Ainsi, l’approche d’un ensemble varié de techniques d’observation et de mesure ouvriront au candidat le monde de la physique expérimentale et complèteront son profil. Le candidat aura également à animer des collaborations dans le but d’analyser les données expérimentales, de développer l’outil de calcul ou de spécifier des calculs atomistiques complémentaires. Il pourra aussi bénéficier d’un environnement de collaboration académique.
Ce travail offre une position centrale et un point de vue synthétique sur la physique du combustible en irradiation. Il vous permettra de contribuer au développement de la physique numérique appliquée à une démarche multiéchelle de modélisation. Vous découvrirez en quoi des outils de simulation basés sur les données microscopiques les plus fondamentales obtenues par le calcul atomistique permettent de traiter et expliquer des situations pratiques.
Pour aller plus loin :
Skorek (2013). Étude par Dynamique d’Amas de l’influence des défauts d’irradiation sur la migration des gaz de fission dans le dioxyde d’uranium. Univ. Aix-Marseille. http://www.theses.fr/2013AIXM4376
Bertolus et al. (2015). Linking atomic and mesoscopic scales for the modelling of the transport properties of uranium dioxide under irradiation. Journal of Nuclear Materials, 462, 475–495.
Génération assistée de noyaux de calculs complexes en mécanique du solide
Les lois de comportement utilisées dans les simulations numériques décrivent les caractéristiques physiques des matériaux simulés. À mesure que notre compréhension de ces matériaux évolue, la complexité de ces lois augmente.L'intégration de ces lois constitue une étape critique pour la performance et la robustesse des calculs scientifiques. De ce fait, cette étape peut conduire à des développements intrusifs et complexes dans le code.
De nombreuses plateformes numériques telles que FEniCS, FireDrake, FreeFEM, Comsol, proposent des techniques de génération de code à la volée (JIT, pour Just In Time) pour gérer différentes physiques. Cette approche JIT réduit considérablement les temps de mise en oeuvre de nouvelles simulations, offrant ainsi une grande versatilité à l'utilisateur. De plus, elle permet une optimisation spécifique aux cas traités et facilite le portage sur diverses architectures (CPU ou GPU). Enfin, cette approche permet de masquer les détails d'implémentation: une évolution de ces derniers est invisible pour l'utilisateur et est absorbée par la couche de génération de code.
Cependant, ces techniques sont généralement limitées aux étapes d'assemblage des systèmes linéaires à résoudre et n'incluent pas l'étape cruciale d'intégration des lois de comportement.
S'inspirant de l'expérience réussie du projet open-source mgis.fenics [1], cette thèse vise à développer une solution de génération de code à la volée dédiée au code de mécanique des structures de nouvelle génération Manta [2] développé par le CEA. L'objectif est de permettre un couplage fort avec les lois de comportement générées par MFront [3], améliorant ainsi la flexibilité, les performances et la robustesse des simulations numériques.
Le doctorant bénéficiera d'un encadrement de la part des développeurs des codes MFront et Manta (CEA), ainsi que des développeurs du code A-Set (collaboration entre Mines-Paris Tech, Onera, et Safran). Cette collaboration au sein d'une équipe multidisciplinaire offrira un environnement stimulant et enrichissant pour le candidat.
De plus, le travail de thèse sera valorisé par la possibilité de participer à des conférences et de publier des articles dans des revues scientifiques à comité de lecture, offrant une visibilité nationale et internationale aux résultats de la thèse.
Le doctorat se déroulera au CEA Cadarache, dans le sud est de la France, au sein du département d'études des combustibles nucléaires de l'institut IRESNE [4]. Le laboratoire d'accueil est le LMPC dont le rôle est de contribuer au développement des composants physiques de la plateforme numérique PLEIADES [5], co-développée par le CEA et EDF.
[1] https://thelfer.github.io/mgis/web/mgis_fenics.html
[2] MANTA : un code HPC généraliste pour la simulation de problèmes complexes en mécanique. https://hal.science/hal-03688160
[3] https://thelfer.github.io/tfel/web/index.html
[4] https://www.cea.fr/energies/iresne/Pages/Accueil.aspx
[5] PLEIADES: A numerical framework dedicated to the multiphysics and multiscale nuclear fuel behavior simulation https://www.sciencedirect.com/science/article/pii/S0306454924002408