Search

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case [Oliveira D. S. et al., Nucl. Fusion 62, 096001 (2022)] and compare these results with previously validated grid simulations [Ulbl P. et al., Phys. Plasmas 30, 107986 (2023)]. The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. Thus, the spectral approach reduces the computational cost by at least an order of magnitude, achieving a speed-up of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge - despite this being the shot with the highest fusion output - is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization.

We perform 3D3V hybrid-Vlasov simulations of turbulence with quasi-isotropic, compressible injection near ion scales to mimic the Earth's magnetosheath plasma, and investigate the novel electron-only reconnection, recently observed by the NASA's MMS mission, and its impact on ion heating. Retaining electron inertia in the generalized Ohm s law enables collisionless magnetic reconnection. Spectral analysis shows a shift from kinetic Alfv\'en waves (KAW) to inertial kinetic Alfv\'en (IKAW) and inertial whistler waves (IWW) near electron scales. To distinguish the roles of inertial scale and gyroradius ($d_{\rm{i}}$ and $\rho_{\rm{i}}$), three ion beta ($\beta_{\rm{i}} = 0.25, 1, 4$) values are studied. Ion-electron decoupling increases with $\beta_{\rm{i}}$, as ions become less mobile when the injection scale is closer to $\rho_{\rm{i}}$ than $d_{\rm{i}}$, highlighting the role of $\rho_{\rm{i}}$ in achieving an electron magnetohydrodynamic (EMHD) regime at sub-ion scales. This regime promotes electron-only reconnection in turbulence with small-scale injection at $\beta_{\rm{i}} \gtrsim 1$. We observe significant ion heating even at large $\beta_{\rm{i}}$, with $Q_{\rm{i}}/\epsilon \approx 69\%, 91\%, 96\%$ at $\beta_{\rm{i}} = 0.25, 1, 4$ respectively. While ion heating is anisotropic at $\beta_{\rm{i}} \leq 1$ ($T_{\rm i,\perp} > T_{\rm i,\parallel}$), it is marginally anisotropic at $\beta_{\rm{i}} > 1$ ($T_{\rm i,\perp} \gtrsim T_{\rm i,\parallel}$). These findings have implications for other collisionless astrophysical environments, like high-$\beta$ plasmas in intracluster medium, where processes such as micro-instabilities or shocks may inject energy near ion-kinetic scales., Comment: Accepted

Turbulent transport is a decisive factor in determining the pedestal structure of H-modes. Here, we present the first comprehensive characterization of gyrokinetic turbulent transport in a JET hybrid H-mode pedestal. Local, linear simulations are performed to identify instabilities and global, nonlinear electromagnetic simulations reveal the turbulent heat and particle flux structure of the pedestal. Our analysis focuses on the Deuterium reference discharge \#97781 performed in the scenario development for the Deuterium-Tritium campaign. We find the pedestal top transport to be dominated by ion temperature gradient (ITG) modes. In the pedestal center turbulent ion transport is suppressed and electron transport is driven by multi-faceted electron temperature gradient (ETG) modes, which extend down to ion-gyroradius scales. A strong impact of $E\times B$ shear on the absolute turbulence level is confirmed by the global, nonlinear simulations. Furthermore, impurities are shown to reduce the main ion transport. Dedicated density and ion temperature profile variations test the sensitivity of the results and do not find strong differences in the turbulent transport in more reactor-like conditions.

We investigate the confinement properties of two recently devised quasi-helically symmetric stellarator configurations, HSK and QSTK. Both have been optimized for large critical gradients of the ion temperature gradient mode, which is an important driver of turbulent transport in magnetic confinement fusion devices. To predict the resulting core plasma profiles, we utilize an advanced theoretical framework based on the gyrokinetic codes GENE and GENE-3D, coupled to the transport code TANGO. Compared to the HSX stellarator, both HSK and QSTK achieve significantly higher core-to-edge temperature ratios, partly thanks to their smaller aspect ratios, with the other part due to more detailed shaping of the magnetic geometry achieved during optimization. The computed confinement time, however, is less sensitive to core temperature than edge temperature, simply due to the disproportionate influence the edge has on stored plasma energy. We therefore emphasize the possible benefits of further optimizing turbulence in the outer core region, and the need to include accurate modelling of confinement in the edge region in order to assess overall plasma performance of turbulence optimized stellarators.

