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

Dit artikel vergelijkt systematisch de VMEC- en HINT-evenwichtsberekeningen voor Large Helical Device-plasma's, waarbij wordt onthuld dat hoewel beide codes overeenstemmen bij een lage bèta, ze uiteenlopen bij hogere bètawaarden, aangezien HINT randstochastiek en het breken van fluxoppervlakken vastlegt die de aanname van geneste fluxoppervlakken van VMEC niet kan weergeven.

De kern

Stel je voor dat je probeert een perfect ronde taart te bakken in een zeer vreemde, gedraaide oven. In de wereld van fusie-energie gebruiken wetenschappers machines die stellaratoren worden genoemd (zoals de Large Helical Device, of LHD) om superheet plasma vast te houden. Om dit plasma stabiel te houden, moeten ze precies berekenen hoe de magnetische "muren" die het plasma vasthouden eruit moeten zien.

Dit artikel vergelijkt twee verschillende "bakkers" (computerprogramma's) die proberen de vorm van deze magnetische muren te bepalen wanneer het plasma erg heet en onder hoge druk staat.

De Twee Bakkers: VMEC en HINT

1. VMEC (De Strikte Architect): Dit programma is als een architect die erop staat dat elke laag van de taart een perfecte, gladde, geneste ui moet zijn. Het gaat ervan uit dat de magnetische muren nooit breken of elkaar raken. Het is geweldig voor eenvoudige, situaties met lage druk, maar het heeft een blinde vlek: het weigert te geloven dat de muren ooit rommelig of gebroken kunnen raken.
2. HINT (De Realistische Waarnemer): Dit programma is als een wetenschapper die de taart daadwerkelijk ziet bakken. Het gaat er niet vanuit dat de lagen perfect zijn. In plaats daarvan laat het de natuurkunde natuurlijk gebeuren. Als de hitte te hoog wordt, staat het toe dat de magnetische muren wankel worden, uit elkaar vallen of een chaotische bende worden.

Het Experiment: De Warmte Opvoeren

De onderzoekers testten deze twee programma's op de LHD-machine met drie verschillende vormen van de magnetische "oven" (sommigen naar binnen verschoven, anderen naar buiten). Ze verhoogden langzaam de druk van het plasma (de "hitte" van de taart) van 0% naar 5%.

Wat gebeurde er bij lage druk?
Toen het plasma koel en kalm was, waren beide bakkers het eens. De magnetische muren bleven glad en genest, precies zoals de Strikte Architect (VMEC) voorspelde. Alles was in orde.

Wat gebeurde er toen de hitte werd opgevoerd?
Zodra de druk een bepaald "kritiek punt" passeerde, begonnen de twee bakkers van mening te verschillen.

• VMEC bleef perfecte, gladde, uitdijende uienlagen tekenen. Het dacht dat het plasma alleen maar groter en ronder werd.
• HINT zag iets anders. Het merkte op dat de magnetische muren "stochastisch" begonnen te worden.

De "Stochastische" Bende: Een Creatieve Analogie

Denk aan de magnetische veldlijnen als een bos spaghetti.

• In een perfecte staat (lage druk) zijn de spaghettestengels netjes gebundeld en lopen ze parallel aan elkaar.
• Naarmate de druk stijgt, werkt de Pfirsch-Schlüter-stroom (een type elektrische stroom dat van nature in het plasma ontstaat) als een chaotische hand die de spaghetti door elkaar mengt.
• Uiteindelijk beginnen de strengen elkaar te overlappen en te verstrengelen. Dit wordt magnetische eilanden en stochastiek genoemd. De nette "uienlagen" vallen uiteen.

Omdat HINT dit verstrengeling toestaat, ziet het de "magnetische kooi" krimpen. De chaotische menging aan de rand van het plasma maakt het effectieve volume kleiner. VMEC mist dit echter volledig en tekent nog steeds de perfecte, uitdijende ui, waardoor het denkt dat het volume groter wordt.

De Belangrijkste Bevindingen

1. Het "Kantelpunt": Er is een specifiek drukgehalte waarbij de nette uienlagen breken. Zodra je dit punt passeert, is VMEC niet langer nauwkeurig omdat het de gebroken muren niet kan zien.
2. Vorm Is Belangrijk: Het "kantelpunt" treedt eerder op (bij lagere druk) als de machine naar buiten is verschoven.
   • Analogie: Stel je voor dat de naar buiten verschoven machine als een wankele tafel is. Het is makkelijker om deze omver te stoten (chaos te creëren) dan een stevige, naar binnen verschoven tafel. De naar buiten verschoven vorm creëert meer "rimpelingen" in het magnetische veld, waardoor de spaghetti sneller in de knoop raakt.
3. Volumeverlies: In de naar buiten verschoven en standaard configuraties, naarmate de druk erg hoog wordt, begint het werkelijke volume van het plasma (volgens het realistische HINT-model) te krimpen omdat de magnetische muren afbreken. VMEC mist dit volledig en denkt dat het volume blijft groeien.

De Kern van de Zaak

Dit artikel laat zien dat we voor hoogdruk fusieplasma's niet alleen kunnen vertrouwen op het "perfecte uienmodel" (VMEC). We hebben de "realistische waarnemer" (HINT) nodig om te zien wanneer de magnetische muren uit elkaar breken en chaotisch worden. Dit is vooral waar, bij machines die naar buiten zijn verschoven, waar het magnetische veld gevoeliger is voor deze rommelige, 3D-effecten. De studie bevestigt dat naarmate we streven naar hogere energie, de aanname van perfecte, gladde magnetische lagen steeds minder geldig wordt.

Technische samenvatting

Technische Samenvatting: Systematische Vergelijking van VMEC en HINT Evenwichtsberekeningen voor Eindige-$\beta$ LHD-plasma's

Probleemstelling
In stellaratorfysica wordt de berekening van driedimensionaal magnetohydrodynamisch (MHD) evenwicht veelal uitgevoerd met de VMEC-solver. VMEC vertrouwt op het variatiemodel en gaat strikt uit van het bestaan van perfect geneste fluxoppervlakken. Echter, in eindige-$\beta$ plasma's kan verhoogde plasmadruk het magnetische veldtopologie veranderen, wat leidt tot de vorming van magnetische eilanden en stochastische regio's. Deze fenomenen, die in apparaten zonder netto stroom met name worden gedreven door Pfirsch-Schlüter (PS) stromen, kunnen de geldigheid van de aanname van geneste fluxoppervlakken in gevaar brengen. Wanneer magnetische eilanden overlappen, breekt de laatste gesloten fluxoppervlak (LCFS), waardoor het omsloten plasma-volume krimpt—een beperking die eenzijdiger is dan de Shafranov-shift. Hoewel de HINT-code dynamische vergelijkingen voor magnetische velden en druk oplost zonder de aanname van geneste oppervlakken, is een systematische vergelijking tussen VMEC en HINT over diverse configuraties en bèta-waarden in de Large Helical Device (LHD) noodzakelijk geweest om de invloed van deze aannames te verhelderen.

Methodologie
De studie voerde een systematische vergelijking uit van evenwichtsberekeningen tussen VMEC en HINT voor LHD-plasma's. De investigatie richtte zich op drie specifieke vacuüm magnetische-as configuraties: $R_{ax}^V = 3.53$ m (naar binnen verschoven), $3.60$ m (standaard), en $3.85$ m (naar buiten verschoven). De analyse besloeg on-as bèta-waarden ($\beta_0$) variërend van $0.0\%$ tot $5.0\%$ in stappen van $0.5\%$.

Beide codes maakten gebruik van een initiële drukprofiel gedefinieerd als $p = p_0(1-s)$, waarbij $s$ de genormaliseerde toroidale flux is. De studie vergeleek drie primaire metrieken:

1. De positie van de magnetische as ($\langle R_{ax} \rangle$).
2. De rotationele transformatie op de as ($\iota/2\pi$).
3. Het volume omsloten door de heldere LCFS.

De HINT-code, die een relaxatiemethode gebruikt om magnetische eilanden en stochastische regio's op te lossen, werd gecontrasteerd met VMEC, dat de restrictie van geneste fluxoppervlakken afdwingt.

Belangrijkste Resultaten

• Laag-$\beta$ Regime: Voor lage $\beta_0$-waarden produceerden VMEC en HINT consistente evenwichtoplossingen. Dit geeft aan dat de aanname van geneste fluxoppervlakken in dit regime grotendeels geldig blijft en dat de 3D MHD-responsen tussen beide codes vergelijkbaar zijn.
• Kritieke $\beta_0$ en Divergentie: Een configuratie-afhankelijke kritieke $\beta_0$ werd geïdentificeerd. Voorbij deze drempelwaarde begonnen de oplossingen van de twee codes te divergeren. In de HINT-berekeningen vertoonde de magnetische as een grotere Shafranov-shift en een hogere rotationele transformatie vergeleken met VMEC.
• Breken van Fluxoppervlakken en Volume: In HINT nam het volume omsloten door de heldere LCFS af bij hogere bèta-waarden voor alle drie de configuraties. Deze afname wordt toegeschreven aan de evolutie van stochastische magnetische velden nabij de plasma-rand, wat de fluxoppervlakken breekt. Daarentegen vertoonde VMEC, beperkt door zijn geneste oppervlak-aanname, een monotone toename van het omsloten volume naarmate $\beta_0$ toenam, waardoor het falen om het breken van fluxoppervlakken weer te geven.
• Configuratieafhankelijkheid: De kritieke $\beta_0$-waarde nam af van naar binnen verschoven naar naar buiten verschoven configuraties. De naar buiten verschoven configuratie ($R_{ax}^V = 3.85$ m) bleek het meest gevoelig voor 3D evenwichtsresponsen. Deze gevoeligheid ontstaat doordat de helicale rimpeling (ripple) toeneemt met $R_{ax}^V$, wat sterkere PS-stromen drijft, terwijl de magnetische shear afneemt. De combinatie van sterkere PS-stromen en lagere shear leidt tot grotere magnetische eilanden en eerdere stochastisiteit bij lagere bèta-waarden.

Betekenis en Claims
Het artikel stelt vast dat 3D evenwichtsresponsen in LHD-plasma's steeds significanter worden naarmate de configuratie verschuift van naar binnen naar naar buiten. De studie demonstreert dat de aanname van geneste fluxoppervlakken, centraal in VMEC, bij eindige bèta wordt aangetast door door PS-stroom gedreven perturbaties en de resulterende rand-stochastisiteit.

De auteurs beweren dat hoewel VMEC en HINT vergelijkbaar zijn voor laag-$\beta_0$ configuraties, HINT een fysisch completere beschrijving biedt bij hogere bèta door het breken van fluxoppervlakken en de daarmee gepaard gaande reductie van het plasma-volume op te lossen. De divergentie tussen de codes benadrukt dat de beperking van de geneste fluxoppervlak-aanname niet uniform is, maar sterk afhangt van de vacuüm magnetische configuratie, specifelijk de interactie tussen toroidiciteit (die PS-stromen drijft) en magnetische shear. De resultaten onderstrepen het belang van het gebruik van codes die in staat zijn om stochastisiteit op te lossen, zoals HINT, bij het analyseren van high-beta evenwichten in stellaratoren, met name in naar buiten verschoven configuraties waar de kritieke bèta-drempel lager is.