Radiation magnetohydrodynamics modeling of an impulsively driven chromospheric jet in the solar atmosphere

Sintesi Tecnica: Modellazione Magnetoidrodinamica Radiativa di un Getto Cromosferico Guidato Impulsivamente

Enunciato del Problema
I getti solari, in particolare le spicole e le macrospicole, sono espulsioni dinamiche di plasma osservate su più lunghezze d'onda. Sebbene la riconnessione magnetica e le perturbazioni impulsive (come impulsi di pressione o velocità) siano ampiamente accettati come meccanismi innescanti, il ruolo dei processi radiativi nell'evoluzione di questi getti rimane oggetto di indagine. Studi numerici precedenti hanno modellato getti guidati impulsivamente utilizzando magnetoidrodinamica (MHD) ideale o non ideale, spesso accoppiata a funzioni di raffreddamento radiativo empiriche. Tuttavia, esiste un vuoto nella simulazione di tali getti mediante equazioni MHD non ideali pienamente accoppiate a un sistema di equazioni di trasporto della radiazione mediato in frequenza e angolo. Nello specifico, mancano simulazioni che evolvano indipendentemente sia la densità di energia radiativa sia le componenti del flusso radiativo per determinarne l'impatto specifico sulla morfologia e sulla dissipazione del getto nella corona solare otticamente sottile.

Metodologia
Lo studio impiega una simulazione numerica bidimensionale utilizzando il codice PLUTO (v4.4-patch3) per modellare un getto cromosferico guidato impulsivamente. Il quadro fisico è costituito da:

• Equazioni Governative: Il sistema risolve equazioni MHD non ideali accoppiate a uno schema di idrodinamica radiativa a due momenti nell'approssimazione grigia. Ciò include l'evoluzione della densità del plasma, della quantità di moto, dell'energia totale e del campo magnetico, insieme alla densità di energia radiativa ($E_r$) e al flusso radiativo ($F_r$).
• Trasporto Radiativo: Viene utilizzata la chiusura M1 per chiudere le equazioni dei momenti radiativi, consentendo la descrizione sia dei limiti di diffusione sia di quelli di flusso libero. I termini di interazione radiazione-materia sono derivati tramite trasformazioni di Lorentz al primo ordine in $\beta = v/c$, tenendo conto delle opacità di assorbimento e di scattering.
• Modello dell'Atmosfera Solare: Lo stato iniziale è un equilibrio idrostatico basato sul modello semi-empirico C7, esteso alla corona. La composizione del plasma assume un peso molecolare medio di $\mu = 0.6$ (idrogeno ed elio completamente ionizzati).
• Parametri Fisici: La simulazione utilizza una bassa opacità di assorbimento media di Planck ($\kappa = 1 \times 10^{-6}$ cm$^2$ g$^{-1}$) e un'opacità di scattering media di Rosseland ($\sigma = 0.34$ cm$^2$ g$^{-1}$), caratteristiche della corona solare otticamente sottile. La conduzione termica è modellata con un approccio a flusso limitato, e una funzione empirica di riscaldamento coronale mantiene l'atmosfera di fondo.
• Perturbazione: La formazione del getto è innescata da un impulso gaussiano localizzato nella componente verticale della velocità ($v_y$) a $y = 1.75$ Mm con un'ampiezza di 100 km s$^{-1}$.
• Schema Numerico: Le equazioni sono risolte utilizzando un integratore temporale Runge-Kutta TVD del secondo ordine. Le equazioni MHD utilizzano il risolutore Riemann HLLD, mentre il trasporto radiativo utilizza il risolutore HLLC. Schemi Runge-Kutta Impliciti-Espliciti (IMEX) gestiscono i termini rigidi di interazione radiazione-materia, e viene impiegata l'approssimazione della velocità ridotta della luce (RSLA) con $\hat{c} = 0.1c$ per garantire la fattibilità computazionale.

