Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

Technische Samenvatting: Foutanalyse van een divergentiebehoudend gemengd eindig-elementenschema voor de incompressibele Hall-MHD-vergelijkingen

Probleemstelling
Het incompressibele Hall-magnetohydrodynamica (Hall-MHD)-systeem modelleert de interactie tussen stroming en magnetische velden in elektrisch geleidende vloeistoffen, waarbij het Hall-effect wordt opgenomen dat significant wordt op schaalvergelijkbaar met de ionen-huiddiepte. Dit effect introduceert een stijve, sterk niet-lineaire term ($\eta \text{curl}(\text{curl} B \times B)$) in de inductievergelijking, wat leidt tot dispersief gedrag, steile gradiënten en oscillaties op hoge frequentie (fluitgolven). Deze kenmerken vormen aanzienlijke analytische en computationele uitdagingen. Specifiek blijft het globale bestaan van sterke oplossingen in drie dimensies een open probleem, en kampen standaard numerieke methoden vaak met stabiliteit en de strikte vereiste dat het magnetische veld solenoidaal blijft ($\text{div } B = 0$). Schendingen van deze divergentievrije voorwaarde kunnen niet-fysische magnetische monopolen induceren en ernstige numerieke instabiliteiten veroorzaken. Bovendien dissiperen bestaande numerieke benaderingen, zoals die gebaseerd zijn op Scalar Auxiliary Variable (SAV)-methodes of spectrale technieken, vaak gemodificeerde energiefunctionalen in plaats van de ware fysische energie, slagen ze er niet in de divergentievrije voorwaarde exact op het discrete niveau af te dwingen, of missen ze een rigoureuze foutanalyse en geometrische flexibiliteit.

Methodologie
Om deze uitdagingen aan te pakken, stellen de auteurs een structuurbehoudend, lineair, volledig discreet gemengd eindig-elementenschema (FEM) voor voor een Voigt-geregulariseerd Hall-MHD-systeem. De Voigt-regularisatie vervangt de tijdsafgeleide in de impulsvergelijking door een ruimtelijk gefilterde versie $(I - \alpha_1 \Delta)\partial_t u$ en behoudt een term voor elektroneninertie $\alpha_2 \partial_t J$ in de gegeneraliseerde wet van Ohm. Deze wijziging verbetert de welgesteldheid en numerieke stabiliteit, terwijl de essentiële fysische structuur en energiewetten van het oorspronkelijke systeem behouden blijven.

Het numerieke schema is opgebouwd uit de volgende componenten:

1. Gemengde Formulering: Het probleem wordt herschreven met het elektrische veld $E$ en de stroomdichtheid $J$ als primaire variabelen, wat leidt tot een gekoppeld systeem van vergelijkingen voor impuls, inductie, de wet van Ohm en de wet van Ampère.
2. Ruimtelijke Discretisatie: Het schema maakt gebruik van Finite Element Exterior Calculus (FEEC) om ervoor te zorgen dat aan de divergentievrije voorwaarde exact wordt voldaan.
   • De snelheid $u$ en druk $p$ worden benaderd met het MINI-eindig-elementenpaar (laagste-orde inf-sup stabiel).
   • Het magnetische veld $B$ en de stroomdichtheid $J$ worden benaderd met respectievelijk de ruimten van Raviart–Thomas ($RT_0$) en Nédélec ($X_0$) van laagste orde.
   • Deze ruimten vormen een discrete de Rham-complex, wat garandeert dat $\text{div } B_h = 0$ exact op het discrete niveau behouden blijft.
3. Temporele Discretisatie: Er wordt een schuifsymmetrisch, lineair impliciet tijdstapschema gebruikt. De convectieve en Lorentz-termen worden semi-impliciet behandeld met waarden uit de vorige tijdstap om lineariteit en onvoorwaardelijke energiestabiliteit te handhaven.
4. Foutanalyse: De auteurs stellen optimale convergentiesnelheden vast door de fout te ontleden in benaderingsfouten (met behulp van Stokes-projecties en orthogonale projecties) en discrete fouten. Zij maken gebruik van discrete Gaffney-ongelijkheden en energie-schattingen om de niet-lineaire koppeltermen te behandelen.

