Entity -- Hardware-agnostic Particle-in-Cell Code for Plasma Astrophysics. III: Higher-order shape functions & generalized field stencils

Technische Samenvatting: Entiteit – Hardware-onafhankelijke Particle-in-Cell-code voor Plasma-Astrofysica. III: Hogere-orde vormfuncties en gegeneraliseerde veldstencils

Probleemstelling
Particle-in-Cell (PIC)-simulaties zijn essentieel voor het modelleren van astrofysische plasma-processen, variërend van turbulentie en magnetische reconnectie tot versnelling van kosmische straling. Deze simulaties staan echter voor aanzienlijke uitdagingen met betrekking tot numerieke ruis, kunstmatige verwarming en de numerieke Cherenkov-instabiliteit. Deze artefacten ontstaan vaak door interpolatieschema's van lage orde en standaard eindige-differentiestencils, die fysische effecten kunnen verduisteren of leiden tot spurious versnelling van deeltjes. Hoewel het verhogen van de numerieke resolutie (roostercellen of deeltjes per cel) een gebruikelijke mitigatiestrategie is, wordt deze beperkt door computercapaciteit. De auteurs identificeren een behoefte aan methoden van hogere orde om ruis te onderdrukken en stabiliteit te verbeteren zonder noodzakelijkerwijs prohibitieve computerkosten te veroorzaken, met name in de context van de Entity PIC-code, die is ontworpen voor exascale GPU-clusters.

Methodologie
Dit werk (Paper III van de Entity-serie) implementeert en test twee primaire vooruitgangen binnen het Entity-kader:

1. Hogere-orde stroomdepositie en veldinterpolatie:
   De auteurs implementeren het ladingsbehoudende stroomdepositieschema voorgesteld door Esirkepov (2001), en breiden dit uit om vormfuncties van willekeurige orde tot de elfde orde ($O=11$) te ondersteunen.

   • Vormfuncties: De code maakt gebruik van van B-splines afgeleide vormfuncties die zijn opgebouwd via opeenvolgende convolutie. Deze functies voldoen aan de voorwaarden voor symmetrie, normalisatie, positiviteit, eindige uitgestrektheid en gladheid.
   • Consistentie: Om consistentie met de stroomdepositie te waarborgen, wordt de veldinterpolatie van het rooster naar deeltjesposities eveneens bijgewerkt om dezelfde hogere-orde vormfuncties te gebruiken, ter vervanging van eerdere indirecte interpolatieschema's.
   • Implementatienoot: Ordes $O=10$ en $O=11$ vereisen dubbele precisie-aritmetiek vanwege kleine numerieke waarden aan de randen van de stencil, wat hun gebruik beperkt op bepaalde GPU-architecturen.

2. Gegeneraliseerde veldstencils voor de Maxwell-oplosser:
   De Faraday-oplosser (die het magnetische veld evolueert via $\partial_t B = -c\nabla \times E$) wordt bijgewerkt om gegeneraliseerde eindige-differentiestencils te gebruiken.

   • Optimalisatie: In plaats van de standaard Yee-stencil van tweede orde, implementeren de auteurs stencils die gebruikmaken van aanvullende aangrenzende roosterpunten. De coëfficiënten van deze stencils zijn geoptimaliseerd om de afwijking tussen de numerieke fase-snelheid en de fysieke lichtsnelheid in de Fourier-ruimte te minimaliseren.
   • Dispersiecontrole: Deze aanpak beoogt de numerieke Cherenkov-instabiliteit te onderdrukken, die optreedt wanneer relativistische deeltjes sneller reizen dan de numerieke lichtsnelheid op het rooster.
   • Behoud van divergentie: De auteurs waarborgen dat specifieke stencilconfiguraties voldoen aan de voorwaarde $\beta_{ij} = \delta_j$ om $\nabla \cdot B = 0$ automatisch te handhaven gedurende de hele simulatie.

Belangrijkste bijdragen en resultaten
De auteurs voeren uitgebreide convergentietests uit om de nieuwe schema's te valideren tegen standaardimplementaties (Zigzag-depositie en Yee-stencil) en diverse filterconfiguraties.