A fast and accurate turbulence transport model based on quasilinear gyrokinetics is developed. The model consists of a set of neural networks trained on a bespoke quasilinear GENE dataset, with a saturation rule calibrated to dedicated nonlinear simulations. The resultant neural network is approximately eight orders of magnitude faster than the original GENE quasilinear calculations. ITER predictions with the new model project a fusion gain in line with ITER targets. While the dataset is currently limited to the ITER baseline regime, this approach illustrates a pathway to develop reduced-order turbulence models both faster and more accurate than the current state-of-the-art.

Flux-tube gyrokinetic codes are widely used to simulate drift-wave turbulence in magnetic confinement devices. While a large number of studies show that flux-tube codes provide an excellent approximation for turbulent transport in medium-large devices, it still needs to be determined whether they are sufficient for modeling supra-thermal particle effects on core turbulence. This is called into question given the large temperature of energetic particles (EPs), which makes them hardly confined on a single flux-surface, but also due to the radially broad mode structure of energetic-particle-driven modes. The primary focus of this manuscript is to assess the range of validity of flux-tube codes in modeling fast ion effects by comparing radially global turbulence simulations with flux-tube results at different radial locations for realistic JET parameters using the gyrokinetic code GENE. To extend our study to a broad range of different plasma scenarios, this comparison is made for four different plasma regimes, which differ only by the profile of the ratio between the plasma kinetic and magnetic pressure. The latter is artificially rescaled to address the electrostatic limit and regimes with marginally stable, weakly unstable and strongly unstable fast ion modes. These energetic-particle-driven modes is identified as an AITG/KBAE via linear ORB5 and LIGKA simulations. It is found that the local flux-tube simulations can recover well the global results only in the electrostatic and marginally stable cases. When the AITG/KBAE becomes linearly unstable, the local approximation fails to correctly model the radially broad fast ion mode structure and the consequent global zonal patterns. According to this study, global turbulence simulations are likely required in regimes with linearly unstable AITG/KBAEs. In conditions with different fast ion-driven modes, these results might change.

The theoretical investigation of relevant turbulent transport mechanisms in H-mode pedestals is a great scientific and numerical challenge. In this study, we address this challenge by global, nonlinear gyrokinetic simulations of a full pedestal up to the separatrix, supported by a detailed characterisation of gyrokinetic instabilities from just inside the pedestal top to pedestal centre and foot. We present ASDEX Upgrade pedestal simulations using an upgraded version of the gyrokinetic, Eulerian, delta-f code GENE (genecode.org) that enables stable global simulations at experimental plasma beta values. The turbulent transport is found to exhibit a multi-channel, multi-scale character throughout the pedestal with the dominant contribution transitioning from ion scale Trapped Electron Modes (TEMs)/Micro Tearing Modes (MTMs) at the pedestal top to electron scale Electron Temperature Gradient modes (ETG) in the steep gradient region. Consequently, the turbulent electron heat flux changes from ion to electron scales and the ion heat flux reduces to almost neoclassic values in the pedestal centre. ExB shear is found to strongly reduce heat flux levels in all channels (electron, ion, electrostatic, electromagnetic) and the interplay of magnetic shear and pressure gradient is found to locally stabilise ion scale instabilities.

The recent progress in high-temperature superconductor technologies has led to the design and construction of SPARC, a compact tokamak device expected to reach plasma breakeven with up to $25$MW of external ion cyclotron resonant heating (ICRH) power. This manuscript presents local (flux-tube) and radially global gyrokinetic GENE (Jenko et al 2000 Phys. Plasmas {\bf 7} 1904) simulations for a reduced-field and current H-mode SPARC scenario showing that supra-thermal particles - generated via ICRH - strongly suppress ion-scale turbulent transport by triggering a fast ion-induced anomalous transport barrier (F-ATB). The trigger mechanism is identified as a wave-particle resonant interaction between the fast particle population and plasma micro-instabilities (Di Siena et al 2021 Phys. Rev. Lett. {\bf 125} 025002). By performing a series of global simulations employing different profiles for the thermal ions, we show that the fusion gain of this SPARC scenario could be substantially enhanced up to $\sim 80\%$ by exploiting this fast ion stabilizing mechanism. A study is also presented to further optimize the energetic particle profiles, thus possibly leading experimentally to an even more significant fusion gain.

