Numerical model for pellet rocket acceleration in PELOTON

Cet article présente un modèle numérique validé au sein du code PELOTON qui simule l'accélération des roquettes de granulés dans les dispositifs de fusion thermonucléaire en tenant compte de l'asymétrie du nuage d'ablation et des gradients de plasma, démontrant une cohérence avec les trajectoires expérimentales de JET et révélant une déviation réduite pour les granulés composites deutérium-néon.

L'essentiel

Imaginez une minuscule boule de neige super froide (un « granulé ») projetée dans un four géant et tourbillonnant de gaz super chaud (plasma) à l'intérieur d'un réacteur de fusion. Ce n'est pas une simple collision ; c'est une danse à haute vitesse où la boule de neige tente de survivre tandis que le four tente de la faire fondre.

Ce document décrit une nouvelle simulation informatique appelée PELOTON qui agit comme un réalisateur de film en haute définition pour cette danse. Sa mission principale est de comprendre pourquoi ces boules de neige ne fondent pas simplement en ligne droite, mais sont plutôt poussées sur le côté, accélérant comme une fusée.

Voici la décomposition de ce que le document a découvert, en utilisant des analogies simples :

1. L'effet « Fusée » : Pourquoi la boule de neige se déplace-t-elle sur le côté ?

Habituellement, si vous soufflez de l'air sur un ballon, cela le pousse loin de vous. Mais ici, l'« air » est en réalité un flux d'électrons invisibles et super rapides provenant du plasma chaud.

• L'installation : À mesure que la boule de neige entre dans le four, elle commence à fondre, créant un épais nuage de gaz froid autour d'elle.
• Le rebondissement : Le four possède un champ magnétique qui est plus fort d'un côté (le « côté haute intensité de champ » ou HFS) et plus faible de l'autre (le « côté basse intensité de champ » ou LFS).
• L'analogie : Imaginez la boule de neige comme une personne debout dans une foule. D'un côté (HFS), la foule est dense et chaotique, ce qui rend difficile pour la « chaleur » (les électrons) d'atteindre la personne. De l'autre côté (LFS), la foule est plus clairsemée, donc la chaleur frappe la personne plus fort.
• Le résultat : Parce que la chaleur frappe le LFS plus fort, le nuage de gaz de ce côté devient plus chaud et pousse plus fort en retour. Cela crée une différence de pression. La boule de neige est compressée du côté chaud et poussée vers le côté froid. C'est comme une fusée poussée par son échappement, mais à l'envers : la pression derrière elle (sur le LFS) est plus élevée que la pression devant elle, la poussant sur le côté.

2. Le modèle informatique : PELOTON

Les auteurs ont construit une simulation 3D pour suivre cela. Considérez PELOTON comme une prévision météorologique ultra-précise pour l'intérieur du réacteur.

• Il suit la boule de neige à mesure qu'elle fond.
• Il calcule comment le nuage de gaz froid se forme et se déplace.
• Il tient compte du fait que le nuage n'est pas uniforme ; il est « chargé » différemment d'un côté à l'autre, ce qui modifie la façon dont les électrons chauds le frappent.
• Ils ont testé ce modèle par rapport à de vraies expériences au JET (un célèbre laboratoire de fusion au Royaume-Uni) et ont constaté que les prédictions de l'ordinateur correspondaient presque parfaitement aux trajectoires réelles des boules de neige.

3. La boule de neige « éclatée » (SPI)

Parfois, au lieu d'une seule grosse boule de neige, on projette une « boule de neige éclatée » (SPI). Imaginez lancer une poignée de copeaux de glace au lieu d'un bloc entier.

• Le chevauchement des nuages : Si deux copeaux de glace sont proches, leurs nuages de gaz peuvent s'entrechoquer. Le document a découvert que si deux copeaux sont côte à côte, celui du bas est poussé plus fort. S'ils sont alignés l'un derrière l'autre sur le même chemin magnétique, ils se tirent l'un vers l'autre car le premier bloque la chaleur qui frapperait l'arrière du second.
• Le mélange de Néon : Ils ont essayé d'ajouter un peu de gaz néon (comme une saveur différente de glace) à la boule de neige. Cela a rendu le nuage de gaz plus frais et plus lent. Bien que la « poussée de la fusée » ait toujours eu lieu, elle était plus faible. Curieusement, dans les expériences réelles, cela ne semblait pas changer beaucoup la trajectoire, probablement parce que le néon a provoqué d'autres changements majeurs dans le plasma qui ont masqué l'effet.

4. La « Loi d'échelle » : Une recette pour la prédiction

