Search

In the context of global gyrokinetic simulations of turbulence using a Particle-In-Cell framework, verifying the delta-f assumption with a fixed background distribution becomes challenging when determining quasi-steady state profiles corresponding to given sources over long time scales, where plasma profiles can evolve significantly. The advantage of low relative sampling noise afforded by the delta-f scheme is shown to be retained by considering the background as a time-evolving Maxwellian with time-dependent density and temperature profiles. Implementation of this adaptive scheme to simulate electrostatic collisionless flux-driven turbulence in tokamak plasmas show small and non-increasing sampling noise levels, which would otherwise increase indefinitely with a stationary background scheme. The adaptive scheme furthermore allows one to reach numerically converged results of quasi-steady state with much lower marker numbers.

In this work, we explore the triangularity effects on turbulent transport employing global gyrokinetic simulations performed with the ORB5 code. Numerous experiments on the Tokamak \`a Configuration Variable (TCV) and, more recently, on the DIII-D machine, have demonstrated superior confinement properties in L-mode of negative triangularity (NT) over positive triangularity (PT) configuration. This presents a particularly attractive scenario, as L-mode operation eliminates or significantly mitigates the presence of hazardous edge-localized modes (ELMs). However, a full theoretical understanding of all these observations remains elusive. Specifically, questions remain about how NT improvements can extend to the core where triangularity is very small, and whether these improvements can scale to larger devices. This paper addresses these two questions. Our analysis is divided into two parts: we first demonstrate that the confinement enhancement in NT configurations arises from the interdependent edge-core dynamics, and then we present the results of a system size scan. Crucially, we find that the relative turbulent transport reduction of NT over PT appears not to be contingent on machine dimensions or fluctuation scales and is moreover robust with respect to variations in plasma profiles. This insight underscores the fundamental nature of the NT confinement advantage and paves the way for its potential application in future fusion devices, regardless of their size.

In this work, we use the global electromagnetic and electrostatic gyro kinetic approaches to investigate the effects of zonal flows forced-driven by Alfv\'en modes due to their excitation by energetic particles (EPs), on the dynamics of ITG (Ion temperature gradient) instabilities. The equilibrium of the 92416 JET tokamak shot is considered. The linear and nonlinear Alfv\'en modes dynamics, as well as the zonal flow dynamics, are investigated and their respective radial structures and saturation levels are reported. ITG dynamics in the presence of the zonal flows excited by these Alfv\'en modes are also investigated. We find that, the zonal flows forced-driven by Alfv\'en modes can significantly impact the ITG dynamics. A zonal flow amplitude scan reveals the existence of an inverse relation between the zonal flow amplitude and the ITG growth rate. These results show that, forced-driven zonal flows can be an important indirect part of turbulence mitigation due to the injection of energetic particles., Comment: 13 pages, 10 figures

Analysis of Lagrangian Coherent Structures (LCSs) has been showed to be a valid mathematical approach to explain the formation of transport barriers in magnetized plasmas. Such LCSs, borrowed from fluid dynamics theory, can be considered as the hidden skeleton of the system and can be used for studying a wide spectrum of transport mechanisms even in plasmas. In this paper, we demonstrate that such structures can be particularly useful for underlying the hidden paths governing the motion of magnetic field lines in chaotic magnetic fields. To perform such an analysis, we developed a numerical tool able to detect LCSs for a general dynamical system. The tool is able to deal with general coordinate systems and it is shown to match with other techniques already used to analyse chaotic magnetic fields, e.g. connection length. After the description of the computational tool, we focus on the heat transport equation and the comparison between temperature profile and topology of the LCSs. We provide evidence that numerical simulations are able to reproduce the temperature profiles similar to those observed in RFP experiments and that our tool successfully predicts the location of temperature gradients. The results suggest that, inside the chaotic region, the field-lines motion is far from stochastic and that the presence of hidden patterns allows the development of high temperature gradients.