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

Dit artikel presenteert een nieuw, vereenvoudigd multi-fluid evenwichtsmodel voor proton-boor fusie in Spherical Tokamaks, waarmee wordt aangetoond dat sterke rotatie en centrifugale krachten leiden tot significante ionenscheiding en elektrostatische potentiaal die met standaard MHD-modellen niet kunnen worden gevangen.

De kern

De Dans van de Atomen: Hoe we een 'Schone' Ster in een Doos Proberen te Temmen

Stel je voor dat je een mini-zon in een laboratorium wilt bouwen. Dat is precies wat wetenschappers proberen met fusie: het samensmelten van atomen om een enorme hoeveelheid energie op te wekken. De meeste wetenschappers werken met een mix van Deuterium en Tritium, maar er is een "heilige graal": de Proton-Boor (p-11B) reactie.

Waarom is die zo speciaal? Omdat deze reactie bijna geen schadelijke neutronen (straling) produceert. Het is de "schone" versie van kernenergie. Maar er is een probleem: om deze reactie te laten werken, moet je de atomen extreem hard laten botsen. Je hebt een soort magnetische "gevangenis" nodig die ongelooflijk sterk is.

Het probleem: De Draaimolen-chaos

De onderzoekers in dit artikel kijken naar een specifiek ontwerp voor die magnetische gevangenis: de Spherical Tokamak. Denk aan een donutvormige machine die de brandstof vasthoudt.

Om de boel stabiel te houden, moet de brandstof in die donut heel hard ronddraaien. Maar hier ontstaat de chaos. In de brandstof zitten twee soorten deeltjes: lichte Protonen en zware Boor-ionen.

De Metafoor: De Draaimolen op de Kermis
Stel je een draaimolen voor die razendsnel rondjes draait. Op de draaimolen zitten twee soorten mensen: kleine kinderen (de lichte protonen) en zware volwassenen met rugzakken (de zware boor-ionen).

Als de draaimolen heel hard gaat draaien, gebeurt er iets natuurkundigs: de zware volwassenen worden door de middelpuntvliedende kracht (centrifugaalkracht) naar de buitenste rand van de draaimolen geslingerd. De kleine kinderen blijven relatief meer in het midden.

In de fusiemachine gebeurt precies hetzelfde: de zware boor-atomen worden naar de buitenkant van de "donut" gedrukt, terwijl de lichte protonen in het midden blijven.

Waarom is dit een probleem voor de energie?

Voor een goede fusie moeten de protonen en de boor-atomen elkaar juist in het midden constant raken om te botsen en energie te maken. Als de boor naar de buitenkant wordt geslingerd en de protonen in het midden blijven, missen ze elkaar! Het is alsof je probeert te dansen met iemand, maar de een staat in het midden van de zaal en de ander wordt tegen de muur gedrukt. De dans (de fusie) stopt.

Bovendien zorgt dit verschil in positie voor een soort "elektrische spanning" in de machine, vergelijkbaar met een batterij die spontaan wordt opgeladen. Dit kan de hele machine verstoren.

De oplossing: Een nieuwe "Verkeersregelaar"

De wetenschappers in dit artikel hebben een nieuw wiskundig model ontwikkeld. Tot nu toe gebruikten ze modellen die ervan uitgingen dat alle deeltjes netjes samen één groep vormden (alsof iedereen op de draaimolen precies hetzelfde gewicht had). Dat werkte niet voor de p-11B reactie.

Het nieuwe model is als een super-geavanceerde verkeersregelaar. Het houdt rekening met:

1. Het gewichtverschil: Het weet dat de "zware boor-volwassenen" naar de buitenkant gaan.
2. De elektrische spanning: Het berekent precies hoeveel "spanning" er ontstaat door dat verschil.
3. De snelheid: Het kijkt hoe hard de draaimolen moet draaien zonder dat de boel uit elkaar spat.

Waarom is dit belangrijk?

Dankzij dit nieuwe model kunnen ingenieurs nu veel nauwkeuriger tekenen hoe een toekomstige energiecentrale eruit moet zien. Ze kunnen berekenen: "Als we de boor naar de buitenkant zien glippen, hoe moeten we dan de magneten instellen om ze toch bij elkaar te houden?"

Het is de eerste stap naar een wereld met bijna onbeperkte, schone energie, zonder de gevaarlijke afvalstoffen van de huidige kernenergie. Ze hebben simpelweg de "handleiding" geschreven voor hoe je de chaos van de draaimolen kunt beheersen.

Technische samenvatting

Technische Samenvatting: Ontwikkeling van een Gereduceerd Multi-Fluid Evenwichtsmodel voor Proton-Boor Sferische Tokamaks

1. Probleemstelling

