Search

The efficacy of fast – slow MHD mode conversion in the surface layers of sunspots has been demonstrated over recent years using a number of modelling techniques, including ray theory, perturbation theory, differential eigensystem analysis, and direct numerical simulation. These show that significant energy may be transferred between the fast and slow modes in the neighbourhood of the equipartition layer where the Alfvén and sound speeds coincide. However, most of the models so far have been two dimensional. In three dimensions the Alfvén wave may couple to the magnetoacoustic waves with important implications for energy loss from helioseismic modes and for oscillations in the atmosphere above the spot. In this paper, we carry out a numerical “scattering experiment,” placing an acoustic driver 4 Mm below the solar surface and monitoring the acoustic and Alfvénic wave energy flux high in an isothermal atmosphere placed above it. These calculations indeed show that energy conversion to upward travelling Alfvén waves can be substantial, in many cases exceeding loss to slow (acoustic) waves. Typically, at penumbral magnetic field strengths, the strongest Alfvén fluxes are produced when the field is inclined 30° – 40° from the vertical, with the vertical plane of wave propagation offset from the vertical plane containing field lines by some 60° – 80°. [ABSTRACT FROM AUTHOR]

The solar corona, modeled by a low- β, resistive plasma slab, sustains MHD wave propagations due to footpoint motions in the photosphere. Simple test cases are undertaken to verify the code. Uniform, smooth and steep density, magnetic profile and driver are considered. The numerical simulations presented here focus on the evolution and properties of the Alfvén, fast and slow waves in coronal loops. The plasma responds to the footpoint motion by kink or sausage waves depending on the amount of shear in the magnetic field. The larger twist in the magnetic field of the loop introduces more fast-wave trapping and destroys initially developed sausage-like wave modes. The transition from sausage to kink waves does not depend much on the steep or smooth profile. The slow waves develop more complex fine structures, thus accounting for several local extrema in the perturbed velocity profiles in the loop. Appearance of the remnants of the ideal singularities characteristic of ideal plasma is the prominent feature of this study. The Alfvén wave which produces remnants of the ideal x singularity, reminiscent of Alfvén resonance at the loop edges, becomes less pronounced for larger twist. Larger shear in the magnetic field makes the development of pseudo-singularity less prominent in case of a steep profile than that in case of a smooth profile. The twist also causes heating at the edges, associated with the resonance and the phase mixing of the Alfvén and slow waves, to slowly shift to layers inside the slab corresponding to peaks in the magnetic field strength. In addition, increasing the twist leads to a higher heating rate of the loop. Remnants of the ideal log ¦ x¦ singularity are observed for fast waves for larger twist. For slow waves they are absent when the plasma experiences large twist in a short time. The steep profiles do not favour the creation of pseudo-singularities as easily as in the smooth case. [ABSTRACT FROM AUTHOR]

A WKB approach, based on the method of Connor, Hastie, and Taylor (1979), is used to obtain simple estimates of the critical conditions for the onset of ideal MHD instabilities in line-tied solar coronal loops. The method is illustrated for the constant twist, Gold-Hoyle (1960) field, and the critical conditions are compared with previous and new numerical results. For the force-free case, the WKB estimate for the critical loop length reduces to $$2\pi m + \sqrt 2 \pi$$ . For the sufficiently non-force-free case the critical length can be expressed in the form l + l / m. The results confirm the findings of De Bruyne and Hood (1992) that for force-free fields the m = 1 mode is the first mode to become unstable but for the sufficiently strong non-force-free case this reverses with the m → ∞ mode being excited first. [ABSTRACT FROM AUTHOR]

