Development of a Reduced Multi-Fluid Equilibrium Model and Its Application to Proton-Boron Spherical Tokamaks

Cette étude développe un modèle multi-fluide réduit pour simuler les effets de la rotation toroïdale et de la séparation centrifuge sur l'équilibre des tokamaks sphériques destinés à la fusion proton-bore, démontrant que ces phénomènes modifient significativement la distribution des ions et le potentiel électrostatique.

L'essentiel

Le Grand Mélange de l'Énergie Propre : Pourquoi les Particules ne font pas ce qu'on leur demande

Imaginez que vous essayez de créer une "étoile miniature" sur Terre pour produire de l'énergie propre. Pour cela, on utilise une réaction appelée Proton-Bore (p-11B). C'est le "Saint Graal" de l'énergie : contrairement aux réacteurs actuels, elle ne produit presque pas de déchets radioactifs.

Mais il y a un problème : pour que cette réaction fonctionne, il faut que les particules (les protons et le bore) soient chauffées à des températures délirantes et maintenues dans un champ magnétique ultra-puissant, un peu comme dans un tourbillon géant appelé Tokamak.

Le problème : La "Danse des Particules" est trop compliquée

Jusqu'à présent, les scientifiques utilisaient des modèles mathématiques simplifiés (appelés modèles "MHD") pour prédire comment le plasma se comporte. C'est un peu comme si, pour prévoir la météo, on disait simplement : "Le vent souffle vers le nord". C'est utile, mais c'est trop simple.

Dans un réacteur p-11B, les particules tournent si vite qu'elles commencent à se comporter de manière très différente selon leur poids :

1. Les Protons sont légers, comme des balles de ping-pong.
2. Le Bore est beaucoup plus lourd, comme des boules de pétanque.

Quand le réacteur tourne à toute vitesse, la force centrifuge (celle qui vous plaque contre le siège quand une voiture prend un virage serré) agit différemment sur eux. Les "boules de pétanque" (le Bore) sont projetées vers l'extérieur du tourbillon, tandis que les "balles de ping-pong" (les Protons) restent plus au centre.

Le souci ? Pour que la fusion fonctionne, il faut que les protons et le bore se rentrent dedans ! Si le Bore s'échappe vers les bords, la réaction s'arrête. C'est comme essayer de faire une recette de cuisine où la farine s'envole d'un côté et les œufs de l'autre : vous n'aurez jamais de gâteau.

La solution de l'équipe : Le modèle "Équilibre Réduit"

Les chercheurs (notamment de l'entreprise ENN) ont créé un nouvel outil mathématique. Ils ont trouvé le juste milieu entre deux extrêmes :

• Le modèle trop simple : Qui ignore que les particules ont des poids différents (et qui se trompe sur la recette).
• Le modèle trop complexe : Qui est tellement difficile à calculer qu'il prendrait des années à un ordinateur pour donner un résultat.

Ils ont inventé un modèle "intelligent et rapide". Il est capable de prédire précisément comment le Bore va se déplacer et, surtout, comment cela crée une "barrière électrique" (un potentiel électrostatique) qui tente de rééquilibrer les choses.

Ce qu'ils ont découvert (L'effet "Crescent")

En testant leur modèle sur des machines virtuelles (les modèles EHL-2 et EHL-3B), ils ont vu quelque chose de fascinant :
Quand le réacteur tourne très vite, le Bore ne se contente pas de bouger, il forme une sorte de "croissant de lune" sur les bords du réacteur. Il quitte le cœur du réacteur, laissant un vide là où la fusion devrait avoir lieu.

