Search

The ALLIANCE code is developed to solve a new set of four-dimensional electromagnetic drift-kinetic equations in slab geometry [J. Plasma Phys. 88 905880117 (2022)]. The nonlinear equations are useful for the study of magnetized plasma systems at scales comparable to, or larger than the ion gyroradius. In particular, it is suited for the study the kinetic turbulent cascade in astrophysical plasma, while preserving finite Larmor radius effects at the fluid-kinetic transition. The equations solved are in spectral Fourier-Laguerre-Hermite form, a pseudo-spectral approach is used for the nonlinear terms, and the code is parallelised over multiple directions. After a presentation of the code, validation runs are shown, and benchmarks for serial and parallel computations are presented.

We present the TRIDK code, which is written in Python, that solves Triad Reduced Interactions for a Drift Kinetic system of equations. The four-dimensional drift kinetic system of equations captures nonlinear electromagnetic (or electrostatic) interactions for a magnetized plasma at scales larger than the ion gyroradius. In the TRIDK code, the nonlinear interaction are modelled using a compound of nested icosahedra-dodecohedra shells in wave space, which allows to perform computationally efficient simulations over a large span of scales. This reduced approach captures essential nonlinear dynamics, while allowing parallel mixing in velocity space to take place. The intent is to offer a cost efficient way to probe the balance between linear and nonlinear dynamics, and analyse the interplay between linear Landau damping and turbulence.

For a magnetised plasma in a straight magnetic guide field, we derive a set of four-dimensional kinetic equations, which can capture electromagnetic turbulence in the drift kinetic limit. To do so, we start from the gyrokinetic equations, employ a Laguerre decomposition in the perpendicular velocity direction, retain only the dominant gyroaverage contributions and only the first two Laguerre moments that source the electromagnetic fluctuations. The model conserves free energy, and can describe electromagnetic turbulence for a plasma at the transition between fluid and gyrokinetic regimes ($k_\perp \rho_i\approx 1$ range of scales), as dominant finite-Larmor-radius (FLR) effects are considered. In addition to the three dimensions in positions space, we retain the parallel velocity dependence, which we describe via a Hermite representation. Employing this system, but without any other physics based assumptions for the plasma species that can bias results, will allow us to investigate how fluid effects transition into the kinetic range, and analyse the interplay between spatial and velocity space mixing for electromagnetic plasma turbulence.

In the present paper, we use a coarse-graining approach to investigate the nonlinear redistribution of free energy in both position and scale space for weakly collisional magnetised plasma turbulence. For this purpose, we use high-resolution numerical simulations of gyrokinetic (GK) turbulence that span the proton-electron range of scales, in a straight magnetic guide field geometry. Accounting for the averaged effect of the particles' fast gyro-motion on the slow plasma fluctuations, the GK approximation captures the dominant energy redistribution mechanisms in strongly magnetised plasma turbulence. Here, the GK system is coarse-grained with respect to a cut-off scale, separating in real space the contributions to the nonlinear interactions from the coarse-grid-scales and the sub-grid-scales (SGS). We concentrate on the analysis of nonlinear SGS effects. Not only that this allows us to investigate the flux of free energy across the scales, but also to now analyse its spatial density. We find that the net value of scale flux is an order of magnitude smaller than both the positive and negative flux density contributions. The dependence of the results on the filter type is also analysed. Moreover, we investigate the advection of energy in position space. This rather novel approach for GK turbulence can help in the development of SGS models that account for advective unstable structures for space and fusion plasmas, and with the analysis of the turbulent transport saturation., Comment: 15 figures Accepted for publication by Journal of Plasma Physics

In tokamak plasmas, sheared flows perpendicular to the driving temperature gradients can strongly stabilize linear modes. While the system is linearly stable, regimes with persistent nonlinear turbulence may develop, i.e. the system is subcritical. A perturbation with small but finite amplitude may be sufficient to push the plasma into a regime where nonlinear effects are dominant and thus allow sustained turbulence. The minimum threshold for nonlinear instability to be triggered provides a criterion for assessing whether a tokamak is likely to stay in the quiescent (laminar) regime. At the critical amplitude, instead of transitioning to the turbulent regime or decaying to a laminar state, the trajectory will map out the edge of chaos. Surprisingly, a quasi-traveling-wave solution is found as an attractor on this edge manifold. This simple advecting solution is qualitatively similar to, but simpler than, the avalanche-like bursts seen in earlier turbulent simulations and provides an insight into how turbulence is sustained in subcritical plasma systems. For large flow shearing rate, the system is only convectively unstable, and given a localised initial perturbation, will eventually return to a laminar state at a fixed spatial location.

In many plasma systems, introducing a small background shear flow is enough to stabilize the system linearly. The nonlinear dynamics are much less sensitive to sheared flows than the average linear growthrates, and very small amplitude perturbations can lead to sustained turbulence. We explore the general problem of characterizing how and when the transition from near-laminar states to sustained turbulence occurs; a model of the interchange instability being used as a concrete example. These questions are fundamentally nonlinear, and the answers must go beyond the linear transient amplification of small perturbations. Two methods that account for nonlinear interactions are therefore explored here. The first method explored is edge tracking, which identifies the boundary between the basins of attraction of the laminar and turbulent states. Here, the edge is found to be structured around an exact, localized, traveling wave solution; a solution that is qualitatively similar to avalanche-like bursts seen in the turbulent regime. The second method is an application of nonlinear, non-modal stability theory which allows us to identify the smallest disturbances which can trigger turbulence (the minimal seed for the problem) and hence to quantify how stable the laminar regime is. The results obtained from these fully nonlinear methods provides confidence in the derivation of a semi-analytic approximation for the minimal seed.