The normal mode spectrum for the linearized MHD equations is investigated for a cylindrical equilibrium. This spectrum is examined for zero perpendicular thermal conduction, with both zero and non-zero scalar resistivity. Particular attention is paid to the continuous branches of this spectrum, or continuous spectra. For zero resistivity there are three types of continuous spectra present, namely the Alfvén, slow, and thermal continua. It is shown that when dissipation due to resistivity is included, the slow and Alfvén continua are removed and that the thermal continuum is shifted to a different position (where the shift is independent of the exact value of resistivity). The 'old' location of the thermal continuum is covered by a dense set of nearly singular discrete modes called a quasi-continuum. The quasi-continuum is investigated numerically, and the eigenfunctions are shown to have rapid spatial oscillating behaviour. These oscillations are confined to the most unstable part of the equilibrium based on the Field criterion, and may be the cause of fine structure in prominences. [ABSTRACT FROM AUTHOR]

Approximate solutions of the linearized non-adiabatic MHD equations, obtained using the ballooning method, are compared with 'exact' numerical solutions of the full equations (including the effects of optically thin plasma radiation). It is shown that the standard ballooning method, developed within the framework of ideal linear MHD, can be generalized to non-ideal linear MHD. The localized (ballooning) spectrum has to be used with caution, but can give valuable (though limited) information on non-ideal stability. The numerical analysis also confirms and quantifies the interesting connection between magnetic and thermal instabilities. The existence of such a coupling is inherent in many qualitative discussions of magnetic disruptions. Finally, the hitherto unrecognized role of the thermal continuum in the unstable part of the 'magnetothermal' spectrum is investigated. [ABSTRACT FROM AUTHOR]

The linear MHD spectrum is investigated for cylindrical equilibrium models under typical coronal conditions. Non-ideal effects are included and attention is focussed on the thermal instability and the influence of perpendicular thermal conduction. It is shown that, when thermal conduction across magnetic field lines is neglected, the 'classic' Alfvén and slow continua are supplemented by a new 'thermal' continuum. Surprisingly, the existence of this non-ideal continuous spectrum appears to have been overlooked for a very long time. Unlike the (still purely oscillatory) Alfvén continuum modes and the slow continuum modes (overstable or damped), the thermal continuum modes are exponentially growing or decaying in time. As with the Alfvén and slow continua, discrete modes may be present above or below the thermal continuum, depending upon the choice of equilibrium parameters. These modes are localized using a simple WKB approach. The knowledge of the thermal subspectrum is then exploited to find necessary and sufficient conditions for instability. The inclusion of perpendicular thermal conduction eliminates this thermal continuum, but replaces it by a dense set of eigenmodes located at the same values of the growth rate. It is shown that perpendicular thermal conduction has no significant influence on the thermal stability of the equilibrium, though individual modes may be strongly influenced. The corresponding eigenfunctions still have a nearly-singular behaviour, but in addition they may contain a rapid spatial oscillation on a scale proportional to some power of the perpendicular thermal conduction coefficient. This rapid oscillation is confined to the most unstable part of the equilibrium and may be of relevance for the formation of fine-scale structure (threads) in prominences. [ABSTRACT FROM AUTHOR]

The stability of ballooning modes in coronal arcades is studied using linear visco-resistive MHD. Rigid wall conditions are adopted for modelling the photospheric line-tying of the magnetic field. The full Braginskii viscosity stress tensor is used and particular attention is given to the effect of the viscosity coefficient η which was left out of an earlier investigation by Van der Linden, Goossens, and Hood (1987, 1988). The numerical results for shearless arcades show that the coefficient η has a stabilizing effect. However, for realistic values of the equilibrium quantities the stabilizing effect by η can be neglected in comparison with the strong stabilizing effect of the perpendicular viscosity. The effect of magnetic field strength and mode number on stability are determined. In particular it is found that there exists a critical field strength for every mode number such that the mode is stable for weaker fields and unstable for stronger fields. [ABSTRACT FROM AUTHOR]