In this work, a comparison of the global gyrokinetic codes EUTERPE and GENE-3D in stellarator configurations of LHD and W7-X is carried out. In linear simulations with adiabatic electrons, excellent agreement is found in the mode numbers, growth rate and frequency, mode structure, and spatial localization of the most unstable mode in LHD. In W7-X, the dependence of the growth rate and frequency with the mode number is well reproduced by both codes. The codes are also compared in linear simulations with kinetic ions and electrons in W7-X using model profiles, and reasonable agreement is found in the wavenumber of the most unstable modes. A stabilization of small-scale modes in kinetic-electron simulations with respect to the adiabatic-electron case is consistently found in both codes. Nonlinear simulations using adiabatic electrons and model profiles are also studied and the heat fluxes are compared. Very good agreement is found in the turbulent ion heat fluxes in both LHD and W7-X. Two problems that cannot be properly accounted for in local flux tube codes are studied: the localization of instabilities and turbulence over the flux surface and the influence of a background long-wavelength electric field. Good agreement between codes is found with respect to the spatial localization of instabilities and turbulence over the flux surface. The localization of saturated turbulence is found in both codes to be much smaller than that of the linear instabilities and smaller than previously reported in full-surface radially-local simulations. The influence of the electric field on the localization is also found to be smaller in the developed turbulent state than in the linear phase, and smaller than in previous works. A stabilizing effect of a constant electric field on the linearly unstable modes is found in both codes. A moderate reduction of turbulent transport by the radial electric field...

In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators -- while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources -- are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neoclassical transport code KNOSOS, and the 1D transport solver TANGO. This framework is used to analyze the recently observed degradation of energy confinement in electron-heated plasmas in the Wendelstein 7-X stellarator, where the central ion temperature was "clamped" to $T_i \approx 1.5$ keV regardless of the external heating power. By performing first-principles based simulations, we provide key evidence to understand this effect, namely the inefficient thermal coupling between electrons and ions in a turbulence-dominated regime, which is exacerbated by the large $T_e/T_i$ ratios, and show that a more efficient ion heat source, such as direct ion heating, will increase the on-axis ion temperature. This work paves the way towards the use of high-fidelity models for the development of the next generation of stellarators, in which neoclassical and turbulent transport are optimized simultaneously.

An accurate description of turbulence up to the transport time scale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport time scales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE-Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF-ASTRA. We show that GENE-Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we demonstrate that the ion temperature profile's peaking is due to electromagnetic effects of submarginal MHD instability. Based on these results, the expected GENE-Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport time scale, possibly making first-principles ITER simulations feasible on current computing resources.

A novel type of internal transport barrier (ITB) called F-ATB (fast ion-induced anomalous transport barrier) has been recently observed in state-of-the-art global gyrokinetic simulations on a properly optimized ASDEX Upgrade experiment [A. Di Siena et al. Phys. Rev. Lett. {\bf 127} 025002 (2021)]. Unlike the transport barriers previously reported in literature, the trigger mechanism for the F-ATB is a basically electrostatic wave-particle resonant interaction between supra-thermal particles - generated via ion cyclotron resonance heating (ICRH) - and ion scale plasma turbulence. This resonant effect strongly depends on the particular shape of the fast ion temperature and density profiles. Therefore, to further improve our theoretical understanding of this transport barrier, we present results exploring the parameter space and physical conditions for the F-ATB generation by performing a systematic study with global GENE simulations. Particular emphasis is given to the transport barrier width and its localization by scanning over different energetic particle temperature profiles. The latter are varied in amplitude, half-width, and radial localization of an ad-hoc Gaussian-like energetic particle logarithmic temperature gradient profile. For the reference parameters at hand, a threshold in the amplitude of the fast ion logarithmic temperature gradient is identified to trigger the transport barrier effectively.

Self-consistent full-size turbulent-transport simulations of the divertor and SOL of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted Ohmic L-mode discharges on the TCV tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from 5 diagnostic systems -- giving a total of 45 1- and 2-D comparison observables in two toroidal magnetic field directions. The dataset is used to validate three flux-driven 3D fluid-turbulence models: GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP -- both in terms of profile shape and absolute magnitude -- while a poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

