Ajustement d'un modèle d’interaction nucléaire effective et propagation des erreurs statistiques
Au cœur de chaque approche « many-body » utilisée pour décrire les propriétés fondamentales d’un noyau atomique, on retrouve l’interaction effective nucléon-nucléon. Une telle interaction effective doit prendre en compte les effets du milieu nucléaire. Pour l’obtenir, il faut utiliser un protocole d’ajustement complexe qui prend en compte une variété d’observables nucléaires comme les rayons, les masses, les centroïdes des résonances géantes ou encore l’équation d’état de la matière nucléaire autour de la densité de saturation.
Un modèle d’interaction forte très utilisé est celui de Gogny, qui est formé par une combinaison linéaire des constantes de couplages et d’opérateurs avec un facteur de forme radial de type Gaussien [1]. Les constantes de couplages sont déterminées via un protocole d’ajustement sur les propriétés d’un nombre restreint de noyaux, typiquement les noyaux sphériques comme 40-48Ca, 56Ni, 120Sn et 208Pb.
L’objectif premier de cette thèse consiste à développer un protocole d’ajustement de l’interaction nucléaire qui puisse donner accès à la matrice de covariance des paramètres du modèle pour ensuite effectuer une analyse de la propagation des erreurs statistiques sur les observables nucléaires [2].
Après avoir analysé les relations entre paramètres et leurs poids relatifs sur les différentes observables, le doctorant explorera la possibilité de modifier certains termes de l’interaction comme le terme à trois corps ou les effets au-delà du champ moyen.
Le doctorant sera positionné dans une équipe de physiciens nucléaires au sein d’un laboratoire d’étude de physique de l'institut CEA IRESNE situé à Cadarache. Le travail s’effectuera en équipe avec le CEA/DIF. Les principaux débouchés professionnels sont la recherche académique et les organismes de R&D dans le domaine nucléaire.
[1] D. Davesne et al. "Infinite matter properties and zero-range limit of non-relativistic finite-range interactions." Annals of Physics 375 (2016): 288-312.
[2] T. Haverinen and M. Kortelainen. "Uncertainty propagation within the UNEDF models." Journal of Physics G: Nuclear and Particle Physics 44.4 (2017): 044008.
Analyse en ligne de simulants d'actinides en solution par LIBS et IA pour les procédés de retraitement du combustible
La construction de nouveaux réacteurs nucléaires dans les années à venir implique une augmentation des capacités de retraitement du combustible. Cette évolution requiert des développements scientifiques et technologiques pour mettre à jour notamment les équipements de surveillance du procédé. L’un des paramètres à suivre en continu est la teneur en actinides en solution, donnée essentielle au pilotage du procédé, actuellement mesurée par des technologies obsolètes. On se propose donc de développer la LIBS (laser-induced breakdown spectroscopy) pour cette application, une technique bien adaptée à l’analyse élémentaire quantitative en ligne. Les spectres des actinides étant particulièrement complexes, on souhaite recourir à des approches multivariées de traitement des données, comme certaines techniques d’intelligence artificielle (IA), pour extraire l’information quantitative des données LIBS et caractériser l’incertitude de mesure.
L’objectif de la thèse est donc d’évaluer les performances de l’analyse en ligne d’actinides en solution par LIBS et IA. On visera en particulier à améliorer la caractérisation des incertitudes à travers des méthodes de machine learning, et à les minimiser fortement pour répondre aux besoins de surveillance de l’usine de retraitement du futur.
Le travail expérimental sera réalisé sur des simulants non radioactifs des actinides, et au moyen d’un équipement LIBS commercial. Les données spectroscopiques alimenteront le volet de la thèse sur le traitement des données, et sur la détermination de l’incertitude obtenue par différents modèles de quantification.
Les résultats obtenus permettront de publier au moins 2 à 3 articles dans des revues à comité de lecture, voire de déposer des brevets. Les perspectives de la thèse sont la montée en maturité de la méthode et de l’instrumentation, pour aller progressivement vers une mise en œuvre sur une installation représentative d’un procédé de retraitement.
Modèles microscopiques de structure nucléaire pour étudier le processus de désexcitation dans la fission nucléaire
Le code FIFRELIN est développé au CEA/IRESNE Cadarache afin de fournir une description détaillée du processus de fission et de calculer avec précision toutes les observables de fission pertinentes. Le code repose en grande partie sur la connaissance détaillée de la structure sous-jacente des noyaux impliqués dans le processus de désexcitation post-fission. Dans la mesure du possible, le code s'appuie sur des bases de données de structures nucléaires telles que RIPL-3, qui fournissent des informations précieuses sur les schémas de niveaux nucléaires, les rapports de branchement et d'autres propriétés nucléaires essentielles. Malheureusement, toutes ces quantités n'ont pas été mesurées, des modèles nucléaires sont donc utilisés.
Le développement de modèles nucléaires avancés est la tâche du groupe de théorie nucléaire nouvellement formé à Cadarache, dont l'expertise principale est l'implémentation de solveurs du problème nucléaire à A corps basés sur des interactions nucléon-nucléon effectives.
Le but de cette thèse est de quantifier l'impact de la fonction de force E1/M1 et E2/M2 sur les observables de fission. Actuellement, cette quantité est principalement estimée à l'aide de modèles simples tels que la Lorentzienne généralisée. Le doctorant devra remplacer ces modèles par des théories entièrement microscopiques basées sur l'interaction effective entre les nucléons via les techniques de type QRPA. Une étude préliminaire a démontré que l'utilisation de modèles macroscopiques (Lorentzienne généralisé) ou microscopiques (QRPA) a un impact non négligeable sur les observables de fission.
Les débouchés de la thèse incluent la recherche académique et les labos de R&D nucléaire théorique et appliquée.
Compréhension et modélisation des mécanismes de découpe laser pour le démantèlement
Depuis plus de 30 ans, le Laboratoire des Technologies d'Assemblage (LTA) du CEA Saclay mène des recherches pour développer des outils innovants pour le démantèlement des installations nucléaires, en mettant au point des procédés de découpe laser adaptés à des environnements hostiles. Cette technologie permet de découper des matériaux épais, en air ou sous l’eau, et s’est révélée particulièrement efficace pour le démantèlement grâce à sa précision et à sa capacité à limiter la génération d’aérosols. Aujourd’hui, cette technologie est considérée comme sûre et éprouvée, notamment grâce au projet européen "LD-SAFE".
Cependant, des verrous techniques persistent, notamment la gestion de l’énergie laser résiduelle qui, en se propageant au-delà de la pièce découpée, peut endommager les structures environnantes. Des premiers travaux, parmi lesquels une thèse, ont permis de développer des modèles numériques pour prédire et maîtriser cette énergie, apportant des avancées significatives. En revanche, des défis technologiques demeurent, comme la gestion des fortes épaisseurs (>10 mm), la découpe de multi-plaques ou encore la prise en compte de l’ajout d’oxygène pour améliorer l'efficacité de la découpe.
L'objectif de la thèse sera de faire face à ces défis afin de mieux comprendre le procédé de découpe laser et la propagation de l'énergie laser résiduelle. Le doctorant perfectionnera un modèle numérique pour prédire son impact sur les structures en arrière-plan, notamment pour des matériaux épais et des configurations multi-plaques. Le travail inclura le développement d'un modèle multiphysique, validé par des expériences, avec un focus particulier sur l'effet de l'oxygène, la création de modèles simples et l'adaptation pour une utilisation par des opérateurs externes.
La thèse sera réalisée en collaboration entre le Laboratoire des Technologies d’Assemblage (LTA) du CEA Saclay et l’Institut de Recherche Dupuy de Lôme (IRDL - UMR CNRS 6027) de l’Université de Bretagne-Sud (Lorient).
Analyse sismique de l’interface sol-fondation : Modélisation physique et numérique du basculement global et du décollement local
Les fondations basculantes offrent un mécanisme potentiel pour améliorer la performance sismique en permettant un soulèvement et un tassement contrôlés, mais les incertitudes dans les interactions sol-fondation limitent leur utilisation généralisée. Les modèles actuels nécessitent des simulations numériques complexes, qui ne représentent pas de manière précise l'interface sol-fondation.
L'objectif principal de cette thèse est de modéliser la transition des effets locaux (friction, soulèvement) à la réponse globale de la structure (basculement, tassement et glissement) sous des charges sismiques, en utilisant une approche expérimentale et numérique combinée. Il s'agit donc d'assurer une modélisation numérique fiable des structures basculantes. Les objectifs clés incluent :
• Étudier la sensibilité des paramètres physiques dans la réponse sismique des systèmes sol-structure basculants en utilisant l'apprentissage automatique et des analyses numériques.
• Développer et réaliser des tests expérimentaux sous charges monotones puis dynamiques pour mesurer les réponses sol-fondation-structure en condition de basculement.
• Implémenter des simulations numériques pour tenir compte des effets d'interaction locaux et valider les résultats avec des résultats expérimentaux.
Enfin, cette recherche vise à proposer un cadre expérimental et numérique fiable pour améliorer la résilience sismique dans la conception en ingénierie. Cette thèse fournira à l'étudiant des compétences pratiques en ingénierie, ainsi qu'une expertise dans les tests en laboratoire et la modélisation numérique. Les résultats seront publiés dans des revues internationales et nationales et présentés lors de conférences, faisant avancer la recherche dans le domaine de la dynamique des sols et des structures.
Validation d'une approche d'identification sans modèle "data-driven" pour la modélisation de la rupture ductile
Ces travaux s'intéressent au passage des modèles constitutifs traditionnels vers un cadre de Mécanique Computationnelle Basée sur des Données (Data-Driven Computational Mechanics, DDCM), discipline introduite il y a peu [1]. Au lieu de s'appuyer sur des équations constitutives complexes, cette approche utilise une base de données d'états contrainte-déformation pour modéliser le comportement des matériaux. L'algorithme minimise la distance entre les états mécaniques calculés et les entrées de la base de données, garantissant le respect des équilibres et des conditions de compatibilité. Ce nouveau paradigme vise à surmonter les incertitudes et les défis empiriques associés aux méthodes conventionnelles.
Outil corollaire de la DDCM, l'Identification Basée sur des Données (Data-Driven Identification, DDI) a émergé comme une méthode puissante pour identifier les réponses en contrainte des matériaux [2, 3]. Elle fonctionne avec peu d'hypothèses et sans modèles constitutifs, ce qui la rend particulièrement adaptée à un potentiel étalonnage de modèles complexes couramment utilisés dans l'industrie.
Les objectifs clés de cette recherche incluent l'adaptation des stratégies de la DDCM vers la DDI pour la plasticité [4] et la rupture [5] en perfectionnant l'existant, l'amélioration de la DDI pour le calcul haute performance, et l'évaluation des équations constitutives à partir d'une base de données uniquement basée sur la donnée expérimentale. La méthodologie proposée consiste à collecter des cartes de mesures de champ à partir d'un essai hétérogène, en utilisant une caméra à haute vitesse et la corrélation d'image numérique. Elle adaptera la DDCM aux scénarios de fracture ductile pour la DDI, mettra en œuvre un solveur DDI dans un cadre de calcul haute performance et réalisera une évaluation d'un modèle constitutif. L'accent sera mis sur l'acier 316L, un matériau largement utilisé dans l'industrie nucléaire.
Cette thèse représente une collaboration des laboratoires CEA et Centrale Nantes qui sont spécialisés en mécanique computationnelle, mathématiques appliquées, l'ingénierie logicielle et traitement de signal.
[1] Kirchdoerfer, Trenton, and Michael Ortiz. "Data-driven computational mechanics." Computer Methods in Applied Mechanics and Engineering 304 (2016): 81-101.
[2] Leygue, Adrien, et al. "Data-based derivation of material response." Computer Methods in Applied Mechanics and Engineering 331 (2018): 184-196.
[3] Dalémat, Marie, et al. "Measuring stress field without constitutive equation." Mechanics of Materials 136 (2019): 103087.
[4] Pham D. et al, Tangent space Data Driven framework for elasto-plastic material behaviors, Finite Elements in Analysis and Design, Volume 216, 2023, https://doi.org/10.1016/j.finel.2022.103895.
[5] P. Carrara, L. De Lorenzis, L. Stainier, M. Ortiz, Data-driven fracture mechanics, Computer Methods in Applied Mechanics and Engineering, Volume 372, 2020, https://doi.org/10.1016/j.cma.2020.113390.
Rhéologie de suspensions concentrées en charges minérales
En tant qu’organisme de recherche dans le domaine du nucléaire et des énergies alternatives, le CEA participe à des études fondamentales mettant en jeu également des suspensions denses. En effet, des particules inorganiques (verre, zéolite, boues, sels ou ciment/sable) suspendues dans un fluide, parfois de très haute viscosité comme le bitume, font parties des systèmes à l’étude pour diverses applications telles que l’optimisation du remplissage de colis de verre (procédé Dem N’ Melt) ou encore de colis cimentaires, où les propriétés d’écoulement sont à optimiser pour garantir une homogénéité des futs de déchets. Outre, la problématique de la reprise (boues historiques), du traitement et du conditionnement de déchets en matrice vitreuse ou bitumineuse, des suspensions concentrées de grains de verre sont à l’étude pour la production de dihydrogène par électrolyse à haute température.
Il s'agira dans un premier de travailler sur des suspensions concentrées modèles et de caractériser leurs propriétés d'écoulement en cisaillement et en compression. Ce mode de sollicitation peut déclencher l’apparition de régimes frictionnels, de séparation de phases liquide/solide, et diverses réponses non-linéaires qu'il s'agira de modéliser. Apres cette première étape, la topologie, la granulométrie et la polydispersité des particules solides seront modifiées pour se rapprocher au plus prêt des suspensions rencontrées dans l'industrie.
Matériaux activés conducteurs pour la conversion énergétique et le stockage de l’énergie par effet capacitif
La production d’énergie à partir de sources renouvelables nécessite des systèmes de stockage efficaces pour pallier les déséquilibres entre offre et demande. Le projet propose de développer des super-condensateurs économiques en utilisant des électrodes composites issues de sous-produits industriels. Les liants minéraux, comme les géopolymères ou les Matériaux Alcali Activés (MAA), rendus conducteurs par dispersion de noir de carbone, sont étudiés pour des applications de stockage d'énergie ou de génération de chaleur. En se basant sur un brevet déposé récemment, nous nous proposons de réaliser une étude approfondie de ces composites conducteurs. Leurs performances seront évaluées en fonction des paramètres de formulation et de mise en forme. Il s’agira également de caractériser finement le réseau poreux et la dispersion des charges conductrices dans le matériau. Enfin, des essais de mise en forme du matériau seront menés et des super-condensateurs seront assemblés, permettant l’étude de l’impact du procédé (impression 3D) et des géométries.
Étude de la synthèse et des propriétés thermodynamiques des composés (An,Zr)O2 et (Zr,An)SiO4
Dans le cas d’un accident nucléaire grave, le combustible présent dans le cœur du réacteur peut entrer en fusion, entrainant la formation d’un composé désigné sous le nom de corium. Les cas d’accidents majeurs, ainsi que les expériences de formation de corium prototypiques, ont permis d’identifier la formation de composés clés que sont les oxydes mixtes (U,Zr)O2 formés par interaction du combustible avec la gaine en zircaloy et les silicates (Zr,U)SiO4 formés par interaction du corium avec des matériaux de structure. Dans le cas de combustibles MOx, (U,Pu)O2, la formation d’un corium pourrait conduire à la formation de phases équivalentes avec des teneurs significatives en plutonium. Cependant, les données thermodynamiques expérimentales sur de tels composés, qui permettraient d’évaluer leur comportement, sont à ce jour inexistantes. Dans ce cadre, la détermination des conditions de synthèse de tels composés avec un bon degré de pureté est essentielle à l’acquisition de telles données. La synthèse de solutions solides (Zr,Pu)O2 et (Zr,Pu)SiO4 constitue alors un premier pas essentiel avant d’étudier les systèmes (Zr,U,Pu)O2 et (Zr,U,Pu)SiO4.
L’objet de cette thèse sera de déterminer les conditions adaptées à la synthèse de ces composés, de réaliser un ensemble de caractérisation permettant d’évaluer leur pureté et d'établir leurs propriétés thermodynamiques. Pour cela des expériences seront conduites sur l’installation ATALANTE et une approche de caractérisation multi techniques sera choisie (DRX, spectroscopie Raman et infrarouge, MEB, caractérisation sur grand instruments…). Des essais de solubilité en milieu contrôlé seront ensuite mis en place et des mesures calorimétriques menées dans le cadre de collaborations internationales.
Mesure et évaluation de la dépendance énergétique des données de neutrons retardés du 239Pu
Cette proposition de thèse vise à mesurer et à caractériser l’émission des neutrons retardés émis par la fission du 239Pu. Cet actinide est impliqué dans divers concepts de réacteurs et la connaissance des données nucléaires qui le caractérisent reste actuellement insuffisante, en particulier en spectre rapide. Ce projet comprend une forte composante expérimentale, avec plusieurs campagnes de mesures sur l'accélérateur électrostatique MONNET au JRC Geel, auxquelles le doctorant prendra activement part.
La première étape de cette thèse consistera à intercomparer les méthodes de mesure du flux neutronique (dosimétrie, chambre à fission, détecteur long-counter et scintillateur à protons de recul) puis de les confronter à des calculs Monte-Carlo simulant l’émission des neutrons par interaction de particules chargées (D+T, D+D, p+T). Ce travail permettra d’assurer la bonne caractérisation du flux neutronique, une étape essentielle pour la suite du projet.
Dans un second temps, le doctorant devra reproduire des mesures de neutrons retardés du 238U, à l’aide d’une cible préexistante, dans une logique d’inter-comparaison par rapport à une campagne expérimentale menée en 2023.
Dans un troisième temps, le doctorant réalisera la mesure des rendements en neutrons retardés et des abondances par groupe du 239Pu, sur une gamme d’énergie de neutrons comprise entre 1 et 8 MeV. In fine, il produira une évaluation dépendante de l’énergie et l’intégrera dans un fichier ENDF pour être testée sur différents calculs de réacteur (beta-eff, transitoires de puissance, calibration d’efficacité d’absorbants…). Ces mesures complèteront une étude en spectre thermique menée à l’ILL en 2022 dans le but de former un modèle cohérent pour le 239Pu sur une gamme d’énergie de 0 à 8 MeV.
Ce projet contribuera au fichier de données nucléaires JEFF-4 de l’OCDE/AEN. Il répond à une forte demande de l’industrie nucléaire (soulignée par l’AIEA) pour améliorer la précision des mesures de multiplicité et des paramètres cinétiques des neutrons retardés, contribuant ainsi à une meilleure maîtrise de la sûreté des réacteurs nucléaires ainsi qu’à la réduction des marges de sûreté.