Improved n=1 Empirical Error Field Penetration Threshold Scaling with Ohmic and L-Mode Conventional Tokamak Plasma Discharges

Dit artikel presenteert een verbeterde n=1 empirische schaalwet voor de penetratiedrempel van foutvelden, afgeleid van een uitgebreide database van conventionele Ohmse en L-modus tokamakontladingen, die een hogere fitkwaliteit en een verminderde onzekerheid biedt om de engineeringstoleranties en het ontwerp van toekomstige conventionele tokamaks beter te onderbouwen.

De kern

Stel je een tokamak (een machine ontworpen om fusie-energie te creëren) voor als een gigantische, onzichtbare draaikolk van superheet gas, op zijn plaats gehouden door krachtige magnetische velden. Idealiter is deze draaikolk perfect symmetrisch, zoals een gladde, draaiende tol. In de echte wereld zijn de magneten die hem bij elkaar houden echter niet perfect. Ze hebben kleine kantelingen, verschuivingen en onvolkomenheden. Deze onvolkomenheden creëren "foutvelden" - kleine, ongewenste magnetische rimpelingen die de gladde draaiing van het gas kunnen verstoren.

Als deze rimpelingen te sterk worden, kunnen ze ervoor zorgen dat de draaikolk een "knoop" ontwikkelt (een magnetisch eiland) dat vast komt te zitten. Eenmaal vastzittend kan het hele systeem instorten, wat leidt tot een plotselinge stillegging die een "disruptie" wordt genoemd. Dit is slecht nieuws voor de machine en de mensen die hem bouwen.

Het Probleem: Hoe sterk is te sterk?
Ingenieurs moeten de exacte limiet kennen: Hoe groot mogen deze magnetische onvolkomenheden worden voordat de machine bezwijkt? Als ze de limiet te laag stellen, moeten ze de machine met onmogelijke precisie bouwen, waardoor deze astronomisch duur wordt en langzaam te construeren is. Als ze de limiet te hoog stellen, kan de machine crasht.

Jarenlang hebben wetenschappers geprobeerd een "vuistregel" (een schaalwet) te creëren om deze limiet te voorspellen op basis van de grootte van de machine en het gedrag van het gas. Maar de oude regels waren wat wankel, zoals een kaart met onscherpe randen.

De Oplossing: Een Scherpere Kaart
Dit artikel presenteert een nieuwe, verbeterde "kaart" (een empirische schaalwet) die veel duidelijker en betrouwbaarder is. Hier is hoe ze dit deden, met behulp van eenvoudige analogieën:

• Schoonmaken van de Data: De onderzoekers gingen door een enorme database van eerdere experimenten uit tokamaks over de hele wereld (zoals DIII-D, JET en KSTAR). Ze besloten zich te richten op specifieke soorten "weer" in de machine: de "Ohmische" en "L-modus" condities. Ze lieten de "H-modus" buiten beschouwing, omdat die toestand als een stevig fort is - het is zeer moeilijk te breken, dus het helpt ons niet om het zwakste punt van de machine te begrijpen. Door zich te richten op de kwetsbare toestanden, vonden ze de ware gevarenzone.
• Toevoegen van Nieuwe Ingrediënten: Ze voegden nieuwe data toe van twee specifieke machines: J-TEXT (die kleiner is en op lagere stromen draait) en meer data van JET (die enorm is, vergelijkbaar met de toekomstige ITER-machine). Denk hierbij aan het toevoegen van nieuwe proefritten aan een auto-veiligheidsdatabase. Je hebt kleine auto's en gigantische vrachtwagens in de data nodig om te weten hoe de veiligheidsregels van toepassing zijn op elk voertuig dat je in de toekomst zou kunnen bouwen.
• Betere Wiskunde: Ze gebruikten een meer geavanceerde wiskundige methode om de relatie te vinden tussen de grootte van de machine, de sterkte van het magnetische veld, de gasdichtheid en de elektrische stroom. Ze ontdekten dat plasma-stroom (hoeveelheid elektriciteit die door het gas stroomt) een kritieke factor is die ze eerder niet volledig hadden meegerekend.

