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.
Impact de la pollution sur la dynamique des écoulements à bulles
En condition d'accident, si le coeur d'un réacteur nucléaire entre en ébullition, la pollution de l'eau qui le constitue peut avoir un rôle important sur les échanges de chaleur. L'enjeu de cette thèse est de comprendre cet impact et d'apprendre à le simuler, le but étant à terme de fournir des données de référence pour l'ébullition en conditions réacteur. Pour y parvenir, cette thèse s'attachera à simuler le transport d'une concentration de polluant au sein d'un écoulement à bulles. L'étudiant simulera la pollution des interfaces par des molécules tensio-actives, un cas particulier de polluant que l'on retrouve dans la plupart des systèmes hydrauliques. Cette étude sera réalisée à partir de Simulations Numériques Directes réalisées avec le code open-source TRUST/TrioCFD. L'étudiant sera accueilli au Laboratoire de Modélisation et de Simulation en mécanique des Fluides (LMSF) au sein d'un groupe de chercheurs et de nombreux doctorants. En collaboration avec le monde académique, l'étudiant publiera ses travaux et participera à des conférences internationales. Nous recherchons donc un étudiant ayant suivi une formation complète en mécanique des fluides numérique (M2 ou équivalent). Une connaissance du langage C++ moderne serait un avantage notable. La réalisation d'un stage en amont de la thèse est possible.
Modélisation sous-maille des transferts interfaciaux de masse et de chaleur appliqués à la condensation des essaims de bulles
Pour évaluer la sûreté des centrales nucléaires, le CEA développe et utilise des outils de simulation multi-échelles en thermohydraulique. L’application de la CFD aux écoulements diphasiques est limitée car elle nécessite de nombreux modèles difficiles à déterminer. Parmi nos autres ces outils, les simulations numériques directes (DNS) à interfaces résolues fournissent des données de référence inaccessibles par des moyens expérimentaux. C'est par exemple le cas des essaims de bulles, où les transferts de chaleur et de masse sont influencés par des effets collectifs complexes.
Afin de réduire le coût de ces simulations DNS, nous avons récemment développé une approche [1] qui montre des résultats prometteurs : elle consiste à coupler une résolution fine des transferts thermiques aux interfaces liquide-vapeur à un champ lointain calculé sur un maillage moins résolu. Pour élargir l'application de cette méthode à des cas plus industriels, il est nécessaire de prendre en compte les collisions entre bulles et d’adapter le modèle au changement de phase.
Nous proposons au cours de cette thèse de commencer par ce travail de modélisation physique et son implémentation en C++ dans notre code open-source de simulation TRUST/TrioCFD [2]. Ensuite, nous utiliserons cette nouvelle capacité pour réaliser une étude paramétrique et une analyse physique approfondie des phénomènes qui mèneraitmèneront, à terme, à une amélioration des modèles de transfert de chaleur dans les codes industriels.
[1] M. Grosso, G. Bois, A. Toutant, Thermal boundary layer modelling for heat flux prediction of bubbles at saturation: A priori analysis based on fully-resolved simulations, International Journal of Heat and Mass Transfer, Vol 222, 2024, https://doi.org/10.1016/j.ijheatmasstransfer.2023.124980
[2] Trio_CFD webpage : http://triocfd.cea.fr/recherche/modelisation-physique/two-phase-flows
Développement d'une approche macroscopique pour la dégradation à long terme des structures en béton sous irradiation
Dans les centrales nucléaires, la protection biologique en béton (CBS) est conçue à proximité de la cuve du réacteur. Cet élément, qui joue également le rôle de structure porteuse, absorbe donc les radiations. Il est ainsi exposé pendant la durée de fonctionnement de la centrale à des niveaux élevés de radiations qui peuvent avoir des conséquences à long terme. Ces radiations peuvent notamment entraîner une diminution des propriétés mécaniques des matériaux et de la structure. Etant donné son rôle clé, il est donc nécessaire de développer des outils et des modèles, pour prédire les comportements de telles structures à l'échelle macroscopique.
Sur la base des résultats existants obtenus à une échelle inférieure - simulations mésoscopiques, à partir desquelles une meilleure compréhension de l'effet de l'irradiation peut être obtenue, et des résultats expérimentaux qui viendront alimentés la simulation (propriétés des matériaux en particulier), il est proposé de développer une méthodologie macroscopique pour le comportement de la protection biologique en béton. Cette approche inclura différents phénomènes, parmi lesquels l'expansion volumétrique induite par le rayonnement, le fluage induit, les déformations thermiques et le chargement mécanique.
Les outils seront développés dans le cadre de la mécanique de l'endommagement. Les principaux défis numériques concernent la proposition et l'implémentation de lois d'évolution adaptées, et en particulier le couplage entre l'endommagement microstructural et l'endommagement au niveau structurel dû aux contraintes appliquées sur la structure.
Ce travail numérique pourra être réalisé dans un contexte de collaboration internationale. Il permettra au candidat retenu de développer un ensemble de compétences autour de la simulation de structures en béton armé en environnement complexe.
Quantification des incertitudes et analyse de sensibilité pour les vibrations de structures minces sous écoulement axial
Les phénomènes d'interaction fluide-structure (IFS) sont omniprésents dans les installations industrielles où des structures sont en contact avec un fluide sous écoulement qui exerce un chargement mécanique. Pour des structures élancées et souples, l’IFS peut induire des phénomènes vibratoires et des instabilités mécaniques à l’origine de grandes amplitudes de déplacements. L'industrie nucléaire est confrontée à cette problématique, notamment au niveau des tuyauteries, des assemblages de combustible, des grappes de commande ou encore des générateurs de vapeur. Dans le cadre de la prévention et de la maîtrise des risques vibratoires, les codes de calcul sont des outils essentiels, qui, à partir de plusieurs paramètres d’entrée, permettent d’accéder à des quantités d’intérêt (variables de sortie) souvent inaccessibles expérimentalement. Cependant, la connaissance des paramètres d’entrée est parfois limitée à cause d’un manque de caractérisation (erreur de mesure ou manque de données) ou tout simplement par la nature intrinsèquement aléatoire de ces paramètres. Il est alors nécessaire de prendre en compte des incertitudes sur les paramètres d’entrée et d’effectuer une analyse de sensibilité pour les variables de sortie.
Dans ce cadre, l'objectif de cette thèse est d’analyser la réponse vibratoire d'une structure mince dont les caractéristiques géométriques sont incertaines (structure présentant un défaut de courbure, localisé ou global). Nous chercherons en particulier à comprendre comment les incertitudes géométriques affectent la stabilité de la structure flexible. Cette caractérisation se fera de manière théorique et numérique. Les calculs numériques seront réalisés avec le code TrioCFD, développé par le CEA. En fonction de l’avancée des travaux, l’effet de différentes incertitudes (liées par exemple aux caractéristiques matérielles de la structure ou aux propriétés de l’écoulement incident) pourront être envisagées. In fine, le travail réalisé dans le cadre de cette thèse permettra d’améliorer la prédiction et le contrôle des vibrations de structures minces sous écoulement axial.
Les interactions fluides-structures et les instabilités associées sont présentes dans de nombreux domaines, que soit en aéronautique avec les phénomènes de flottement d’ailes, dans le nucléaire avec les vibrations de composants sous écoulements, en biologie pour la compréhension de la locomotion animale sous-marine, dans la botanique pour la compréhension des croissances végétales, dans le sport pour l’optimisation des performances, dans la récupération d’énergie des structures flexibles excitées par un fluide. La thèse permettra à l’étudiant d’acquérir de nombreuses compétences en mathématiques, en simulation numérique, en mécanique des fluides, en mécanique des solides, de se former à la recherche dans le domaine de la mécanique des fluides et des solides, pour in-fine une carrière dans ce domaine, que ce soit dans le monde académique ou de la R&D appliquée dans de nombreux domaines d’intérêt pour les scientifiques et plus généralement pour la société. Un sujet de stage de 6 mois est également proposé en préambule de la thèse (optionnel).
Niveau exigé : Master 2 / Dernière année d’École d’Ingénieur.
Formation exigée : mécanique des milieux continus, résistance des matériaux (théorie des
poutres et des plaques), mécanique des fluides, interaction fluides-structures, simulation
numérique (éléments finis).
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.
Etude des transitions de régime d’écoulement en post-assèchement
Les écoulements diphasiques interviennent dans de nombreux systèmes fluides, notamment pour le refroidissement des réacteurs nucléaires. Selon le flux thermique échangé dans le cœur du réacteur, le débit, la sous-saturation ou la pression, on peut constater des écoulements purement monophasiques, des écoulements à bulles ou annulaires (avec un film liquide au contact de la paroi et un cœur de vapeur).
Lors d’un accident de perte de réfrigérant primaire, le cœur du réacteur qui contient les crayons combustibles s’échauffe jusqu’à la crise d’ébullition lorsque le flux thermique est suffisamment important. Une illustration des régimes d’écoulements diphasiques lors de cet accident est présenté en figure 1. Un film de vapeur se forme rapidement et isole thermiquement les crayons, tandis que du liquide subsiste dans le cœur de l’écoulement. Les crayons du cœur sont asséchés, leur surface n’est refroidie que par de la vapeur et l’échange thermique à la paroi est ainsi dégradé [1]. Cet écoulement est du type "inverted annular film boiling". Au fur et à mesure que le liquide se vaporise, le film de vapeur s’épaissit et la turbulence induite aura tendance former des vagues a l’interface liquide-vapeur et à déstabiliser l’interface jusqu’à la formation de poches de liquides (inverted slug film boiling). Puis, l’évaporation et la fragmentation de ces poches vont mener à la formation d’un écoulement dispersé à gouttes (dispersed film boiling).
Actuellement, les transitions de régime d’écoulement dans cette configuration sont très mal identifiées [1], [2] bien que la connaissance de celles-ci soit importante pour l’étude de refroidissement du cœur du réacteur nucléaire. Une des principales difficultés expérimentales réside dans la nécessité de chauffer fortement les parois pour établir un film de vapeur et maintenir celui-ci, rendant de ce fait les sections d’essai opaques. Il est donc particulièrement complexe d’accéder à une visualisation directe et plus encore à des mesures de paramètres locaux comme les champs de températures et vitesses fluides. Les résultats expérimentaux disponibles dans la littérature sur ce sujet sont donc très limités et insuffisants pour développer un modèle physique [1], [3], [4], [5].
Cette thèse, qui constitue une première étape vers l’identification précise des transitions de régime, porte sur l’étude de l’effet purement hydrodynamique, en couplant des approches expérimentale et analytique. Afin d’obtenir une compréhension de la physique des différents phénomènes, la configuration d’un écoulement de liquide au cœur d’un écoulement gazeux est proposée. Dans celle-ci, la déformation de l’interface, la vitesse du gaz et la vitesse du liquide peuvent jouer sur la transition d’un régime à l’autre [6], [7] : l’interface lisse devient perturbée par des vagues (instabilités de Kelvin-Helmholtz), des gouttes sont arrachées de l’interface. Une analyse paramétrique sera effectuée en faisant varier les débits liquides et gazeux et ainsi l’épaisseur du film gazeux pour observer ces différents phénomènes et comprendre les effets de chaque paramètre sur les transitions de régimes. Expérimentalement, un banc a récemment été conçu au DM2S/STMF/LE2H afin d’étudier plus particulièrement ces transitions grâce à une visualisation des déformations de l’interface. Des adaptations pourront être apportées avec de nouvelles mesures ou éventuellement une nouvelle méthodologie si nécessaire.
A partir des résultats expérimentaux, il sera nécessaire d’identifier, voire de définir, les nombres adimensionnels pertinents pour décrire les phénomènes observés. L’analyse portera ensuite sur la caractérisation des transitions de régimes sur la base de ces nombres adimensionnels, afin de proposer une carte des transitions de régimes d’écoulements.
La combinaison des résultats obtenus permettra de renforcer les modèles utilisés dans les codes de calcul comme le code de thermohydraulique CATHARE, développé au CEA en particulier pour les études de sureté des réacteurs nucléaires. Cette thèse présente donc un fort intérêt académique par l’exploitation d’une installation expérimentale innovante et la production de résultats nouveaux qui confirmeront également son intérêt sur le plan industriel pour l’amélioration de la connaissance des phénomènes importants dans la démonstration de sûreté des réacteurs nucléaires.
Références :
[1] M. Ishii et G. De Jarlais, « Flow visualization study of inverted annular flow of post-dryout heat transfer region », Nuclear Engineering and Design, 1987.
[2] G. De jarlais, M. Ishii, et J. Linehan, « Hydrodynamic stability of inverted annular flow in an adiabatic simulation », Argonne National Laboratory, CONF-830702-9, 1983.
[3] T. G. Theofanous, « The boiling crisis in nuclear reactor safety and performance », International Journal of Multiphase Flow, vol. 6, no 1, p. 69-95, févr. 1980, doi: 10.1016/0301-9322(80)90040-3.
[4] N. Takenaka, T. Fujii, et others, « Flow pattern transition and heat transfer of inverted annular flow », Int. J. Multiphase Flow, 1989.
[5] M. A. El Nakla, D. C. Groeneveld, et S. C. Cheng, « Experimental study of inverted annular film boiling in a vertical tube cooled by R-134a », International Journal of Multiphase Flow, vol. 37, p. 37-75, 2011.
[6] Q. Liu, J. Kelly, et X. Sun, « Study on interfacial friction in the inverted annular film boiling regime », Nuclear Engineering and Design, vol. 375, 2021.
[7] K. K. Fung, « Subcooled and low quality film boiling of water in vertical flow at atmospheric pressure », PhD Thesis, Argonne National Laboratory, 1981.
Combustion d'hydrogène et d'ammoniac en milieux poreux : expériences et modélisation
- Contexte
Les perspectives énergétiques actuelles suggèrent l'utilisation de l'hydrogène (H2) et de l'ammoniac (NH3) comme vecteurs d'énergie décarbonés. La combustion du NH3 offre des avantages tels qu'une densité énergétique élevée et un stockage sûr, mais présente une plage d’inflammabilité étroite et des émissions élevées de NOx. Il est possible d'obtenir de l'hydrogène par craquage partiel d’ammoniac pour créer des mélanges ayant des propriétés de combustion plus favorables, mais il reste des questions ouvertes concernant les émissions de polluants et la teneur en NH3 imbrûlé.
- Défis
Les brûleurs poreux sont des candidats prometteurs pour la combustion de mélanges NH3/H2 à faibles émissions polluantes. Malheureusement, les problèmes de durabilité des matériaux et la complexité de la stabilisation des flammes constituent encore des obstacles importants à leur industrialisation. Toutefois, les récentes avancées dans le domaine de la fabrication additive permettent un design avancé de matrices poreuses, leur caractérisation expérimentale restant difficile en raison de l'opacité de la matrice solide.
- Objectifs de recherche
Le doctorant exploitera un banc expérimental au CEA Saclay pour mener des expériences de combustion avec des mélanges NH3/H2/N2+air dans différents brûleurs poreux. Les tâches principales incluront la conception de nouvelles géométries poreuses, la comparaison des résultats expérimentaux avec les simulations numériques, un travail de modélisation 1D par moyennes volumiques et théorie asymptotique. Les mesures expérimentales comprendront : l'anémométrie à fil chaud, la thermométrie infrarouge, l'analyse de la composition des gaz de sortie, la chimiluminescence et les diagnostics laser. Les brûleurs poreux seront fabriqués à l'aide de techniques d'impression 3D avec des matériaux tels que l'acier inoxydable, l'inconel, l'alumine, la zircone et le carbure de silicium.
La recherche vise à développer des brûleurs poreux plus robustes et plus efficaces pour la combustion de mélanges NH3/H2, améliorant ainsi leur application pratique pour atteindre la neutralité carbone. Le candidat contribuera à faire progresser le domaine grâce à des données expérimentales, des conceptions innovantes et des techniques de modélisation améliorées.
Modélisation de la chute de gouttes dans un volume libre, en support au code système CATHARE
Cette thèse porte sur l'étude de la chute de gouttes dans des volumes libres, dans le cadre de l'amélioration continue des modèles physiques du code CATHARE, utilisé pour les études de sûreté des Réacteurs à Eau Pressurisée. Les modèles actuels reposent sur les travaux d'Ishii et Zuber, qui modélisent la vitesse de chute des gouttes dans un fluide diphasique. L'objectif de la thèse est de raffiner la précision de ce modèle en y intégrant des paramètres supplémentaires et en le validant grâce à des expériences telles que celles de Dampierre et CARAYDAS. Le doctorant devra concevoir un modèle mécaniste plus représentatif, fondé sur des données expérimentales ou des simulations CFD si nécessaire. L'innovation réside dans le développement d'une modélisation plus fidèle des processus de chute de gouttes, ouvrant la voie à des applications spécifiques, telles que la modélisation des sprays, et contribuant ainsi à la validation du code CATHARE dans des domaines supplémentaire.