From pulsar scintillations to coronal heating: discontinuities in magnetohydrodynamics

TL;DR

This paper investigates the formation of magnetic current sheets in magnetohydrodynamics through topological analysis and numerical simulations across various conditions, revealing their ubiquitous formation and dependence on boundary conditions and parameters.

Contribution

It demonstrates that current sheets form under a wide range of conditions in MHD simulations, expanding understanding of their formation beyond specific boundary setups.

Findings

Current sheets form under all simulated conditions.

Thinner and more numerous sheets appear at higher resolution.

Magnetic field relaxes into a minimum energy state influenced by helicity and boundaries.

Abstract

From pulsar scintillations we infer the presence of sheet-like structures in the ISM; it has been suggested that these are current sheets. Current sheets probably play an important role in heating the solar corona, and there is evidence for their presence in the solar wind. Such magnetic discontinuities have been found in numerical simulations with particular boundary conditions, as well as in simulations using an incompressible equation of state. Here, I investigate their formation under more general circumstances by means of topological considerations as well as numerical simulations of the relaxation of an arbitrary smoothly-varying magnetic field. The simulations are performed with a variety of parameters and boundary conditions: in low, high and of-order-unity plasma-$\beta$ regimes, with periodic and fixed boundaries, with and without a friction force, at various resolutions and with various diffusivities. Current sheets form, over a dynamical timescale, under {\it all} conditions explored. At higher resolution they are thinner, and there is a greater number of weaker current sheets. The magnetic field eventually relaxes into a smooth minimum energy state, the energy of which depends on the magnetic helicity, as well as on the nature of the boundaries.