In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermodynamic effects on the energy release process. Specifically, we investigate convergence of the magnetic field at the loop footpoints, atmospheric stratification, as well as thermal conduction. In all cases, the application of vortical driving at the footpoints of an initally potential field leads to an internal kink instability. With the exception of the curved loop with high footpoint convergence, the global geometry of the loop change little during the simulation. Footpoint convergence, curvature and atmospheric structure clearly influences the rapidity with which a loop achieves instability as well as the size of the subsequent energy release. Footpoint convergence has a stabilising influence and thus the loop requires more energy for instability, which means that the subsequent relaxation has a larger heating effect. Large-scale curvature has the opposite result: less energy is needed for instability and so the amount of energy released from the field is reduced. Introducing a stratified atmosphere gives rise to decaying wave phenomena during the driving phase, and also results in a loop that is less stable.

Context: Relaxation theory offers a straightforward method for estimating the energy that is released when a magnetic field becomes unstable, as a result of continual convective driving. Aims: We present new results obtained from nonlinear magnetohydrodynamic (MHD) simulations of idealised coronal loops. The purpose of this work is to determine whether or not the simulation results agree with Taylor relaxation, which will require a modified version of relaxation theory applicable to unbounded field configurations. Methods: A three-dimensional (3D) MHD Lagrangian-remap code is used to simulate the evolution of a line-tied cylindrical coronal loop model. This model comprises three concentric layers surrounded by a potential envelope; hence, being twisted locally, each loop configuration is distinguished by a piecewise-constant current profile. Initially, all configurations carry zero-net-current fields and are in ideally unstable equilibrium. The simulation results are compared with the predictions of helicity conserving relaxation theory. Results: For all simulations, the change in helicity is no more than 2% of the initial value; also, the numerical helicities match the analytically-determined values. Magnetic energy dissipation predominantly occurs via shock heating associated with magnetic reconnection in distributed current sheets. The energy release and final field profiles produced by the numerical simulations are in agreement with the predictions given by a new model of partial relaxation theory: the relaxed field is close to a linear force free state; however, the extent of the relaxation region is limited, while the loop undergoes some radial expansion. Conclusions: The results presented here support the use of partial relaxation theory, specifically, when calculating the heating-event distributions produced by ensembles of kink-unstable loops.

It has been proposed that the million degree temperature of the corona is due to the combined effect of barely-detectable energy releases, so called nanoflares, that occur throughout the solar atmosphere. Alas, the nanoflare density and brightness implied by this hypothesis means that conclusive verification is beyond present observational abilities. Nevertheless, we investigate the plausibility of the nanoflare hypothesis by constructing a magnetohydrodynamic (MHD) model that can derive the energy of a nanoflare from the nature of an ideal kink instability. The set of energy-releasing instabilities is captured by an instability threshold for linear kink modes. Each point on the threshold is associated with a unique energy release and so we can predict a distribution of nanoflare energies. When the linear instability threshold is crossed, the instability enters a nonlinear phase as it is driven by current sheet reconnection. As the ensuing flare erupts and declines, the field transitions to a lower energy state, which is modelled by relaxation theory, i.e., helicity is conserved and the ratio of current to field becomes invariant within the loop. We apply the model so that all the loops within an ensemble achieve instability followed by energy-releasing relaxation. The result is a nanoflare energy distribution. Furthermore, we produce different distributions by varying the loop aspect ratio, the nature of the path to instability taken by each loop and also the level of radial expansion that may accompany loop relaxation. The heating rate obtained is just sufficient for coronal heating. In addition, we also show that kink instability cannot be associated with a critical magnetic twist value for every point along the instability threshold.

Context: It is thought likely that vast numbers of nanoflares are responsible for the corona having a temperature of millions of degrees. Current observational technologies lack the resolving power to confirm the nanoflare hypothesis. An alternative approach is to construct a magnetohydrodynamic coronal loop model that has the ability to predict nanoflare energy distributions. Aims: This paper presents the initial results generated by such a model. It predicts heating events with a range of sizes, depending on where the instability threshold for linear kink modes is encountered. The aims are to calculate the distribution of event energies and to investigate whether kink instability can be predicted from a single parameter. Methods: The loop is represented as a straight line-tied cylinder. The twisting caused by random photospheric motions is captured by two parameters, representing the ratio of current density to field strength for specific regions of the loop. Dissipation of the loop's magnetic energy begins during the nonlinear stage of the instability, which develops as a consequence of current sheet reconnection. After flaring, the loop evolves to the state of lowest energy where, in accordance with relaxation theory, the ratio of current to field is constant throughout the loop and helicity is conserved. Results: The results suggest that instability cannot be predicted by any simple twist-derived property reaching a critical value. The model is applied such that the loop undergoes repeated episodes of instability followed by energy-releasing relaxation. Hence, an energy distribution of the nanoflares produced is collated. Conclusions: The final energy distribution features two nanoflare populations that follow different power laws. The power law index for the higher energy population is more than sufficient for coronal heating.