Tearing Instability in Gyrotropic MHD: Effects of Equilibrium Pressure Anisotropy

Technische Samenvatting: Tearing-instabiliteit in Gyrotrope MHD: Effecten van Evenwichtsdruk-anisotropie

Probleemstelling
De tearing-instabiliteit is een fundamenteel mechanisme voor magnetische reconnectie in dunne stroomsheets. Klassieke resistente Magnetohydrodynamica (MHD)-theorie, gevestigd door Furth et al. (1963) en uitgebreid door Coppi et al. (1976), gaat uit van een isotrope evenwichtstensor van de druk. In dit kader zijn de leidende orde schalingen voor de tearing-groeisnelheid en golfgetal onafhankelijk van de plasma-$\beta$ (de verhouding tussen thermische en magnetische druk). Veel astrofysische en ruimteplasma's zijn echter zwak collisioneel, waardoor ze een significante druk-anisotropie kunnen handhaven ($p_\parallel \neq p_\perp$). Hoewel kinetische en fluïdumstudies hebben gesuggereerd dat anisotropie de modustructuur en groeisnelheden beïnvloedt, ontbrak een systematisch analytisch kader dat de voorgeschreven evenwichtsanisotropie verbindt met de standaard outer-inner asymptotische constructie van de tearing-theorie. Specifiek waren de afzonderlijke effecten van anisotropie op de ideale outer matching-parameter ($\Delta'$) en de resistente inner-layer krachtbalans nog niet rigoureus afgeleid.

Methodologie
De auteurs ontwikkelen een lineaire theorie van resistente tearing binnen een niet-ideale gyrotrope MHD-context. De studie richt zich op een krachtloze Harris-type stroomsheet die gekenmerkt wordt door uniforme parallelle en loodrechte evenwichtspredden, gedefinieerd door plasma-$\beta$ ($\beta_0$) en de anisotropie-verschil ($\Delta\beta_0 = \beta_{\parallel,0} - \beta_{\perp,0}$).

1. Beheersende Vergelijkingen: De auteurs maken gebruik van incompressibele, niet-ideale gyrotrope MHD-vergelijkingen die worden gesloten door polytrope relaties voor parallelle en loodrechte drukken ($\gamma_\parallel$ en $\gamma_\perp$). Dit omvat Ohmse en viskeuze verwarmingstermen om het energiebalans te handhaven, volgens de formulering van Hau (2002).
2. Linearisatie: Kleine perturbaties worden toegepast op een "bevroren" evenwichtsprofiel. Het systeem wordt gelineariseerd om een set gekoppelde vergelijkingen af te leiden voor snelheid, magnetisch veld en druk-anisotropie-perturbaties.
3. Asymptotische Analyse: De standaard outer-inner laagconstructie wordt toegepast:
   • Outer Regio: Een ideale regio waar resistiviteit wordt verwaarloosd. De auteurs leiden een tweede-orde differentiaalvergelijking af voor de magnetische flux-perturbatie, waarbij zij identificeren hoe anisotropie de verre-veld vervaldatum en de tearing-stabiliteitsindex $\Delta'$ wijzigt.
   • Inner Regio: Een resistente laag waar de momentumbalans wordt gedomineerd door de anisotrope spanning.
   • Matching: De FKR (constant-$\psi$) en Coppi (niet-constante-$\psi$) takken worden afgeleid en gematcht om de snelst groeiende modus te bepalen.
4. Numerieke Validatie: De analytische voorspellingen worden getest tegen numerieke eigenwaarde-berekeningen met behulp van het PSECAS (Pseudo-Spectral Eigenvalue Calculator with an Automated Solver) framework. Drie specifieke afsluitingen worden onderzocht: dubbel-adiabatisch ($\gamma_\parallel=3, \gamma_\perp=2$), dubbel-polytroop ($\gamma_\parallel=0.5, \gamma_\perp=2$), en dubbel-isotherm ($\gamma_\parallel=1, \gamma_\perp=1$).

Belangrijkste Bijdragen en Resultaten