De Nieuwe Bevindingen
De nieuwe "vuistregel" vertelt ons dat:

1. Hogere dichtheid is je vriend: Het in de machine persen van meer gas maakt het moeilijker voor foutvelden om een crash te veroorzaken.
2. Grotere machines zijn verrassend veerkrachtig: Grotere machines (zoals de toekomstige ITER) kunnen grotere magnetische onvolkomenheden aan dan we eerder dachten.
3. Stroom is van belang: De hoeveelheid stroom die door het plasma stroomt, verandert hoe de machine reageert op deze fouten.

Waarom Dit Belangrijk Is voor de Toekomst
Het artikel kijkt specifiek naar het ITER-project, een massaal internationaal fusie-experiment dat momenteel in aanbouw is. Met hun nieuwe, scherpere kaart lieten de onderzoekers miljoenen simulaties draaien (zoals het uitvoeren van een weersvoorspelling een miljoen keer met licht verschillende startvoorwaarden).

Het Resultaat: Ze ontdekten dat ITER veel beter in orde is dan we dachten. De "gevarenzone" voor magnetische fouten ligt veel verder weg dan de werkelijke onvolkomenheden die van ITER worden verwacht.

• De Oude Kaart: Suggereerde dat er een aardige kans was dat ITER over zijn eigen veters zou struikelen (vastgelopen modi).
• De Nieuwe Kaart: Toont aan dat de kans hierop ongelooflijk klein is (minder dan 1 op een miljoen voor het meest waarschijnlijke scenario).

De Conclusie
Dit artikel zegt niet alleen "fusie is moeilijk". Het geeft ingenieurs een veel zelfverzekerder, preciezer liniaal om de toleranties van hun machines te meten. Omdat de nieuwe regels aantonen dat de machines robuuster zijn tegen magnetische fouten, hoeven ingenieurs de magneten misschien niet met zulke extreme, dure precisie te bouwen. Dit kan tijd en geld besparen terwijl de machine veilig blijft.

Kortom: Ze namen een wazige, verwarrende kaart van magnetische veiligheidslimieten, maakten deze schoon met betere data en slimmere wiskunde, en ontdekten dat de toekomst van fusiecentrales veiliger en haalbaarder is dan we eerder geloofden.

Technische samenvatting

1. Probleemstelling

