Lithium Droplet Transport in Tokamak Edge Plasmas

Ce papier présente un modèle de transport et d'évaporation de gouttelettes de lithium intégré au code OpenEdge, validé avec une grande précision et appliqué au concept de réacteur CAT, démontrant que la trajectoire et la survie des gouttelettes dépendent de leur taille, vitesse et position de lancement, tout en permettant une évaluation auto-cohérente de leur impact sur le plasma de bord via un couplage bidirectionnel avec SOLPS-ITER.

L'essentiel

🌟 Le Grand Voyage des Gouttes de Lithium dans le Soleil Artificiel

Imaginez que vous essayez de construire une centrale électrique capable de reproduire l'énergie du Soleil : c'est le but d'un tokamak (un réacteur à fusion). Le problème ? La chaleur est si intense qu'elle pourrait fondre n'importe quel matériau solide.

Pour résoudre ce problème, les scientifiques ont une idée brillante : utiliser du lithium liquide. C'est comme avoir un "bouclier liquide" qui se répare tout seul. Mais attention, sous l'effet de la chaleur et des champs magnétiques, ce liquide ne reste pas calme : il se met à éclabousser, projetant des milliers de gouttelettes dans le plasma (le gaz surchauffé).

Ce papier raconte l'histoire de ces gouttelettes. Les chercheurs ont créé un simulateur informatique (appelé OpenEdge) pour comprendre ce qu'elles deviennent une fois lancées dans la fournaise.

---

🎈 Les Gouttelettes : Des Ballons de Savon dans un Cyclone

Pour comprendre le comportement de ces gouttes, imaginez-les comme de minuscules ballons de baudruche traversant un cyclone violent. Trois forces principales les guident :