Global gradient driven GENE gyrokinetic simulations are used to investigate TCV plasmas with negative triangularity. Considering a limited L-mode plasma, corresponding to an experimental triangularity scan, numerical results are able to reproduce the actual transport level over a major fraction of the plasma minor radius for a plasma with $\delta_{\rm LCFS}=-0.3$ and its equivalent with standard positive triangularity $\delta$. For the same heat flux, a larger electron temperature gradient is sustained by $\delta<0$, in turn resulting in an improved electron energy confinement. Consistently with the experiments, a reduction of the electron density fluctuations is also seen. Local flux-tube simulations are used to gauge the magnitude of nonlocal effects. Surprisingly, very little differences are found between local and global approaches for $\delta>0$, while local results yield a strong overestimation of the heat fluxes when $\delta<0$. Despite the high sensitivity of the turbulence level with respect to the input parameters, global effects appear to play a crucial role in the negative triangularity plasma and must be retained to reconcile simulations and experiments. Finally, a general stabilizing effect of negative triangularity, reducing fluxes and fluctuations by a factor dependent on the actual profiles, is recovered.

The formation of a substantial post-disruption runaway electron current in ASDEX Upgrade material injection experiments is determined by avalanche multiplication of a small seed population of runaway electrons. For the investigation of these scenarios, the runaway electron description of the coupled 1.5D transport solvers ASTRA-STRAHL is amended by a fluid-model describing electron runaway caused by the hot-tail mechanism. Applied in simulations of combined background plasma evolution, material injection, and runaway electron generation in ASDEX Upgrade discharge #33108, both the Dreicer and hot-tail mechanism for electron runaway produce only $\sim$ 3$~$kA of runaway current. In colder plasmas with core electron temperatures $T_\mathrm{e,c}$ below 9$~$keV, the post-disruption runaway current is predicted to be insensitive to the initial temperature, in agreement with experimental observations. Yet in hotter plasmas with $T_\mathrm{e,c} > 10~\mathrm{keV}$, hot-tail runaway can be increased by up to an order of magnitude, contributing considerably to the total post-disruption runaway current. In ASDEX Upgrade high temperature runaway experiments, however, no runaway current is observed at the end of the disruption, despite favourable conditions for both primary and secondary runaway., Comment: 26 pages, 11 figures; published in the Journal of Plasma Physics

The key result of the present work is the theoretical prediction and observation of the formation of a new type of transport barrier in fusion plasmas, called F-ATB (fast ion-induced anomalous transport barrier). As demonstrated through state-of-the-art global electrostatic and electromagnetic simulations, the F-ATB is characterized by a full suppression of the turbulent transport - caused by strongly sheared, axisymmetric $E \times B$ flows - and an increase of the neoclassical counterpart, albeit keeping the overall fluxes at significantly reduced levels. The trigger mechanism is shown to be a mainly electrostatic resonant interaction between supra-thermal particles, generated via ion-cyclotron-resonance heating, and plasma micro-turbulence. These findings are obtained by realistic simulations of the ASDEX Upgrade discharge $\#36637$ - properly designed to maximized the beneficial role of the wave-particle resonance interaction - which exhibits the expected properties of improved confinement produced by energetic particles., Comment: Version 2

The nonlinear dynamics of beta-induced Alfv\'en Eigenmodes (BAE) driven by energetic particles (EP) in the presence of ion-temperature-gradient (ITG) turbulence is investigated, by means of selfconsistent global gyrokinetic simulations and analytical theory. A tokamak magnetic equilibrium with large aspect ratio and reversed shear is considered. A previous study of this configuration has shown that the electron species plays an important role in determining the nonlinear saturation level of a BAE in the absence of turbulence [A. Biancalani, et al., J. Plasma Phys. (2020)]. Here, we extend the study to a turbulent plasma. The EPs are found modify the heat fluxes by introducing energy at the large spatial scales, mainly at the toroidal mode number of the dominant BAE and its harmonics. In this regime, BAEs are found to carry a strong electron heat flux. The feed-back of the global relaxation of the temperature profiles induced by the BAE, and on the turbulence dynamics, is also discussed.

We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data, to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands with corresponding toroidal mode numbers n=4 and 8. Such disparate fluctuation scales, with substantial gaps between toroidal mode numbers, are commonly observed in pedestal fluctuations. Here we provide a clear explanation, namely the alignment of the relevant rational surfaces (and not others) with the peak in the omega star profile, which is localized in the steep gradient region of the pedestal. We demonstrate that a global treatment is required to capture this effect. Nonlinear simulations suggest that the MTM fluctuations produce experimentally-relevant transport levels and saturate by relaxing the background electron temperature gradient, slightly downshifting the fluctuation frequencies from the linear predictions. Scans in collisionality are compared with simple MTM dispersion relations. At the experimental points considered, MTM growth rates can either increase or decrease with collision frequency depending on the parameters thus defying any simple characterization of collisionality dependence.