Risultati Chiave
La simulazione riproduce con successo la formazione e l'evoluzione di una struttura di plasma collimata simile a una macrospicola:

• Morfologia e Dinamica del Getto: L'impulso impulsivo genera un'onda d'urto che propelle il plasma verso l'alto, formando un getto che raggiunge un'altezza di circa 8 Mm. Il getto mostra una velocità verticale massima di $\sim$60 km s$^{-1}$ al fronte d'urto e di $\sim$30 km s$^{-1}$ all'interno del corpo del getto. La struttura mostra una morfologia a doppio filo nelle fasi successive, una caratteristica comunemente associata alle spicole solari.
• Evoluzione Termica: La testa del getto raggiunge temperature di $\sim$2 $\times$ 10$^5$ K, mentre la massa del getto si stabilizza intorno a 2 $\times$ 10$^5$ K. La regione d'urto inizialmente raggiunge $\sim$9 $\times$ 10$^5$ K prima di raffreddarsi.
• Effetti Radiativi: L'inclusione di un trasporto radiativo auto-consistente rivela che la radiazione agisce come un meccanismo di raffreddamento efficiente nella corona otticamente sottile. La densità di energia radiativa aumenta linearmente, mentre il flusso radiativo diminuisce, coerentemente con le perdite radiative otticamente sottili.
• Meccanismo di Dissipazione: Crucialmente, i termini di trasporto radiativo contribuiscono alla dissipazione del getto. L'accoppiamento tra il plasma e il campo radiativo induce velocità di flusso negative (moto discendente) nella corona, causando la diffusione e la dissipazione del getto di nuovo nell'atmosfera inferiore. Le variabili radiative non dettano la forma o la morfologia iniziale del getto, ma ne governano significativamente l'evoluzione e la durata.
• Stabilità di Fondo: L'atmosfera solare di fondo rimane globalmente stabile rispetto ai termini sorgente (conduzione termica, riscaldamento e radiazione) in assenza della perturbazione del getto, con osservate solo lievi variazioni di temperatura.

Significato e Affermazioni
L'articolo afferma che questo lavoro rappresenta la prima simulazione numerica di un getto cromosferico guidato impulsivamente in cui la densità di energia radiativa e il flusso sono evoluti in modo auto-consistente attraverso equazioni di trasporto a due momenti, piuttosto che affidarsi a funzioni di raffreddamento empiriche. Il significato principale risiede nel dimostrare che:

1. Il Raffreddamento Radiativo è Critico: Anche nel regime otticamente sottile, dove le opacità sono basse, il trattamento esplicito del trasporto radiativo agisce come un meccanismo di raffreddamento dominante che limita l'altezza, la durata e il contenuto energetico di eventi simili a macrospicole.
2. Dinamica della Dissipazione: I termini radiativi inducono flussi discendenti responsabili della rapida attenuazione e dissipazione del getto, un processo coerente con studi precedenti che utilizzavano funzioni di perdita empiriche ma ora derivato dai primi principi all'interno del quadro del trasporto radiativo.
3. Avanzamento Metodologico: Lo studio convalida l'uso della chiusura M1 accoppiata alla MHD non ideale per la modellazione dei getti solari, mostrando che, sebbene la radiazione non alteri la morfologia iniziale del getto, è essenziale per catturare la corretta evoluzione termodinamica e le scale temporali di dissipazione.

Gli autori concludono che, sebbene il getto simulato riproduca le caratteristiche dinamiche chiave delle macrospicole osservate (come la formazione di onde d'urto e le strutture a doppio filo), i suoi parametri specifici (altezza, velocità, temperatura) rientrano nei limiti inferiori delle gamme osservate, suggerendo che l'evento modellato sia un evento a relativamente bassa energia. Il lavoro sottolinea che un trattamento auto-consistente del trasporto radiativo è necessario per modellare accuratamente la dissipazione di tali fenomeni atmosferici solari.