Développements des outils multiphysiques dédiés à la modélisation des réacteurs RNR-Na et études associées.
Le groupe sodium du DM2S (département du CEA Saclay) développe des outils numériques de couplage afin de réaliser des études de cas accidentels (transitoires rapides). Les domaines physiques concernés sont la neutronique, la thermo-hydraulique et la mécanique. Le sujet de ce post-doc s’inscrit dans ce cadre.
Il s’agit de mener plusieurs travaux : l’intégration d’un couplage au sein de la plateforme CORPUS, réaliser des études dans le but de tester les effets et introduire dans le couplage l’impact, sur l’écoulement du sodium, de la déformation des assemblages par la température, l’utilisation des sections efficaces neutroniques générées par le code APOLLO3, l’étude d’autres cas accidentels, et étendre la modélisation à l’échelle sous-canal et aiguille.
Développement d’un framework de calcul dédié à la réduction de modèles par la méthode des bases réduites certifiées.
De nombreux domaines de l’ingénierie requièrent de pouvoir résoudre numériquement des équations aux dérivées partielles (EDP) modélisant des phénomènes physiques.
Lorsque nous nous intéressons à un modèle mathématique qui décrit le comportement physique d’un système en s’appuyant sur une ou plusieurs EDPs paramétrées (paramètres géométriques ou physiques), il peut être souhaitable de pouvoir évaluer rapidement et de manière fiable la sortie du modèle (quantité d’intérêt) pour différentes valeurs des paramètres.
Le contexte temps-réel, nécessaire pour faire du contrôle commande, ainsi que les contextes demandant beaucoup d’évaluations des sorties du modèle (typiquement pour des méthodes d’optimisation ou d’analyse d’incertitudes et de sensibilité) s’y prêtent parfaitement.
La méthode des bases réduites est une méthode de réduction de modèle dite intrusive car, à la différence des méthodes de type non-intrusives, la réduction est basée sur la projection des opérateurs des EDPs du modèle physique.
Cette méthode permet d’obtenir de manière rapide, pour un ensemble de valeurs de paramètres donné, une approximation de l’évaluation de la sortie du modèle.
Un des points forts de la méthode est l’aspect "certifié" qui permet d’estimer l’erreur d’approximation de l’évaluation de la sortie du modèle.
L’objectif du post-doctorat est de développer un framework de calcul pour la méthode des bases réduites certifiées. Ce framework devra être basé sur la plateforme TRUST (https://sourceforge.net/projects/trust-platform/) développée au CEA et devra être suffisamment générique pour permettre de traiter différents types de problèmes (linéaires ou non, stationnaires ou non, coercifs ou non...).
Le framework devra pouvoir être utilisé dans le cas d’un modèle de mélange de deux fluides.
Stratégie d’identification de modèles continus à partir d’une approche discrète 3D
Afin de développer une stratégie d’identification de modèles de comportements pour matériaux quasi-fragiles, adaptés au calcul de structure, très souvent réalisée arbitrairement, un modèle fondé sur la méthode des éléments discrets a été élaboré. Le modèle discret est utilisé pour compenser le manque de données expérimentales nécessaires à l’identification du modèle continu. Grâce à certaines prédispositions vis-à-vis du mécanisme de fissuration, la mise en oeuvre d’un modèle discret est extrêment aisée, et son efficacité a été démontrée, mais se limite pour le moment à un cadre 2D, pour des raisons liées au temps de calculs.
Un cadre 2D, réduit fortement les possibilités d’analyses pouvant être conduites avec un tel modèle, notamment pour des structures renforcées, où les effets 3D sont prépondérents. L’objectif des travaux proposés dans le cadre de ce travail post-doctoral est donc d’étendre en 3D l’approche discrète initialement développée. Les développements se feront dans le code de calcul CAST3M-CEA développé au DEN/DANS/DM2S/SEMT. Le code de calcul discret sera, par la même occasion, optimisé à l’aide d’outils déjà présents dans l’environnement CAST3M-CEA. En fonction des gains en efficacité, le calcul de structures complètes par la méthode des éléments discrets pourra être envisagée.
A l’issue de ce travail, un outil numérique sera disponible permettant d’étendre la stratégie d’identification à des modèles de comportement intégrant des effets 3D, tels que les modèles d’interface acier/béton (confinement) et de béton (dilatance).
Analyse dynamique non-linéaire d’une structure en béton armé sous sollicitations sismiques : Etude déterministe et probabiliste des spectres de plancher
L’étude qui est proposée s’appuie sur la campagne ENISTAT et comporte trois axes de travail :
1. Calibration (5 mois)
Le candidat devra, dans un premier temps, s’approprier le modèle aux éléments finis de maquette qui a été développé en phase d’avant projet. Sur la base des résultats expérimentaux, les éventuelles sources d’écart entre les résultats numériques initiaux et expérimentaux devront être identifiées en vue de mieux calibrer le modèle.
2. Étude déterministe et probabiliste des spectres de plancher (5 mois)
Sur la base du modèle calibré, les spectres de planchers pourront être déterminés. Ces derniers seront confrontés aux spectres de demande issus des textes normatifs, tels que l’Eurocode 8. Enfin, en s’appuyant sur des méthodes issues de la théorie de la fiabilité, les incertitudes qui entachent non seulement les paramètres d’entrée du modèle mais aussi les signaux d’excitation seront prises en compte et la variabilité des spectres de planchers pourra être quantifiée. Cette donnée est particulièrement essentielle dans un contexte de dimensionnement et de retombées au niveau de l’ingénierie parasismique.
3. Étude de la tenue au séisme des rupteurs thermiques (2 mois)
Sur la base des résultats, non seulement expérimentaux mais aussi numériques, une discussion sur l’effet des rupteurs thermiques sur la tenue globale de la structure au séisme sera réalisée.
Il peut être noté que l’ensemble de ces recherches s’inscrit pleinement dans le cadre de la feuille de route de l’institut SEERA (ex. IMRS) dont le CEA est membre fondateur.
Réduction de modèle en dynamique vibratoire : application au génie parasismique
La complexité et la finesse des modèles numériques mis en oeuvre pour prédire le comportement des structures soumises à des sollicitations sismiques imposent bien souvent des temps de calculs de plusieurs jours pour la résolution des équations aux dérivées partielles du problème de référence.
En outre, l’étude des marges, les méconnaissances sur les paramètres constitutifs, les analyses stochastiques ou encore le recalage des modèles impose généralement de fournir cet effort numérique, non plus pour la simulation d’un seul modèle, mais d’une famille de modèles.
Afin de diminuer les temps de calculs, des techniques de réduction de modèle (Proper Orthogonal/Generalized Decomposition) peuvent être envisagées. La présente étude post-doctorale propose de définir et mettre en oeuvre (notamment dans CAST3M) une technique adaptée pour la réduction de modèles de génie civil (type béton armé) soumis à des chargements sismiques.
Calculs d'évolution à grande échelle avec code Monte-Carlo de transport de neutrons
L'un des principaux objectifs de la physique des réacteurs modernes est d'effectuer des simulations multi-physiques fidèles du comportement d'un cœur de réacteur nucléaire, avec une description détaillée de la géométrie à l'échelle des crayons de combustible. Les calculs multi-physiques en conditions nominales impliquent un couplage entre un solveur de transport pour les neutrons et les précurseurs, des solveurs thermiques et thermohydrauliques pour le transfert de chaleur, et un solveur Bateman pour calculer l'évolution isotopique du combustible nucléaire au cours d'un cycle du réacteur. L'objectif de ce post-doc est de réaliser un tel calcul entièrement couplé, à l'aide de PATMOS, une mini-app Monte-Carlo pour le transport de neutrons, et de la plateforme de couplage C3PO, toutes deux développées au CEA. Le système cible est un cœur de la taille d'un réacteur commercial.
Simulation des écoulements multiphasiques réactifs gaz-liquide
L'objectif de ce travail postdoctoral est de développer et mettre en œuvre une méthode de simulation pour la modélisation d'incendies provoqués par une pulvérisation de sodium. Pour cela, deux aspects du problème sont à considérer: tout d'abord, il est nécessaire de proposer une représentation adéquate des gouttelettes de sodium (phase dispersée) depuis leur génération par un jet (phase séparée) jusqu'à leur comportement (mouvement, oxydation, combustion) dans l'atmosphère. Cela nécessite de dériver un modèle d'écoulement' prenant en compte plusieurs composants avec plusieurs régimes topologiques d'interface (dispersé et séparé). Deuxièmement, il est nécessaire de développer une stratégie de discrétisation robuste pour ce modèle d'écoulement complexe.
Le travail numérique sera mis en œuvre dans un nouvel outil numérique pour effectuer des simulations d'incendies provoqués par la pulvérisation de sodium développé au CEA. Cet outil est basé sur la bibliothèque CanoP. Canop est une bibliothèque conçue pour résoudre des problèmes de dynamique des fluides numériques en utilisant une approche de raffinement de maillage adaptatif basée sur des cellules et un calcul parallèle.
Méthodes robustes de pilotage indirect du chargement pour la simulation de structures en béton armé
Les algorithmes de pilotage indirect (« path-following ») sont généralement employées pour décrire des réponses structurales instables caractérisées par des « snap-backs » et/ou des « snap-troughs ». Dans ces formulations, l’évolution des actions extérieures (efforts/déplacements) est calculée pendant la simulation pour satisfaire un critère de pilotage donné. Adapter le chargement externe pendant le calcul est utile pour obtenir la solution du problème, mais également pour réduire le nombre d’itérations à convergence. Ce second aspect est d’une importance primordiale, notamment pour les calculs à grande échelle (c.-à-d., à l’échelle de la structure). Différentes formulations « path-following » ont été proposées dans la littérature. Malheureusement, un critère objectif pour choisir une formulation plutôt qu’une autre pour la simulation de la réponse de structures en béton armé (BA) (en présence de mécanismes dissipatifs différents et complexes) n’est pas encore disponible. Le travail proposé portera sur la formulation d’algorithmes de pilotage indirect du chargement adaptés pour simuler des structures BA.
Accélération GPU d’un code de transport déterministe DSN 3D en neutronique.
Dans le cadre des Programmes Transversaux de Compétences (PTC), les équipes du DES/ISAS/DM2S et celles du CEA-DIF montent une collaboration sur le portage GPU de leurs codes de transport déterministe en neutronique.
D’un côté, les équipe du DES/ISAS/DM2S sont en charge du développement et de la prospective sur les codes de calcul de neutronique déterministe pour la physique des réacteurs, en particulier le code APOLLO3®. De l’autre, le laboratoire de neutronique du CEA-DIF est impliqué dans le développement des codes de neutronique déterministe utilisés dans le cadre du programme Simulation.
Les deux unités cherchent aujourd’hui à préparer l’arrivée d’une nouvelle génération de supercalculateurs massivement dotés en GPU. Elles entament simultanément des études de faisabilité et de prospective sur cette thématique. Parce que les problématiques à étudier, la démarche mise en œuvre et les conclusions qui en seront tirées peuvent largement être mutualisées, il a été identifié un fort intérêt pour créer des synergies entre les deux unités. Ces synergies passent par la mise en place d’échanges réguliers et ont abouti à l’ouverture d’un poste de post-doctorat commun. Ce poste sera accueilli par les équipes du SERMA au CEA Saclay, mais constituera l’interface privilégiée entre les deux unités.
Le sujet de post-doctorat a pour objectif d’étudier le portage d’un code-jouet de neutronique déterministe DSN3D sur GPU.
Il s’appuie sur les expériences de portage déjà réalisées dans les deux unités d’accueil sur la base d’approches complémentaires : une approche haut-niveau est choisie côté DES au travers de la plate-forme Kokkos, tandis qu’une approche bas-niveau en langage Cuda est retenue côté CEA-DIF.
Etude du couplage neutronique-thermohydraulique dans les coeurs RNR-Na hétérogènes
Dans le cadre du développement du prototype ASTRID (Advanced Sodium Test Reactor for Industrial Demonstration) un processus de mise à jour des méthodes de calcul du cœur est initié. Ces méthodes bénéficieront du développement des calculs à hautes performances ainsi que des capacités de couplage avancé des codes de nouvelle génération. L’objectif est d’intégrer ces méthodes dans les démonstrations de sûreté à venir avec en particulier les méthodes couplées neutronique, thermo-hydraulique et physique du combustible.
Parmi ces méthodes, les couplages neutronique/thermo-hydrauliques diphasiques sont envisagés en particulier pour les études de sûreté de type :
• Perte de réfrigérant (Loss Of Flow transients)
• Passage d’incondensables gazeux
Ces couplages sont particulièrement pertinents pour l’analyse de réacteurs hétérogènes (exemple du réacteur CFV) dont la démonstration de sûreté repose en partie sur un calcul fidèle des fuites neutroniques.