We study the effect of turbulent transport in different magnetic configurations of the Weldenstein 7-X stellarator. In particular, we performed direct numerical simulations with the global gyrokinetic code GENE-3D, modeling the behavior of Ion Temperature Gradient turbulence in the Standard, High-Mirror, and Low-Mirror configurations of W7-X. We found that the Low-Mirror configuration produces more transport than both the High-Mirror and the Standard configurations. By comparison with radially local simulations, we have demonstrated the importance of performing global nonlinear simulations to predict the turbulent fluxes quantitatively.

We present the first successful simulation of a induced disruption in ASDEX Upgrade from massive material injection (MMI) up to established runaway electron (RE) beam, thus covering pre-thermal quench, thermal quench and current quench (CQ) of the discharge. For future high-current fusion devices such as ITER, the successful suppression of REs through MMI is of critical importance to ensure the structural integrity of the vessel. To computationally study the interplay between MMI, background plasma response, and RE generation, a toolkit based on the 1.5D transport code coupling ASTRA-STRAHL is developed. Electron runaway is described by state-of-the-art reduced kinetic models in the presence of partially ionized impurities. Applied to argon MMI in ASDEX Upgrade discharge #33108, key plasma parameters measured experimentally, such as temporal evolution of the line averaged electron density, plasma current decay rate and post-CQ RE current, are well reproduced by the simulation presented. Impurity ions are transported into the central plasma by the combined effect of neoclassical processes and additional effects prescribed inside the $q = 2$ rational surface to explain experimental time scales. Thus, a thermal collapse is induced through strong impurity radiation, giving rise to a substantial RE population as observed experimentally., Comment: 15 pages, 11 figures; submitted to Nuclear Fusion

The global and electromagnetic gyrokinetic particle-in-cell code ORB5 is employed to investigate the self-consistent interactions between Alfv\'en modes (AM) and ion temperature gradient (ITG) turbulence in a magnetically confined plasma. Here, an axisymmetric magnetic equilibrium with reversed shear and large aspect ratio is considered. An AM with toroidal mode number n=5 is driven unstable by introducing a population of suprathermal ions. Once the AM saturates in the presence of the fully developed turbulence, the ion heat flux is dominated by the AM and its main harmonics. ITG-induced transport is found to also be enhanced in the presence of the unstable AM.

During a tokamak discharge, the plasma can vary between different confinement regimes: Low (L), High (H) and, in some cases, a temporary (intermediate state), called Dithering (D). In addition, while the plasma is in H mode, Edge Localized Modes (ELMs) can occur. The automatic detection of changes between these states, and of ELMs, is important for tokamak operation. Motivated by this, and by recent developments in Deep Learning (DL), we developed and compared two methods for automatic detection of the occurrence of L-D-H transitions and ELMs, applied on data from the TCV tokamak. These methods consist in a Convolutional Neural Network (CNN) and a Convolutional Long Short Term Memory Neural Network (Conv-LSTM). We measured our results with regards to ELMs using ROC curves and Youden's score index, and regarding state detection using Cohen's Kappa Index.

Turbulent dynamics in the scrape-off layer (SOL) of magnetic fusion devices is intermittent with large fluctuations in density and pressure. Therefore, a model is required that allows perturbations of similar or even larger magnitude to the time-averaged background value. The fluid-turbulence code GRILLIX is extended to such a global model, which consistently accounts for large variation in plasma parameters. Derived from the drift reduced Braginskii equations, the new GRILLIX model includes electromagnetic and electron-thermal dynamics, retains global parametric dependencies and the Boussinesq approximation is not applied. The penalisation technique is combined with the flux-coordinate independent (FCI) approach [F. Hariri and M. Ottaviani, Comput.Phys.Commun. 184:2419, (2013); A. Stegmeir et al., Comput.Phys.Commun. 198:139, (2016)], which allows to study realistic diverted geometries with X-point(s) and general boundary contours. We characterise results from turbulence simulations and investigate the effect of geometry by comparing simulations in circular geometry with toroidal limiter against realistic diverted geometry at otherwise comparable parameters. Turbulence is found to be intermittent with relative fluctuation levels of up to 40% showing that a global description is indeed important. At the same time via direct comparison, we find that the Boussinesq approximation has only a small quantitative impact in a turbulent environment. In comparison to circular geometry the fluctuations are reduced in diverted geometry, which is related to a different zonal flow structure. Moreover, the fluctuation level has a more complex spatial distribution in diverted geometry. Due to local magnetic shear, which differs fundamentally in circular and diverted geometry, turbulent structures become strongly distorted in the perpendicular direction and are eventually damped away towards the X-point.