L'équipe a analysé des centaines de simulations pour créer une « recette » simple (une loi d'échelle).

• La recette : La force de la poussée latérale dépend principalement de la température et de la densité du plasma.
• La surprise : La taille de la boule de neige (le rayon du granulé) importe très peu ! Un minuscule copeau et un gros morceau sont poussés avec environ la même force par unité de masse. C'est une simplification énorme pour les scientifiques qui tentent de prédire le comportement de ces granulés.

5. Ce que cela signifie pour l'avenir

Le document conclut que ce modèle est prêt à être utilisé pour la prochaine machine de fusion géante, ITER.

• Ils prévoient d'utiliser cette « physique de la fusée » pour prédire comment les granulés éclatés se comporteront dans le plasma massif d'ITER.
• Ils veulent affiner le modèle pour inclure la façon dont les particules du plasma se dispersent (diffusion) afin de rendre les prédictions encore plus précises.

En résumé : Le document explique que lorsque les granulés froids fondent dans un réacteur de fusion, ils sont poussés sur le côté par un « vent » invisible de chaleur qui les frappe de manière inégale. Les auteurs ont construit un modèle informatique qui prédit cette poussée parfaitement, montant que la taille du granulé importe peu, mais que la température et la densité du plasma sont cruciales. Cela aide les scientifiques à comprendre comment injecter du carburant en toute sécurité dans les futures centrales à fusion.

Résumé technique

Résumé technique : Modèle numérique de l'accélération par effet fusée des pellets dans PELOTON

Énoncé du problème
L'injection de pellets cryogéniques ou de fragments d'injection de pellets fragmentés (SPI) dans les dispositifs de fusion thermonucléaire est sujette à une « accélération par effet fusée », un phénomène où les fragments dévient de leurs trajectoires prévues. Cette accélération provient du chauffage asymétrique du nuage d'ablation entourant le pellet. Alors que les électrons du plasma chaud circulent le long des lignes de champ magnétique, le nuage d'ablation s'ionise et dérive à travers les lignes de champ en raison de la courbure du champ magnétique (dérive grad-B). Cette dérive crée une asymétrie : le côté haut champ (HFS) du nuage offre un blindage électrostatique moindre et des chemins plus courts pour les électrons entrants par rapport au côté bas champ (LFS). Par conséquent, le HFS subit un dépôt d'énergie électronique et une pression plus élevés, accélérant le pellet vers le LFS. Bien que observé expérimentalement, les modèles existants traitent souvent la différence de pression motrice comme un paramètre ajustable ou reposent sur des approximations semi-analytiques qui peuvent manquer de réalisme dans les termes d'ordre supérieur.

Méthodologie
Les auteurs ont développé un modèle de simulation numérique directe pour l'accélération par effet fusée des pellets, implémenté dans PELOTON, un code de particules lagrangiennes en 3D. La méthodologie intègre plusieurs composantes physiques clés :

1. Charge non uniforme du nuage : Contrairement aux versions précédentes de PELOTON qui supposaient un potentiel de blindage électrostatique constant, ce modèle résout un potentiel de blindage non uniforme ($\Phi$) à travers le nuage d'ablation. Ceci est réalisé en résolvant une équation non linéaire pour le potentiel de blindage le long des lignes de champ magnétique, en tenant compte de l'effet d'albédo et de la rétrodiffusion collisionnelle. Cette variation conduit à des densités d'électrons effectives ($n_{eff}$) et des dépôts de chaleur différents sur le HFS par rapport au LFS.
2. Physique MHD et d'ablation : L'évolution du nuage d'ablation froid et partiellement ionisé est régie par des équations de magnétohydrodynamique (MHD) à faible nombre de Reynolds magnétique dans des coordonnées lagrangiennes. Le modèle inclut :
   • Des conditions limites d'ablation de surface où le flux de chaleur pilote la sublimation.
   • Une dynamique de dérive grad-B, modélisée via une équation d'accélération de dérive transversale simplifiée pilotée par l'intégrale de la différence de pression.
   • Un modèle avancé de dépôt de chaleur électronique basé sur des solutions de l'équation cinétique de Fokker-Planck linéarisée en 3D, incorporant l'opacité et les effets de blindage.
   • Un solveur d'équation d'état multi-espèces pour les mélanges deutérium-néon, utilisant les équations de Saha pour déterminer les fractions d'ionisation, la température et la vitesse du son.
