Systematic comparison of VMEC and HINT equilibrium calculations for finite-beta LHD plasmas

Cet article compare systématiquement les calculs d'équilibre VMEC et HINT pour les plasmas du Large Helical Device, révélant que, bien que les deux codes concordent à faible bêta, ils divergent à des valeurs de bêta plus élevées car HINT capture la stochasticité de bord et la rupture des surfaces de flux que l'hypothèse de surfaces de flux emboîtées de VMEC ne peut représenter.

L'essentiel

Imaginez essayer de cuire un gâteau parfaitement rond à l'intérieur d'un four très étrange et tordu. Dans le monde de la fusion nucléaire, les scientifiques utilisent des machines appelées stellarators (comme le Large Helical Device, ou LHD) pour piéger un plasma surchauffé. Pour maintenir la stabilité de ce plasma, ils doivent calculer exactement l'apparence des « murs » magnétiques qui le retiennent.

Cet article compare deux « boulangers » différents (programmes informatiques) essayant de déterminer la forme de ces murs magnétiques lorsque le plasma devient très chaud et sous pression.

Les deux boulangers : VMEC et HINT

1. VMEC (L'Architecte Rigide) : Ce programme est comme un architecte qui insiste pour que chaque couche du gâteau soit un oignon parfait, lisse et emboîté. Il suppose que les murs magnétiques ne se brisent jamais et ne se touchent jamais. Il est excellent pour les situations simples à basse pression, mais il a un angle mort : il refuse de croire que les parois puissent devenir désordonnées ou brisées.
2. HINT (L'Observateur Réaliste) : Ce programme est comme un scientifique qui regarde le gâteau en train de cuire réellement. Il ne suppose pas que les couches sont parfaites. Au contraire, il laisse la physique se produire naturellement. Si la chaleur devient trop élevée, il permet aux murs magnétiques de devenir vacillants, de se briser ou de se transformer en un désordre chaotique.

L'expérience : Monter la température

Les chercheurs ont testé ces deux programmes sur la machine LHD avec trois formes différentes de l'« four » magnétique (certaines décalées vers l'intérieur, d'autres vers l'extérieur). Ils ont augmenté lentement la pression du plasma (la « chaleur » du gâteau) de 0 % à 5 %.

Qu'est-ce qui s'est passé à basse pression ?
Lorsque le plasma était frais et calme, les deux boulangers étaient d'accord. Les murs magnétiques restaient lisses et emboîtés, tout comme le prédisait l'Architecte Rigide (VMEC). Tout allait bien.

Qu'est-ce qui s'est passé quand la chaleur a monté ?
Une fois que la pression a franchi un certain « point critique », les deux boulangers ont commencé à être en désaccord.

• VMEC continuait de dessiner des couches d'oignon parfaites, lisses et en expansion. Il pensait que le plasma devenait simplement plus grand et plus rond.
• HINT voyait quelque chose de différent. Il a remarqué que les murs magnétiques commençaient à devenir « stochastiques ».

Le désordre « stochastique » : Une analogie créative

Imaginez les lignes de champ magnétique comme un faisceau de spaghetti.

• Dans un état parfait (basse pression), les brins de spaghetti sont proprement regroupés et courent parallèlement les uns aux autres.
• À mesure que la pression augmente, le courant de Pfirsch-Schlüter (un type de courant électrique qui se forme naturellement dans le plasma) agit comme une main chaotique mélangeant les spaghettis.
• Finalement, les brins commencent à se chevaucher et à s'emmêler. C'est ce qu'on appelle les îlots magnétiques et la stochasticité. Les « couches d'oignon » ordonnées se brisent.

Parce que HINT permet ce mélange, il voit la « cage magnétique » rétrécir. Le mélange chaotique à la périphérie du plasma rend le volume effectif plus petit. VMEC, quant à lui, continue de dessiner l'oignon parfait en expansion, et pense donc que le volume augmente.

Les principales conclusions

1. Le « point de bascule » : Il existe un niveau de pression spécifique où les couches d'oignon ordonnées se brisent. Une fois ce point dépassé, VMEC n'est plus précis car il ne peut pas voir les parois brisées.
2. La forme compte : Le « point de bascule » survient plus tôt (à une pression plus faible) si la machine est décalée vers l'extérieur.
   • Analogie : Imaginez que la machine décalée vers l'extérieur est comme une table bancale. Elle est plus facile à renverser (créer du chaos) qu'une table décalée vers l'intérieur qui est robuste. La forme vers l'extérieur crée plus de « rides » dans le champ magnétique, ce qui fait que les spaghettis s'emmêlent plus vite.
3. Perte de volume : Dans les configurations décalées vers l'extérieur et standards, à mesure que la pression devient très élevée, le volume réel du plasma (selon le modèle réaliste de HINT) commence à rétrécir car les murs magnétiques se brisent. VMEC manque totalement cela, pensant que le volume continue de croître.

L'essentiel

Cet article montre que pour les plasmas de fusion à haute pression, nous ne pouvons pas nous contenter de compter sur le modèle de « l'oignon parfait » (VMEC). Nous avons besoin de « l'observateur réaliste » (HINT) pour voir quand les murs magnétiques se brisent et deviennent chaotiques. Cela est particulièrement vrai pour les machines décalées vers l'extérieur, où le champ magnétique est plus sensible à ces effets 3D désordonnés. L'étude confirme qu'en poussant vers des énergies plus hautes, l'hypothèse de couches magnétiques parfaites et lisses devient de moins en moins valide.

Résumé technique

Résumé technique : Comparaison systématique des calculs d'équilibre VMEC et HINT pour les plasmas LHD à $\beta$ fini

Énoncé du problème
Dans la physique des stellarators, le calcul de l'équilibre magnétohydrodynamique (MHD) tridimensionnel est fréquemment réalisé à l'aide du solveur VMEC. VMEC repose sur le principe variationnel et suppose strictement l'existence de surfaces de flux parfaitement imbriquées. Cependant, dans les plasmas à $\beta$ fini, l'augmentation de la pression du plasma peut altérer la topologie du champ magnétique, générant des îles magnétiques et des régions stochastiques. Ces phénomènes, particulièrement dans les dispositifs à courant net nul où ils sont pilotés par les courants de Pfirsch-Schlüter (PS), peuvent compromettre la validité de l'hypothèse des surfaces de flux imbriquées. Lorsque les îles magnétiques se chevauchent, la dernière surface de flux fermée (LCFS) se brise, réduisant le volume de plasma enfermé — une limitation observée comme étant plus sévère que le décalage de Shafranov. Alors que le code HINT résout les équations dynamiques des champs magnétiques et de la pression sans supposer de surfaces imbriquées, une comparaison systématique entre VMEC et HINT à travers diverses configurations et valeurs de bêta dans le dispositif LHD a été nécessaire pour clarifier l'influence de ces hypothèses.

Méthodologie
L'étude a mené une comparaison systématique des calculs d'équilibre entre VMEC et HINT pour les plasmas LHD. L'investigation s'est concentrée sur trois configurations spécifiques d'axe magnétique sous vide : $R_{ax}^V = 3,53$ m (décalé vers l'intérieur), $3,60$ m (standard), et $3,85$ m (décalé vers l'extérieur). L'analyse a couvert des valeurs de bêta sur l'axe ($\beta_0$) allant de $0,0\,\%$ à $5,0\,\%$ par incréments de $0,5\,\%$.

Les deux codes ont utilisé un profil de pression initial défini par $p = p_0(1-s)$, où $s$ est le flux toroïdal normalisé. L'étude a comparé trois métriques primaires :

1. La position de l'axe magnétique ($\langle R_{ax} \rangle$).
2. Le transform de rotation sur l'axe ($\iota/2\pi$).
3. Le volume délimité par la LCFS claire.

Le code HINT, qui emploie une méthode de relaxation pour résoudre les îles magnétiques et les régions stochastiques, a été contrasté avec VMEC, qui impose la contrainte de surfaces de flux imbriquées.

Résultats clés

• Régime de faible $\beta$ : Pour de faibles valeurs de $\beta_0$, VMEC et HINT ont produit des solutions d'équilibre cohérentes. Cela indique que l'hypothèse des surfaces de flux imbriquées reste largement valide dans ce régime, et que les réponses MHD 3D sont comparables entre les deux codes.
• $\beta_0$ critique et divergence : Un $\beta_0$ critique dépendant de la configuration a été identifié. Au-delà de ce seuil, les solutions des deux codes ont commencé à diverger. Dans les calculs HINT, l'axe magnétique présentait un décalage de Shafranov plus important et un transform de rotation plus élevé par rapport à VMEC.
• Rupture des surfaces de flux et volume : Dans HINT, le volume délimité par la LCFS claire a diminué à des valeurs de bêta plus élevées pour les trois configurations. Cette réduction est attribuée à l'évolution des champs magnétiques stochastiques près du bord du plasma, ce qui brise les surfaces de flux. Inversement, VMEC, contraint par son hypothèse de surfaces imbriquées, a montré une augmentation monotone du volume délimité à mesure que $\beta_0$ augmentait, échouant à représenter la rupture des surfaces de flux.
• Dépendance à la configuration : La valeur de $\beta_0$ critique a diminué en passant de la configuration décalée vers l'intérieur à la configuration décalée vers l'extérieur. La configuration décalée vers l'extérieur ($R_{ax}^V = 3,85$ m) a été jugée la plus sensible aux réponses d'équilibre 3D. Cette sensibilité provient du fait que l'ondulation hélicoïdale augmente avec $R_{ax}^V$, pilotant des courants PS plus forts, tandis que le cisaillement magnétique diminue. La combinaison de courants PS plus forts et d'un cisaillement plus faible conduit à des îles magnétiques plus grandes et à une stochasticité plus précoce à des valeurs de bêta plus basses.

Signification et revendications
L'article établit que les réponses d'équilibre 3D dans les plasmas LHD deviennent de plus en plus significatives à mesure que la configuration passe d'un décalage vers l'intérieur à un décalage vers l'extérieur. L'étude démontre que l'hypothèse des surfaces de flux imbriquées, centrale à VMEC, est compromise à $\beta$ fini en raison des perturbations pilotées par le courant PS et de la stochasticité de borden résultante.

Les auteurs affirment que bien que VMEC et HINT soient comparables pour les configurations de faible $\beta_0$, HINT fournit une description physiquement plus complète à $\beta$ plus élevé en résolvant la rupture des surfaces de flux et la réduction associée du volume de plasma. La divergence entre les codes souligne que la limitation de l'hypothèse des surfaces de flux imbriquées n'est pas uniforme mais dépend fortement de la configuration magnétique sous vide, spécifiquement de l'interaction entre la toroidicité (pilotant les courants PS) et le cisaillement magnétique. Les résultats soulignent l'importance d'utiliser des codes capables de résoudre la stochasticité, tels que HINT, lors de l'analyse des équilibres à haut $\beta$ dans les stellarators, particulièrement dans les configurations décalées vers l'extérieur où le seuil de $\beta$ critique est plus bas.