In recent years, a strong reduction of plasma turbulence in the presence of energetic particles has been reported in a number of magnetic confinement experiments and corresponding gyrokinetic simulations. While highly relevant to performance predictions for burning plasmas, an explanation for this primarily nonlinear effect has remained elusive so far. A thorough analysis finds that linearly marginally stable energetic particle driven modes are excited nonlinearly, depleting the energy content of the turbulence and acting as an additional catalyst for energy transfer to zonal modes (the dominant turbulence saturation channel). Respective signatures are found in a number of simulations for different JET and ASDEX Upgrade discharges with reduced transport levels attributed to energetic ion effects.

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position $\rho=0.80$, simulations with the ion temperature gradient increased by $40\%$ above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty $\sigma$ of the Charge Exchange Recombination Spectroscopy (CER) measurements at the $1.6 \sigma$ level. Multi-scale simulations are carried out with realistic mass ratio and geometry for the first time in the near-edge. These multi-scale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position $\rho=0.90$, nonlinear simulations show a hybrid state of ITG and trapped electron modes~(TEMs), which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of $\mathbf{E} \times \mathbf{B}$ shear effects and an increase in the electron temperature gradient by $\sim 23\%$. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the $1.3 \sigma$ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region., Comment: 14 pages, 17 figures, published in Physics of Plasmas

The question of the relative importance of coherent structures and waves has for a long time attracted a great deal of interest in astrophysical plasma turbulence research, with a more recent focus on kinetic scale dynamics. Here we utilize high-resolution observational and simulation data to investigate the nature of waves and structures emerging in a weakly collisional, turbulent kinetic plasma. Observational results are based on in situ solar wind measurements from the Cluster and MMS spacecraft, and the simulation results are obtained from an externally driven, three-dimensional fully kinetic simulation. Using a set of novel diagnostic measures we show that both the large-amplitude structures and the lower-amplitude background fluctuations preserve linear features of kinetic Alfven waves to order unity. This quantitative evidence suggests that the kinetic turbulence cannot be described as a mixture of mutually exclusive waves and structures but may instead be pictured as an ensemble of localized, anisotropic wave packets or "eddies" of varying amplitudes, which preserve certain linear wave properties during their nonlinear evolution., Comment: Revised and extended version; accepted for publication in Phys. Rev. X

We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation $k_{\parallel}(k_{\perp})$, where $k_{\parallel}$ is a characteristic wave number along the local mean magnetic field at perpendicular scale $l_{\perp}\sim 1/k_{\perp}$. The subion range anisotropy is scale dependent with $k_{\parallel}

A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-law energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ~-2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of -5/3 in the MHD range and of ~-3 at sub-ion scales., Comment: 6 pages, 5 figures

The observation of distinct peaks in tokamak core reflectometry measurements - named quasi-coherent-modes (QCMs) - are identified as a signature of Trapped-Electron-Mode (TEM) turbulence [H. Arnichand et al. 2016 Plasma Phys. Control. Fusion 58 014037]. This phenomenon is investigated with detailed linear and nonlinear gyrokinetic simulations using the \gene code. A Tore-Supra density scan is studied, which traverses through a Linear (LOC) to Saturated (SOC) Ohmic Confinement transition. The LOC and SOC phases are both simulated separately. In the LOC phase, where QCMs are observed, TEMs are robustly predicted unstable in linear studies. In the later SOC phase, where QCMs are no longer observed, ITG modes are identified. In nonlinear simulations, in the ITG (SOC) phase, a broadband spectrum is seen. In the TEM (LOC) phase, a clear emergence of a peak at the TEM frequencies is seen. This is due to reduced nonlinear frequency broadening of the underlying linear modes in the TEM regime compared with the ITG regime. A synthetic diagnostic of the nonlinearly simulated frequency spectra reproduces the features observed in the reflectometry measurements. These results support the identification of core QCMs as an experimental marker for TEM turbulence

