No Arabic abstract
We study the propagation and dissipation of magnetohydrodynamic waves in a set of numerical models that each include a solar--like stratified atmosphere and a magnetic field with a null point. All simulations have the same magnetic field configuration but different transition region heights. Compressive wave packets introduced in the photospheric portion of the simulations refract towards the null and collapse it into a current sheet, which then undergoes reconnection. The collapsed null forms a current sheet due to a strong magnetic pressure gradient caused by the inability of magnetic perturbations to cross the null. Although the null current sheet undergoes multiple reconnection episodes due to repeated reflections off the lower boundary, we find no evidence of oscillatory reconnection arising from the dynamics of the null itself. Wave mode conversion around the null generates a series of slow mode shocks localized near each separatrix. The shock strength is asymmetric across each separatrix, and subsequent shock damping therefore creates a tangential discontinuity across each separatrix, with long--lived current densities. A parameter study of the injected wave energy to reach the null confirms our previous WKB estimates. Finally, using current estimates of the photospheric acoustic power, we estimate that the shock and Ohmic heating we describe may account for $approx1-10%$ of the radiative losses from coronal bright points with similar topologies, and are similarly insufficient to account for losses from larger structures such as ephemeral regions. At the same time, the dynamics are comparable to proposed mechanisms for generating type--II spicules.
We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that give rise to a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into multiple secondary wavepackets. A portion of the packet refracts towards the null due to the varying Alfven speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfven speeds coincide, partially transmit, reflect, and mode convert between branches of the local dispersion relation. Approximately $15.5%$ of the wavepackets initial energy ($E_{input}$) converges on the null, mostly as a fast magnetoacoustic wave. Conversion is very efficient: $70%$ of the energy incident on the null is converted to slow modes propagating away from the null, $7%$ leaves as a fast wave, and the remaining $23%$ (0.036$E_{input}$) is locally dissipated. The acoustic energy leaving the null is strongly concentrated along field lines near each of the nulls four separatrices. The portion of the wavepacket that refracts towards the null, and the amount of current accumulation, depends on the vertical and horizontal wavenumbers and the centroid position of the wavepacket as it crosses the photosphere. Regions that refract towards or away from the null do not simply coincide with regions of open versus closed magnetic field or the local field orientation. We also modeled wavepacket propagation using a WKB method and found that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation.
The plasma of the lower solar atmosphere consists of mostly neutral particles, whereas the upper solar atmosphere is mostly ionised particles and electrons. A shock that propagates upwards in the solar atmosphere therefore undergoes a transition where the dominant fluid is either neutral or ionised. An upwards propagating shock also passes a point where the sound and Alfven speed are equal. At this point the energy of the acoustic shock can separated into fast and slow components. How the energy is distributed between the two modes depends on the angle of magnetic field. The separation of neutral and ionised species in a gravitationally stratified atmosphere is investigated. The role of two-fluid effects on the structure of the shocks post-mode-conversion and the frictional heating is quantified for different levels of collisional coupling. Two-fluid numerical simulations are performed using the (Punderline{I}P) code of a wave steepening into a shock in an isothermal, partially-ionised atmosphere. The collisional coefficient is varied to investigate the regimes where the plasma and neutral species are weakly, strongly and finitely coupled. The propagation speeds of the compressional waves hosted by neutral and ionised species vary, therefore velocity drift between the two species is produced as the plasma attempts to propagate faster than the neutrals. This is most extreme for a fast-mode shock. We find that the collisional coefficient drastically changes the features present in the system, specifically the mode conversion height, type of shocks present, and the finite shock widths created by the two-fluid effects. In the finitely-coupled regime fast-mode shock widths can exceed the pressure scale height leading to a new potential observable of two-fluid effects in the lower solar atmosphere.
We investigate the dynamical evolution of magnetic fields in closed regions of solar and stellar coronae. To understand under which conditions current sheets form, we examine dissipative and ideal reduced magnetohydrodynamic models in cartesian geometry, where two magnetic field components are present: the strong guide field $B_0$, extended along the axial direction, and the dynamical orthogonal field $mathbf{b}$. Magnetic field lines thread the system along the axial direction, that spans the length $L$, and are line-tied at the top and bottom plates. The magnetic field $b$ initially has only large scales, with its gradient (current) length-scale of order $ell_b$. We identify the magnetic intensity threshold $b/B_0 sim ell_b/L$. For values of $b$ below this threshold, field-line tension inhibits the formation of current sheets, while above the threshold they form quickly on fast ideal timescales. In the ideal case, above the magnetic threshold, we show that current sheets thickness decreases in time until it becomes smaller than the grid resolution, with the analyticity strip width $delta$ decreasing at least exponentially, after which the simulations become under-resolved.
This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer.
We propose a mechanism for the fast dissipation of magnetic fields which is effective in a stratified medium where ion motions can be neglected. In such a medium, the field is frozen into the electrons and Hall currents prevail. Although Hall currents conserve magnetic energy, in the presence of density gradients, they are able to create current sheets which can be the sites for efficient dissipation of magnetic fields. We recover the frequency, $omega_{MH}$, for Hall oscillations modified by the presence of density gradients. We show that these oscillations can lead to the exchange of energy between different components of the field. We calculate the time evolution and show that magnetic fields can dissipate on a timescale of order $1/omega_{MH}$. This mechanism can play an important role for magnetic dissipation in systems with very steep density gradients where the ions are static such as those found in the solid crust of neutron stars.