• Ladingsbehoud: Het bijgewerkte Esirkepov-schema toont uitstekend ladingsbehoud over alle geteste ordes ($O=1$ tot $O=11$), waarbij de nauwkeurigheid voornamelijk wordt beperkt door afrondingsfouten in plaats van de orde van de vormfunctie.
• Remvermogen: In tests die de elektrostatische remkracht op relativistische deeltjes meten, verminderen vormfuncties van hogere orde de numerieke energieverlies aanzienlijk. Een overstap van 1e naar 9e orde verlaagt het effectieve remvermogen, wat overeenkomt met het verhogen van de numerieke resolutie met een factor 3–4.
• Numerieke verwarming en stabiliteit:
  • Koude plasma's: Voor koude plasma's met onder-resolvene Debye-lengten verbeteren schema's van hogere orde de stabiliteit drastisch. Een overstap van 1e naar 2e orde verlaagt de vereiste roosterresolutie met een factor $\sim4$; een overstap naar 3e orde verlaagt deze met $\sim8$.
  • Oblique modi: In tests met kosmische straalbundels onderdrukken vormen van hogere orde de groei van numerieke aliasing-ruis. Hoewel 6e orde lijkt voldoende voor convergentie in deze specifieke tests, vertragen hogere ordes het begin van instabiliteit verder.
  • Vervalende turbulentie: In 2D vervalende turbulentie verbeteren vormfuncties van hogere orde het energiebehoud, waardoor de totale energiefout met maximaal twee ordes van grootte wordt verminderd in vergelijking met de standaard Zigzag-depositie.
• Prestaties van veldstencils:
  • Enkeldeltje: Tests met een enkel relativistisch elektron tonen aan dat geoptimaliseerde stencils (bijv. Cowan, min3, min6) numerieke Cherenkov-artefacten aanzienlijk verminderen in vergelijking met de standaard Yee-stencil.
  • Drijvende plasma's: In warme, drijvende plasmasimulaties onderdrukken de min3- en min6-stencils effectief de groei van de Cherenkov-instabiliteit en de bijbehorende numerieke verwarming.
  • Relativistische shocks: In een relativistische paar-shocksimulatie stabiliseert de combinatie van de min3-stencil en vormfuncties van 3e orde het stroomopwaartse plasma, waardoor de kunstmatige verwarming en stroomfilamentatie die worden gezien met de Yee-stencil, worden voorkomen. Dit leidt tot nauwkeurigere shock-rampstructuren en deeltjesversnellingsspectra.
• Prestatieafwegingen:
  • 1D: De prestatieimpact van schema's van hogere orde is verwaarloosbaar.
  • 2D: De kosten nemen gematigd toe (bijv. $\sim2,5\times$ voor 3e orde, $\sim24\times$ voor 11e orde).
  • 3D: De kosten worden prohibitief voor ordes $O > 3$ (bijv. $\sim5\times$ voor 3e orde).
  • Mitigatie: De auteurs tonen aan dat de toegenomen computerkosten van interpolatie van hogere orde kunnen worden gecompenseerd door het aantal stroomfilterpasses of het aantal deeltjes per cel te verminderen, aangezien schema's van hogere orde inherent meer stabiliteit bieden.

Betekenis en claims
Het artikel stelt dat hoewel geen enkele methode alle numerieke uitdagingen in PIC-simulaties oplost, de integratie van vormfuncties van hogere orde en gegeneraliseerde veldstencils in Entity een robuust toolkit biedt voor het verbeteren van de simulatiefideliteit. De auteurs benadrukken dat:

1. Stabiliteit versus resolutie: Methoden van hogere orde stabiele simulaties van koude plasma's en relativistische shocks mogelijk maken bij lagere resoluties dan voorheen mogelijk was, wat potentieel de totale computerkosten kan verminderen ondanks de overhead per stap.
2. Onderdrukking van ruis: Deze methoden onderdrukken effectief numerieke Cherenkov-effecten en aliasing-ruis, wat cruciaal is voor het bestuderen van fysische instabiliteiten zoals de Bell-instabiliteit of stroming van kosmische straling.
3. Flexibiliteit: De gegeneraliseerde stencils stellen gebruikers in staat de oplosser af te stemmen op specifieke dispersierelaties of directionele symmetrieën (bijv. het gebruik van de min3-stencil voor algemene stabiliteit of Cowan/min6 voor grote CFL-getallen).
4. Praktische aanbevelingen: De auteurs raden het gebruik van vormfuncties van 2e of 3e orde aan als een praktische balans tussen nauwkeurigheid en prestaties, en het gebruik van de min3 (2D) of min3D2 (3D) stencils voor simulaties voor algemeen gebruik.

Het werk concludeert dat deze verbeteringen nu beschikbaar zijn in Entity v1.3.0, waardoor nauwkeurigere en stabielere modellering van complexe astrofysische plasma-verschijnselen op exascale-hardware mogelijk wordt.