Tokamaks vereisen een nauwkeurige axiale symmetrie om stabiel te opereren. Echter, fabricageonvolkomenheden (spoelkanten, verschuivingen en misaligneringen) creëren niet-axiale foutvelden (EF's). Als deze velden een kritieke penetratiedrempel overschrijden, kunnen ze magnetische reconnectie triggeren, magnetische eilanden vormen en ervoor zorgen dat deze eilanden "vastlopen" aan de vesselwand. Dit leidt tot vastgelopen modi en daaropvolgende plasma-disrupties.

Het voorspellen van deze drempel is cruciaal voor:

• Technische toleranties: Bepalen van toelaatbare spoelmisaligneringen tijdens de constructie.
• Kosten en planning: Te strenge toleranties verhogen de bouwkosten en -tijd; te losse toleranties riskeren apparaatfalen.
• Toekomstige apparaten: Projecten zoals ITER en SPARC hebben een hoge opgeslagen energie, waardoor disrupties bijzonder gevaarlijk zijn.

Vorige empirische schaalwetten leden aan:

• Lage aanpassingskwaliteit: Bepalingscoëfficiënten ($R^2$) zo laag als 0,38.
• Hoge onzekerheid: Tot 48% onzekerheid in projecties voor ITER.
• Inconsistente data: Gemengde opname van H-mode, L-mode en Ohmische plasma's, sferische tokamaks (NSTX), en variërende niveaus van intrinsieke foutveldmodellering over verschillende apparaten.
• Beperkte parameterruimte: Gebrek aan data voor apparaten met grote hoofdstraal (zoals ITER) en specifieke dichtheidsregimes.

2. Methodologie

De auteurs ontwikkelden een bijgewerkte empirische schaalwet door de database en de fitmethode te verfijnen.

A. Database verfijning

• Regime-selectie: De studie beperkte de dataset tot conventionele (niet-sferische) tokamaks die opereren in Ohmische en L-mode regimes. H-mode werd uitgesloten omdat het resistenter is tegen foutveldpenetratie, en de auteurs wilden het "meest gevaarlijke" regime modelleren. Shot met verzadigde Ohmische opsluiting (SOC) werden uitgesloten om verstorende dichtheidsafhankelijkheden te vermijden.
• Apparaatopname:
  • Toegevoegd: Nieuwe data van J-TEXT (leverend lage stroom/dichtheidspunten en een nieuwe hoofdstraalgrootte) en een uitgebreide dataset van JET (21 shots versus de vorige 4, dekkend grote hoofdstraal en hoge plasmastroom).
  • Uitgesloten: NSTX (sferische tokamak) en H-mode shots.
• Standaardisatie:
  • DIII-D: Overgeschakeld van een gegroepeerde evenwichtsbenadering naar individuele evenwichtsreconstructies voor elke shot met de SURFMN-code om intrinsieke foutvelden nauwkeuriger te modelleren.
  • J-TEXT: Gebruik van de TokaMaker-code om evenwichten te reconstrueren (aangezien routinematige magnetische reconstructies niet beschikbaar waren), zodat compatibiliteit met de GPEC (General Perturbed Equilibrium Code) die voor alle andere apparaten wordt gebruikt, gewaarborgd was.
  • Maatstaf: Gebruik van de Dominante Mode-overlap ($\delta$), een coördinatenstelsel-onafhankelijke maatstaf die rekening houdt met de plasmarespons, berekend via GPEC.

B. Schaalmethodologie

• Regressietechniek: Toepassing van forward selection least-squares regressie om de optimale log-lineaire machtswet te identificeren.
• Parameterselectie: Uitbreiding van de parameterruimte boven eerdere studies uit om expliciete afhankelijkheden op te nemen van Plasmastroom ($I_p$), Toriaal Veld ($B_T$), Hoofdstraal ($R_0$), Dichtheid ($n_e$) en Genormaliseerde Druk ($\beta_n/\ell_i$).
• Validatie: Gebruik van Conditiestallen om multicollineariteit te detecteren (stoppen bij 5 variabelen om het conditiestal < 20 te houden) en vergelijking van Ordinary Least Squares (OLS) met Weighted Least Squares (WLS) met Kernel Density Estimation (KDE) om oversampling van data-rijke gebieden te vermijden.
• Onzekerheidskwantificering: Aannemen van een basale 10% onzekerheid in de identificatietijd van vastgelopen modi, doorgegeven via Monte Carlo-simulaties om waarschijnlijkheidsdichtheidsfuncties (PDF's) en cumulatieve verdelingsfuncties (CDF's) voor voorspellingen te genereren.

3. Belangrijkste bijdragen

1. Nieuwe schaalwetten: Afgeleid twee nieuwe schaalvergelijkingen (Vergelijking 9 voor OLS en Vergelijking 10 voor WLS) met significant verbeterde statistische maatstaven.
   • Vergelijking 9 (OLS): $\delta \propto (\beta_n/\ell_i)^{0.25} I_p^{-0.97} R_0^{1.88} n_e^{0.77} B_T^{0.20}$
   • Vergelijking 10 (WLS): $\delta \propto (\beta_n/\ell_i)^{0.13} I_p^{-1.01} R_0^{1.57} n_e^{0.56} B_T^{0.30}$
2. Opname van expliciete $I_p$-afhankelijkheid: In tegenstelling tot eerdere schalingen die vertrouwden op impliciete afhankelijkheden, omvatten de nieuwe modellen expliciet de plasmastroom, wat aantoont dat lagere plasmastroom de weerstand tegen foutvelden verhoogt.
3. Verschuiving in toriaal veld-afhankelijkheid: De schaling voor $B_T$ verschuift van een sterke negatieve afhankelijkheid (in eerdere modellen) naar een kleine positieve afhankelijkheid, gedreven door de opname van $I_p$ en de verwijdering van sferische tokamak-data.
4. Robuustheidstesten: Aangetoond dat de schaling robuust is tegen verwijdering van individuele apparaten (bijv. DIII-D, J-TEXT, C-MOD), terwijl "minimale" apparaat-subsets (bijv. JET+KSTAR+C-MOD) leidden tot ernstig overfitting.

4. Resultaten

• Verbeterde aanpassingskwaliteit: De nieuwe schaling bereikte een $R^2$ van 0,63 (OLS) en 0,66 (WLS), een significante verbetering ten opzichte van eerdere waarden (0,38–0,58).
• Verminderde onzekerheid:
  • Verminderde Chi-kwadraat: Verbeterd van 17,0 (2020-model) naar 12,99.
  • Gemiddelde absolute log-fout: Verbeterd van 1,336 naar 1,306 (equivalent aan het verminderen van de voorspellingonzekerheid van $\pm33,6\%$ naar $\pm30,6\%$).
• Parameterafhankelijkheden:
  • Dichtheid ($n_e$): Bevestigde een positieve afhankelijkheid ($n_e^{0.56-0.77}$), consistent met veel individuele apparaten maar verschillend van sommige theoretische modellen.
  • Hoofdstraal ($R_0$): Een sterke positieve afhankelijkheid ($R_0^{1.57-1.88}$), wat aangeeft dat grotere apparaten resistenter zijn.
  • Plasmastroom ($I_p$): Een sterke negatieve afhankelijkheid ($I_p^{-0.97}$ tot $-1.01$), wat aangeeft dat hogere stroom het plasma vatbaarder maakt voor vastlopen.
• Apparaatspecifieke uitbijters:
  • NSTX: De nieuwe schaling onderschat de foutveldtolerantie van NSTX aanzienlijk (voorspellend $\delta \sim 10^{-5}$ versus werkelijk $\sim 10^{-4}$), waarschijnlijk vanwege zijn sferische aard en hoge $q_{95}$-profiel, wat bevestigt dat het moet worden uitgesloten van conventionele tokamak-schalingen.
  • C-MOD: Hoge-dichtheidsdata van C-MOD was cruciaal voor het vastleggen van niet-monotoon gedrag van de dichtheidsschaling.

5. Betekenis en implicaties voor ITER

• ITER-projectie:
  • De nieuwe schaling voorspelt een veel hogere foutveldpenetratiedrempel voor ITER in vergelijking met eerdere modellen.
  • Meest waarschijnlijke intrinsieke $\delta$ voor ITER: $0,38 \times 10^{-4}$.
  • Voorspelde drempel (Vergelijking 10): $4,01 \times 10^{-4}$ (met een 1-sigma onzekerheidsbereik).
  • Vastloopkans: De kans dat ITER vastloopt onder zijn basisscenario wordt berekend als extreem laag (0,0001% voor het meest waarschijnlijke intrinsieke veld, en 22,3% zelfs voor het maximale mogelijke intrinsieke veld).
• Technische impact:
  • De resultaten suggereren dat de huidige technische toleranties van ITER voor spoelalignering conservatief en toereikend zijn.
  • De schaling biedt een betrouwbaarder hulpmiddel voor het ontwerpen van toekomstige apparaten (zoals SPARC) en het optimaliseren van operaties op bestaande machines door de waarschijnlijkheid van modus-vastlopen te voorspellen op basis van Monte Carlo-simulaties van intrinsieke foutvelden.
• Toekomstig werk: De auteurs suggereren dat het toevoegen van data van JT-60SA (grote straal) en ASDEX-U (hoge $\beta_n$) de schaling verder zou verfijnen. Bovendien zou het ontwikkelen van een schaling gebaseerd op shots met expliciete rotatiemetingen de voorspellende kracht verder kunnen verbeteren, aangezien rotatie een dominante factor is in de fysica van foutvelden.

Concluderend stelt dit artikel een robuustere, fysisch-consistente en statistisch betrouwbare empirische schaling vast voor n=1 foutveldpenetratie, wat het vertrouwen in het ontwerp en de operatie van volgende-generatie conventionele tokamaks aanzienlijk verhoogt.