Thermal stability of magnetohydrodynamic modes in homogeneous plasmas

Abstract

Context. Thermal instabilities give rise to condensations in the solar corona, and are the most probable scenario for coronal rain and prominence formation. We revisit the original theoretical treatment done by Field (1965, ApJ, 142, 531) in a homogeneous plasma with heat-loss effects and combine this with state-of-the-art numerical simulations to verify growth-rate predictions and address the long-term non-linear regime. We especially investigate interaction between multiple magnetohydrodynamic (MHD) wave modes and how they in turn trigger thermal mode development. Aims. We assess how well the numerical MHD simulations retrieve the analytically predicted growth rates. We complete the original theory with quantifications of the eigenfunctions, calculated to consistently excite each wave mode. Thermal growth rates are fitted also in the non-linear regime of multiple wave-wave interaction setups, at the onset of thermal instability formation. Methods. We performed 2D numerical MHD simulations, including an additional (radiative) heat-loss term and a constant heating term to the energy equation. We mainly focus on the thermal (i.e. entropy) and slow MHD wave modes and use the wave amplitude as a function of time to make a comparison to predicted growth rates. Results. It is shown that the numerical MHD simulations retrieve analytically predicted growth rates for all modes, of thermal and slow or fast MHD type. In typical coronal conditions, the latter are damped due to radiative losses, but their interaction can cause slowly changing equilibrium conditions which ultimately trigger thermal mode development. Even in these non-linear wave-wave interaction setups, the growth rate of the thermal instability agrees with the exponential profile predicted by linear theory. The nonlinear evolutions show systematic field-guided motions of condensations with fairly complex morphologies, resulting from thermal modes excited through damped slow MHD waves. These results are of direct interest to the study of solar coronal rain and prominence fine structure. Our wave-wave interaction setups are relevant for coronal loop sections which are known to host slow wave modes, and hence provide a new route to explain the sudden onset of coronal condensation.

Thermal instabilities: Fragmentation and field misalignment of filament fine structure

“…Finally, such a backto-basics approach allows us to omit additional physical effects such as gravity or evaporation-condensation processes and associated flow patterns, which further complicate detailed analysis. This extends our previous work with analogous local 2D simulations (Claes & Keppens 2019) by including anisotropic thermal conduction, three-dimensional simulations, and much higher resolutions.…”

The effect of the magnetic field on the damping of slow waves in the solar corona

“…It is important to remark explicitly that, following from the previous works on thermal misbalance, we allow both the heating and cooling functions to be perturbed. This is is contrast to several previous works in which the heating term is held constant, which is to say remains unperturbed by the wave, such as Claes & Keppens (2019); Kaneko & Yokoyama (2017);De Moortel & Hood (2003) do when setting up their simulations.…”

Seismological constraints on the solar coronal heating function

“…In other words, the wave is subject to a back-reaction of the modification of the cooling and heating rates by the density and temperature perturbations in the wave, i.e., a wave-induced heating/cooling misbalance. This effect can cause enhanced damping or amplification of slow waves (Nakariakov et al 2000;Carbonell et al 2006;Kumar et al 2016;Nakariakov et al 2017;Perelomova 2018;Claes & Keppens 2019;Kolotkov et al 2019). Characteristic time scales of the thermal misbalance lead to the wave dispersion (Ibanez S. & Escalona T. 1993;Zavershinskii et al 2019).…”