The Alfvenic nature of chromospheric swirls

Context. Observations show that small-scale vortical plasma motions are ubiquitous in the quiet solar atmosphere. They have received increasing attention in recent years because they are a viable candidate mechanism for the heating of the outer solar atmospheric layers. However, the true nature and the origin of these swirls, and their effective role in the energy transport, are still unclear.Aims. We investigate the evolution and origin of chromospheric swirls by analyzing numerical simulations of the quiet solar atmosphere. In particular, we are interested in finding their relation with magnetic field perturbations and in the processes driving their evolution.Methods. The radiative magnetohydrodynamic code CO5BOLD is used to perform realistic numerical simulations of a small portion of the solar atmosphere, ranging from the top layers of the convection zone to the middle chromosphere. For the analysis, the swirling strength criterion and its evolution equation are applied in order to identify vortical motions and to study their dynamics. As a new criterion, we introduce the magnetic swirling strength, which allows us to recognize torsional perturbations in the magnetic field.Results. We find a strong correlation between swirling strength and magnetic swirling strength, in particular in intense magnetic flux concentrations, which suggests a tight relation between vortical motions and torsional magnetic field perturbations. Furthermore, we find that swirls propagate upward with the local Alfven speed as unidirectional swirls driven by magnetic tension forces alone. In the photosphere and low chromosphere, the rotation of the plasma co-occurs with a twist in the upwardly directed magnetic field that is in the opposite direction of the plasma flow. All together, these are clear characteristics of torsional Alfven waves. Yet, the Alfven wave is not oscillatory but takes the form of a unidirectional pulse. The novelty of the present work is that these Alfven pulses naturally emerge from realistic numerical simulations of the solar atmosphere. We also find indications of an imbalance between the hydrodynamic and magnetohydrodynamic baroclinic effects being at the origin of the swirls. At the base of the chromosphere, we find a mean net upwardly directed Poynting flux of 12.86.5 kW m(-2), which is mainly due to swirling motions. This energy flux is mostly associated with large and complex swirling structures, which we interpret as the superposition of various small-scale vortices.Conclusions. We conclude that the ubiquitous swirling events observed in numerical simulations are tightly correlated with perturbations of the magnetic field. At photospheric and chromospheric levels, they form Alfven pulses that propagate upward and may contribute to chromospheric heating.

A comprehensive study of electrostatic turbulence and transport in the laboratory basic plasma physics device TORPEX

Fabio Avino, Alexandre Dominique Bovet, Ambrogio Fasoli, Ivo Furno, Kyle Gustafson, Davoud Iraji, Benoît Labit, Joaquim Loizu Cisquella, Paolo Ricci, Christian Gabriel Theiler

TORPEX is a toroidal device located at the CRPP-EPFL in Lausanne. In TORPEX, a vertical magnetic field superposed on a toroidal field creates helicoidal field lines with both ends terminating on the torus vessel. The turbulence driven by magnetic curvature and plasma gradients causes plasma transport in the radial direction while at the same time plasma is progressively lost along the field lines. The relatively simple magnetic geometry and diagnostic access of the TORPEX configuration facilitate the experimental study of low frequency instabilities and related turbulent transport, and make an accurate comparison between simulations and experiments possible. We first present a detailed investigation of electrostatic interchange turbulence, associated structures and their effect on plasma using high-resolution diagnostics of plasma parameters and wave fields throughout the whole device cross-section, fluid models and numerical simulations. Interchange modes nonlinearly develop blobs, radially propagating filaments of enhanced plasma pressure. Blob velocities and sizes are obtained from probe measurements using pattern recognition and are described by an analytical expression that includes ion polarization currents, parallel sheath currents and ion-neutral collisions. Then, we describe recent advances of a non-perturbative Li 6+ miniaturized ion source and a detector for the investigation of the interaction between supra thermal ions and interchange–driven turbulence. We present first measurements of the spatial and energy space distribution of the fast ion beam in different plasma scenarios, in which the plasma turbulence is fully characterized. The experiments are interpreted using two-dimensional fluid simulations describing the low-frequency interchange turbulence, taking into account the plasma source and plasma losses at the torus vessel. By treating fast ions as test particles, we integrate their equations of motion in the simulated electromagnetic fields, and we compare their time-averaged and statistical properties with experimental data. Finally, we discuss future developments including the possibility of closing the magnetic field lines and of performing magnetic reconnection experiments.

