Développement et optimisation de techniques de rafinement de maillage adaptatif (AMR) pour des problèmes d'intéraction fluide/structure dans un contexte de calcul haute performance
Le CEA développe actuellement un nouveau code de simulation pour la mécanique des structures et des fluides compressibles : Manta. Ce code a pour double objectif d'unifier les fonctionnalités des codes historiques implicite et explicite du CEA et d'être nativement orienté vers le calcul intensif. Grâce à de nombreuses méthodes numériques (éléments finis, volumes finis, résolutions de problèmes implicites ou explicites, ...), Manta permet de simuler différents types de problèmes mécaniques dynamiques ou statiques pour la structure et le fluide, ainsi que l'interaction fluide-structure.
Dans le cadre de la recherche d'optimisation et de gain en temps de calcul, une des techniques incontournables pour améliorer la précision des solutions tout en maîtrisant les coûts de calcul est l'adaptation dynamique du maillage (ou AMR pour « Adaptive Mesh Refinement »).
Ce postdoc s'attache à la définition et à la mise en œuvre d'algorithmes d'AMR dans un contexte de calcul haute performance pour des problèmes faisant intervenir des fluides et des structures en intéraction.
Une tâche préliminaire consistera à implémenter des fonctionnalités de raffinement de maillage hiérarchique dans Manta (sous-découpage/fusion de cellules, transferts des champs, critères de raffinement, création de liaisons pour les « hanging-nodes »). Ces travaux se feront si possible en s'appuyant sur des librairies externes.
Dans un second temps, il s'agira d'optimiser les performances des calculs parallèles à mémoire distribuée. En particulier, il sera essentiel de définir une stratégie d'équilibrage de charge entre les processus MPI, en particulier dans le cadre de problèmes d'intéraction fluide/structure.
Enfin, en particulier pour des calculs explicites, il faudra définir et mettre en œuvre des techniques d'adaptation du pas de temps en fonction du niveau de raffinement.
Ces deux derniers points donneront lieu à une ou plusieurs publications dans des revues spécialisées.
Mise en œuvre de capteurs permettant le suivi en ligne de la corrosion des aciers inoxydables en milieu acide nitrique chaud et concentré
La maîtrise du vieillissement des matériaux des équipements (principalement en acier inoxydable) de l'usine de recyclage du combustible nucléaire usé, fait l'objet d'une attention permanente. Certaines installations de l’usine de la Hague devront d’ailleurs être remplacées très prochainement. Dans ce contexte, il est important pour l’industriel de développer des capteurs, résistants à l’acide nitrique concentré (˜ 2,5 mol/L) et à la température (de l’ambiante à 130 °C), permettant de suivre la corrosion en ligne.
L’objectif de ce travail est de fabriquer un capteur permettant de détecter la corrosion de l’acier.
Les challenges de ce sujet de post-doc sont essentiellement technologiques puisqu’il s’agira de développer ou d’utiliser des matériaux adaptés à des milieux acides nitriques concentrés et chauds.
Le laboratoire est spécialisé dans l'étude de la corrosion dans des conditions extrêmes. Il est composé d'une équipe scientifique très dynamique et motivée.
Description de l’évolution de la taille de grain et des densités de dislocations lors de la consolidation des aciers ODS
Les aciers renforcés par une dispersion d’oxydes nanométriques (généralement désignés aciers ODS) sont envisagés notamment comme matériau de gainage combustible des réacteurs de 4ème génération. Ces matériaux sont à l’heure actuelle élaborée classiquement par métallurgie des poudres. L‘évolution de la microstructure lors de l’élaboration n’est pas encore bien décrite. Des travaux récents menés au laboratoire ont porté sur l’évolution de la nano-précipitation lors de l’élaboration. L’objectif du post-doctorat est donc d’affiner la description de cette évolution, plus particulièrement vis-à-vis de la taille de grain et de la densité de dislocations. Ce sujet couple une approche expérimentale, aux travers d’analyses en microscopie électronique et en diffraction des rayons X, et une approche numérique, visant à définir une méthode optimisée pour le traitement de l’évolution des dislocations.
Modélisation de l’altération du combustible irradié en milieu saturé avec effet de la température
La modélisation de l'altération des combustibles irradiés dans l'éventualité d'un entreposage de longue durée en piscine ou d'un stockage géologique profond est essentiel pour prédire leur comportement à long terme. Dans l'éventualité d'un assemblage défectueux et d'un contact direct entre le combustible et l'eau, l'altération par l'eau peut conduire à une dégradation des crayons et au relâchements des radionucléides en solution. Un modèle géochimique couplant la chimie au transport (transport réactif) a fait l'objet de premiers développements en lien avec le comportement des combustibles en situation de stockage géologique profond. La plate-forme HYTEC développée par l'Ecole des Mines de Paris a été utilisée pour ces premiers développements de simulation. Ces simulations menées à 25°C prennent en compte les mécanismes d'altération des combustibles, les cinétiques réactionnelles associées et des bases de données thermodynamiques robustes. Il est aujourd'hui important dans le cadre de ce post-doctorat de poursuivre ces développements dans une gamme de température allant jusqu'à 70°C. le modèle existant devra également être adapté à d'autres conditions que celles d'un stockage et notamment à la situation d'un entreposage de longue durée dans des piscines dédiées.
Stabilité sous irradiation de l'interface oxyde / métal d'alliage d'aluminium 6061-T6 traitée par oxydation anodique
L’alliage d’aluminium 6061-T6 a été utilisé par la fabrication des composants principaux du cœur du réacteur Jules Horowitz (RJH) en raison de sa bonne résistance à la corrosion, de ses bonnes propriétés mécaniques et de sa haute transparence aux neutrons. Afin d’améliorer la résistance à l’usure et à l’oxydation, certaines pièces subissent un traitement de surface appelé Oxydation Anodique dure (OAd) qui permet de forcer la croissance d’une couche d’oxyde qui renforce les propriétés tribologiques.
Sous irradiation, les structures du cœur du RJH sont soumises à un fort flux de neutrons qui modifie la microstructure et les propriétés de l’alliage, les neutrons rapides induisent des cascades de déplacements d’atomes et créent des défauts ponctuels qui se regroupent sous forme d’amas pour former des boucles de dislocation et des amas lacunaires tridimensionnels. Ces amas durcissent l’alliage et induisent un gonflement macroscopique. Dans l’alliage métallique d’aluminium 6061-T6, on observe également une dissolution des nano-précipités préexistants et la re-précipitation des éléments sous une forme hors équilibre. L’oxyde (Oad) est également sujet à des évolutions microstructurales sous irradiation (densité de boucles après irradiation, gonflement, amorphisation) mais il a néanmoins fait l’objet d’un nombre plus limité d’études.
Cette étude postdoctorale vise donc à caractériser la microstructure de la couche OAd, son évolution sous irradiation et sa résistance mécanique avant et après irradiation.
Implantation et extension de la bibliothèque de solveur Alien dans la proto-application Hélix
Le travail du post-doctorant sera dans un premier temps d’intégrer la solution Alien dans Helix, de réaliser des évaluations de performances et d’usage en configuration solveur direct ou itératif. Ces évaluations seront réalisées sur divers architectures allant de la machine de bureau jusqu'aux supercalculateurs nationaux sur plusieurs centaines voire milliers de cœurs de calcul.
Dans un second temps, le post-doctorant traitera la possibilité d’ajouter de nouvelles fonctionnalités dans la bibliothèque Alien pour résoudre des systèmes non-linéaires composés d'équations et d'inéquations afin de pouvoir résoudre dans un cadre HPC certains de problèmes de mécanique comme les méthodes de champs de phase ou bien les problèmes de contacts, sujet encore ouvert dans la communauté. Les résultats seront comparés aux cas tests et benchmark classiques de l’état de l’art du domaine.
Le candidat intégrera l'équipe de développement d'Helix, constituée pour le moment de 3/4 personnes au sein du laboratoire LM2S d'une quinzaine de personnes. Le post-doc se réalise dans le cadre d'un projet transverse entre les différentes directions du CEA. Le candidat collaborera donc également avec les auteurs de la bibliothèque Alien, à la DAM du CEA.
Modélisation des vents de vallées par descente d'échelle statistique
Pour modéliser et suivre les émissions atmosphériques dans une zone ayant un relief significatif, il est essentiel de représenter les vents à l'échelle de ce relief. Le modèle météorologique opérationnel de Cadarache ne dispose que d’une résolution horizontale de 1km, ce qui ne lui permet donc pas de résoudre les effets orographiques de la vallée.
Cependant, obtenir des résultats de simulation avec un modèle haute résolution nécessite des temps de calculs encore incompatibles avec les contraintes de la prévision météorologique opérationnelle (6h de calcul sur nos serveurs pour 1h de prévision pour Cadarache en 2020). Ceci contraint la résolution horizontale des calculs et ne permet pas de résoudre les effets orographiques de vallée.
L’objet du post-doc est donc de mettre au point un modèle de descente d'échelle appliqué à un maillage 3D de la vallée, avec une résolution suffisante pour, à la fois, modéliser l’aérologie de la vallée et y suivre un panache de pollution à l’aide d’un modèle de dispersion atmosphérique. Elle sera mise en œuvre via l'utilisation d'un réseau de neurones artificiels dont l'apprentissage sera fondé sur des mesures de climatologie et d'aérologie locales, complétées par des données synthétiques utilisant un modèle local à haute résolution.
Le/la candidat(e) évoluera au sein d’une petite équipe CEA attentive et bienveillante tout en restant connecté à la recherche universitaire via le Laboratoire d’Aérologie de Toulouse. Il pourra à la fois devenir un spécialiste en recherche appliquée dans le domaine météorologique et acquérir des compétences numériques et scientifiques valorisables en entreprise.
Optimisation des scénarios de transition énergétique par une approche dynamique de l’Analyse de Cycle de Vie
La modélisation de la transition énergétique, avec une projection jusqu’en 2050 et adaptable à différents pays ou stratégies, est complexe en termes d’ACV car elle fait intervenir de nombreux paramètres :
- une dizaine d’énergies possibles, avec des inventaires évolutifs de construction des infrastructures de production/stockage d’électricité
- une évolution délicate à estimer des technologies pour une filière donnée
- une production d’électricité en adéquation avec la consommation nationale, pouvant prendre des trajectoires très variées
- des scénarios très contrastés possibles, incluant des montées en puissance plus ou moins rapides des renouvelables et une baisse temporaire du nucléaire, compensée ou non par des centrales à cycle combiné au gaz
- une nécessité de prévoir plusieurs formes de stockage de l’électricité selon l’importance du parc d’énergies non pilotables, avec des puissances dépendant de la durée de stockage
- la corrélation ou non de la puissance de stockage avec le niveau d’interconnexion des réseaux électriques européens.
Le travail consistera à analyser les inventaires disponibles dans la base Ecoinvent couplée à SimaPro, à les modifier selon les technologies prévisibles pour le moyen terme, à compléter la modélisation en langage Python pour inclure la totalité des paramètres envisagés.
L’objectif est de déterminer les meilleurs trajectoires possibles d’un point de vue environnemental, en s’appuyant sur les « midpoints » puis les « endpoints » pour réduire le nombre d’indicateurs et fournir in fine un outil flexible d’aide à la décision.
modélisation de la cinétique de précipitation de l’uranium en fonction du pH. Application à un réacteur à lit fluidisé
L’usine Orano au Niger (Somaïr) précipite son concentré uranium dans un réacteur à lit fluidisé par ajout d’hydroxyde de sodium. Le concentré obtenu contient environ 6% de sodium qui entraine des pénalités du convertisseur. Orano a effectué en fin d’année 2019 des essais sur un lit fluidisé au laboratoire pour changer le point de fonctionnement de la précipitation et former préférentiellement de l’UO3 via un changement de pH. Pour affiner le pilotage de l’unité industrielle, une modélisation des réactions de précipitation de l’uranium s’avère nécessaire. Le candidat devra proposer et calibrer un modèle de la précipitation compétitive de Na2U2O7 et UO3 basé sur les constantes d’équilibre et des cinétiques des réactions, en fonction du pH au sein du réacteur. En particulier, le modèle devra permettre de comprendre l’impact du pH sur la répartition des deux espèces principales identifiées dans le concentré : Na2U2O7 et UO3. Ce modèle chimique devra servir de donnée d’entrée à un modèle physique existant du réacteur à lit fluidisé. Un élargissement du modèle à d’autres réactifs de précipitation, notamment la magnésie pourrait également être étudié.
Synthèse et caractérisation structurale de phases minérales d’uranium de référence pour l’identification de phases porteuses d’uranium en environnement minier par SLRT.
Dans le cadre d’un projet de recherche collaboratif entre le CEA et Orano, une étude est menée pour détecter et identifier les phases porteuses d’uranium dans les résidus des traitements miniers et les stériles par spectrofluorimétrie laser à résolution temporelle (SLRT). La spectroscopie SLRT est un outil de spéciation permettant de sonder l’environnement chimique local de l’uranyle. Néanmoins, l’identification des phases uranifères au sein d’un échantillon prélevé sur site par SLRT repose sur l’utilisation d’un modèle chimiométrique construit à partir d’une base de spectres de référence. Afin d’enrichir et diversifier la base de spectres de référence la préparation et la caractérisation structurale d’échantillons uranifères synthétiques au sein de trois familles principales s’avèrent nécessaire. Ces familles phases minérales sont les oxyhydroxydes, les sulfates et silicates d’uranyle. L’acquisition des spectres SLRT des phases synthétisées et la mise en évidence de l’impact de la chimie de coordination de l’ion uranyle sur les spectres enregistrés constitue le second volet de ce travail.