Prominences and filaments are thought to arise as a consequence of a magnetized plasma undergoing thermal instability. Therefore, the thermal stability of a magnetized plasma is investigated under coronal conditions. The equilibrium structure of the plasma is approximated by a 1-D slab configuration. This is investigated in thermal instability taking into account optically thin plasma radiation and anisotropic thermal conduction. The thermal conduction perpendicular to the magnetic field is taken to be small but non-zero. The classic rigid wall boundary conditions which are often applied in the literature, either directly on the plasma or indirectly through some other medium, are replaced by a more physical situation in which the plasma column is placed in a low-density background stretching towards infinity. Results for a uniform equilibrium structure indicate the major effect of this change is on the eigenfunctions rather than on the growth rate. Essentially, perpendicular thermal conduction introduces field-aligned fine structure. It is also shown that in the presence of perpendicular thermal conduction, thermal instability in a slab model is only possible if the inner plasma has the shortest thermal instability time scale. [ABSTRACT FROM AUTHOR]

The stability equations for localized (or ballooning) modes in the solar atmosphere are formulated. Dissipation due to viscosity, resistivity, and thermal conduction are included using the general forms due to Braginskii (1965). In addition, the effect of gravity, plasma radiation, and coronal heating are included. The resulting equations are one-dimensional and only involve derivatives along the equilibrium magnetic field. Thus, the stabilising influence of photospheric line-tying, which is normally neglected in most numerical simulations, can be studied in a simple manner. Two applications to sound wave propagation and thermal instabilities in a low-beta plasma are considered with a view to determining realistic coronal boundary conditions that model the lower, denser levels of the solar atmosphere in a simple manner. [ABSTRACT FROM AUTHOR]

The study of resistive ballooning instabilities in line-tied coronal magnetic fields is extended by including viscosity in the stability analysis. The equations that govern the resistive ballooning instabilities are derived and the effects of parallel and perpendicular viscosity are included using Braginskii's stress tensor. Numerical solutions to these equations are obtained under the rigid wall boundary conditions for arcades with cylindrically-symmetric magnetic fields. It is found that viscosity has a stabilizing effect on the resistive ballooning instabilities with perpendicular viscosity being more important by far than parallel viscosity. The strong stabilizing effect of perpendicular viscosity can lead to complete stabilization for realistic values of the equilibrium quantities. [ABSTRACT FROM AUTHOR]

A first attempt is made to study the continuous spectrum of linear ideal MHD for 2D solar loops and to understand how 2D effects change the continuum eigenfrequencies and continuum eigenfunctions. The continuous spectrum is computed for 2D solar loops with purely poloidal magnetic fields and it is investigated how non-circularity of the cross-sections of the poloidal magnetic surfaces and variations of density along the poloidal magnetic field lines change the continuous spectrum and induce poloidal wave number coupling in the eigenfunctions. Approximate analytical results and numerical results are obtained for the eigenfrequencies and the eigenfunctions and the poloidal wave number coupling is clearly illustrated. It is found that the continuum frequencies are substantially increased, that the ranges of the continuum frequencies are considerably enlarged and that the derivatives of the continuum frequencies normal to the magnetic surfaces are substantially increased. The eigenfunctions are strongly influenced by poloidal wave number coupling. Implications of these findings for the heating mechanisms of resonant absorption and phase mixing are briefly considered. [ABSTRACT FROM AUTHOR]

The continuous spectrum of a 2D magnetostatic equilibrium with y-invariance is derived. It is shown that the continuous spectrum is given by an eigenvalue problem on each magnetic surface and is related to the different behaviour of the equilibrium quantities in different magnetic surfaces. The special case of a uniform poloidal magnetic field in a 1D equilibrium that is stratified with height, has been considered in detail and it is found that there is no continuous spectrum. It is shown that this result relies completely on the artificial property that the behaviour of the equilibrium quantities along a magnetic field line is independent of the field line considered. As a consequence the non-existence of a continuous spectrum in a 1D equilibrium with a uniform magnetic field cannot be used to argue that the continuous spectrum has no physical relevance. [ABSTRACT FROM AUTHOR]