3. Module de refroidissement du plasma : Un nouveau module a été développé pour fournir des états de plasma de fond évolutifs lors des simulations SPI. Ce module discrétise le plasma en volumes de flux, distribue le matériau ablaté et prend en compte les pertes d'ionisation, le chauffage ohmique et l'échange d'énergie ion-électron. Il a été validé par rapport aux simulations JOREK et INDEX.
4. Validation et simulation : Le modèle a été appliqué pour simuler des expériences SPI de deutérium dans le tokamak JET. Les simulations comprenaient des mélanges de deutérium pur et de deutérium-néon (0,5 % de néon). Le code a suivi les trajectoires des fragments, les taux d'ablation et les forces d'accélération par effet fusée, comparant les résultats aux données expérimentales (spécifiquement JET #96874).

Contributions clés et résultats

• Simulation numérique directe de l'accélération par effet fusée : L'article présente un modèle numérique de premier principe qui calcule directement la force d'accélération par effet fusée à partir de la différence de pression ($P_{HFS} - P_{LFS}$) générée par le chauffage électronique asymétrique, plutôt que de s'appuyer sur des paramètres d'ajustement empiriques.
• Loi d'échelle pour JET : Sur la base de 20 points de données de simulation couvrant les paramètres de plasma pertinents pour JET ($T_e \in [450, 2000]$ eV, $n_e \in [6\times10^{19}, 3\times10^{20}]$ m$^{-3}$), les auteurs ont dérivé une loi d'échelle pour la différence de pression pilotant l'accélération :
  $$\Delta P [\text{bar}] = 1.0 \left( \frac{T_e [\text{keV}]}{2} \right)^{1.26} \left( \frac{n_e [1/m^3]}{10^{20}} \right)^{0.4}$$
  L'ajustement montre une erreur absolue relative moyenne d'environ 13 %. L'étude a constaté que $\Delta P$ est pratiquement indépendant du rayon du pellet dans la plage testée.
• Validation avec les expériences JET : Les simulations PELOTON des trajectoires SPI de deutérium, intégrant l'accélération par effet fusée calculée, ont montré une cohérence avec les trajectoires mesurées expérimentalement dans JET, particulièrement concernant la déviation vers le côté bas champ.
• Effets du mélange de néon : Les simulations de pellets deutérium-néon (0,5 % de néon) ont révélé que la composante de néon réduit la température du nuage d'ablation, la vitesse d'expansion longitudinale et la vitesse de dérive grad-B. Bien que cela se traduise par une accélération par effet fusée plus faible mais toujours mesurable par rapport au deutérium pur, les auteurs notent que dans les expériences réelles de JET, l'effet de fusée sur les trajectoires de néon est probablement masqué par les événements MHD globaux déclenchés par le rayonnement du néon.
• Interactions entre fragments : L'étude a quantifié l'impact du chevauchement des nuages. Les fragments séparés transversalement de 2 à 4 cm ont montré une accélération modifiée due à la redistribution de la pression. Les fragments alignés le long de la même ligne de champ magnétique ont subi une forte attraction mutuelle due au blindage partiel du flux de chaleur électronique, bien que les auteurs notent que la probabilité d'un tel alignement dans les panaches SPI réels soit faible.
• Gradients de plasma : Le modèle a démontré que des gradients de plasma locaux forts peuvent amplifier l'accélération par effet fusée d'un facteur d'environ 3, mais que ces effets sont transitoires, n'agissant que pendant de courts intervalles lorsque les fragments traversent les limites de dépôt.

Signification et revendications
L'article affirme fournir un cadre numérique complet qui capture la physique essentielle de l'accélération par effet fusée des pellets, spécifiquement le rôle de la charge électrostatique non uniforme et de la dérive grad-B. Les auteurs soulignent que leur modèle va au-delà des approches semi-empiriques en résolvant directement l'asymétrie de pression.

La signification de ce travail réside dans sa capacité à :

1. Reproduire les trajectoires SPI expérimentales dans JET avec une fidélité quantitative.
2. Fournir une loi d'échelle prédictive pour la différence de pression pilotant l'accélération par effet fusée, laquelle peut être utilisée dans d'autres codes.
3. Démontrer que, bien que l'accélération par effet fusée soit un mécanisme physique significatif, son impact observable sur les trajectoires peut être modulé par la composition du pellet (par exemple, l'ajout de néon) et les gradients de plasma de fond.

Les auteurs restent modestes quant à l'application immédiate à ITER, affirmant que des travaux futurs sont nécessaires pour développer des lois d'échelle spécifiques pour les plasmas d'ITER et pour affiner le module de refroidissement afin d'inclure la diffusion du plasma et les effets MHD. Ils précisent explicitement que la loi d'échelle rapportée est spécifique à l'espace des paramètres et à la géométrie magnétique de JET.