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

Titel: Afleiding van horizontale snelheidsvelden uit de inductievergelijking in de zonne-atmosfeer: I. Analytische en numerieke oplossingen in 2D

Auteurs: H. Vila Crespo et al.
Publicatie: Astronomy & Astrophysics (2026)

---

1. Probleemstelling

Traditionele spectroscopische en spectropolarimetrische waarnemingen van de zon maken gebruik van het Dopplereffect om uitsluitend de snelheidscomponent langs de zichtlijn ($v_z$) te meten. Veel dynamische processen in de zonne-atmosfeer (zoals subgranulaire stromingen, wervelingen en magnetische tornado's) worden echter gedomineerd door sterke horizontale bewegingen ($v_x, v_y$).

Bestaande methoden om deze horizontale snelheden te schatten, zoals Local Correlation Tracking (LCT), zijn over het algemeen ongevoelig voor variaties in hoogte ($z$). Ze leveren snelheidsvelden op die beperkt zijn tot één enkel vlak loodrecht op de zichtlijn ($v_x(x,y), v_y(x,y)$). Het doel van dit onderzoek is om de afhankelijkheid van de hoogte voor het horizontale snelheidsveld te kunnen afleiden, dus $v_h(x, y, z)$, door gebruik te maken van de volledige inductievergelijking.

2. Methodologie

De auteurs ontwikkelen een nieuwe numerieke methode om de horizontale snelheidscomponenten af te leiden uit de ideale magnetohydrodynamische (MHD) inductievergelijking:

$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B})$$

Aannames en Input:
De methode gaat ervan uit dat de volgende grootheden bekend zijn (via moderne Stokes-inversietechnieken toegepast op spectropolarimetrische data):

• Het magnetische veld $\mathbf{B}(x, y, z, t)$.
• De tijdsafgeleide van het magnetische veld $\dot{\mathbf{B}}$.
• De snelheidscomponent langs de zichtlijn $v_z(x, y, z, t)$.

Numerieke Implementatie (2D):
In dit eerste deel van het onderzoek wordt de methode getest in een tweedimensionaal domein $(y, z)$, waarbij $z$ de richting van de zichtlijn is en $v_x = B_x = 0$.

1. Discretisatie: De inductievergelijking wordt opgelost met behulp van het eindige-differentiemethode (finite differences) met een 5-punts centraal stencil voor de ruimtelijke afgeleiden.
2. Overbepaald Systeem: Voor elk punt in het rooster leveren de twee niet-triviale componenten van de inductievergelijking ($\dot{B}_y$ en $\dot{B}_z$) twee vergelijkingen op voor één onbekende: de horizontale snelheid $v_y(y, z)$.
3. Lineair Systeem: Dit leidt tot een overbepaald lineair stelsel van vergelijkingen ($\mathbf{A} \cdot \mathbf{x} = \mathbf{c}$), waarbij $\mathbf{x}$ de vector is van de onbekende snelheden.
4. Oplossing: Het systeem wordt opgelost in de zin van de kleinste kwadraten (least squares) via LU-decompositie (LAPACK).
5. Randvoorwaarden: Er worden Neumann-randvoorwaarden (symmetrisch, afgeleide nul) toegepast aan de randen van het domein, wat geschikter is voor praktische toepassingen waar de exacte snelheid aan de rand niet bekend is.

3. Key Contributions (Belangrijkste Bijdragen)

• Nieuwe Methodiek: De introductie van een methode die de volledige inductievergelijking gebruikt om de hoogte-afhankelijkheid van horizontale snelheden te reconstrueren, in plaats van ze te beperken tot één vlak.
• 2D Validatie: Een robuuste validatie van de methode in een 2D-domein, wat de basis legt voor toekomstige 3D-reconstructies.
• Analytische en Numerieke Tests: De methode wordt getest op zowel exact opgeloste analytische functies als op realistische data uit een CO5BOLD MHD-simulatie van zonne-magnetoconectie.

4. Resultaten

De auteurs testen hun methode in twee scenario's:

A. Analytisch geval:

• Gebruikmakend van exact bekende magnetische en snelheidsvelden.
• Resultaat: De horizontale snelheid $v_y(y, z)$ wordt met een gemiddelde fout van ongeveer 1% en een mediaan van 1,5% gereconstrueerd.
• De hoek tussen het oorspronkelijke en het afgeleide snelheidsveld is gemiddeld slechts 0,5°.
• De correlatiecoëfficiënt (Pearson) is bijna 1,0.

B. CO5BOLD MHD-simulatie:

• Toepassing op een realistische simulatie van de zonneoppervlakte.
• Resultaat: De methode levert een visueel identiek resultaat op ten opzichte van de originele simulatie.
• De fouten zijn zelfs lager dan in het analytische geval (gemiddelde fout ~1%, mediaan ~0,11%), mede omdat de snelheid in de simulatie nergens exact nul wordt (wat de hoekfouten beperkt).
• Beperkingen: De fouten nemen toe in gebieden waar de horizontale snelheid $v_y$ dicht bij nul ligt (waar de snelheid voornamelijk langs de $z$-as is gericht). In deze gebieden wordt de richting en grootte minder nauwkeurig bepaald.

5. Betekenis en Toekomstperspectief

• Fundamentele Vooruitgang: Deze studie bewijst dat het mogelijk is om de verticale structuur van horizontale stromingen in de zonne-atmosfeer te reconstrueren, wat essentieel is voor het begrijpen van energietransport (magnetische energie en helicaliteit) van de fotosfeer naar de corona.
• Toekomstige Uitbreiding: De methode is direct uitbreidbaar naar 3D door de derde component van de inductievergelijking toe te voegen. Dit zal echter leiden tot een aanzienlijk groter lineair stelsel (van $2NM \times NM$ naar $3NML \times 2NML$), wat meer rekenkracht vereist.
• Realistische Scenarios: De auteurs benadrukken dat toekomstig werk zich moet richten op de invloed van meetfouten in de input-data (uit Stokes-inversie), eindige tijdsbemonstering voor $\dot{\mathbf{B}}$, en magnetische diffusie, aangezien deze factoren de ideale foutmarges van 1% waarschijnlijk zullen overschaduwen.

Conclusie: De voorgestelde methode slaagt erin om het horizontale snelheidsveld in het $(y, z)$-vlak succesvol te reconstrueren met hoge nauwkeurigheid, en legt zo de fundering voor toekomstige 3D-reconstructies van zonne-stromingen.