Ils ont aussi découvert que cela crée une tension électrique énorme (jusqu'à 10 000 volts !), ce qui change complètement la structure interne du réacteur.

Pourquoi c'est important ?

C'est comme si, avant de construire une Formule 1, on avait enfin trouvé le logiciel capable de simuler précisément comment les pneus et le carburant se comportent dans un virage à 300 km/h.

Grâce à ce travail, les ingénieurs ne construisent plus "à l'aveugle". Ils savent maintenant qu'ils doivent concevoir des réacteurs capables de gérer ce "mélange inégal" des particules. C'est une étape cruciale pour passer du rêve de l'énergie propre infinie à la réalité d'une centrale électrique qui alimente nos maisons.

Résumé technique

Résumé Technique : Développement d'un modèle d'équilibre multi-fluides réduit pour les tokamaks sphériques Proton-Bore

1. Problématique

La fusion Proton-Bore ($p\text{-}^{11}\text{B}$) est une voie prometteuse pour l'énergie propre car elle est aneutronique (minimisant les déchets radioactifs). Cependant, elle nécessite des températures ioniques extrêmement élevées ($> 100\text{ keV}$) et un confinement magnétique robuste, des conditions que les Tokamaks Sphériques (ST) alimentés par injection de faisceaux neutres (NBI) peuvent offrir.

Le problème majeur réside dans la physique du plasma :

• Disparité de masse : La différence de masse entre les protons ($1\text{ m}_p$) et les ions bore ($\sim 11\text{ m}_p$) est extrême.
• Effets centrifuges : Sous une forte rotation toroïdale, la force centrifuge déplace les ions lourds (bore) vers le côté champ faible (LFS), créant une séparation des espèces.
• Polarisation électrostatique : Cette séparation induit un potentiel électrostatique ($\Phi$) nécessaire pour maintenir la quasi-neutralité.
• Limites des modèles actuels : Les modèles MHD à fluide unique échouent à capturer ces effets. À l'inverse, les modèles multi-fluides complets (incluant les flux poloïdaux) sont numériquement instables et trop complexes pour une utilisation en ingénierie de routine.

2. Méthodologie

Les auteurs ont développé un modèle multi-fluides réduit conçu pour équilibrer fidélité physique et robustesse numérique.

• Simplifications stratégiques : Le modèle privilégie l'inertie toroïdale et le couplage électrostatique, tout en négligeant l'inertie des flux poloïdaux (pour éviter les singularités mathématiques de type transsonique) et l'anisotropie de pression.
• Cadre mathématique : Le modèle est formulé comme une équation de Grad-Shafranov généralisée, couplée à des relations de Bernoulli spécifiques à chaque espèce et à une contrainte de quasi-neutralité.
• Algorithme numérique (NFEQ) : Un schéma itératif imbriqué a été développé :
  • Une boucle externe de Picard pour le flux poloidal ($\psi$).
  • Une boucle interne de Newton-Raphson pour résoudre le potentiel électrostatique ($\Phi$) à chaque point de la grille.
  • Un mécanisme de contrôle par rétroaction (feedback) pour maintenir un courant total ($I_p$) fixe lors des simulations.

3. Contributions Clés

• Nouveau modèle théorique : Établissement d'un cadre mathématique capable de décrire la séparation des espèces et la formation de potentiels électrostatiques sans la complexité des modèles à deux fluides complets.
• Outil de conception robuste : Fourniture d'un outil de simulation efficace pour l'optimisation de l'ingénierie des réacteurs de fusion.
• Validation sur des dispositifs réels : Application du modèle aux configurations de l'entreprise ENN : le dispositif expérimental EHL-2 et le réacteur à grande échelle EHL-3B.

4. Résultats Principaux

L'étude démontre que l'impact des effets multi-fluides est régi par le nombre de Mach ionique ($M$) :

• Régime de faible rotation ($M < 0,5$) : Les effets multi-fluides sont faibles ; les solutions convergent vers les limites de la MHD à fluide unique.
• Régime de forte rotation ($M > 1$) :
  • Séparation centrifuge : Les ions bore s'accumulent massivement sur le côté champ faible (LFS), créant une distribution en forme de "croissant" et appauvrissant le cœur du plasma. Cela peut réduire la densité de puissance de fusion locale.
  • Potentiel électrostatique : Dans le réacteur EHL-3B, un potentiel interne massif de l'ordre de 10 à 15 kV est généré pour compenser la séparation des charges.
  • Modification de l'équilibre : La redistribution de la pression modifie le facteur de sécurité ($q$) et le profil de courant, ce qui est crucial pour la stabilité du plasma.

5. Signification et Portée

Ce travail démontre que la modélisation multi-fluides n'est pas une simple option théorique, mais une nécessité pratique pour la conception des réacteurs $p\text{-}^{11}\text{B}$.

Au-delà de la fusion Bore, le modèle est applicable à :

• Le transport des impuretés de haut Z (ex: Tungstène) : Prédiction de l'accumulation des impuretés dans les tokamaks D-T (comme ITER).
• La séparation isotopique : Analyse de la distribution du combustible dans les plasmas rotatifs.
• Les régimes de Tokamak Avancé : Étude des barrières de transport interne (ITB) liées à la rotation.