Numerical thermalization in $n$-D particle-in-cell simulations

Technische Samenvatting: Numerieke Thermalisatie in n-D Particle-in-Cell Simulaties

Probleemstelling
Particle-in-cell (PIC) simulaties worden breed gebruikt om kinetisch plasma-gedrag te modelleren door fysieke deeltjes te representeren als macrodeeltjes met eindige vormfuncties. Hoewel vaak wordt aangenomen dat ze de botsingsloze Vlasov-vergelijking oplossen door de effecten van grid-interpolatie, lijden PIC-simulaties vaak aan een kunstmatig hoge collisionaliteit. Deze "numerieke thermalisatie" ontstaat wanneer het gewicht van het macrodeeltje (de ratio tussen fysieke deeltjes en macrodeeltjes) groot is, wat de interactiekrachten tussen deeltjes versterkt ondanks de grid-smoothing. In hogere dimensies (2D en 3D) is het vaak onhaalbaar om een numerieke thermalisatietijd ($\tau_{num}$) te bereiken die langer is dan de fysieke thermalisatietijd ($\tau_{phys}$), zelfs wanneer het grid de Debye-lengte oplost. Deze kunstmatige collisionaliteit vervormt de snelheidsverdelingen en beïnvloedt cruciale grootheden (bijv. transport, emissie, doorslag) over simulatietijdschalen $t \gtrsim \tau_{num}$, wat potentieel kan leiden tot valse convergentie of incorrecte fysische voorspellingen.

Methodologie
De auteurs valideren een kinetische theorie voor numerieke botsingstijden tegenover elektrostatische PIC-simulaties in 1D, 2D en 3D dimensionaliteiten.

1. Simulaties: Single-species elektrostatische PIC-simulaties werden uitgevoerd met uniforme initiële condities. De snelheid-autocorrelatiefunctie (VACF), $Z(t)$, werd berekend uit de trajecten van de deeltjes. De numerieke botsingstijd werd gedefinieerd als de exponentiële vervaltijd van $Z(t)$.
2. Parameters: De simulaties varieerden het gewicht van het macrodeeltje ($W$), de grid-afstand ($\Delta$) en de straal van de macrodeeltjesvorm ($R$). Zowel ongefilterde simulaties (waarbij $R = \Delta/2$) als gefilterde simulaties (waarbij $R > \Delta/2$ via $k$-ruimte filtering) werden uitgevoerd om effecten van de vorm te isoleren van grid-heating.
3. Theorie: De studie maakt gebruik van de kinetische theorie ontwikkeld door Okuda, Birdsall en Langdon. Dit raamwerk incorporeert de vormfunctie en het gewicht van het macrodeeltje in een Lenard-Balescu-achtige kinetische vergelijking. De theorie leidt botsingsoperatoren af voor gevormde macrodeeltjes, waarbij rekening wordt gehouden met de modificatie van de inter-deeltjes krachtkern door de vormfunctie $S(k)$ en het gewicht $W$.
4. Vergelijking: Theoretische voorspellingen voor de botsingstijd werden afgeleid met behulp van dynamische, statische en traag-bewegende benaderingen van de drag-integraal. Deze werden direct vergeleken met de vervalsnelheden waargenomen in de PIC-simulaties.

Belangrijkste Bijdragen en Resultaten

