NUMERICAL SIMULATIONS OF CORONAL HEATING THROUGH FOOTPOINT BRAIDING

Acoustic-gravity wave propagation characteristics in three-dimensional radiation hydrodynamic simulations of the solar atmosphere

There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including ±180° phase jumps. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

On the generation of solar spicules and Alfvénic waves

In the lower solar atmosphere, the chromosphere is permeated by jets known as spicules, in which plasma is propelled at speeds of 50 to 150 kilometers per second into the corona. The origin of the spicules is poorly understood, although they are expected to play a role in heating the million-degree corona and are associated with Alfvénic waves that help drive the solar wind. We compare magnetohydrodynamic simulations of spicules with observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules are shown to occur when magnetic tension is amplified and transported upward through interactions between ions and neutrals or ambipolar diffusion. The tension is impulsively released to drive flows, heat plasma (through ambipolar diffusion), and generate Alfvénic waves.

Quantifying the tangling of trajectories using the topological entropy

We present a simple method to efficiently compute a lower limit of the topological entropy and its spatial distribution for two-dimensional mappings. These mappings could represent either two-dimensional time-periodic fluid flows or three-dimensional magnetic fields, which are periodic in one direction. This method is based on measuring the length of a material line in the flow. Depending on the nature of the flow, the fluid can be mixed very efficiently which causes the line to stretch. Here, we study a method that adaptively increases the resolution at locations along the line where folds lead to a high curvature. This reduces the computational cost greatly which allows us to study unprecedented parameter regimes. We demonstrate how this efficient implementation allows the computation of the variation of the finite-time topological entropy in the mapping. This measure quantifies spatial variations of the braiding efficiency, important in many practical applications.