L'instabilité de dérive non-résonnante dans les plasmas turbulents
La turbulence magnétique présente dans de nombreux systèmes astrophysiques tels que le vent solaire et les rémanents de supernova joue un rôle clé dans l'accélération de particules à haute énergie, notamment dans le contexte des ondes de choc non collisionnelles. En confinant certaines particules près du front de choc, elle leur permet de gagner de l'énergie par réflexions successives entre les milieux amont et aval, un mécanisme connu sous le nom de processus de Fermi et considéré comme étant à l'origine des rayons cosmiques.
La turbulence autour des rémanents de supernova est probablement déclenchée par les rayons cosmiques eux-mêmes, via des instabilités plasma lors de leur propagation en amont du choc. Dans le cas d'une onde de choc se déplaçant parallèlement au champ magnétique ambiant, l'instabilité dominante est celle de "dérive non-résonnante" - ou instabilité de Bell -, qui contribue à renforcer la turbulence préexistante.
L'objectif de cette thèse est de développer un modèle analytique complet de cette instabilité dans un plasma turbulent, puis d'en comparer les résultats à des simulations numériques avancées.
Effets électroniques dans les cascades de collisions dans le GaN
Dans les environnements radiatifs tels que l'espace et les installations nucléaires, les composants microélectroniques sont soumis à des
flux intenses de particules qui détériorent leur fonctionnement en dégradant les matériaux les constituant. Les particules entrent en
collision avec des atomes dans les matériaux semi-conducteurs, leur cèdent une partie de leur énergie cinétique et les éjectent de
leur site cristallin. Les atomes éjectés vont à leur tour générer des collisions, formant une cascade de collisions qui conduira à la
création de défauts de déplacements. De plus, les particules chargées primaires ou secondaires (issues de l’interaction avec un
neutron par exemple) vont aussi interagir spécifiquement avec les électrons du réseau et leur céder une partie de leur énergie en générant des paires électron-trou. On parle de freinage électronique. Une simulation complète de cascade de collisions se doit donc d’intégrer
ces deux éléments : collisions avec les noyaux des atomes et effets électroniques.
La méthode de prédilection pour la simulation de cascades de collisions à l’échelle atomique est la dynamique moléculaire (DM).
Cependant, les effets électroniques ne sont pas inclus car la méthode ne traite pas explicitement les électrons. Pour pallier à ce
problème, des modules additionnels à la DM imitant le plus fidèlement possible les effets des électrons doivent être utilisés. L’état de
l’art en ce qui concerne la simulation du freinage électronique d’un projectile dans un solide est la méthode de la théorie de la
fonctionnelle de la densité dépendante du temps (TDDFT). L’objectif de cette thèse est de combiner DM et TDDFT pour réaliser des
simulations de cascades de collisions dans le GaN et étudier l’influence des effets électroniques. En plus de compétences transverses
communes à toute thèse, le/la candidat.e sera amené.e à développer des compétences dans plusieurs méthodes de modélisation
à l’échelle atomique, en physique du solide, en interactions particules-matière, en environnement linux ainsi qu’en programmation.
Mesure de la vitesse du son dans H2, He constitutifs des intérieurs des géantes gazeuses
L'objectif de la thèse est d'étudier les mélanges hydrogène-hélium en phase fluide à haute pression et haute température par spectroscopie Raman et Brillouin. Les expériences seront effectuées en cellule a enclumes diamant sous chauffage laser permettant d’explorer un vaste domaine de pression et de température représentatif des intérieurs planétaires des géantes de gaz (1-300 GPa, 300-4000 K). La spectroscopie Raman sera utiliser pour sonder les changements chimiques susceptibles d’apparaître en conditions extrêmes. La spectroscopie Brillouin donnera accès à la vitesse du son adiabatique et aux équations d’état de ces systèmes en phase fluide. Ces données seront particulièrement utiles pour améliorer la modélisation des intérieurs de Jupiter et Saturne.
Étude numérique de la sensibilité de la transition à la turbulence des couches-limites hypersoniques aux modèles de gaz pour en conditions de vol.
Cette thèse vise à améliorer la compréhension des écoulements hypersoniques et notamment le phénomène de transition laminaire-turbulent. Pour assurer une conception correcte des véhicules hypersoniques, il est important de prédire précisément les flux thermiques à la paroi. Ces flux sont très fortement contraints par la nature de la couche limite (laminaire/ transitionnelle/ turbulente). Or, les mécanismes à l’origine de la transition laminaire-turbulent sont complexes et souvent mal compris.
Ce constat est d’autant plus vrai lorsque la complexité thermochimique de l’écoulement augmente, notamment pour des hauts nombres de Mach, où les effets dits de gaz réel tels que la dissociation ou le non-équilibre thermochimique ne sont plus négligeables. Cependant, devant certaines difficultés de simulation de ces effets, la plupart des études sont aujourd’hui encore effectuée en gaz parfait, même pour des nombres de Mach intermédiaire (7 par exemple). Ainsi, le rôle des effets thermochimiques sur les étapes successives de la transition reste encore mal connu.
Dans ce contexte, cette thèse vise à étudier numériquement l’impact de la montée en complexité des modèles de gaz sur la transition laminaire-turbulent. L’objectif final de la thèse est d’être en capacité d’effectuer des simulations numériques haute-fidélité de transition laminaire-turbulente en conditions de vol, et de mieux maitriser les modèles physiques et les paramètres de similitudes nécessaires pour restituer les essais en vol.
Modélisation multiéchelle de la réponse magnétique de matériaux hétérogènes
La dépendance spectrale de la perméabilité des matériaux magnétiques, que ce soit dans les matériaux composites ou massifs, reste un sujet complexe, en raison des différentes échelles des phénomènes impliqués. Des modèles analytiques approximatifs sont souvent utilisés pour décrire la réponse en fréquence des matériaux magnétiques, notamment pour améliorer leurs performances dans des domaines comme l’électronique de puissance. Des résultats récents ont montré que des codes de micro-magnétisme permettent maintenant de prédire la réponse d’un ensemble de nanoparticules couplées, ou d’une particule d’un volume représentatif des matériaux en question. Cette thèse vise à utiliser ces outils pour améliorer les modèles analytiques existants. Une inclusion baignant dans un champ effectif sera le paradigme à partir duquel la structure en domaine et la réponse spectrale de la particule seront calculées en utilisant un code de micro-magnétisme. Les matériaux étudiés incluent des particules sphériques ou à fort rapport de forme (oxydes magnétiques, pétales ferromagnétiques) à concentration variable, allant des milieux dilués aux matériaux massifs. Des pistes seront ainsi dégagées pour optimiser la microstructure des matériaux, en vue de meilleures performances dans des applications comme l’électronique de puissance et les composants hyperfréquences. A cet effet, le CEA offre un environnement de calcul scientifique avec des ressources HPC, ainsi que des capacités pour l’élaboration d'échantillons et les caractérisations magnétiques statiques et dynamiques. A l’issue de ce travail, le candidat aura acquis une très bonne maîtrise des relations microstructure-propriétés décrites par une approche numérique appliquées aux matériaux magnétiques. Plus généralement, cette démarche est en pleine expansion dans le domaine des matériaux ("materials par design", ou conception numérique des matériaux).
Description de phénomènes collectifs au sein des noyaux d’atome au-delà de l’approche de la fonctionnelle de la densité dépendant du temps
Contexte :
Prédire l'organisation et la dynamique des neutrons et protons au sein des noyaux d'atome est un véritable
défi scientifique, crucial pour concevoir les technologiques nucléaires du futur mais aussi pour répondre à
des questions fondamentales comme l'origine des atomes lourds dans notre univers. Dans ce cadre, le CEA,
DAM, DIF développe des approches théoriques visant à simuler la dynamique des constituants élémentaires
du noyau atomique. Les équations du mouvement obtenues, dans le cadre de la mécanique quantique, sont
résolues sur nos supercalculateurs. Les années 2010 ont vu une montée en puissance de l’approche dite de la
fonctionnelle de la densité dépendant du temps (TDDFT) pour traiter ce problème. Malgré sa description
révolutionnaire de certains phénomènes tels que les résonances géantes observées dans les noyaux d’atome
ou encore la fission nucléaire, cette approximation possède des limites intrinsèques.
Objectifs :
Cette thèse vise à développer et explorer une nouvelle approche théorique permettant de décrire les
phénomènes de mouvements collectifs des protons et neutrons au sein du noyau atomique. L’idée est de
généraliser l’approche TDDFT afin de mieux prédire certaines propriétés des réactions nucléaires comme la
répartition de l’énergie entre les fragments issus d’une fission nucléaire. Partant de premier travaux allant
dans cette direction, le/la doctorante devra dériver les équations du mouvement de cette nouvelle approche
puis les implémenter sous forme d’une librairie C++ optimisée et tirant profit des supercalculateurs du CEA.
Le but sera enfin d’étudier comment cette nouvelle approche améliore les prédictions de certains
phénomènes comme l’amortissement des résonances géantes dans les noyaux d’atomes ou encore la
naissance des fragments générés lors d’une fission nucléaire.
Assimilation de données hétérogènes dans des simulations de dispersion atmosphérique de radionucléides à échelle régionale
La modélisation et la simulation apportent des connaissances essentielles à la dispersion aérienne de gaz et de particules et au marquage de l’environnement qui en résulte. Ceci s’applique notamment aux rejets qui ont été engendrés par les essais nucléaires atmosphériques effectués dans le passé par la France en Polynésie. Si les calculs météorologiques et de dispersion à l’échelle régionale sont raisonnablement fiables, leurs résultats ont une part d’incertitude et présentent des écarts aux mesures hétérogènes des activités ou débits de doses dans l’air, sur le sol et dans les compartiments biologiques. La thèse visera à développer des méthodes d’inversion, basées sur l’assimilation de données, afin de réduire les erreurs et incertitudes des simulations de dispersion régionale de radionucléides. L’application concernera certains essais nucléaires dans l’atmosphère. Toutefois, les méthodes développées au cours de la thèse, telles que l’échantillonnage de Monte-Carlo par des chaînes de Markov, auront un domaine de mise en œuvre plus général. Après une revue bibliographique portant sur les essais nucléaires et les méthodes d’assimilation de données, des algorithmes originaux de modélisation inverse seront programmés, testés et appliqués à la simulation de la dispersion des rejets aériens issus d’essais. Ceci permettra d’estimer le rôle pressenti important de l’assimilation des mesures pour améliorer les simulations.
Etude du comportement d'un composite CMC en température par essais in situ en tomographie X
Le sujet proposé concerne l’étude du comportement mécanique d’un matériau composites à matrice céramique de type oxyde/oxyde en température (jusqu’à 1000°C). L’originalité du sujet est l’utilisation de la tomographie X in situ pour accéder d’une part à la déformation macroscopique des éprouvettes testées et d’autre part aux mécanismes d’endommagement à l’échelle microscopique qui caractérisent ce type de matériaux dit « endommageables ».
Cette technique a été développée à température ambiante lors d'une thèse précédente : il s'agit ici de l'appliquer en température et sur des sollicitations plus complexes (eg traction-torsion). Il s'agira également de proposer des développements au protocole d'analyse par
corrélation d'image volumique existant.
Potentiels réactifs par réseaux de neurones : optimisation de l'acquisition des données d'entraînement et application aux réactions mécanochimiques
La décomposition spontanée de molécules organiques lors de leur synthèse, manipulation, ou stockage, pose de sérieux problème de sécurité dans le cas de matériaux énergétiques. En plus de l'activation thermique, de récentes études montrent que les déformations intramoléculaires telles que celles induites par le passage d'une onde de choc influencent la réactivité chimique et peuvent modifier les mécanismes de décomposition. Les études à l'échelle moléculaire de ces phénomènes représentent un défi car elles nécessitent à la fois des calculs de chimie quantique pour décrire la formation et la rupture de liaison chimique, mais aussi des effets de phase condensée.
Pour remédier à cela, nous proposons le développement et l'application d'un potentiel réactif par réseaux de neurones (MLIP) pouvant combiner ces deux aspects. En particulier, nous visons à avancer significativement la méthodologie de construction de base de données contenant des structures moléculaires hors-équilibre, telles que celles pouvant être rencontrées lors de réactions chimiques complexes activées thermiquement ou mécaniquement. Ce potentiel réactif sera ensuite utilisé pour étudier la décomposition d'un matériaux modèle dans différentes conditions de pression et de température. Les outils et connaissances développés serviront à la fois aux études de mécanismes de décomposition de molécules énergétiques mais aussi à la communauté de mécanochimie.
Description microscopique des propriétés des fragments de fission évaluées à la scission
La fission est l’une des réactions nucléaires les plus difficiles à décrire, reflétant la diversité des aspects dynamiques du problème à N-corps. Au cours de ce processus, le noyau explore des états de déformations extrêmes aboutissant à la formation de deux fragments. Alors que la nombre de degrés de liberté (DLC) mis en jeu est extrêmement grand, l’approximation de champ moyen est un bon point de départ qui opère une réduction drastique des DLC, l’élongation et l’asymétrie étant incontournables. Cette réduction introduit des discontinuités dans la génération successive des états par lesquels le noyau transite, la continuité en énergie n’assurant pas la continuité des états issus d’un principe variationnel. Récemment, une nouvelle méthode basée sur des contraintes associées aux recouvrements de fonctions d’onde, a été mise en oeuvre afin d’assurer cette continuité jusqu’à la scission et au-delà (vallée de Coulomb). Cette continuité est capitale pour décrire la dynamique du processus.
L’objectif de la thèse que nous proposons est de réaliser pour la première fois la mise en oeuvre à deux dimensions de cette nouvelle approche afin de prendre en compte l’ensemble de la collectivité générée par les DLC d’élongation et d’asymétrie. Les développements théoriques et numériques s’inscrivent dans le cadre de la méthode de la coordonnée génératice dépendant du temps. Ce type d’approche contient une première étape statique qui consiste en la génération de surfaces d’énergie potentielle (PES) obtenues par des calculs Hartree-Fock-Bogoliubov sous contraintes et une seconde étape dynamique qui décrit la propagation dynamique d’un paquet d’onde sur ces surfaces via la résolution de l’équation de Schrödinger dépendant du temps. C’est à partir de cette deuxième étape que les observables sont en général extraites.
Dans le cadre de cette thèse, le doctorant devra:
- en première étape, construire des PES à deux dimensions continues pour l’état adiabatique et les premiers états excités. Cela mettra en oeuvre les trois algorithmes Link, drop et Deflation
- en deuxième étape, extraire des observables accessibles par ce type d’approches: les rendements, le bilan énergétique à la scission, la déformation des fragments, le nombre moyen de neutrons émis. Nous souhaitons en particulier étudier l’impact des excitations intrinsèques sur les observables de fission, qui se manifestent essentiellement dans la descente du point scelle vers la scission.
Enfin, ces résulats seront confrontés aux données expérimentales, dans des actinides et des pré-actinides d’intérêt. Notamment, les mesures récentes très précises obtenues par les expériences SOFIA pour des noyaux peu à très exotiques devraient contribuer à tester la précision et la prédictivité de nos approches, et à guider les futurs développements des approches à N-corps et de l’interaction nucléaire dans le cadre de la fission.