We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully-kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta ($\beta_i$) regimes are considered: 0.1 and 0.5. For $\beta_i=0.5$, good agreement between the GK and FK models is found at scales ranging from the ion to the electron gyroradius, thus providing firm evidence for a kinetic Alfv\'en cascade scenario. In the same range, the HK model produces shallower spectral slopes, presumably due to the lack of electron Landau damping. For $\beta_i=0.1$, a detailed analysis of spectral ratios reveals a slight disagreement between the GK and FK descriptions at kinetic scales, even though kinetic Alfv\'en fluctuations likely still play a significant role. The discrepancy can be traced back to scales above the ion gyroradius, where the FK and HK results seem to suggest the presence of fast magnetosonic and ion Bernstein modes in both plasma beta regimes, but with a more notable deviation from GK in the low-beta case. The identified practical limits and strengths of reduced-kinetic approximations, compared here against the fully-kinetic model on a case-by-case basis, may provide valuable insight into the main kinetic effects at play in turbulent collisionless plasmas, such as the solar wind.

We present two-dimensional hybrid kinetic/magnetohydrodynamic simulations of planned laser-ablation experiments in the Large Plasma Device (LAPD). Our results, based on parameters which have been validated in previous experiments, show that a parallel collisionless shock can begin forming within the available space. Carbon-debris ions that stream along the magnetic-field direction with a blow-off speed of four times the Alfven velocity excite strong magnetic fluctuations, eventually transfering part of their kinetic energy to the surrounding hydrogen ions. This acceleration and compression of the background plasma creates a shock front, which satisfies the Rankine-Hugoniot conditions and can therefore propagate on its own. Furthermore, we analyze the upstream turbulence and show that it is dominated by the right-hand resonant instability., Comment: 13 pages, 17 figures. Accepted at Physics of Plasmas

We analyze plasma heating in weakly collisional kinetic Alfv\'en wave (KAW) turbulence using high resolution gyrokinetic simulations spanning the range of scales between the ion and the electron gyroradii. Real space structures that have a higher than average heating rate are shown not to be confined to current sheets. This novel result is at odds with previous studies, which use the electromagnetic work in the local electron fluid frame, i.e. $\mathbf{J} \!\cdot\! (\mathbf{E} + \mathbf{v}_e\times\mathbf{B})$, as a proxy for turbulent dissipation to argue that heating follows the intermittent spatial structure of the electric current. Furthermore, we show that electrons are dominated by parallel heating while the ions prefer the perpendicular heating route. We comment on the implications of the results presented here., Comment: 5 pages, 3 figures

The efficient production of electricity from nuclear fusion in magnetically confined plasmas relies on a good confinement of the thermal energy. For more than thirty years, the observation that such confinement depends on the mass of the plasma isotope and its interaction with apparently unrelated plasma conditions has remained largely unexplained and it has become one of the main unsolved issues. By means of numerical studies based on the gyrokinetic theory, we quantitatively show how the plasma microturbulence depends on the isotope mass through nonlinear multiscale microturbulence effects involving the interplay between zonal flows, electromagnetic effects and the torque applied. This finding has crucial consequences for the design of future reactors since, in spite of the fact that they will be composed by multiple ion species, their extrapolation from present day experiments heavily relies on the knowledge obtained from a long experimental tradition based in single isotope plasmas.

A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfv\'en wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at $\beta\gtrsim1$ and by magnetosonic/whistler fluctuations at lower $\beta$. The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma $\beta$, and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions., Comment: 6 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

Nonlinear energy transfer and dissipation in Alfv\'en wave turbulence are analyzed in the first gyrokinetic simulation spanning all scales from the tail of the MHD range to the electron gyroradius scale. For typical solar wind parameters at 1 AU, about 30% of the nonlinear energy transfer close to the electron gyroradius scale is mediated by modes in the tail of the MHD cascade. Collisional dissipation occurs across the entire kinetic range $k_\perp\rho_i\gtrsim 1$. Both mechanisms thus act on multiple coupled scales, which have to be retained for a comprehensive picture of the dissipation range in Alfv\'enic turbulence., Comment: Made several improvements to figures and text suggested by referees

In this work, we compare gyrokinetic (GK) and fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields ($b_g$). Here we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma ($\beta_i=0.01$), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field ($b_g\gtrsim 30$) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive ($b_g\gtrsim 5$). Kinetic PIC simulations using guide fields $b_g \lesssim 30$ reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma ($\beta_i=1.0$) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields ($b_g\lesssim 3$)., Comment: 18 pages, 13 figures. Revised to match with the published version in Physics of Plasmas

Gyrokinetic validation studies are crucial in developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

