Decoupled iterative schemes for solving stationary MHD problems

Technische Samenvatting: Ontkoppelde Iteratieve Schema's voor Stationaire MHD-problemen

Probleemstelling
Het artikel behandelt de numerieke oplossing van het stationaire, incompressibele magnetohydrodynamica (MHD) systeem. Dit systeem modelleert de interactie tussen elektrisch geleidende vloeistoffen en elektromagnetische velden, die worden beheerst door de koppeling van de Navier-Stokes-vergelijkingen en de Maxwell-vergelijkingen. De wiskundige formulering behelst het vinden van het snelheidveld ($u$), het magnetische veld ($B$), de gemodificeerde druk ($p$) en een Lagrange-multiplier ($\lambda$) die voldoen aan een reeks gekoppelde partiële differentiaalvergelijkingen. De primaire uitdaging ligt in de sterke nietlineariteit en de nauwe koppeling tussen de fluïdum- en magnetische variabelen, wat de constructie van efficiënte en stabiele numerieke solvers bemoeilijkt.

Methodologie
De auteurs stellen twee nieuwe ontkoppelde iteratieve schema's voor en analyseren deze, geïnspireerd door de Incremental Viscosity Splitting (IVS) strategie, die oorspronkelijk is ontwikkeld voor onstabiele Navier-Stokes-stromingen.

1. Standaard IVS-schema (Algoritme 5.1):

   • Splitsingsstrategie: De methode splitst de impuls- en inductievergelijkingen met betrekking tot hun diffusieve termen.
   • Stap 1: Lost voor intermediaire variabelen ($\tilde{u}, \tilde{B}$) door de convectieve en koppelende termen expliciet te behandelen met behulp van waarden uit de vorige iteratie, terwijl de druk en de Lagrange-multiplier van de vorige stap behouden blijven. Dit resulteert in een gekoppeld lineair systeem voor de intermediaire velden.
   • Stappen 2 & 3: Corrigeert de intermediaire velden om de uiteindelijke snelheid ($u$) en het magnetische veld ($B$) te verkrijgen door twee afzonderlijke, ontkoppelde Stokes-type systemen op te lossen. Deze systemen handhaven de divergentievrije restricties ($\nabla \cdot u = 0, \nabla \cdot B = 0$) en updaten de druk en de Lagrange-multiplier.
   • Lineair Algebraïsch Voordeel: Een belangrijk kenmerk is dat de Schur-complementmatrices voor de Stokes-correcties onafhankelijk zijn van het iteratieniveau ($k$). Bijgevolg hoeven ze slechts één keer geassembleerd en gepreconditioneerd te worden, wat de computationele overhead vermindert in vergelijking met methoden waarbij de systeemmatrix bij elke stap verandert.

2. Elsässer IVS-schema (Algoritme 6.1):

   • Herformulering: Om het systeem verder te ontkoppelen, herformuleren de auteurs de MHD-vergelijkingen met behulp van Elsässer-variabelen ($w = u + \sqrt{\kappa}B$ en $z = u - \sqrt{\kappa}B$).
   • Ontkoppeling: In deze formulering vindt de koppeling tussen de nieuwe variabelen plaats via één nietlineaire en één lineaire term. Het voorgestelde algoritme lost $w$ en $z$ in afzonderlijke stappen op, waarbij het grote gekoppelde systeem van de standaard IVS wordt vervangen door twee kleinere, volledig ontkoppelde lineaire systemen.
   • Toepasbaarheid: Deze benadering is theoretisch voordelig onder homogene Dirichlet-randvoorwaarden, hoewel de auteurs opmerken dat er uitdagingen ontstaan bij gemengde randvoorwaarden.

Theoretische Bijdragen
Het artikel biedt een rigoureuze wiskundige analyse voor beide schema's:

• Wel-gesteltheid: De auteurs vestigen resultaten voor existentie en uniciteit voor het continue MHD-probleem onder een "small data" aanname (specifiek, een begrenzing op de Reynolds- en magnetische Reynolds-getallen ten opzichte van de viscositeit).
• Convergentiebewijzen:
  • Het klassieke Picard-iteratieve schema wordt geanalyseerd en bewezen te convergeren onder dezelfde "small data" aanname.
  • Zowel het Standaard IVS- als het Elsässer IVS-schema worden bewezen uniform begrensd en geometrisch convergerend naar de unieke oplossing.
  • De bewijzen maken gebruik van energie-schattingen, de inf-sup (LBB) conditie en specifieke lemma's betreffende de convergentie van niet-negatieve sequenties.
• Foutanalyse: Theoretische grenzen worden afgeleid die aantonen dat de fouten in snelheid, magnetisch veld, druk en Lagrange-multipliers naar nul dalen naarmate het aantal iteraties toeneemt.

Numerieke Resultaten
De auteurs valideren hun theoretische bevindingen met behulp van de FreeFem++ software met $P_2-P_1$ eindele elementparen.

• Convergentietests: Gemanufactureerde oplossingen in 2D demonstreren optimale convergentiecijfers ($O(h^2)$ in $L^2$ en $O(h)$ in $H^1$) voor zowel de snelheid als de magnetische velden, wat de nauwkeurigheid van de schema's bevestigt.
• Benchmark Problemen:
  • Hartmann-stroom: Zowel 2D als 3D Hartmann-stroomproblemen werden opgelost. De numerieke resultaten kwamen met hoge nauwkeurigheid overeen met analytische oplossingen voor diverse Hartmann-getallen ($Ha$).
  • Lid-Driven Cavity: De schema's simuleerden succesvol complexe stromingspatronen voor diverse Reynolds ($Re$), magnetische Reynolds ($Rm$) en koppelingsgetallen ($\kappa$), inclus\n bij hoge koppelingsgevallen ($\kappa=5000$).
  • Stroming over een Trede: De methode legde stromingsseparatie en drukval accuraat vast, wat de robuustheid demonstreert bij het hanteren van geometrische singulariteiten en variërende koppelingssterktes.

Betekenis en Claims
Het artikel claimt dat de voorgestelde ontkoppelde iteratieve schema's een praktisch en efficiënt alternatief bieden voor klassieke monolithische of volledig gekoppelde iteratieve methoden (zoals Newton- of standaard Picard-methoden).

• Efficiëntie: Door de diffusieve termen te splitsen, vermindert de methode de complexiteit van de lineaire systemen die in elke stap worden opgelost. Specifiek staat de Elsässer-formulering het oplossen van twee kleinere, ontkoppelde systemen toe in plaats van één groot gekoppeld systeem.
• Herbruikbaarheid: De IVS-benadering zorgt ervoor dat de Schur-complementmatrices voor de Stokes-correcties statisch zijn, wat eenmalige assemblage en preconditioning mogelijk maakt, wat computationeel voordelig is voor 3D-problemen.
• Robuustheid: De schema's behouden de divergentievrije eigenschap van de snelheid en de magnetische velden op het discrete niveau en blijven stabiel onder de vastgestelde "small data" condities.

De auteurs concluderen dat, hoewel het huidige werk zich richt op stationaire problemen, de methodologie een fundament biedt voor toekomstige uitbreidingen naar tijd-afhankelijke (onstabiele) MHD-problemen, met een focus op temporele en ruimtelijke stabiliteitsanalyse.