De fusie van protonen en boor ($p\text{-}^{11}\text{B}$) wordt beschouwd als een veelbelovende route naar schone, aneutronische energie (minimale neutronenstraling). Echter, dit proces vereist extreem hoge ionentemperaturen ($>100\text{ keV}$) en robuuste magnetische opsluiting. Sferische Tokamaks (ST), aangedreven door krachtige Neutral Beam Injection (NBI), zijn hiervoor een ideale kandidaat.

In dergelijke apparaten treden twee complexe fenomenen op die standaardmodellen problematisch maken:

• Centrifugale scheiding: Door de enorme massaverhouding tussen protonen en boor-ionen zorgt de sterke toroidale rotatie ervoor dat zwaardere boor-ionen naar de buitenkant (low-field side) van de magnetische oppervlakken worden gedreven.
• Elektrostatische polarisatie: Deze scheiding van ladingen vereist een zelfconsistent elektrostatisch potentiaal om de quasi-neutraliteit van het plasma te behouden.

Het tekortschieten van bestaande modellen:

• Single-fluid MHD-modellen: Kunnen de differentiële verplaatsing van verschillende soorten deeltjes (species) niet beschrijven.
• Volledige multi-fluid modellen: Zijn wiskundig te complex en numeriek instabiel (stijf) vanwege de aanwezigheid van poloidale stromingen, wat ze ongeschikt maakt voor routinematige engineering-analyses.

2. Methodologie

De auteurs hebben een "gereduceerd multi-fluid evenwichtsmodel" ontwikkeld dat een balans zoekt tussen fysische nauwkeurigheid en computationele robuustheid.

Kernconcepten van het model:

• Strategische vereenvoudiging: Het model behoudt de dominante toroidale rotatie en de elektrostatische koppeling, maar verwaarloost secundaire effecten zoals de traagheid van poloidale stromingen en druk-anisotropie. Hierdoor blijft de onderliggende Grad-Shafranov (GS) vergelijking elliptisch en numeriek stabiel.
• Wiskundige formulering: Het model is gebaseerd op een gegeneraliseerde GS-vergelijking, gekoppeld aan soortspecifieke Bernoulli-relaties en een quasi-neutraliteitsvoorwaarde.
• Numerieke implementatie: Er is een nieuwe code ontwikkeld, genaamd NFEQ (N-Fluid EQuilibrium). Deze maakt gebruik van een genest iteratief schema: een Picard-iteratie voor de magnetische flux ($\psi$) en een Newton-Raphson methode voor het oplossen van het elektrostatische potentiaal ($\Phi$).

3. Belangrijkste Bijdragen

• Ontwikkeling van een hanteerbaar model: Een theoretisch kader dat specifiek is afgestemd op de unieke fysica van $p\text{-}^{11}\text{B}$ sferische tokamaks.
• NFEQ-solver: Een efficiënt numeriek algoritme dat in staat is om de complexe interactie tussen rotatie, species-scheiding en magnetische geometrie op te lossen zonder numerieke singulariteiten.
• Validatie op reactor-schaal: Het model is toegepast op twee specifieke ontwerpen van de ENN Group: de experimentele EHL-2 en de reactor-schaal EHL-3B.

4. Resultaten

De simulaties tonen aan dat de impact van multi-fluid effecten wordt bepaald door het ion-Machgetal ($M$):

• Lage rotatie ($M < 0.5$): De effecten zijn verwaarloosbaar en het model convergeert naar de standaard single-fluid MHD-oplossing.
• Hoge rotatie ($M > 1$):
  • Centrifugale scheiding: Er treedt een significante accumulatie van boor-ionen op aan de buitenkant (low-field side), wat resulteert in een "halve maan"-vormige dichtheidsverdeling. Dit veroorzaakt een mismatch in de brandstofdichtheid tussen protonen en boor in de kern.
  • Elektrostatisch potentiaal: In de reactor-schaal EHL-3B wordt een aanzienlijk elektrostatisch potentiaal gegenereerd van ongeveer 10–15 kV. Dit beïnvloedt de interne elektrische velden en daarmee de turbulentie-onderdrukking ($E \times B$ shear).
  • Geometrische wijziging: De drukgradiënten en de Shafranov-verschuiving worden fundamenteel beïnvloed door de herverdeling van de ionen.

5. Betekenis en Toepassing

Dit onderzoek bewijst dat multi-fluid modellering essentieel is voor het ontwerp van $p\text{-}^{11}\text{B}$ reactoren; single-fluid modellen geven een onnauwkeurig beeld van de fusie-output en de stabiliteit.

Brede relevantie:
Hoewel gericht op $p\text{-}^{11}\text{B}$, is het model generiek toepasbaar op:

• Zware onzuiverheden (bijv. Wolfraam): Voorspellen van de accumulatie van zware ionen in standaard D-T reactoren (zoals ITER).
• Isotopen-scheiding: Analyse van brandstofverdeling in D-T of D-$^3\text{He}$ plasma's.
• Advanced Tokamak regimes: Onderzoek naar interne transportbarrières (ITB's) in sterk roterende plasma's.