• Gemodificeerde Outer Regio Dynamica: De auteurs leiden af dat evenwichtsdruk-anisotropie de verre-veld vervaldatum van de tearing-modus verandert. Dit introduceert een nieuwe marginale stabiliteitsafkapwaarde voor gelokaliseerde modi: $\alpha_c = \sqrt{A/R_0}$, waarbij $A = 1 - \Delta\beta_0/2$ de effectieve anisotrope spanning vertegenwoordigt en $R_0$ de drukrespons incorporeert. Modi met golfgetallen $\alpha \geq \alpha_c$ worden gestabiliseerd tegen gelokaliseerde tearing.
• Gemodificeerde Inner Layer Dynamica: In tegen van het isotrope geval, waar anisotropie primair de outer matching beïnvloedt, komt de voorgeschreven evenwichtsanisotropie via de factor $A$ in de leidende orde momentumbalans van de inner layer terecht. Dit wijzigt direct de krachtbalans die zowel de FKR als de Coppi regimes beheerst.
• Herziene Schalingswetten: Door de gemodificeerde FKR en Coppi takken te matchen, leiden de auteurs schalingswetten af voor de snelst groeiende modus. Hoewel de klassieke Lundquist-getal ($S$) exponenten behouden blijven ($\gamma_{max} \sim S^{-1/2}$, $k_{max} \sim S^{-1/4}$), zijn de prefactoren niet langer universeel. Ze zijn afhankelijk van de evenwichtsanisotropie en de afsluiting:
  • $\gamma_{max}\tau_A \sim A^{1/2} R_0^{-1/4} S^{-1/2}$
  • $k_{max}a \sim A^{1/4} R_0^{-3/8} S^{-1/4}$
• Numerieke Bevestiging: PSECAS-berekeningen bevestigen de Coppi-tak schaling en de FKR-tak wanneer een passende eindige-golflengte benadering voor $\Delta'$ wordt gebruikt. De numerieke resultaten valideren de analytische voorspelling dat negatieve $\Delta\beta_0$ (perpendicular druk dominantie) de groeisnelheden versterkt en de snelste modus naar hogere golfgetallen verschuift, terwijl positieve $\Delta\beta_0$ tearing onderdrukt en de inner layer verbredt.
• Afhankelijkheid van Afsluiting: De studie benadrukt dat de omvang van het anisotropie-effect afhangt van de polytrope afsluiting. In de dubbel-isotherme limiet ($\gamma_\parallel=\gamma_\perp=1$) verdwijnt de plasma-$\beta$ afhankelijke drukrespons, waardoor het directe effect van de evenwichtsanisotropie via $A$ en $R_0$ geïsoleerd wordt.

Betekenis en Claims
Het artikel claimt een "gecontroleerde uitbreiding" van de klassieke tearing-theorie naar het gyrotrope-MHD-regime te bieden. De primaire betekenis ligt in het aantonen dat een voorgeschreven evenwichtsdruk-anisotropie de lineaire aanloop van tearing fundamenteel verandert door zowel de ideale outer matching-condities als de resistente inner-layer dynamica te modificeren.

De auteurs stellen expliciet dat hun werk een theoretisch model is dat bedoeld is om het lineaire resistente gyrotrope-MHD regime te verhelderen, in plaats van een direct, volledig kinetisch model van astrofysische reconnectie te bieden. Zij merken op dat de resultaten conditioneel zijn: zij tonen hoe een voorgeschreven anisotropie de stabiliteitsproblematiek wijzigt, maar houden geen rekening met kinetische relaxatie, Hall-fysica, of de zelfconsistente evolutie van de anisotropie. Bijgevolg beargumenteert het artikel dat in hoog-$\beta$ of zwak collisionele omgevingen, het aannemen van druk-isotropie kan leiden tot fouten van orde één bij het schatten van tearing-groeisnelheden, voorkeursgolflengten en inner-layer diktes. Het werk dient als een basis voor toekomstige uitbreidingen die compressibiliteit, zelfconsistente drukprofielen en twee-vloeistof effecten zullen incorporeren.