SPRAY: A smoothed particle radiation hydrodynamics code for modeling high intensity laser-plasma interactions

In dit artikel wordt SPRAY voorgesteld, een nieuw, op GPU's gebaseerd en mesh-vrij SPH-code voor stralingshydrodynamica dat specifiek is ontwikkeld voor het nauwkeurig simuleren van complexe laser-plasma-interacties bij hoge energie-dichtheden.

De kern

Wat is SPRAY eigenlijk?

Stel je voor dat je een heel krachtige laserstraal op een stukje metaal schijnt. Het metaal wordt niet gewoon warm; het smelt, verdampt en explodeert in een razendsnel dans van plasma (een superheet gas van geladen deeltjes). Dit noemen we laser-plasma interactie.

Het simuleren van dit proces op een computer is als proberen een explosie in slow-motion te filmen, maar dan met wiskunde. Het probleem is dat het materiaal zich op een heel rare manier gedraagt: het rekt uit, krimpt, vormt golven en breekt open.

SPRAY is een nieuwe computercode (een programma) die deze explosies simuleert. Maar in plaats van de ruimte op te delen in een vast raster van bakjes (zoals een schaakbord), gebruikt SPRAY duizenden kleine "deeltjes" die vrij door de ruimte kunnen bewegen.

De grote uitdaging: Het "Vast Raster" vs. "Vrije Deeltjes"

Om dit te begrijpen, kun je twee manieren van kijken vergelijken:

1. De oude methode (Eulerisch): Stel je voor dat je een zwembad hebt dat je hebt ingedeeld in een rooster van bakjes. Als het water stroomt, moet je berekenen hoeveel water van het ene bakje naar het andere stroomt. Als het water echter een enorme golf vormt of een gat in het water valt, raken de bakjes in de war. De randen van het water worden onnauwkeurig, en de computer moet steeds nieuwe bakjes toevoegen of verwijderen. Dit is traag en lastig.
2. De SPRAY-methode (SPH - Smoothed Particle Hydrodynamics): Stel je nu voor dat het water bestaat uit miljoenen kleine, zwevende balletjes. Als je een golf maakt, bewegen de balletjes gewoon mee. Er zijn geen bakjes. Als het water uit elkaar valt of in een kring draait, volgen de balletjes de beweging moeiteloos. SPRAY gebruikt deze "balletjes" om de fysica te berekenen.

De kern van het artikel: SPRAY is de eerste code die deze "balletjes-methode" succesvol toepast op de extreme wereld van hoge-energie-fysica (zoals kernfusie en sterreninterieurs), waar lasers en plasma samenkomen.

De magische trucs van SPRAY

De auteurs hebben SPRAY uitgerust met een paar slimme trucs om de problemen van de oude methoden op te lossen:

• De "Spiegel"-truc (Vrije oppervlakken):
  Als een laser op een doelwit schijnt, verdampt de buitenkant en ontstaat er een wolk van plasma die zich razendsnel uitbreidt. In de oude methoden verdwijnt de "rand" van het materiaal vaak uit beeld, waardoor de berekening fout gaat.
  De oplossing: SPRAY gebruikt een slimme truc. Als een balletje aan de rand staat en er zijn geen buren meer aan de buitenkant, "verzonnt" de computer een spiegelbeeld-balletje aan de andere kant van de rand. Dit zorgt ervoor dat de berekening aan de rand net zo nauwkeurig blijft als in het midden. Het is alsof je een spiegel voor een schilderij zet om te zien wat er achter de rand gebeurt.

• De "Laser-straal" die niet vastzit:
  Normaal gesproken volgen laserstralen in computersimulaties een vast rooster. Maar als het plasma beweegt, zit de straal vast in een bakje waar hij niet meer hoort.
  De oplossing: SPRAY behandelt elke laserstraal als een onafhankelijk, zwevend deeltje (een "virtueel deeltje"). Deze deeltjes kunnen vrij door het plasma vliegen, buigen waar het plasma dikker is (breking) en worden geabsorbeerd waar het nodig is. Het is alsof je duizenden muggen hebt die door een mist vliegen en precies meten waar ze botsen, in plaats van ze te forceren om op een ruit te blijven zitten.

• De "Versneller" (GPU):
  Omdat er miljoenen balletjes zijn, moet de computer enorm veel berekeningen doen. SPRAY is gebouwd om te werken op GPU's (de krachtige videokaarten die ook in gaming-computers zitten).
  De analogie: Stel je voor dat je 10.000 mensen hebt die een muur moeten bouwen. Een oude computer is als één bouwvakker die alles één voor één doet. SPRAY gebruikt duizenden GPU-kernen die allemaal tegelijkertijd een steen leggen. Hierdoor is de simulatie honderden keren sneller.

Wat hebben ze bewezen?

De auteurs hebben SPRAY getest met verschillende "proefballonnen":