Using large resolution numerical simulations of GK turbulence, spanning an interval ranging from the end of the fluid scales to the electron gyroradius, we study the energy transfers in the perpendicular direction for a proton-electron plasma in a slab magnetic geometry. In addition, to aid our understanding of the nonlinear cascade, we use an idealized test representation for the energy transfers between two scales, mimicking the dynamics of turbulence in an infinite inertial range. For GK turbulence, a detailed analysis of nonlinear energy transfers that account for the separation of energy exchanging scales is performed. We show that locality functions associated with the energy cascade across dyadic (i.e. multiple of two) separated scales achieve an asymptotic state, recovering clear values for the locality exponents. We relate these exponents to the energy exchange between two scales, diagnostics that are less computationally intensive than the locality functions. It is the first time asymptotic locality is shown to exist for GK turbulence and the contributions made by highly non-local interactions, previously reported in the literature, are explained as very local transfers of energy that occur between wavenumbers within the same dyadic signal. The results presented here and their implications are discussed from the perspective of previous findings reported in the literature and the idea of universality of GK turbulence.

We employ magnetized plasma turbulence, described by a gyrokinetic formalism in an interval ranging from the end of the fluid scales to the electron gyroradius, to introduce the first study of kinetic intermittency, in which nonlinear structures formed directly in the distribution functions are analyzed by accounting for velocity space correlations generated by linear (Landau resonance) and nonlinear phase mixing. Electron structures are found to be strongly intermittent and dominated by linear phase mixing, while nonlinear phase mixing dominates the weakly intermittent ions. This is the first time spatial intermittency and linear phase mixing are shown to be self-consistently linked for the electrons and, as the magnetic field follows the intermittency of the electrons at small scales, explain why magnetic islands are places dominated by Landau damping in steady state turbulence.

For a Z-pinch geometry, we report on the nonlinear redistribution of free energy across scales perpendicular to the magnetic guide field, for a turbulent plasma described in the framework of gyrokinetics. The analysis is performed using a local flux-surface approximation, in a regime dominated by electrostatic fluctuations driven by the entropy mode, with both ion and electron species being treated kinetically. To explore the anisotropic nature of the free energy redistribution caused by the emergence of zonal flows, we use a polar coordinate representation for the field-perpendicular directions and define an angular density for the scale flux. Positive values for the classically defined (angle integrated) scale flux, which denote a direct energy cascade, are shown to be also composed of negative angular sections, a fact that impacts our understanding of the backscatter of energy and the way in which it enters the modeling of sub-grid scales for turbulence. A definition for the flux of free energy across each perpendicular direction is introduced as well, which shows that the redistribution of energy in the presence of zonal flows is highly anisotropic.

Pitch-angle scattering rates for cosmic-ray particles in magnetohydrodynamic (MHD) simulations with imbalanced turbulence are calculated for fully evolving electromagnetic turbulence. We compare with theoretical predictions derived from the quasilinear theory of cosmic-ray diffusion for an idealized slab spectrum and demonstrate how cross helicity affects the shape of the pitch-angle diffusion coefficient. Additional simulations in evolving magnetic fields or static field configurations provide evidence that the scattering anisotropy in imbalanced turbulence is not primarily due to coherence with propagating Alfven waves, but an effect of the spatial structure of electric fields in cross-helical MHD turbulence., Comment: 13 pages, 15 figures. Accepted by ApJ

In magnetized plasma turbulence, the couplings of perpendicular spatial scales that arise due to the nonlinear interactions are analyzed from the perspective of the free-energy exchanges. The plasmas considered here, with appropriate ion or electron adiabatic electro-neutrality responses, are described by the gyrokinetic formalism in a toroidal magnetic geometry. Turbulence develops due to the electrostatic fluctuations driven by temperature gradient instabilities, either ion temperature gradient (ITG) or electron temperature gradient (ETG). The analysis consists in decomposing the system into a series of scale structures, while accounting separately for contributions made by modes possessing special symmetries (e.g., the zonal flow modes). The interaction of these scales is analyzed using the energy transfer functions, including a forward and backward decomposition, scale fluxes and locality functions. The comparison between the ITG and ETG cases shows that ETG turbulence has a more pronounced classical turbulent behavior, exhibiting a stronger energy cascade, with implications for gyrokinetic turbulence modeling., Comment: 15 pages, 14 figures

The present work investigates the acceleration of test particles in balanced and imbalanced Alfv\'{e}nic turbulence, relevant to the solar-wind problem. These turbulent states, obtained numerically by prescribing the injection rates for the ideal invariants, are evolved dynamically with the particles. While the energy spectrum for balanced and imbalanced states is known, the impact made on particle heating is a matter of debate, with different considerations giving different results. By performing direct numerical simulations, resonant and non-resonant particle accelerations are automatically considered and the correct turbulent phases are taken into account. For imbalanced turbulence, it is found that the acceleration rate of charged particles is reduced and the heating rate diminished. This behaviour is independent of the particle gyroradius, although particles that have a stronger adiabatic motion (smaller gyroradius) tend to experience a larger heating., Comment: 5 pages 6 figures

The nature of the nonlinear interactions in gyrokinetic (GK) turbulence, driven by an ion-temperature gradient instability, is investigated using numerical simulations of single ion species plasma in three-dimensional flux tube geometry. To account for the level of separation existing between scales involved in an energetic interaction, the degree of locality of the free energy scale flux is analyzed employing Kraichnan's infrared (IR) and ultraviolet locality functions. Due to the nontrivial dissipative nature of GK turbulence, an asymptotic level for the locality exponents, indicative of a universal dynamical regime for GK's, is not recovered and an accentuated non-local behavior of the IR interactions is found instead, in spite of the local energy cascade observed., Comment: 4 pages, 3 figures