• Validatie van Kinetische Theorie: Het artikel toont aan dat de kinetische theorie van Okuda, Birdsall en Langdon de numerieke botsingstijd in 1D, 2D en 3D PIC-simulaties nauwkeurig voorspelt. De theorie vangt correct de afhankelijkheid van $\tau_{num}$ van het gewicht van het macrodeeltje, de vormgrootte en dimensionaliteit.
• Interpretatie van de Krachtkern: De studie verheldert dat numerieke botsingen worden gedreven door een effectieve inter-deeltjes kracht waarbij de eindige vormfunctie de kracht verzwakt en gladstrijkt, terwijl het gewicht van het macrodeeltje deze juist versterkt. De theorie toont aan dat de ratio van numerieke naar fysieke botsingstijden onafhankelijk is van de massa- en ladingratio's van de soort, en primair afhangt van $W$ en $R$.
• Dimensionaliteitseffecten: De resultaten benadrukken dat het bereiken van de fysieke thermalisatietijd aanzienlijk uitdagender is in 3D dan in 1D of 2D voor macrodeeltjesgroottes die de Debye-lengte oplossen. In 3D schaalt de botsingstijd logaritmisch met de plasma-parameter $\Lambda$ in het limiet van een kleine deeltjesstraal ($R \ll \lambda_D$).
• Trade-off Analyse: Het artikel kwantificeert de afweging tussen het verminderen van numerieke thermalisatie en het behouden van fysieke nauwkeurigheid. Het vergroten van de macrodeeltjesgrootte ($R$) om botsingen te verminderen (door krachten te verzachten) verzwakt simultaan het gemiddelde Vlasov-veld en verandert botsingsloze processen (bijv. golfdispersie en demping). Er bestaat een "vervormd veld" limiet waarbij het verder vergroten van $R$ de fysische beperkingen op ruimtelijke resolutie schendt.
• Hogere-orde Vormen: De analyse breidt zich uit naar hogere-orde vormfuncties (bijv. B-splines). Er wordt getoond dat het verhogen van de spline-orde wiskundig equivalent is aan het vergroten van de effectieve deeltjesbreedte. Hieruit volgt dat het gebruik van hogere-orde vormen de schaling van numerieke thermalisatie niet fundamenteel verandert, maar de effectieve deeltjesgrootte effectief vergroot.
• Toepassing op Reële Systemen: De auteurs passen de afgeleide restricties toe op verschillende bestaande PIC-studies (bijv. ICF hotspots, magnetische reconnectie, Penning-ontladingen, Hall-thrusters, streamer-ontladingen en ITER neutrale bundel injectoren). In veel 3D-gevallen wordt vastgesteld dat de numerieke botsingstijd korter is dan de fysieke botsingstijd, wat suggereert dat numerieke thermalisatie dominant kan zijn of de resultaten significant kan beïnvloeden, tenzij specifieke mitigatiestrategieën (zoals Monte Carlo-botsingen of extreme grid-verfijning) worden toegepast.

Betekenis
De primaire betekenis van dit werk is het bieden van een gevalideerd, voorspellend raamwerk voor het schatten van de tijdschaal van zelfconsistente Coulomb-interacties in PIC-simulaties. Door vast te stellen dat de kinetische theorie $\tau_{num}$ nauwkeurig voorspelt als functie van de grid-afstand en het gewicht van het macrodeeltje, biedt het artikel PIC-gebruikers de middelen om:

1. Relevantie te Karakteriseren: A priori bepalen of numerieke thermalisatie de specifieke fysica van belang zal beïnvloeden (bijv. of $\tau_{num}$ korter is dan de simulatietijd of de fysieke botsingstijd).
2. Restricties Vast te Stellen: Toelaatbare operationele regimes definiëren voor het gewicht en de vorm van het macrodeeltje om te waarborgen dat numerieke artefacten de fysische fenomenen niet maskeren.
3. Valse Convergentie te Vermijden: Erkennen dat convergentie met betrekking tot deeltjes per cel niet garandeert dat numerieke thermalisatie wordt geëlimineerd; er zijn voldoende macrodeeltjes nodig om de fysieke botsingssnelheid op te lossen, niet alleen om ruiskenmerken op te lossen.

Het artikel concludeert dat hoewel PIC-methoden krachtig zijn, de "botsingsloze" aanname vaak ongeldig is in 3D-simulaties met standaardparameters. Gebruikers moeten expliciet rekening houden met numerieke thermalisatie om de getrouwheid van langdurige transport- en collisionele processen te waarborgen.