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

Sintesi Tecnica: Analisi dell'Errore di uno Schema agli Elementi Finiti Misti che Preserva la Divergenza per le Equazioni Hall–MHD Incomprimibili

Enunciato del Problema
Il sistema Hall-magnetoidrodinamico (Hall–MHD) incomprimibile modella l'interazione tra il moto dei fluidi e i campi magnetici nei fluidi elettricamente conduttori, incorporando l'effetto Hall che diventa significativo a scale comparabili alla profondità della pelle ionica. Questo effetto introduce un termine rigido e altamente non lineare ($\eta \text{curl}(\text{curl} B \times B)$) nell'equazione di induzione, portando a comportamenti dispersivi, gradienti ripidi e oscillazioni ad alta frequenza (onde fischianti). Queste caratteristiche pongono sfide analitiche e computazionali sostanziali. In particolare, l'esistenza globale di soluzioni forti in tre dimensioni rimane un problema aperto, e i metodi numerici standard spesso faticano con la stabilità e il requisito rigoroso che il campo magnetico rimanga solenoidale ($\text{div } B = 0$). Le violazioni di questa condizione di divergenza nulla possono indurre monopoli magnetici non fisici e innescare gravi instabilità numeriche. Inoltre, gli approcci numerici esistenti, come quelli basati su metodi a Variabile Ausiliaria Scalare (SAV) o tecniche spettrali, spesso dissipano funzionali di energia modificati anziché la vera energia fisica, non riescono a imporre la condizione di divergenza nulla esattamente a livello discreto, o mancano di un'analisi rigorosa dell'errore e di flessibilità geometrica.

Metodologia
Per affrontare queste sfide, gli autori propongono uno schema agli elementi finiti misti (FEM), lineare e completamente discreto, che preserva la struttura, per un sistema Hall–MHD regolarizzato di Voigt. La regolarizzazione di Voigt sostituisce la derivata temporale nell'equazione della quantità di moto con una versione filtrata spazialmente $(I - \alpha_1 \Delta)\partial_t u$ e mantiene un termine di inerzia elettronica $\alpha_2 \partial_t J$ nella legge di Ohm generalizzata. Questa modifica migliora la benposta e la stabilità numerica, preservando al contempo la struttura fisica essenziale e le leggi energetiche del sistema originale.

Lo schema numerico è costruito utilizzando i seguenti componenti:

1. Formulazione Mista: Il problema è riformulato utilizzando il campo elettrico $E$ e la densità di corrente $J$ come variabili primarie, portando a un sistema accoppiato di equazioni della quantità di moto, induzione, legge di Ohm e legge di Ampère.
2. Discretizzazione Spaziale: Lo schema impiega il Calcolo Esterno agli Elementi Finiti (FEEC) per garantire che la condizione di divergenza nulla sia soddisfatta esattamente.
   • La velocità $u$ e la pressione $p$ sono approssimate utilizzando la coppia di elementi finiti MINI (di ordine più basso, stabile inf-sup).
   • Il campo magnetico $B$ e la densità di corrente $J$ sono approssimati utilizzando rispettivamente gli spazi di Raviart–Thomas ($RT_0$) e Nédélec ($X_0$) di ordine più basso.
   • Questi spazi formano un complesso di de Rham discreto, garantendo che $\text{div } B_h = 0$ sia preservato esattamente a livello discreto.
3. Discretizzazione Temporale: Viene utilizzato uno schema di integrazione temporale lineare implicito e antisimmetrico. I termini convettivi e di Lorentz sono trattati semi-implicitamente utilizzando valori dal passo temporale precedente per mantenere la linearità e la stabilità energetica incondizionata.
4. Analisi dell'Errore: Gli autori stabiliscono tassi di convergenza ottimali decomponendo l'errore in errori di approssimazione (utilizzando proiezioni di Stokes e proiezioni ortogonali) ed errori discreti. Utilizzano disuguaglianze di Gaffney discrete e stime energetiche per gestire i termini di accoppiamento non lineari.

Contributi Chiave

• Schema che Preserva la Struttura: Il documento presenta un FEM lineare e completamente discreto che impone il vincolo solenoidale sul campo magnetico esattamente, evitando la necessità di passaggi di proiezione o metodi di penalizzazione che potrebbero degradare accuratezza o stabilità.
• Stabilità Energetica Incondizionata: Lo schema proposto soddisfa una legge di dissipazione dell'energia discreta incondizionatamente, rispecchiando la legge fisica dell'energia continua, anche in presenza del termine rigido di Hall.
• Stime Rigorose dell'Errore: Gli autori forniscono un'analisi completa dell'errore per il problema regolarizzato di Voigt, stabilendo tassi di convergenza ottimali dell'ordine $O(h^s + \tau)$ sotto opportune ipotesi di regolarità.
• Estensione al Sistema Non Regolarizzato: Crucialmente, l'analisi si estende al sistema Hall–MHD non regolarizzato (dove $\alpha_1 = \alpha_2 = 0$). Gli autori dimostrano che la convergenza può essere ancora provata per il sistema non regolarizzato, purché la resistività $\sigma > 0$. Nel regime non resistivo ($\sigma = 0$), si dimostra che la regolarizzazione di Voigt gioca un ruolo cruciale nelle stime dell'errore.
• Implementazione 2.5D e 3D: Il framework è validato sia in impostazioni 2.5D (invarianza traslazionale) che completamente 3D, dimostrando l'applicabilità del metodo a geometrie complesse e fenomeni fisici.

Risultati
Le simulazioni numeriche in 2.5D e 3D confermano i risultati teorici:

• Convergenza: Lo schema raggiunge i tassi di convergenza spaziale e temporale previsti per la velocità, il campo magnetico e la densità di corrente. Sebbene i tassi di convergenza per velocità e densità di corrente siano leggermente ridotti in alcuni casi (attribuiti alla regolarità limitata della soluzione dovuta al termine di Hall), il campo magnetico mantiene una convergenza ottimale.
• Preservazione della Divergenza: L'errore massimo di divergenza ($\max | \text{div } B_h |$) rimane al livello della tolleranza del risolutore per tutta la durata delle simulazioni, confermando la preservazione esatta del vincolo solenoidale.
• Stabilità Energetica: L'energia totale del sistema decade monotonicamente nel tempo, in coerenza con la legge teorica di dissipazione dell'energia.
• Fedeltà Fisica: Le simulazioni del flusso di Arnold–Beltrami–Childress (ABC), del vortice di Orszag–Tang e della riconnessione guidata da Hall in una lastra di Harris catturano con successo fenomeni fisici complessi, inclusa la riconnessione magnetica, la formazione di isole magnetiche, funi di flusso e i caratteristici pattern di corrente quadrupolari associati alla dinamica Hall–MHD. Lo schema risolve queste caratteristiche accuratamente anche su mesh relativamente grezze.

Significato
Il documento afferma che questo lavoro fornisce un framework matematicamente rigoroso e computazionalmente robusto per la simulazione di sistemi Hall–MHD. Combinando la regolarizzazione di Voigt con una formulazione agli elementi finiti misti che preserva la divergenza, gli autori superano le doppie sfide della non linearità del termine di Hall e del rigoroso requisito di divergenza nulla. Il metodo offre un'approssimazione fisicamente coerente che è stabile, convergente e capace di gestire geometrie complesse, colmando un vuoto lasciato da metodi precedenti privi di un'analisi rigorosa dell'errore o di una preservazione esatta della divergenza. I risultati suggeriscono che questo approccio è uno strumento valido per lo studio di fenomeni di plasma astrofisici e di laboratorio in cui gli effetti di Hall sono dominanti.