Belangrijkste Bijdragen

• Structuurbehoudend Schema: Het artikel presenteert een lineair, volledig discreet FEM dat de solenoidale voorwaarde voor het magnetische veld exact afdwingt, waardoor de noodzaak voor projectiestappen of penalty-methodes die de nauwkeurigheid of stabiliteit kunnen aantasten, wordt vermeden.
• Onvoorwaardelijke Energietabiliteit: Het voorgestelde schema voldoet onvoorwaardelijk aan een discrete energiedissipatiewet, die de continue fysische energiewet weerspiegelt, zelfs in aanwezigheid van de stijve Hall-term.
• Rigoureuze Foutschattingen: De auteurs leveren een volledige foutanalyse voor het Voigt-geregulariseerde probleem, waarbij optimale convergentiesnelheden van de orde $O(h^s + \tau)$ worden vastgesteld onder geschikte regulariteitsaannames.
• Uitbreiding tot Ongeregulariseerd Systeem: Cruciaal is dat de analyse zich uitstrekt tot het ongeregulariseerde Hall-MHD-systeem (waarbij $\alpha_1 = \alpha_2 = 0$). De auteurs tonen aan dat convergentie nog steeds kan worden bewezen voor het ongeregulariseerde systeem, mits de weerstand $\sigma > 0$. In het niet-resistieve regime ($\sigma = 0$) blijkt de Voigt-regularisatie een cruciale rol te spelen in de foutschattingen.
• Implementatie in 2,5D en 3D: Het raamwerk is gevalideerd in zowel 2,5D (translatie-invariantie) als volledig 3D-instellingen, waarbij de toepasbaarheid van de methode op complexe geometrieën en fysische verschijnselen wordt aangetoond.

Resultaten
Numerieke simulaties in 2,5D en 3D bevestigen de theoretische bevindingen:

• Convergentie: Het schema bereikt de voorspelde ruimtelijke en temporele convergentiesnelheden voor de snelheid, het magnetische veld en de stroomdichtheid. Hoewel de convergentiesnelheden voor snelheid en stroomdichtheid in sommige gevallen iets verminderd zijn (toegeschreven aan beperkte oplossingregulariteit door de Hall-term), behoudt het magnetische veld optimale convergentie.
• Divergentiebehoud: De maximale divergentiefout ($\max | \text{div } B_h |$) blijft gedurende de simulaties op het niveau van de tolerantie van de solver, wat de exacte behouding van de solenoidale voorwaarde bevestigt.
• Energietabiliteit: De totale energie van het systeem neemt monotoon af in de tijd, in overeenstemming met de theoretische energiedissipatiewet.
• Fysische Betrouwbaarheid: Simulaties van de Arnold–Beltrami–Childress (ABC)-stroom, de Orszag–Tang-wervel en Hall-gedreven reconnectie in een Harris-sheet slagen erin complexe fysische verschijnselen nauwkeurig vast te leggen, waaronder magnetische reconnectie, de vorming van magnetische eilanden, fluxtouwen en de karakteristieke quadrupolaire stroompatronen die geassocieerd zijn met Hall-MHD-dynamica. Het schema lost deze kenmerken nauwkeurig op, zelfs op relatief grove meshes.

Betekenis
Het artikel stelt dat dit werk een wiskundig rigoureus en computationeel robuust raamwerk biedt voor het simuleren van Hall-MHD-systemen. Door Voigt-regularisatie te combineren met een divergentiebehoudende gemengde eindig-elementenformulering, overwinnen de auteurs de dubbele uitdagingen van de niet-lineariteit van de Hall-term en de strikte divergentievrije vereiste. De methode biedt een fysisch consistente benadering die stabiel, convergent is en complexe geometrieën aankan, waarmee een gat wordt opgevuld dat door eerdere methoden is achtergelaten die geen rigoureuze foutanalyse of exact divergentiebehoud bezaten. De resultaten suggereren dat deze aanpak een levensvatbaar hulpmiddel is voor het bestuderen van astrofysische en laboratoriumplasma-verschijnselen waarbij Hall-effecten dominant zijn.