An Adjoint Projection Formulation for Enforcing the divergence-free Constraint in Smoothed Particle Magnetohydrodynamics

Sintesi Tecnica: Una Formulazione di Proiezione Adjoint per l'Imposizione del Vincolo $\nabla \cdot \mathbf{B} = 0$ nella Magnetoidrodinamica a Particelle Smooth (SPMHD)

Enunciato del Problema
Nella magnetoidrodinamica (MHD), il vincolo solenoidale $\nabla \cdot \mathbf{B} = 0$ è un requisito fisico fondamentale. Sebbene i codici basati su griglia utilizzino spesso il trasporto vincolato (CT) per preservare questo vincolo in modo discreto, i metodi senza griglia come la Magnetoidrodinamica a Particelle Smooth (SPMHD) mancano di una naturale topologia faccia-spigolo per applicare direttamente il CT. Di conseguenza, l'approccio prevalente in SPMHD è stato la "pulizia della divergenza" (divergence cleaning), che evolve un campo scalare ausiliario per propagare e smorzare gli errori di divergenza. Tuttavia, i metodi di pulizia introducono parametri fisici (velocità di propagazione e scala temporale di smorzamento) che possono influenzare i risultati della simulazione se non scelti in modo coerente, particolarmente in sistemi con ampi intervalli dinamici. Precedenti tentativi di proiezione in SPMHD, come quelli di Price e Monaghan, hanno utilizzato somme di tipo gravitazionale per le correzioni, ma soffrivano di inconsistenze tra l'operatore di divergenza usato per misurare gli errori e l'operatore di gradiente usato per calcolare le correzioni, portando a scarse prestazioni sul rumore a breve lunghezza d'onda.

Metodologia
Il lavoro propone un metodo di proiezione che impone il vincolo discreto di divergenza nulla risolvendo un problema ellittico dopo ogni aggiornamento MHD. L'innovazione principale è la costruzione di un operatore di correzione che è strettamente l'adianto (adjoint) dell'operatore di divergenza discreto usato per misurare l'errore, definito all'interno di una metrica pesata sul volume.

1. Formulazione: Il metodo cerca un campo magnetico corretto $\mathbf{B}^{n+1}$ che minimizzi la variaione della norma dell'energia magnetica discreta (definita dai volumi delle particelle $V_i = m_i/\rho_i$) soggetto al vincolo $D\mathbf{B}^{n+1} = 0$, dove $D$ è l'operatore di divergenza SPH discreto.
2. Gradiente Adjoint: Introducendo un campo di moltiplicatori di Lagrange $\pi$, la correzione viene derivata come $\mathbf{B}^{n+1} = \mathbf{B}^* - \mathcal{G}\pi$. L'operatore $\mathcal{G}$ è esplicitamente costruito come l'adianto della metrica di volume di $D$ ($\mathcal{G} = M_V^{-1} D^T$). Ciò assicura che lo spazio di correzione sia ortogonale allo spazio privo di divergenza rispetto al prodotto interno di volume.
3. Discretizzazione: Il documento deriva formule SPH esplicite sia per l'operatore di divergenza $D$ che per il suo gradiente adjanto $\mathcal{G}$. Fondamentalmente, il gradiente adjanto tiene conto dell'asimmetria nei coefficienti delle derivate SPH (derivante dalle variazioni di lunghezza di smoothing e densità) utilizzando la specifica simmetria della lista dei vicini e la ponderazione del volume.
4. Soluzione Numerica: L'equazione di proiezione risultante $L\pi = r$ (dove $L = D\mathcal{G}$) forma un sistema lineare simmetrico semidefinito positivo. Ciò consente di risolvere il sistema efficientemente utilizzando il metodo del Gradiente Coniugato Precondizionato (PCG) senza dover assemblare esplicitamente la matrice globale. L'implementazione distingue tra particelle attive (dove si risolve $\pi$) e particelle di confine fisse.

Contributi Chiave