1. De Sod-schokbuis: Een klassieke test waarbij een schokgolf door een buis schiet. SPRAY kon dit perfect nabootsen, net als de wiskundige theorie voorspelde.
2. Laser op aluminium: Ze schoten een laser op een aluminiumplaat en vergeleken het met een ander, zeer betrouwbaar programma (MULTI-IFE). De resultaten waren bijna identiek.
3. Rayleigh-Taylor instabiliteit: Dit is een fenomeen waar zwaar vloeistof op licht vloeistof ligt (zoals olie op water, maar dan andersom), waardoor er mooie, kronkelende patronen ontstaan. SPRAY kon deze patronen precies zo goed simuleren als de beste bestaande programma's.
4. Kernfusie-implosie: Ze simuleerden hoe een laser een kleine capsule samendrukt (zoals in een kernfusie-reactor). Ook hier werkte het goed.

Waarom is dit belangrijk?

Tot nu toe was het moeilijk om deze soort explosies en vervormingen nauwkeurig te simuleren met de oude "bakjes"-methoden. SPRAY bewijst dat de "balletjes-methode" (SPH) een krachtig alternatief is.

Het opent de deur voor betere ontwerpen van:

• Kernfusie-energie: Om schoner en onbeperkt energie te maken.
• Sterrenkunde: Om te begrijpen wat er in het binnenste van sterren gebeurt.
• Nieuwe materialen: Om te zien hoe materialen reageren op extreme krachten.

Kortom: SPRAY is een nieuwe, supersnelle en slimme manier om de chaos van een laser-explosie op de computer te ordenen, door te vertrouwen op vrij bewegende deeltjes in plaats van starre bakjes. Het is een grote stap voorwaarts voor de wetenschap.

Technische samenvatting

Probleemstelling

De interactie tussen hoge-intensiteit lasers en plasma's is een centraal onderzoeksonderwerp in de fysica van hoge energiedichtheid (HEDP) en inertieel opsluitingfusie (ICF). Wanneer een doelwit wordt bestraald met een intense laser, treden er complexe stromingsvervormingen op door instabiliteiten (zoals de Rayleigh-Taylor-instabiliteit) en extreme druk- en dichtheidsgradiënten.

Bestaande numerieke methoden hebben hierbij beperkingen:

• Eulerische methoden (rooster-based): Moeilijkheden bij het volgen van bewegende vloeistofgrenzen en complexe vervormingen. Het vereist vaak adaptief roosterverfijning (AMR), wat computatief zwaar is.
• Lagranse methoden (cell-based): Hoewel ze grensvlakken beter volgen, zijn ze beperkt in het hanteren van extreme vloeistofvervormingen die voorkomen bij hydrodynamische instabiliteiten.
• Stralingstransport: Veel bestaande SPH-codes (Smoothed Particle Hydrodynamics) gebruiken vereenvoudigde benaderingen voor straling of vereisen dure Monte Carlo-methoden die niet schaalbaar zijn voor laser-plasma simulaties.

Er is behoefte aan een robuuste, roostervrije methode die in staat is om zowel de hydrodynamica als de stralingstransport en laser-energiekoppeling nauwkeurig te modelleren in deze extreme omstandigheden.

Methodologie: De SPRAY Code

De auteurs hebben SPRAY (Smoothed Particle RAdiation hYdrodynamics) ontwikkeld, een massaal parallelle, GPU-versnelde code die gebaseerd is op de SPH-methode.

Kerncomponenten en Numerieke Methoden:

1. SPH Formuleren:

   • SPRAY is volledig roostervrij en Lagranse (deeltjes volgen de vloeistof).
   • Het gebruikt de Wendland C2-kernfunctie voor stabiliteit tegen koppelingsinstabiliteiten.
   • De vergelijkingen voor impulsbehoud en energiebehoud zijn geformuleerd om energiebehoud exact te garanderen, zelfs bij schokgolven, door het gebruik van geconjugeerde paren voor gradiëntbenaderingen.

2. Drie-Temperatuur Model:

   • In plaats van een enkel temperatuurmodel, gebruikt SPRAY een één-vloeistof drie-temperatuur model (ionen, elektronen en straling). Dit is cruciaal omdat tijdens intense laser-plasma interacties de temperatuurverschillen tussen ionen en elektronen groot kunnen zijn, en straling een significante rol speelt.
   • De energieoverdracht tussen de componenten wordt gemodelleerd via elektron-phonon botsingen (voor koude vastestoffen) en Spitzer-Härm botsingen (voor heet plasma).