1. La Gravité (Le poids) : Comme une pierre qui tombe, la gravité tire les gouttes vers le bas.
2. Le Freinage (La résistance de l'air) : Le plasma est très dense. C'est comme essayer de courir dans de l'eau épaisse : les gouttes sont freinées par les particules du plasma.
3. Le "Moteur Fusée" (La force de recul) : C'est le point le plus fascinant ! Quand une goutte de lithium chauffe, elle s'évapore. Mais dans un tokamak, la chaleur arrive surtout d'un seul côté (comme un four à micro-ondes qui chauffe une seule face de votre pop-corn).
   • L'analogie : Imaginez une fusée qui perd du carburant. Si le carburant s'échappe par l'arrière, la fusée avance. Ici, si le lithium s'évapore d'un seul côté de la goutte, cela crée une poussée dans la direction opposée. C'est ce qu'on appelle la force de fusée. Cela peut faire changer de cap à la goutte, la propulsant parfois vers l'extérieur au lieu de tomber.

🎯 Le Jeu de la "Grosse ou la Petite"

Les chercheurs ont lancé des millions de gouttes virtuelles de différentes tailles pour voir où elles atterrissent. C'est ici que la taille fait toute la différence :

• Les petites gouttes (comme des grains de sable) : Elles sont si légères et ont une surface énorme par rapport à leur masse. Elles s'évaporent presque instantanément, comme une goutte d'eau sur un fer à repasser. Elles ne vont pas loin et disparaissent avant d'atteindre le cœur du réacteur.
• Les grosses gouttes (comme des billes) : Elles sont plus robustes. Elles résistent à la chaleur, gardent leur masse et voyagent plus loin. Elles finissent souvent par se déposer sur les parois voisines, comme de la pluie sur un pare-brise.

Le résultat clé : Seules les toutes petites gouttes (ou les débris de gouttes éclatées) réussissent à pénétrer profondément dans le cœur du réacteur. Les grosses gouttes restent près de leur point de départ.

🤝 La Danse à Deux Pas : Le Dialogue entre les Codes

Pour être précis, les scientifiques ont dû faire dialoguer deux programmes informatiques :

1. OpenEdge : Qui suit le destin de chaque gouttelette (où elle va, combien elle perd de poids).
2. SOLPS-ITER : Qui simule le comportement du plasma (le gaz chaud).

C'est comme une conversation en boucle :

• Le programme des gouttes dit au programme du plasma : "Hé, j'ai évaporé tant de lithium ici, ça va refroidir ton gaz !".
• Le programme du plasma répond : "OK, le gaz est maintenant plus froid et plus dense, voici comment cela change ton trajet.".
• Le programme des gouttes ajuste sa trajectoire en conséquence.

Cette boucle permet de voir comment les gouttes modifient le plasma, et comment le plasma modifie le sort des gouttes, en temps réel.

🏁 La Conclusion du Voyage

Ce travail est crucial pour la conception des futurs réacteurs à fusion (comme le projet CAT mentionné). Il nous apprend que :

• Si on veut que le lithium aille au cœur du réacteur pour le refroidir, il faut des très petites gouttes.
• Si on veut que le lithium reste près des parois pour les protéger, il faut des gouttes plus grosses.
• La "force de fusée" due à l'évaporation est un facteur clé qu'on ne peut pas ignorer : elle peut faire rebondir les gouttes vers l'extérieur au lieu de les laisser tomber.

En résumé, ce papier est la carte routière qui permet aux ingénieurs de savoir exactement comment gérer ces gouttes de lithium pour que notre futur réacteur nucléaire fonctionne de manière sûre et efficace. C'est de la physique de haute précision appliquée à un problème de survie thermique ! 🔥⚛️

Résumé technique

1. Problématique et Contexte

Les composants face au plasma (PFC) des tokamaks de réacteurs doivent résister à des flux de chaleur et de particules extrêmes. Les concepts de limiteurs et de diverteurs en lithium liquide offrent des avantages prometteurs, notamment des surfaces auto-cicatrisantes, une réduction des pertes radiatives et un blindage par la vapeur. Cependant, des campagnes expérimentales (CDX-U, NSTX, EAST) ont révélé un mécanisme limitant : l'instabilité de la surface liquide et les éclaboussures.

Sous l'effet de forces transitoires et magnétohydrodynamiques (MHD), des gouttelettes de lithium sont éjectées dans le plasma de bord. Ces gouttelettes se fragmentent, s'évaporent et modifient les conditions locales du plasma. Les modèles de source de plasma de bord conventionnels ne prennent pas explicitement en compte le transport de ces gouttelettes, pourtant essentielles pour prédire si le lithium reste localisé au diverteur, se redépose sur les tuiles adjacentes ou pénètre dans le cœur du plasma.

2. Méthodologie

L'étude développe un modèle de transport et d'évaporation de gouttelettes de lithium intégré dans le code de simulation Monte Carlo Directe (DSMC) OpenEdge.

A. Modèle Physique de la Gouttelette
Chaque gouttelette est modélisée comme une sphère caractérisée par son rayon ($r_d$), sa masse ($m_d$) et sa température ($T_d$). Les équations du mouvement incluent :

• Forces mécaniques : Gravité, traînée collisionnelle ionique (composantes de collecte directe et diffusion de Coulomb), et force de Lorentz (négligeable pour les gouttelettes de taille millimétrique par rapport à la gravité et la traînée).
• Chargement électrique : Modèle de limite orbitale (OML) pour le potentiel flottant, incluant l'émission thermoionique (Richardson-Dushman) pour les métaux liquides chauds.
• Évaporation : Un modèle d'équilibre énergétique couple la perte de masse et l'évolution thermique. Le flux de chaleur incident provient du plasma (électrons, ions, rayonnement, neutres).
• Force de fusée anisotrope : En raison du flux de chaleur aligné avec le champ magnétique, l'évaporation est préférentielle sur l'hémisphère face au plasma, générant une force de recul (« rocket force ») modélisée par un paramètre d'asymétrie $\eta$.

B. Schéma Numérique
Le modèle utilise un intégrateur Strang-split pour résoudre les équations couplées :

1. Évaporation et recul : Appliqués par demi-pas symétriques.
2. Forces non-Lorentziennes (gravité, traînée) : Résolues par une méthode implicite en demi-pas.
3. Forces de Lorentz et advection : Avancées via le pousseur de Boris sur le pas de temps complet.
   Cette approche garantit une précision d'ordre deux en temps.

C. Couplage avec SOLPS-ITER
L'étude met en œuvre deux cadres de couplage avec le code de transport de bord SOLPS-ITER (basé sur le concept de réacteur CAT) :

• Couplage unidirectionnel : OpenEdge calcule les sources de lithium évaporé dans un fond de plasma fixe.
• Couplage bidirectionnel itératif : Un pilote Python échange les sources de lithium et les profils de plasma mis à jour entre les deux codes, utilisant une sous-relaxation ($\alpha = 0.5$) pour assurer la stabilité. Les sources sont mappées de la grille cartésienne d'OpenEdge vers la grille alignée sur les surfaces de flux de SOLPS.

3. Contributions Clés

• Développement d'un modèle complet : Intégration de la dynamique, du chargement, de l'évaporation et de la force de recul dans un code DSMC dédié.
• Validation rigoureuse : Comparaison avec des solutions analytiques (traînée-gravité) et une intégration RK45 indépendante pour l'évaporation, montrant des erreurs relatives inférieures à $10^{-5}$.
• Cadre de couplage flexible : Mise en place d'une interface robuste entre OpenEdge et SOLPS-ITER permettant une évaluation auto-cohérente de l'impact des gouttelettes sur le plasma de bord.

4. Résultats Principaux

Des simulations d'ensembles de $10^5$ gouttelettes lancées depuis les diverteurs interne et externe ont été réalisées pour des rayons initiaux de 1,5 mm, 2,5 mm et 3,5 mm.

• Dépendance à la taille :
  • Les petites gouttelettes (1,5 mm) perdent une grande fraction de leur masse par évaporation avant d'atteindre le cœur. Seule une petite fraction (environ 5 % pour le diverteur interne) pénètre jusqu'au cœur.
  • Les grosses gouttelettes (3,5 mm) conservent presque toute leur masse et se redéposent localement sur les tuiles proches du point de lancement.
• Dépendance à la position de lancement :
  • Les gouttelettes du diverteur externe se redéposent principalement localement.
  • Les gouttelettes du diverteur interne ont tendance à atteindre la paroi du côté à faible champ magnétique (LFS).
• Impact de la force de fusée : L'activation de la force de recul ($\eta > 0$) inverse la dérive radiale des gouttelettes. Au lieu de tomber vers le cœur sous l'effet de la gravité, elles sont éjectées vers l'extérieur, réduisant leur temps de séjour dans les zones à fort flux de chaleur et limitant ainsi l'évaporation.
• Couplage bidirectionnel : Dans une démonstration avec 50 gouttelettes, le couplage itératif montre que la densité de lithium et la température du plasma convergent vers un état stable après environ 10 itérations. La densité de lithium totale atteint des valeurs de l'ordre de $10^{20} \, m^{-3}$ près du point de frappe du diverteur interne.

5. Signification et Perspectives

Cette étude comble une lacune critique dans la modélisation des diverteurs à métal liquide en quantifiant le devenir des gouttelettes éjectées. Elle démontre que la taille initiale est le paramètre déterminant pour la survie et la pénétration du lithium.

Les résultats suggèrent que seuls les produits de fragmentation les plus petits offrent un chemin efficace pour transporter le lithium vers le cœur du plasma, tandis que les gouttelettes plus grosses agissent comme un système de refroidissement et de redéposition local.

Les extensions futures prévues incluent :

• L'inclusion du puits d'énergie d'évaporation dans le couplage SOLPS.
• Un couplage plus serré en mémoire pour éliminer les entrées/sorties de fichiers.
• L'intégration d'un code MHD de surface liquide pour simuler la génération auto-cohérente des gouttelettes à partir des instabilités de surface libre.

En conclusion, ce travail fournit un outil prédictif essentiel pour la conception de diverteurs à lithium liquide, permettant d'optimiser le contrôle des impuretés et la gestion thermique dans les réacteurs de fusion de prochaine génération.