Inference of horizontal velocity fields from the induction equation in the solar atmosphere. I. Analytical and numerical solutions in 2D

Titolo

Inferenza dei campi di velocità orizzontale dall'equazione di induzione nell'atmosfera solare. I. Soluzioni analitiche e numeriche in 2D.

1. Il Problema

Le osservazioni spettroscopiche e spettropolarimetriche del Sole, basate sull'effetto Doppler, permettono di misurare direttamente solo la componente della velocità del plasma lungo la linea di vista (indicata come $v_z$). Tuttavia, molti processi dinamici nell'atmosfera solare (come flussi su scala sub-granulare, moti vorticosi, riconnessione magnetica e tornado magnetici) coinvolgono movimenti significativi nel piano perpendicolare alla linea di vista ($v_x, v_y$).

Le tecniche esistenti per stimare queste velocità orizzontali, come il Local Correlation Tracking (LCT), sono generalmente insensibili alle variazioni di altezza ($z$), fornendo campi di velocità bidimensionali su un singolo piano ($v_x(x,y), v_y(x,y)$). L'obiettivo di questo studio è superare questa limitazione, inferendo la dipendenza dall'altezza del campo di velocità orizzontale ($v_h(x,y,z)$), sfruttando le moderne tecniche di inversione Stokes che permettono già di recuperare la dipendenza da $z$ per il campo magnetico ($B$) e la velocità longitudinale ($v_z$).

2. Metodologia

Gli autori propongono un nuovo metodo basato sull'equazione di induzione magnetica nell'ambito della Magnetoidrodinamica (MHD) ideale (diffusione magnetica nulla):
$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B})$$

Ipotesi e Assunzioni:

• Si assume che il vettore campo magnetico $\mathbf{B}(x,y,z,t)$, la sua derivata temporale $\dot{\mathbf{B}}$ e la componente di velocità lungo la linea di vista $v_z$ siano noti (ottenibili tramite inversione Stokes).
• L'obiettivo è risolvere per le componenti orizzontali della velocità, in particolare $v_y$ nel caso bidimensionale studiato.

Implementazione Numerica (Piano 2D $y, z$):

1. Riduzione 2D: Il problema è inizialmente testato nel piano $(y, z)$, assumendo $v_x = B_x = 0$. In questo piano, le due componenti non banali dell'equazione di induzione forniscono due equazioni per l'unica incognita $v_y(y, z)$.
2. Discretizzazione: Le equazioni sono discretizzate utilizzando il metodo delle differenze finite con uno stencil centrato a 5 punti.
3. Sistema Sovradeterminato: Per ogni punto della griglia $(i, j)$, si ottengono due equazioni lineari per una singola incognita ($v_y$). Questo crea un sistema lineare sovradeterminato della forma $\mathbf{A} \cdot \mathbf{x} = \mathbf{c}$, dove $\mathbf{x}$ è il vettore delle velocità incognite.
4. Risoluzione: Il sistema viene risolto in senso dei minimi quadrati ($\mathbf{x} = (\mathbf{A}^T\mathbf{A})^{-1}\mathbf{A}^T\mathbf{c}$) utilizzando la libreria LAPACK.
5. Condizioni al Contorno: Vengono implementate condizioni al contorno di Neumann (derivata nulla) o Dirichlet (valore fissato) tramite l'uso di "celle fantasma" (ghost cells) per mantenere lo schema di derivata valido ai bordi del dominio.

3. Contributi Chiave

• Estensione della dimensionalità: Per la prima volta, un metodo basato sull'equazione di induzione viene esteso per recuperare la dipendenza verticale ($z$) delle velocità orizzontali, non limitandosi a un singolo piano.
• Sfruttamento completo dell'equazione di induzione: A differenza di metodi precedenti che usavano solo la componente verticale dell'equazione di induzione, questo approccio utilizza tutte le componenti disponibili nel piano 2D per vincolare la soluzione.
• Validazione rigorosa: Il metodo è stato testato sia su soluzioni analitiche prescritte che su simulazioni MHD realistiche (codice CO5BOLD), dimostrando la sua robustezza in scenari ideali.

4. Risultati

Gli autori hanno valutato le prestazioni del metodo confrontando la velocità orizzontale recuperata ($\tilde{v}_y$) con quella reale ($v_y$).

• Caso Analitico:
  • Errore medio sulla magnitudine di $v_y$: ~0.82%.
  • Errore mediano: ~1.41%.
  • Errore angolare medio tra il vettore velocità reale e quello recuperato: ~0.6°.
  • Il coefficiente di correlazione di Pearson è quasi 1.0.
• Simulazione MHD (CO5BOLD):
  • Il metodo ha recuperato con successo il campo di velocità da una simulazione complessa di convezione magnetosolare.
  • L'errore medio è risultato leggermente inferiore rispetto al caso analitico (~1.06%), con una mediana di 0.11%.
  • L'errore angolare medio è stato di 0.1°.

Limitazioni e Osservazioni:

• L'accuratezza diminuisce nelle regioni dove la componente orizzontale della velocità ($v_y$) o il modulo totale della velocità tendono a zero. In queste zone, l'errore relativo e l'errore angolare aumentano significativamente (fino al 50% e 5° nel caso analitico), poiché la direzione del vettore diventa indefinita quando la sua magnitudine è nulla.
• Tuttavia, la correlazione lineare rimane estremamente alta, indicando che il metodo cattura correttamente la struttura del campo di velocità.

5. Significato e Prospettive Future

Questo lavoro stabilisce le fondamenta teoriche e numeriche per la ricostruzione tridimensionale dei campi di velocità orizzontali nell'atmosfera solare.

• Futuro 3D: Il metodo è progettato per essere esteso al caso tridimensionale completo $(x, y, z)$. Sebbene ciò comporti un aumento significativo delle dimensioni della matrice da invertire (da $2NM \times NM$ a $3NML \times 2NML$) e un maggiore costo computazionale, la procedura di discretizzazione rimane analoga.
• Impatto Scientifico: La capacità di conoscere $v_h(x,y,z)$ è cruciale per calcolare con precisione i flussi di energia magnetica ed elicità dalla fotosfera alla corona, nonché per migliorare l'inferenza delle correnti elettriche solari.
• Sfide Future: I prossimi studi dovranno affrontare l'impatto di errori sistematici e casuali nei dati di input (derivanti dall'inversione Stokes), la risoluzione temporale finita per il calcolo di $\dot{B}$ e gli effetti della diffusione magnetica non ideale.

In sintesi, il paper dimostra che, in condizioni ideali, è possibile recuperare con alta precisione (errore medio ~1%) la struttura verticale delle velocità orizzontali solari, aprendo la strada a una comprensione più completa della dinamica magnetosolare 3D.