Local cubic spline interpolation for Vlasov-type equations on a multi-patch geometry

Technische Samenvatting: Lokale Cubische Spline-interpolatie voor Vlasov-type Vergelijkingen op een Multi-patch Geometrie

Probleemstelling
Het artikel behandelt de numerieke resolutie van Vlasov-type vergelijkingen, specifiek guiding-center modellen, op complexe multi-patch geometrieën die typerend zijn voor tokamak-simulaties. In deze simulaties zijn roosterlijnen uitgelijnd met magnetische fluxoppervlakken om variaties langs veldlijnen te minimaliseren. Deze geometrie creëert specifieke uitdagingen: de kern bevat een O-punt waar roosterlijnen samenkomen, en de rand bevat een X-punt (separatrix) waar lijnen uiteenlopen. Hoewel het opsplitsen van het globale domein in subdomeinen (patches) het mogelijk maakt om lokale verfijning en geoptimaliseerde uniforme roosters binnen elke patch toe te passen, introduceert dit problemen bij het behoud van de oplossingsregulariteit op de patch-interfaces. Eerdere studies naar $C^0$-continuïteit (waarbij alleen functiewaarden worden gematcht) bij interfaces bleken potentieel instabiel voor semi-Lagrangiaanse advectie-schema's. Het kernprobleem is hoe men hogere-orde continuïteit ($C^1$) kan afdwingen en interface-afgeleiden nauwkeurig kan berekenen voor lokale cubische B-spline interpolaties op willekeurige (uniforme en niet-uniforme) multi-patch roosters zonder het rendement van lokale berekeningen op te offeren.

Methodologie
De auteurs stellen een backward semi-Lagrangiaanse (BSL) methode voor die gebruikmaakt van lokale cubische B-splines gedefinieerd op elke patch. Om stabiliteit en nauwkeurigheid te garanderen, dwingt de methode $C^1$-regulariteit af op patch-interfaces door zowel functiewaarden als hun eerste afgeleiden te matchen.

1. Reconstructie van Afgeleiden: De centrale innovatie is een gegeneraliseerde formule voor het berekenen van interface-afgeleiden. Bouwend op een formule van Crouseilles et al. [11], ontworpen voor uniforme roosters, breiden de auteurs de aanpak uit naar willekeurige niet-uniforme roosters met behulp van Hermite-polynomen.

   • Exacte Formulering: Door $C^2$-continuïteit op te leggen tussen lokale Hermite-polynomen gedefinieerd op aangrenzende cellen, leiden de auteurs een recursieve relatie af die de afgeleide op een interface koppelt aan de afgeleiden op naburige interfaces en de functiewaarden. Dit resulteert in een klein globaal lineair systeem (omvang $(N_p-1) \times (N_p-1)$ voor $N_p$ patches) dat moet worden opgelost om afgeleiden te verkrijgen die equivalent zijn aan die van een globale spline.
   • Benaderde Formulering: De auteurs analyseren het asymptotische gedrag van de coëfficiënten in het globale systeem. Zij demonstreren dat de invloed van verre interfaces exponentieel afneemt naarmate het aantal cellen per patch toeneemt. Bij voldoende grote patches (bijv. $\ge 30$ cellen) wordt de afhankelijkheid van verre afgeleiden verwaarloosbaar (bereikt machineprecisie). Dit maakt een afgeknotte, lokale benadering mogelijk waarbij afgeleiden worden berekend met behulp van een beperkt aantal naburige cellen, waardoor de globale oplossing wordt vermeden terwijl een hoge precisie behouden blijft.
   • Expliciete Formules: Voor het specifieke geval van uniforme patches bieden de auteurs expliciete gesloten vorm-expressies voor de coëfficiënten, wat de berekening verder optimaliseert.

2. 2D-Extensie: De methode wordt uitgebreid naar 2D multi-patch geometrieën, waarbij drie gevallen worden behandeld:

   • Conforme Roosters: Een directe tensorproduct-toepassing van de 1D-formules om $r$-afgeleiden, $\theta$-afgeleiden en kruis-afgeleiden te berekenen.
   • Niet-conforme Roosters: Waarbij breukpunten niet over patches heen uitgelijnd zijn. De auteurs stellen een algoritme voor dat progressief afgeleiden op grovere interfaces berekent en deze interpoleert op verfijnde patches om $C^1$-continuïteit af te dwingen.
   • T-joints: Verbindingen waarbij drie of meer patches betrokken zijn en een zijde slechts gedeeltelijk verbonden is. De methode past de volgorde van de afgeleideberekening aan op basis van de geometrie om ervoor te zorgen dat alle noodzakelijke afgeleiden beschikbaar zijn.