• Formulazione Adjoint Coerente: Il lavoro fornisce la prima formulazione di proiezione SPMHD in cui il gradiente di correzione è rigorosamente costruito come l'adianto dell'operatore di divergenza discreto scelto sotto una metrica pesata sul volume. Ciò risolve le inconsistenze a livello di discretizzazione presenti in precedenti tentativi di proiezione.
• Proprietà di Minimizzazione dell'Energia: Viene dimostrato che la metrica di volume è la metrica naturale per l'energia magnetica. Di conseguenza, la proiezione produce una correzione che minimizza l'energia e non aumenta l'energia magnetica discreta, una proprietà critica per la fedeltà fisica.
• Solver Efficiente: La formulazione porta a un sistema simmetrico semidefinito positivo risolvibile tramite PCG, evitando il costo $O(N^2)$ dei metodi di sommatoria diretta usati nelle precedenti approcci.
• Criteri di Arresto Pratici: Gli autori introducono criteri di arresto pratici basati sull'accumulo dell'errore di divergenza tra i passi di proiezione, piuttosto che richiedere la convergenza alla precisione di macchina per ogni passo, bilanciando accuratezza e costo computazionale.

Risultati
Il metodo è stato testato utilizzando test di divergenza di tipo Dedner in due dimensioni e simulazioni di collasso gravitazionale magnetizzato in tre dimensioni.

• Test di Divergenza: In test 2D con sia reticoli perturbati che distribuzioni di particelle casuali, il metodo ha ridotto con successo il residuo di divergenza discreta al livello di errore di arrotondamento in virgola mobile ($\sim 10^{-15}$) con un numero sufficiente di iterazioni PCG. Notevolmente, ha corretto gli errori a lunghezze d'onda comparabili alla lunghezza di smoothing ($h$), con cui i metodi precedenti hanno faticato.
• Simulazioni di Collasso: Nelle simulazioni di collasso 3D, il metodo di proiezione ha superato significativamente la standard pulizia della divergenza iperbolica.
  • Riduzione dell'Errore: Con un criterio di arresto pratico (riducendo l'errore accumulato a una frazione dell'errore di aggiornamento MHD), il metodo di proiezione ha soppresso gli errori di divergenza normalizzati ben al di sotto di quelli ottenuti con la pulizia della divergenza.
  • Robustezza: Variando l'intervallo di proiezione ($N_{interval} = 1, 10, 100$), sono state ottenute strutture di densità e plasma-$\beta$ coerenti. Al contrario, la corsa con pulizia della divergenza ha mostrato differenze quantitative nella struttura di outflow e nel plasma $\beta$.
  • Costo: Il costo computazionale della proiezione è stato modesto, variando approssimativamente dall'1% al 10% del costo dell'aggiornamento SPMHD per intervalli pratici ($N_{interval} \approx 10$). Anche in test di stress estremi, lo schema può portare gli errori ai livelli di arrotondamento.

Significatività e Rivendicazioni
Il lavoro sostiene che questa formulazione di proiezione adjanta è un'alternativa robusta e attraente alla pulizia della divergenza per il controllo degli errori di $\nabla \cdot \mathbf{B}$ in SPMHD e in schemi MHD correlati a particelle o a griglia mobile.

I principali vantaggi evidenziati includono:

1. Indipendenza dai Parametri: A differenza della pulizia della divergenza, che dipende da parametri fisici (velocità dell'onda, scala di smorzamento) che possono alterare i risultati astrofisici, i parametri del metodo di proiezione (criterio di arresto, intervallo) sono adimensionali e si è dimostrato che hanno un impatto minimo sulle soluzioni entro gli intervalli testati.
2. Coerenza Fisica: Il metodo garantisce che la correzione non inietti energia magnetica artificiale.
3. Scalabilità: L'approccio non è limitato a SPMHD; gli autori suggeriscono che l'idea della proiezione basata su operatore sia applicabile ad altri schemi a particelle o a griglia mobile (es. GIZMO, AREPO) con le opportune modifiche.

Gli autori concludono che, sebbene il metodo richieda la scelta di parametri liberi per il termine e la frequenza, la sua capacità di produrre risultati coerenti e ad alta fedeltà senza i problemi di sensibilità dei metodi di pulizia lo rende una scelta superiore per il controllo della divergenza numerica nelle particelle MHD.