Improving divergence cleaning in cosmological SPMHD simulations

Sintesi Tecnica: Miglioramento della Pulizia della Divergenza nelle Simulazioni Cosmologiche SPMHD

Problematica
Nelle simulazioni cosmologiche di Magnetoidrodinamica a Particelle Smootizzate (SPMHD), il mantenimento della condizione solenoidale ($\nabla \cdot \mathbf{B} = 0$) rappresenta una sfida numerica persistente. Gli schemi standard in codici come Gadget si affidano tipicamente all'approccio "8-wave" di Powell et al. (1999). Sebbene questo metodo avviti gli errori di divergenza con il flusso del fluido, non li rimuove attivamente. Nei codici basati su SPH, gli errori residui di $\nabla \cdot \mathbf{B}$ possono raggiungere l'ordine unitario, portando a un'amplificazione non fisica del campo magnetico e alla "instabilità tensile", in cui termini sorgente spurii nell'equazione del momento dominano quando la pressione magnetica supera la pressione del gas. Sebbene lo schema di pulizia iperbolico/parabolico proposto da Dedner et al. (2002) offra una soluzione, la sua formulazione originale si è rivelata instabile in SPMHD, in particolare in ambienti dinamici con timestep individuali e forti salti di densità, agendo spesso come una sorgente non fisica di energia magnetica.

Metodologia
Gli autori implementano una versione vincolata e conservativa dell'energia dell'algoritmo di pulizia di tipo Dedner (Tricco & Price 2012; Tricco et al. 2016a) nel codice SPMHD cosmologico OpenGadget3. Le principali adattamenti metodologici includono:

• Conservazione dell'Energia: Le equazioni governanti sono riformulate per garantire la conservazione dell'energia totale. Ciò comporta l'evoluzione della variabile di campo di pulizia come $\psi/c_h$ (dove $\psi$ è la variabile scalare ausiliaria e $c_h$ è la velocità di pulizia) anziché solo $\psi$, e l'inclusione di un termine sorgente specifico $-\frac{1}{2}\psi(\nabla \cdot \mathbf{v})$ nell'equazione di evoluzione. Ciò assicura che lo smorzamento parabolico del campo di pulizia rimuova l'energia magnetica senza introdurre una crescita illimitata.
• Adattamento Cosmologico: Le equazioni sono adattate per coordinate comoventi e timestep individuali delle particelle. La velocità di pulizia $c_h$ è impostata sulla velocità locale magnetosonica veloce ($\sqrt{v_A^2 + c_s^2}$), garantendo che la pulizia sia intrinsecamente stabile rispetto al criterio di Courant entro i vincoli standard del timestep SPH.
• Configurazione della Simulazione: Gli autori eseguono una serie di simulazioni su un ammasso di galassie massiccio ($M_{200c} \approx 10^{15} M_\odot$) seminato con un campo magnetico uniforme ($10^{-14}$ G comoving). Vengono confrontati cinque casi distinti:
  1. NIFTY-noclean: Run di riferimento che utilizza solo lo schema di advezione di Powell.
  2. NIFTY-clean: Run fiduciale con il pieno schema di pulizia vincolato (inclusa la retroazione dell'entropia).
  3. NIFTY-clean2: Identico al precedente ma con un fattore di sovra-pulizia $\zeta=2$ (pulizia a due volte la velocità d'onda locale).
  4. NIFTY-advect: Schema di pulizia con il termine sorgente della divergenza di velocità disattivato.
  5. NIFTY-heat: Schema di pulizia in cui l'energia di pulizia dissipata viene scartata invece di essere reimmessa nell'equazione dell'entropia.

Risultati Chiave

• Riduzione dell'Errore di Divergenza: Il risultato primario è una riduzione degli errori di divergenza relativa ($|\nabla \cdot \mathbf{B}|/B$) di 2–3 ordini di grandezza in tutto il volume dell'ammasso rispetto allo schema solo Powell. Questa riduzione è uniforme tra il nucleo dell'ammasso e le periferie.
• Amplificazione del Campo Magnetico: Nelle periferie dell'ammasso ($r \approx 1\text{--}3 \, h^{-1}$ Mpc), dove il campo magnetico è più debole e la risoluzione è minore, lo schema di pulizia attivo permette un'amplificazione del campo magnetico di un fattore di 5–10 rispetto alla run solo Powell. Gli autori attribuiscono ciò alla rimozione delle forze di Lorentz non fisiche causate dagli errori di divergenza, che precedentemente sopprimevano il dinamo a piccola scala in queste regioni.
• Stabilità del Nucleo: La struttura della densità globale e la geometria del campo magnetico nel nucleo dell'ammasso rimangono ampiamente invariate tra le run con e senza pulizia, indicando che la pulizia non altera artificialmente lo stato saturo del dinamo centrale.
• Consistenza Termodinamica: L'inclusione della retroazione dell'entropia (reimmissione dell'energia dissipata dalla pulizia sotto forma di calore) si dimostra critica per la consistenza termodinamica. Le run senza questa retroazione (NIFTY-heat) risultano in un'entropia centrale inferiore e in campi magnetici centrali saturi più elevati a causa della compressione adiabatica di un nucleo più freddo.
• Costo Computazionale: L'implementazione non introduce alcun sovraccarico misurabile nel tempo di esecuzione; la pulizia migliorata viene ottenuta essenzialmente "gratuitamente" entro i vincoli di timestep di Courant esistenti.

Significatività e Rivendicazioni
Il lavoro sostiene che un controllo accurato degli errori di divergenza è essenziale per catturare l'amplificazione del campo magnetico nelle regioni a bassa densità e scarsamente risolte delle periferie dei cluster galattici. Gli autori sostengono che gli errori di divergenza incontrollati negli schemi Powell standard agiscono come rumore numerico che sopprime il dinamo nelle regioni sub-sature. Implementando lo schema di pulizia vincolato e conservativo dell'energia, la simulazione raggiunge una realizzazione più fedele della fisica MHD senza penalità computazionali.

Gli autori concludono che lo schema di pulità vincolato:

1. Riduce drasticamente gli errori di divergenza numerica (2–3 ordi di grandezza).
2. Permette una rappresentazione più accurata dell'amplificazione del campo nelle regioni sotto-risolte.
3. Produce strutture del campo magnetico più lisce e coerenti.
4. Stabilizza l'evoluzione del campo magnetico attraverso diverse risoluzioni, migliorando la fedeltà fisica dei risultati anche a risoluzioni moderate.

Il lavoro raccomanda che le applicazioni SPMHD cosmologiche adottino questo schema di pulizia di default, utilizzando la velocità di pulizia naturale ($\zeta=1$) e mantenendo attiva la retroazione energetica parabolica per chiudere il bilancio energetico. L'implementazione è ora pubblicamente disponibile in OpenGadget3.