Decoupled iterative schemes for solving stationary MHD problems

Sintesi Tecnica: Schemi Iterativi Disaccoppiati per Problemi MHD Stazionari

Enunciato del Problema
Il documento affronta la soluzione numerica del sistema magnetoidrodinamico (MHD) stazionario e incomprimibile. Questo sistema modella l'interazione tra fluidi elettricamente conduttori e campi elettromagnetici, governata dall'accoppiamento tra le equazioni di Navier-Stokes e le equazioni di Maxwell. La formulazione matematica consiste nel trovare il campo di velocità ($u$), il campo magnetico ($B$), la pressione modificata ($p$) e un moltiplicatore di Lagrange ($\lambda$) che soddisfino un insieme di equazioni differenziali alle derivate parziali accoppiate. La sfida principale risiede nella forte non linearità e nel serrato accoppiamento tra le variabili del fluido e quelle magnetiche, il che complica la costruzione di solver numerici efficienti e stabili.

Metodologia
Gli autori propongono e analizzano due nuovi schemi iterativi disaccoppiati ispirati alla strategia di Incremental Viscosity Splitting (IVS), originariamente sviluppata per i flussi di Navier-Stokes non stazionari.

1. Schema IVS Standard (Algoritmo 5.1):

   • Strategia di Splitting: Il metodo scinde le equazioni di quantità di moto e di induzione rispetto ai loro termini diffusivi.
   • Passaggio 1: Si risolvono le variabili intermedie ($\tilde{u}, \tilde{B}$) trattando i termini convettivi e di accoppiamento esplicitamente utilizzando i valori dell'iterazione precedente, mantenendo la pressione e il moltiplicatore di Lagrange del passo precedente. Ciò produce un sistema lineare accoppiato per i campi intermedi.
   • Passaggi 2 e 3: Si correggono i campi intermedi per ottenere la velocità finale ($u$) e il campo magnetico ($B$) risolvendo due sistemi separati e disaccoppiati di tipo Stokes. Questi sistemi impongono i vincoli di divergenza nulla ($\nabla \cdot u = 0, \nabla \cdot B = 0$) e aggiornano la pressione e il moltiplicatore di Lagrange.
   • Vantaggio dell'Algebra Lineare: Una caratteristica chiave è che le matrici del complemento di Schur per le correzioni di Stokes sono indipendenti dal livello di iterazione ($k$). Di conseguenza, devono essere assemblate e precondizionate una sola volta, riducendo il carico computazionale rispetto ai metodi in cui la matrice del sistema cambia ad ogni passaggio.

2. Schema Elsässer IVS (Algoritmo 6.1):

   • Riformulazione: Per disaccoppiare ulteriormente il sistema, gli autori riformulano le equazioni MHD utilizzando le variabili di Elsässer ($w = u + \sqrt{\kappa}B$ e $z = u - \sqrt{\kappa}B$).
   • Disaccoppiamento: In questa formulazione, l'accoppiamento tra le nuove variabili avviene attraverso un termine non lineare e uno lineare. L'algoritmo proposto risolve $w$ e $z$ in passi separati, sostituendo il grande sistema accoppiato dello standard IVS con due sistemi lineari più piccoli e completamente disaccoppiati.
   • Applicabilità: Questo approccio è teoricamente vantaggioso sotto condizioni di Dirichlet omogenee, sebbene gli autori notino che emergono sfide con condizioni al contorno miste.

Contributi Teorici
Il documento fornisce un'analisi matematica rigorosa per entrambi gli schemi:

• Ben posto (Well-posedness): Gli autori stabiliscono risultati di esistenza e unicità per il problema MHD continuo sotto un'ipotesi di "piccoli dati" (specificamente, un limite sui numeri di Reynolds e di Reynolds magnetico rispetto alla viscosità).
• Dimostrazioni di Convergenza:
  • Lo schema iterativo di Picard classico viene analizzato e dimostrato convergente sotto la stessa ipotesi di piccoli dati.
  • Entrambi gli schemi (Standard IVS ed Elsässer IVS) sono dimostrati essere uniformemente limitati e geometricamente convergenti alla soluzione unica.
  • Le dimostrazioni utilizzano stime di energia, la condizione inf-sup (LBB) e lemmi specifici riguardanti la convergenza di sequenze non negative.
• Analisi dell'Errore: Vengono derivati limiti teorici che mostrano come gli errori in velocità, campo magnetico, pressione e moltiplicatori di Lagrange decadano a zero all'aumentare del numero di iterazioni.

Risultati Numerici
Gli autori validano i loro risultati teorici utilizzando il software FreeFem++ con coppie di elementi finiti $P_2-P_1$.

• Test di Convergenza: Soluzioni prodotte (manufactured solutions) in 2D dimostrano tassi di convergenza ottimali ($O(h^2)$ in $L^2$ e $O(h)$ in $H^1$) sia per la velocità che per il campo magnetico, confermando l'accuratezza degli schemi.
• Problemi di Benchmark:
  • Flusso di Hartmann: Sia i problemi di flusso di Hartmann 2D che 3D sono stati risolti. I risultati numerici hanno corrisponduto alle soluzioni analitiche con alta precisione per vari numeri di Hartmann ($Ha$).
  • Cavità a parete mossa (Lid-Driven Cavity): Gli schemi hanno simulato con successo pattern di flusso complessi per vari numeri di Reynolds ($Re$), Reynolds magnetico ($Rm$) e accoppiamento ($\kappa$), inclusi casi di alto accoppiamento ($\kappa=5000$).
  • Flusso su uno scalino: Il metodo ha catturato accuratamente la separazione del flusso e i cali di pressione, dimostrando robustezza nel gestire singolarità geometriche e diverse intensità di accoppiamento.

Significatività e Rivendicazioni
Il documento afferma che i processi iterativi disaccoppiati proposti offrono un'alternativa pratica ed efficiente ai classici metodi iterativi monolitici o completamente accoppiati (come i metodi di Newton o Picard standard).

• Efficienza: Dividendo i termini diffusivi, il metodo riduce la complessità dei sistemi lineari risolti ad ogni passaggio. Nello specifico, la formulazione di Elsässer permette la soluzione di due sistemi più piccoli e disaccoppiati invece di un unico grande sistema accoppiato.
• Riutilizzabilità: L'approccio IVS garantisce che le matrici del complemento di Schur per le correzioni di Stokes siano statiche, consentendo un singolo passaggio di assemblaggio e precondizionamento, il che è computazionalmente vantaggioso per i problemi 3D.
• Robustezza: Gli schemi mantengono la proprietà di divergenza nulla dei campi di velocità e magnetico a livello discreto e rimangono stabili sotto le condizioni di piccoli dati stabilite.

Gli autori concludono che, sebbene il presente lavoro si concentri su problemi stazionari, la metodologia fornisce una base per futuri estensioni a problemi MHD tempo-dipendenti (non stazionari), con un focus sull'analisi della stabilità temporale e spaziale.