3. Stralingstransport:

   • De code gebruikt de flux-gelimiteerde diffusiebenadering (flux-limited diffusion), wat geschikt is voor optisch dikke plasma's.
   • De stralingstransportvergelijkingen worden impliciet opgelost (via een Jacobi-iteratie met een BiCGSTAB-oplosser op GPU's), terwijl de hydrodynamica expliciet wordt opgelost (operator splitting).
   • Een unieke aanpak behoudt de niet-lineaire kwartische afhankelijkheid van de Planck-energiedichtheid van de elektronentemperatuur in de impliciete solver, wat nauwkeuriger is dan lineaire benaderingen.

4. Laser-Energiekoppeling:

   • Laserabsorptie wordt gemodelleerd via inverse bremsstrahlung in de WKB-benadering (Wentzel-Kramers-Brillouin).
   • Een nieuwe, roostervrije straal-tracing module is ontwikkeld. Laserstralen worden behandeld als hypothetische deeltjes die door het elektronendichtheidsveld bewegen. De stralen worden gebroken (refractie) op basis van de dichtheidgradiënt, berekend via SPH-interpolatie. Dit elimineert de noodzaak van een onderliggend rooster.

5. Randvoorwaarden en Stabiliteit:

   • Vrij oppervlak: Een innovatief schema voor het volgen van het vrije oppervlak en een deeltjes-mirroring-methode worden gebruikt om de dichtheidsgradiënten aan de randen nauwkeurig te schatten en onderbepaling te voorkomen.
   • Numerieke dissipatie: Een energiebehoudend kunstmatige viscositeitsschema (gebaseerd op Price) wordt gebruikt om schokgolven te vangen zonder onfysische oscillaties.
   • Anisotrope Kernel: Om de sterke uitrekking van deeltjes tijdens ablatie te compenseren, wordt een ellipsoïdale kern gebruikt die zich aanpast aan de vervormingstensor van het snelheidsveld.

6. GPU Parallelisatie:

   • De code is geschreven in CUDA C/C++ en maakt gebruik van een Cell-based Nearest Neighbor Particle Search (NNPS) algoritme voor efficiënte $O(N \log N)$ zoekopdrachten.
   • Het ondersteunt multi-GPU uitvoering met geoptimaliseerde geheugenoverdracht tussen GPU's.

Belangrijkste Bijdragen

• Eerste toepassing van SPH in HEDP: SPRAY is, naar weten van de auteurs, de eerste code die de SPH-methode toepast voor simulaties van laser-plasma interacties in de HEDP-wereld, met gedetailleerde fysische modellering en materiaaldata.
• Roostervrije Laser-Raytracing: De ontwikkeling van een volledig mesh-free ray-tracing schema dat gekoppeld is aan de SPH-deeltjes, wat ideaal is voor willekeurige geometrieën en complexe vervormingen.
• Geavanceerde Randbehandeling: Een nieuw schema voor het volgen van vrij oppervlak en het corrigeren van gradiënten aan de randen, essentieel voor het modelleren van ablatie en uitdijend plasma.
• Impliciete Stralingstransport op GPU: Een efficiënte implementatie van een impliciete solver voor stralingstransport met behoud van de niet-lineariteit van de Planck-functie, gekoppeld aan een GPU-versnelde lineaire solver.

Resultaten en Validatie

De code is gevalideerd via een reeks benchmarkproblemen:

1. Sod Shock Tube: Bevestigt de robuustheid van het schema voor het vangen van schokgolven en het behoud van energie. De resultaten komen overeen met de analytische oplossing.
2. Laser-Irradiatie van Aluminium: Vergelijking met de gevestigde code MULTI-IFE. SPRAY levert bijna identieke resultaten op voor dichtheid, snelheid en temperatuurprofielen, wat de nauwkeurigheid van de laser-energiekoppeling en stralingstransport aantoont.
3. Rayleigh-Taylor Instabiliteit: Vergelijking met de ATHENA code (Eulerisch). SPRAY reproduceert succesvol de groeisnelheden van "spikes" en "bubbles" en de kwantitatieve resultaten komen overeen met lineaire theorie en de referentiecode.
4. Implosie Benchmark: Een grote simulatie (2,8 miljoen deeltjes) van een imploderende schaal toont de vorming van instabiliteiten en turbulente menging, wat de schaalbaarheid en robuustheid van de code bewijst.
5. ICF Doelwitcompressie: Een simulatie van een laser-gedreven implosie van een DT-ice schaal. De resultaten tonen een nauwkeurige koppeling van laserenergie en een correct implosiegedrag, vergeleken met MULTI-IFE.

Betekenis en Toekomstperspectief

SPRAY vertegenwoordigt een doorbraak in de numerieke modellering van HEDP. Door de voordelen van SPH (roostervrij, uitstekend in het hanteren van extreme vervormingen) te combineren met geavanceerde stralingstransport en laser-energiekoppeling, biedt het een krachtig alternatief voor traditionele rooster-gebaseerde methoden.

Beperkingen en Toekomstig Werk:

• Huidige beperkingen omvatten de uitdagingen bij het modelleren van faseovergangen (vast naar plasma) en het gebruik van een "gray" (grijze) benadering voor straling in plaats van multi-groep.
• Toekomstige ontwikkelingen omvatten het modelleren van laserbundel-bundel interacties (zoals CBET) en de implementatie van een geavanceerde multi-groep stralingstransportmodule voor nog hogere fysieke nauwkeurigheid.

Samenvattend opent SPRAY nieuwe mogelijkheden voor het bestuderen van complexe laser-plasma fenomenen met behulp van de snel evoluerende rekenkracht van GPU's en de flexibiliteit van de SPH-methode.