Modélisation de frottement interfacial en géométrie grappe dans le code de thermo hydraulique système CATHARE
Le code système de thermo-hydraulique CATHARE, développé au CEA en partenariat avec EDF, Framatome et l’ASNR, est un code permettant de simuler les comportements normaux et accidentels du circuit hydraulique d’un Réacteur à Eau Pressurisée (REP). Il est aujourd’hui une référence en France en matière de simulation des transitoires dans les réacteurs nucléaires, et est notamment un composant central dans les rapports de sûreté produits par EDF et Framatome.
Des études précédentes montrent le besoin d’améliorer la validité du modèle de frottement interfacial en coeur à basse pression ou à grand diamètre hydraulique. De plus, le modèle actuel de frottement interfacial en coeur repose sur de nombreuses simplifications et sur une calibration effectuée sur un nombre restreint de données expérimentales pour des écoulements eau-vapeur à haute température. L’émergence de nouvelles modélisations dans la littérature et la disponibilité de nouvelles données expérimentales pourraient permettre une complexification du modèle et une calibration sur une plus grande base de données.
Cette thèse a donc pour objectif d’améliorer le modèle de frottement interfacial en géométrie grappe en investiguant les phénomènes physiques associés à cette problématique. Ce travail permettra d’implémenter à la suite de la thèse un modèle plus complet dans le code CATHARE et donc d’élargir les conditions de validité du code à de nouvelles applications.
Modélisation de la Condensation en Paroi et des Interactions avec le Film Liquide
Dans cette thèse, on s’intéresse à la modélisation des transferts de masse et d’énergie associés à la condensation pariétale pour un écoulement turbulent d’un mélange de vapeur et de gaz incondensables. L’écoulement est diphasique, en régime turbulent, où les modes de convection forcée, mixte et naturelle peuvent être rencontrés. Le cadre de ce travail est l’approche RANS des équations de Navier–Stokes compressibles, où la condensation en paroi sera décrite par des fonctions de paroi semi-analytiques développées dans un précédent travail de thèse cite{iziquel2023}. Ces fonctions intègrent les différents modes de convection ainsi que les phénomènes d’aspiration et d’interdiffusion des espèces, mais négligent la présence du film liquide.
Dans la littérature, l’effet de la formation et de l’écoulement du film liquide sur les transferts de masse et d’énergie est souvent négligé, car on considère qu’en présence de gaz incondensables, la résistance de la couche gazeuse à la diffusion de la vapeur est nettement plus importante que la résistance thermique du film liquide.
L’objectif de cette thèse est d’améliorer la prédiction des transferts thermique et massique par l’investigation, au-delà de la résistance thermique du condensat, de l’effet dynamique du liquide et de son interaction avec la couche gazeuse de diffusion lors de la condensation pariétale, en considérant d’abord un écoulement de film laminaire, puis en tentant de prendre en compte le régime turbulent.
Dans la phase gazeuse, le modèle de fonctions de paroi développé dans la thèse de A. Iziquel (2023) pour un mélange binaire de vapeur et d’un gaz incondensable sera étendu aux mélanges de vapeur et de $n>1$ gaz incondensables (N2, H2, …) afin de traiter la thématique du risque hydrogène.
La validation des modèles implantés sera effectuée sur la base de résultats d’expériences à effets séparés (SET) et à effets couplés (CET) disponibles dans la littérature (Huhtiniemi, COPAIN, ISP47-MISTRA, ISP47-TOSQAN, RIVA). Les comparaisons à l’échelle CFD, avec des fonctions de paroi pour la condensation négligeant le film sur des cas tests de la littérature et des expériences de condensation (COPAIN), permettront d’évaluer l’impact de cette hypothèse ainsi que l’apport de la nouvelle modélisation en termes de précision et de coût de calcul.
Conception et développement d’algorithmes asynchrones pour la résolution de l’équation du transport des neutrons sur des architectures massivement parallèles et hétérogènes
Cette proposition de thèse s’inscrit dans le cadre de la résolution numérique d’équations aux dérivées partielles par le biais d’une discrétisation des variables. Elle s’intéresse, dans un formalisme d’éléments finis, à travailler sur la conception d’algorithmes au travers de modèles de programmation parallèle et asynchrone pour la résolution de ces équations.
Le cadre industriel applicatif est la résolution de l’équation de Boltzmann appliquée au transport des neutrons dans le cœur d’un réacteur nucléaire. Dans ce contexte, beaucoup de codes modernes de simulations’appuient sur une discrétisation par éléments finis (plus précisément, un schéma Galerkin discontinu décentré amont) pour des maillages cartésiens ou hexagonaux du domaine spatial. L’intérêt de ce travail de thèse prolonge des travaux précédents pour explorer leur extension dans un cadre d’architecture distribuée qui n’ont pas été abordé jusque-là dans notre contexte. Il s’agira de coupler des stratégies algorithmiques et numériques pour la résolution du problème à un modèle de programmation qui expose du parallélisme asynchrone.
Ce sujet s’inscrit dans le cadre de la simulation numérique des réacteurs nucléaires. Ces simulations multiphysiques coûteuses requièrent le calcul du transport des neutrons en cinétique qui peuvent être associées à des transitoires de puissance violents. La stratégie de recherche adopté pour cette thèse permettra de gagner en coût de calcul, et alliée à un modèle massivement parallèle, peut définir les contours d’un solveur neutronique efficace pour ces problèmes multiphysiques.
Un travail réussi dans le cadre de cette thèse permettra à l’étudiant de prétendre à un poste de recherche en simulation et analyse numérique de problèmes physiques complexes, par-delà la seule physique des réacteurs nucléaires.
Mécanismes de communications unidirectionnelles pour la décomposition de données des applications de transport de particules Monte-Carlo
Dans le cadre d’un calcul Monte-Carlo d’évolution d’un cœur de REP (réacteur à eau pressurisée), il est nécessaire de calculer un très grand nombre de taux de réaction neutron-noyau, impliquant un volume de données pouvant dépasser la capacité mémoire d’un nœud de calcul sur les supercalculateurs actuels. Dans le cadre de Tripoli-5, les architectures à mémoire distribuée ont été identifiées comme cible pour le déploiement de calcul à haute performance. Pour exploiter de telles architectures il convient donc d’utiliser des approches de décomposition de données, notamment sur les taux de réaction. Toutefois, avec une méthode de parallélisation classique, les processus n’ont pas d’affinité particulière pour les taux qu’ils hébergent localement ; au contraire, chaque taux reçoit des contributions de manière uniforme de tous les processus. Les accès aux données décomposées peuvent s’avérer coûteux quand ces derniers imposent une utilisation intensive des communications. Toutefois, des mécanismes de communications unidirectionnelles comme par exemple les MPI RMA (Message Passing Interface, Remote Memory Access) permettent de faciliter ses accès aussi bien en termes d’expressions que de performances.
L’objectif de cette thèse est de proposer une méthode de décomposition partielle de données en s’appuyant sur des mécanismes de communications unidirectionnelles pour accéder aux données stockées à distance, telles que les taux de réaction. Une telle approche permettra de réduire considérablement le volume de donnée stocké en mémoire sur chaque nœud de calcul sans engendrer une forte dégradation des performances.
Simulations HPC diphasiques par méthodes de Boltzmann sur réseaux et adaptation automatique de maillage
Le CEA/STMF développe des outils de calcul scientifique en thermohydraulique qui ont pour objectif de quantifier les transferts de masse et d’énergie dans les systèmes nucléaires du cycle tels que les réacteurs et les dispositifs de retraitement ou de confinement des déchets radioactifs. Dans ce travail de thèse, on s'intéresse aux méthodes “Lattice Boltzmann” (LBM) adaptés aux maillages dynamiques (Adaptative Mesh Refinement – AMR) dans un environnement informatique mutualisé et générique sur base Kokkos et exécutable sur les supercalculateurs multi-GPU. Le travail proposé consiste à développer dans le code Kalypsso-lbm les méthodes LB pour simuler des Equations aux Dérivées Partielles (EDPs) couplées qui modélisent les écoulements incompressibles diphasiques et multi-composants comme ceux rencontrées dans les dispositifs de l'aval du cycle. Une fois les développements réalisés, ils seront validés avec des solutions de références qui permettront une comparaison des méthodes d'interpolation entre les blocs de différentes taille du maillage AMR. Une discussion sera réalisée sur les critères de raffinement et de déraffinement qui seront généralisés pour ces nouvelles EDPs. Enfin, des benchmarks de performance quantifieront l'apport de l'AMR sur des simulations 3D lorsque la simulation de référence est réalisée sur un maillage statique et uniforme. Ce travail exploitera les supercalculateurs déjà opérationnels (e.g. Topaze-A100 du CEA-CCRT), ainsi que le supercalculateur exascale Alice Recoque selon l'état d'avancement de son installation.
Architecture innovante et traitement du signal pour des télécommunications optiques mobiles
Les communications optiques en espace libre reposent sur la transmission de données par la lumière entre deux points distants, sans recourir à des fibres ou à des câbles. Cette approche s’avère particulièrement intéressante lorsque les connexions filaires sont difficiles à déployer ou trop coûteuses.
Cependant, ces liaisons sont fortement affectées par les conditions atmosphériques : brouillard, pluie, poussières et turbulences thermiques atténuent ou déforment le faisceau lumineux, entraînant une dégradation notable de la qualité de la communication. Les solutions existantes restent coûteuses et limitées, tant du point de vue des dispositifs optiques de compensation que des algorithmes de traitement du signal.
Dans ce cadre, la thèse vise à concevoir des liaisons optiques mobiles performantes et robustes, capables de s’adapter à des environnements dynamiques et perturbés. L’étude portera notamment sur l’exploitation de dispositifs de type Optical Phased Array (OPA) sur Silicium — une technologie issue des systèmes LiDAR « low cost » — offrant une voie prometteuse vers des architectures compactes, intégrées et à faible coût.
L’orientation principale des travaux concernera le développement d’approches algorithmiques avancées pour le traitement et la compensation du signal. Le ou la doctorant·e sera amené·e à concevoir un environnement de simulation dédié, permettant d’évaluer et de valider les choix architecturaux et les stratégies algorithmiques avant toute expérimentation pratique.
L’objectif global est de proposer une architecture intégrée, flexible et fiable, garantissant la continuité des communications optiques en mouvement, avec des applications potentielles dans les domaines aérien, spatial et terrestre.
Modélisation de l’impact de défauts dans les structures acier-béton. Identification des défauts critiques par méta-modélisation et algorithmes d’optimisation.
Pour faire face à des enjeux de constructibilité grandissant, les structures « acier-béton » (structures « steel – concrete » ou « SC ») deviennent une alternative prometteuse face aux structures classiques en béton armé. Ces éléments sont constitués d’un béton de remplissage, de deux plaques externes métalliques et de goujons en acier permettant d’assurer l’action composite. Si ces structures présentent un intérêt certain lié à leur comportement mécanique d’ensemble, la présence des plaques métalliques empêche un contrôle visuel de la qualité du bétonnage. Il apparaît donc essentiel de caractériser l’impact de défauts potentiels. C’est dans ce contexte que s’inscrit le sujet de thèse. Il s’agira, en s’appuyant sur des résultats récents du laboratoire, de proposer une démarche numérique de prise en compte des défauts dans les structures acier-béton. Le travail de thèse s’articulera en plusieurs étapes : validation d’une stratégie de modélisation du comportement mécanique de structures acier-béton sans défaut, introduction de défauts dans la simulation et évaluation de l’applicabilité de la stratégie numérique, construction d’un métamodèle et analyse de sensibilité et définition de la (ou des) configuration(s) de défauts critiques par algorithmes d’optimisation. L’un des objectifs opérationnels de la thèse est de disposer d’un outil permettant de déterminer les configurations de défauts critiques (taille, position et nombre) en lien avec un objectif fixé sur une quantité d’intérêt donnée (perte de résistance ou de rigidité moyenne par exemple). Ce travail s’appuiera donc sur l’utilisation et le développement d’outils numériques à l’état de l’art, dans les domaines de la modélisation par éléments finis, des techniques d’optimisation, d’analyse de sensibilité et d’optimisation. La thèse sera réalisée dans un cadre collaboratif riche, notamment en partenariat avec EDF.
Fonctions optiques intégrées sur plan focal micro-bolométrique pour l’imagerie infrarouge non refroidie
L’imagerie infrarouge en bande thermique (longueurs d’onde 8-14 µm) est un domaine en forte croissance, particulièrement dans les domaines de l’industrie, du transport, de l’environnement. Elle s’appuie sur une technologie de détection, les microbolomètres, pour laquelle le CEA-Léti est au meilleur niveau de l’état de l’art mondial. L’intégration de fonctions optiques avancées directement sur les détecteurs est une voie très prometteuse pour gagner en performance, en compacité et en coût dans les futures caméras infrarouges.
Les fonctions optiques envisagées comprennent le filtrage spectral, la polarimétrie, la correction de front d’onde, etc. Certaines visent à enrichir l’image par des informations indispensables aux applications telles que la thermographie absolue (mesure de température et d’émissivité), l’identification pour l’interprétation automatique de scène (machine vision), la détection de gaz, etc.
Les travaux proposés comprendront des activités de conception, réalisation et caractérisation électro-optique de matrices de microbolomètres fonctionnalisés. La conception de ces fonctions optiques utilisera des moyens de simulation 3D électromagnétique, elle prendra en compte leur compatibilité avec nos technologies de microbolomètres et les capacités de nos moyens de micro-fabrication. La réalisation se fera dans les salles blanches du CEA-Léti par des personnels dédiés mais le (la) candidat(e) prendra part à la définition et au suivi des travaux. Enfin, les caractérisations optiques et électro-optiques seront faites dans notre laboratoire, si besoin avec le développement de bancs de caractérisation dédiés.
Effets des hétérogénéités structurales sur les écoulements d’air à travers une paroi en béton armé
Le bâtiment réacteur représente la troisième barrière de confinement dans les centrales nucléaires. Il a pour rôle de protéger l’environnement en cas d’accident hypothétique en limitant les rejets vers l’extérieur. Sa fonction est donc étroitement liée à son étanchéité. Classiquement, l’estimation du taux de fuite s’appuie sur une bonne connaissance des propriétés de transfert (comme la perméabilité) associée à une démarche de simulation chaînée (thermo-)hygro-mécanique. Si le comportement mécanique de la structure est aujourd’hui globalement maîtrisé, des progrès restent nécessaires dans la compréhension et la quantification des écoulements. C’est particulièrement le cas en présence d’hétérogénéités (fissures, nid de cailloux, reprise, armatures, câbles, etc.) qui représentent autant de situations pouvant perturber localement la perméabilité. C’est dans ce cadre que se place le sujet de thèse.
Il s’agira, en s’appuyant sur une démarche combinant essais expérimentaux et simulation, d’améliorer la représentation des écoulements en prenant en compte l’effet des hétérogénéités. Une première analyse permettra de définir un plan d’expérience qui sera ensuite mis en œuvre. Les résultats seront analysés afin de caractériser empiriquement l’influence de chaque type d’hétérogénéités testé sur les propriétés de transfert. Une approche de simulation, exploitant les résultats expérimentaux, sera ensuite développée, s’appuyant la méthode des éléments finis et discrets. Enfin, l’applicabilité de la méthode à une structure réelle sera étudiée, en rendant compte des incertitudes concernant la présence et l’impact de ces hétérogénéités (approche probabiliste).
La thèse s’appuie donc sur des outils et des méthodes expérimentaux et numériques à l’état de l’art et sera réalisée dans un contexte collaboratif riche (CEA, ASNR, EDF).
Simulations hydrodynamiques de matériaux poreux pour l'endommagement ductile
Le comportement mécanique des matériaux métalliques sous sollicitation fortement dynamique (choc), et en particulier leur endommagement, est une thématique d'intérêt pour le CEA-DAM. Pour le tantale, l'endommagement est de nature ductile : par germination, croissance et coalescence de pores (vides) au sein du matériau. Les modèles usuels d'endommagement ductiles ont été développés à partir d'hypothèses simplificatrices de pores isolés dans la matière. Cependant des études récentes par simulations directes décrivant explicitement une population de pores répartis dans le matériau (ainsi que des observations expérimentales après rupture) ont montré l'importance de l'interaction entre pores pour la prévision de l'endommagement ductile. Toutefois, les mécanismes microscopiques de cette interaction restent à élucider. De plus, ces études numériques doivent être étendues aux échelles de longueur et de vitesses de sollicitation d'intérêt.
L'objectif de la thèse est d'étudier les phases de croissance et de coalescence de l'endommagement ductile au travers de simulations numériques directes d'un milieux poreux soumis à une sollicitation dynamique. Des simulations hydrodynamiques, dans lesquelles des pores seront maillés explicitement au sein d'une matrice continue, seront utilisées afin de se placer aux échelles d'intérêt de temps et de longueur. Le suivi de la population de pores au cours de la simulation renseignera à différents niveaux sur l'influence de l'interaction entre pores pendant l'endommagement ductile. D'abord, le comportement du massif sera comparé à celui prédit par les modèles classiques à pores isolés, montrant l'effet macroscopique de l'interaction entre pores. On s'intéressera également à l'évolution de la distribution de tailles dans la population de pores. Enfin, un dernier objectif sera de comprendre l'interaction microscopique pore à pore. Afin de tirer parti de la richesse des résultats de simulation, des approches issues de l'intelligence artificielle (réseau de neurones sur le graphe associé à la population de pores) seront utilisées afin d'apprendre le lien entre voisinage d'un pore et croissance de celui-ci.
Le/la doctorant(e) aura l'occasion de développer ses compétences en physique des chocs et en mécanique, en simulations numériques (avec l'accès aux supercalculateurs du CEA-DAM) et en science des données.