Observations and numerical simulations of laboratory and space plasmas in almost collisionless regimes reveal anisotropic and non-gyrotropic particle distribution functions. We investigate how such states can persist in the presence of a sheared flow. We focus our attention on the pressure tensor equation in a magnetized plasma and derive analytical self-consistent plasma equilibria which exhibit a novel asymmetry with respect to the magnetic field direction. These results are relevant for investigating, within fluid models that retain the full pressure tensor dynamics, plasma configurations where a background shear flow is present., Comment: 13 pages, 7 figures

The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$\beta$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive linear and nonlinear gyrokinetic simulations are then carried out for the high-$\beta$ JET hybrid discharge 75225, at two separate locations at inner and outer radii. It is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment. The enhancement of this effect by suprathermal pressure also remains significant. It is also found that flow shear stabilization is not effective at the inner radii. However, at outer radii the situation is reversed. Electromagnetic stabilization is negligible while the flow shear stabilization is significant. These results constitute the high-$\beta$ generalization of comparable observations found at low-$\beta$ at JET. This is encouraging for the extrapolation of electromagnetic ITG stabilization to future devices. An estimation of the impact of this effect on the ITER hybrid scenario leads to a 20% fusion power improvement., Comment: 10 pages, 13 figures. Paper coupled to invited talk at the 41st EPS conference, Berlin, 2014

Nonlinear gyrokinetics provides a suitable framework to describe short-wavelength turbulence in magnetized laboratory and astrophysical plasmas. In the electrostatic limit, this system is known to exhibit a free energy cascade towards small scales in (perpendicular) real and/or velocity space. The dissipation of free energy is always due to collisions (no matter how weak the collisionality), but may be spread out across a wide range of scales. Here, we focus on freely-decaying 2D electrostatic turbulence on sub-ion-gyroradius scales. An existing scaling theory for the turbulent cascade in the weakly collisional limit is generalized to the moderately collisional regime. In this context, non-universal power law scalings due to multiscale dissipation are predicted, and this prediction is confirmed by means of direct numerical simulations., Comment: 7 pages, 5 figures, accepted for publication in Physics of Plasmas

The nonlinear gyrokinetic equations describe plasma turbulence in laboratory and astrophysical plasmas. To solve these equations, massively parallel codes have been developed and run on present-day supercomputers. This paper describes measures to improve the efficiency of such computations, thereby making them more realistic. Explicit Runge-Kutta schemes are considered to be well suited for time-stepping. Although the numerical algorithms are often highly optimized, performance can still be improved by a suitable choice of the time-stepping scheme, based on spectral analysis of the underlying operator. Here, an operator splitting technique is introduced to combine first-order Runge-Kutta-Chebychev schemes for the collision term with fourth-order schemes for the remaining terms. In the nonlinear regime, based on the observation of eigenvalue shifts due to the (generalized) $E\times B$ advection term, an accurate and robust estimate for the nonlinear timestep is developed. The presented techniques can reduce simulation times by factors of up to three in realistic cases. This substantial speedup encourages the use of similar timestep optimized explicit schemes not only for the gyrokinetic equation, but also for other applications with comparable properties., Comment: 11 pages, 5 figures, accepted for publication in Computer Physics Communications

We propose a quantitative model of ion temperature gradient driven turbulence in toroidal magnetized plasmas. In this model, the turbulence is regulated by zonal flows, i.e. mode saturation occurs by a zonal-flow-mediated energy cascade ("shearing"), and zonal flow amplitude is controlled by nonlinear decay. Our model is tested in detail against numerical simulations to confirm that both its assumptions and predictions are satisfied. Key results include (1) a sensitivity of the nonlinear zonal flow response to the energy content of the linear instability, (2) a persistence of zonal-flow-regulated saturation at high temperature gradients, (3) a physical explanation of the nonlinear saturation process in terms of secondary and tertiary instabilities, and (4) dependence of heat flux in terms of dimensionless parameters., Comment: Final journal version. Some clarifications and a new Fig. 4

The Large Eddy Simulation (LES) approach - solving numerically the large scales of a turbulent system and accounting for the small-scale influence through a model - is applied to nonlinear gyrokinetic systems that are driven by a number of different microinstabilities. Comparisons between modeled, lower resolution, and higher resolution simulations are performed for an experimental measurable quantity, the electron density fluctuation spectrum. Moreover, the validation and applicability of LES is demonstrated through a series of diagnostics based on the free energetics of the system., Comment: 14 pages, 9 figures