2012

Numerical Analysis of the Effects of the Magnetic Self-Field on the Transport Properties of a Multilayer HTS Cable

Bertrand Dutoit, Francesco Grilli

In this paper the effects of the magnetic self-field on the transport properties of a multi-layer high-Tc superconducting (HTS) cable are investigated by means of 2D finite element method (FEM) simulations. Analyzed is a 3-layer HTS cable, but the developed methods can be used for a different number of layers. The superconductor is described by the non-linear power-law relation E=Ec(J/Jc)^n, where the parameters Jc and n depend on the magnetic field experienced by the material. This dependence decreases the global transport capacity of the superconductor, enhancing its AC losses. It is shown that, especially at high transport currents, the AC losses are considerably higher than in the case where the dependence on the magnetic field is neglected. A simple electrical model, considering the cable from macroscopic point of view, has been proposed for finding the optimal winding pitches, leading to a uniform current repartition. The use of this electrical model allows to overcome the difficulties of direct 3D FEM computations. In addition, the rapidity of solutions by the electric model gives the possibility of testing quickly many geometrical configurations in order to find the ones leading to an even current repartition. This optimization process would not be possible with detailed FEM simulations.

2004

Numerical modelling of high temperature superconducting tapes and cables

Francesco Grilli

This Ph.D. thesis is focused on the numerical modelling of high-Tc superconductors (HTS) at the operating temperature of 77 K (liquid nitrogen). The purpose of numerical modelling is to precisely calculate the current and field distributions inside HTS devices (tapes, cables) and the corresponding AC losses, which are one of the most important limiting factors for a large-scale application of such materials. From the electrical point of view, superconductors are characterized by a strongly non-linear voltage-current relation, which defines the transition from superconducting to normal state. In the case of HTS, the steepness of this transition is smoother than for low-Tc superconductors (LTS), so that the commonly used critical state model (CSM) gives a too simplified representation of their electromagnetic behaviour and can be used for a qualitative description only. In this thesis the finite element method (FEM) has been used for precisely computing the current and field distributions as well as for evaluating the AC losses in HTS devices. The superconducting transition is modelled with a power-law relation, E(J) = Ec(J/Jc)n, which is derived from the fit of transport measurements. Firstly, the results obtained with the software package FLUX3D on multi-filamentary tapes have been validated by means of a comparison with the ones obtained by another software package (FLUX2D). The results from FLUX2D, having already been successfully compared with experimental measurements within the framework of two previous Ph.D. theses at LANOS, have been used as reference. Secondly, two power-law models, which take into account the spatial variation of the critical current density inside HTS tapes and its strongly anisotropic dependence on the magnetic field, have been implemented in FLUX3D. In most cases, this latter dependence is extremely important, since the transport capacity of the superconductor is considerably reduced (and its power loss sensibly increased) by the presence of a magnetic field. Afterwards, FEM modelling of HTS tapes has been extended to cables. HTS cables are in general composed by different layers of several tapes and have a quite complex three-dimensional structure: in fact, the layers are wound around a central cylindrical support with different pitch lengths and relative winding orientations, in order to obtain a uniform repartition of the transport current among the layers, which minimizes the total AC losses. For overcoming the difficulties of a direct 3D FEM simulation, a simple electrical model, which allows to find the optimal pitch lengths and whose results are the input data for a 2D FEM evaluation of the AC losses, has been developed. FEM computations have also been used to investigate the influence of the non-uniformity of the tape properties (contact resistance, critical current, power index) on the global loss performance of a single-layer HTS cable. As an alternative to FEM computations, an equivalent circuit model of HTS cables has been utilized. It describes the cable from the macroscopic point of view and allows to compute the current repartition among the layers and the corresponding AC losses, without the necessity of having detailed information about individual tapes. In the framework of the European project BIG-POWA, I have collaborated to the assembling process of a HTS power-link at the Pirelli Labs, in the Milan region, Italy. Finally, 3D simulations have been performed to extensively study the coupling effect between two superconducting filaments via the resistive matrix, which is a typical case where analytical solutions exist for very peculiar geometries and physical conditions only. FEM simulations have been utilized to study the dependence of the filament coupling on the physical and geometrical parameters of the conductors.