The kink instability of magnetohydrodynamics is believed to be fundamental to many aspects of the dynamic activity of the solar atmosphere, such as the initiation of flares and the heating of the solar corona. In this work, we investigate the importance of viscosity on the kink instability. In particular, we focus on two forms of viscosity; isotropic viscosity (independent of the magnetic field) and anisotropic viscosity (with a preferred direction following the magnetic field). Through the detailed analysis of magnetohydrodynamic simulations of the kink instability with both types of viscosity, we show that the form of viscosity has a significant effect on the nonlinear dynamics of the instability. The different viscosities allow for different flow and current structures to develop, thus affecting the behaviour of magnetic relaxation, the formation of secondary instabilities and the Ohmic and viscous heating produced. Our results have important consequences for the interpretation of solar observations of the kink instability.
We investigate single-fluid magnetohydrodynamics (MHD) with anisotropic viscosity, often referred to as Braginskii MHD, with a particular eye to solar coronal applications. First, we examine the full Braginskii viscous tensor in the single-fluid limit. We pay particular attention to how the Braginskii tensor behaves as the magnetic field strength vanishes. The solar corona contains a magnetic field with a complex and evolving topology, so the viscosity must revert to its isotropic form when the field strength is zero, e.g. at null points. We highlight that the standard form in which the Braginskii tensor is written is not suitable for inclusion in simulations as singularities in the individual terms can develop. Instead, an altered form, where the parallel and perpendicular tensors are combined, provides the required asymptotic behaviour in the weak-field limit. We implement this combined form of the tensor into the Lare3D code, which is widely used for coronal simulations. Since our main focus is the viscous heating of the solar corona, we drop the drift terms of the Braginskii tensor. In a stressed null point simulation, we discover that small-scale structures, which develop very close to the null, lead to anisotropic viscous heating at the null itself (that is, heating due to the anisotropic terms in the viscosity tensor). The null point simulation we present has a much higher resolution than many other simulations containing null points so this excess heating is a practical concern in coronal simulations. To remedy this unwanted heating at the null point, we develop a model for the viscosity tensor that captures the most important physics of viscosity in the corona: parallel viscosity for strong field and isotropic viscosity at null points. We derive a continuum model of viscosity where momentum transport, described by this viscosity model, has the magnetic field as its preferred orientation. When the field strength is zero, there is no preferred direction for momentum transport and viscosity reverts to the standard isotropic form. The most general viscous stress tensor of a (single-fluid) plasma satisfying these conditions is found. It is shown that the Braginskii model, without the drift terms, is a specialization of the general model. Performing the stressed null point simulation with this simplified model of viscosity reveals very similar heating profiles compared to the full Braginskii model. The new model, however, does not produce anisotropic heating at the null point, as required. Since the vast majority of coronal simulations use only isotropic viscosity, we perform the stressed null point simulation with isotropic viscosity and compare the heating profiles to those of the anisotropic models. It is shown than the fully isotropic viscosity can over-estimate the viscous heating by an order of magnitude.
Context. Magnetic null points are associated with high-energy coronal phenomena such as solar flares and are often sites of reconnection and particle acceleration. Dynamic twisting of a magnetic null point can generate a Kelvin-Helmholtz instability (KHI) within its fan plane and can instigate spine-fan reconnection and an associated collapse of the null point under continued twisting. Aims. This article aims to compare the effects of isotropic and anisotropic viscosity in simulations of the KHI and collapse in a dynamically twisted magnetic null point. Methods. We performed simulations using the 3D magnetohydrodynamics code Lare3d with a custom anisotropic viscosity module. A pair of high-resolution simulations were performed, one using isotropic viscosity and another using anisotropic viscosity, keeping all other factors identical. We analysed the results in detail. A further parameter study was performed over a range of values for viscosity and resistivity. Results. Both viscosity models permit the growth of the KHI and the eventual collapse of the null point. Over all studied parameters, anisotropic viscosity allows a faster growing instability, while isotropic viscosity damps the instability to the extent of stabilisation in some cases. Although the viscous heating associated with anisotropic viscosity is generally smaller, the ohmic heating dominates and is enhanced by the current sheets generated by the instability. This leads to a greater overall heating rate when using anisotropic viscosity. The collapse of the null point occurs significantly sooner when anisotropic viscosity is employed.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.