3. Implementatie: Het algoritme is geïmplementeerd in de Gyselalib++ bibliotheek. De Poisson-vergelijking wordt opgelost met een Finite Element Method (FEM) met cubische B-splines op een equivalente globale spline, terwijl de advectievergelijking wordt opgelost met de voorgestelde BSL-methode op lokale splines.

Belangrijkste Bijdragen

• Generalisatie van Interface-afgeleideformules: Het artikel generaliseert de berekening van de afgeleide-formule van uniforme naar willekeurige niet-uniforme roosters met behulp van Hermite-polynomen, wat het afhandelen van niet-conforme situaties en lokale verfijning mogelijk maakt.
• Stabiliteitsanalyse: De auteurs tonen aan dat $C^0$-regulariteit onvoldoende is voor stabiliteit in BSL-schema's op multi-patch domeinen en dat $C^1$-regulariteit vereist is.
• Trade-off tussen Lokaal en Globaal: Het werk stelt vast dat hoewel een exacte globale oplossing $C^2$-continuïteit garandeert, een lokale benadering met een afgetopte stencil (bijv. 30 cellen) een machineprecisie-equivalent bereikt van de globale oplossing, wat een computationeel efficiënt alternatief biedt.
• Afhandeling van Complexe Geometrieën: De methode handelt T-joints en niet-conforme roosters succesvol af, die noodzakelijk zijn voor het modelleren van complexe tokamak-doorsneden die zowel O-punten als X-punten bevatten.

Numerieke Resultaten
De methode werd gevalideerd met de volgende tests in de Gyselalib++ bibliotheek:

1. Spline-interpolatie: Tests op concentrische patches (uniform en niet-uniform) bevestigden dat de lokale spline-representaties convergeren naar de equivalente globale spline met orde 4. De fout in interface-afgeleiden met de afgetopte benadering kwam overeen met de exacte formule zodra het aantal cellen de $\approx 30$ overschreed, waarbij machineprecisie ($10^{-13}$) werd bereikt.
2. Advectie-tests:
   • Niet-conforme patches: De lokale splines op een niet-conform rooster bereikten een nauwkeurigheid die vergelijkbaar is met een verfijnd globaal rooster, ondanks het gebruik van minder totale cellen.
   • T-joints: In een 2D T-joint configuratie kwamen de lokale spline-resultaten overeen met de globale spline-resultaten tot machineprecisie ($\approx 6 \times 10^{-13}$).
3. Guiding-Center Model (Diocotron Instabiliteit):
   • Simulaties op conforme multi-patch roosters reproduceerden de resultaten van globale spline-simulaties (groeisnelheden en massaconservatie) met fouten die consistent zijn met de interpolatieprecisie.
   • Een simulatie met een niet-optimaal rooster (gedepopuleerde centrale patch en verschoven interfaces) toonde aan dat hoewel de keuze van het rooster de absolute foutgrootte beïnvloedde, de lokale spline-methode het gedrag van de globale oplossing correct volgde, wat de robuustheid van de aanpak valideert.

Betekenis en Claims
Het artikel claimt een robuust kader te bieden voor het oplossen van Vlasov-type vergelijkingen op multi-patch geometrieën, die essentieel zijn voor realistische tokamak-simulaties. De primaire betekenis ligt in het mogelijk maken van lokale cubische spline-interpolatie met Hermite-randvoorwaarden om stabiliteit en hoge-orde nauwkeurigheid over patch-interfaces te handhaven.

De auteurs benadrukken dat hun methode de beperking wegneemt die in eerdere studies werd gevonden met betrekking tot het aantal punten per patch in conforme gevallen, door deze te vervangen door de vereiste om een klein globaal systeem (of een lokale benadering) op te lossen. Zij merken op dat de methode momenteel gevalideerd is voor geometrieën met O-punten. Hoewel de introductie van T-joints een stap is naar het verbinden van O-punt en X-punt geometrieën, geven de auteurs expliciet aan dat het afhandelen van X-punt singulariteiten verdere studie vereist en nog niet volledig is opgelost in dit werk. De methode wordt gepresenteerd als een belangrijke stap naar efficiënte, hoge-orde semi-Lagrangiaanse solvers voor complexe magnetische opsluitingsgeometrieën.