G. Merlo

An arbitrary wavelength solver for global gyrokinetic simulations. Application to the study of fine radial structures on microturbulence due to non-adiabatic passing electron dynamics

The influence of the fine layers of the non-adiabatic passing electron response on electrostatic turbulent transport, previously studied systematically in flux tube geometry [Dominski et al., Phys. Plasmas 22, 062303 (2015)], is pursued in global geometry in conditions relevant for the TCV tokamak with a deuterium plasma (mi/me = 3672). The spectral organization of the passing electron turbulent flux and its dependence on the radial profile of the safety factor are revealed. A radially dependent toroidal spectral analysis of the turbulent fluxes led to the key result that the particle and heat diffusivities of passing-electrons are proportional to the local density of low-order mode rational surfaces. To permit this study of the short radial scales associated with the passing electron dynamics, a new field solver valid at an arbitrary wavelength is implemented in ORB5, for the gyrokinetic quasi-neutrality equation. A benchmark is conducted against the global version of the gyrokinetic code GENE, showing v...

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant ă confinement improvement in plasmas with negative triangularity (delta < ă 0) compared to those with standard positive triangularity has been ă interpreted in terms of different degrees of profile stiffness (Sauter ă et al 2014 Phys. Plasmas 21 055906) and/or different critical gradients. ă Employing the Eulerian gyrokinetic code GENE (Jenko et al 2000 Phys. ă Plasmas 7 1904), profile stiffness and critical gradients are studied ă under TCV relevant conditions. For the considered experimental ă discharges, trapped electron modes (TEMs) and electron temperature ă gradient (ETG) modes are the dominant microinstabilities, with the ă latter providing a significant contribution to the non-linear electron ă heat fluxes near the plasma edge. Two series of simulations with ă different levels of realism are performed, addressing the question of ă profile stiffness at various radial locations. Retaining finite ă collisionality, impurities and electromagnetic effects, as well as the ă physical electron-to-ion mass ratio are all necessary in order to ă approach the experimental flux measurements. However, flux-tube ă simulations are unable to fully reproduce the TCV results, pointing ă towards the need to carry out radially nonlocal (global) simulations, ă i.e. retaining finite machine size effects, in a future study. Some ă conclusions about the effect of triangularity can nevertheless be drawn ă based on the flux-tube results. In particular, the importance of ă considering the sensitivity to both temperature and density gradient is ă shown. The flux tube results show an increase of the critical gradients ă towards the edge, further enhanced when d < 0, and they also appear to ă indicate a reduction of profile stiffness towards plasma edge.

Linear multispecies gyrokinetic flux tube benchmarks in shaped tokamak plasmas

Verification is the fundamental step that any turbulence simulation code has to be submitted in order to assess the proper implementation of the underlying equations. We have carried out a cross comparison of three flux tube gyrokinetic codes, GENE [F. Jenko et al., Phys. Plasmas 7, 1904 (2000)], GKW [A. G. Peeters et al., Comput. Phys. Commun. 180, 2650 (2009)], and GS2 [W. Dorland et al., Phys. Rev. Lett. 85, 5579 (2000)], focusing our attention on the effect of realistic geometries described by a series of MHD equilibria with increasing shaping complexity. To simplify the effort, the benchmark has been limited to the electrostatic collisionless linear behaviour of the system. A fully gyrokinetic model has been used to describe the dynamics of both ions and electrons. Several tests have been carried out looking at linear stability at ion and electron scales, where for the assumed profiles Ion Temperature Gradient (ITG)/Trapped Electron Modes and Electron Temperature Gradient modes are unstable. The capa...

Bringing global gyrokinetic turbulence simulations to the transport timescale using a multiscale approach

The vast separation dividing the characteristic times of energy confinement and turbulence in the core of toroidal plasmas makes first-principles prediction on long timescales extremely challenging. Here we report the demonstration of a multiple-timescale method that enables coupling global gyrokinetic simulations with a transport solver to calculate the evolution of the self-consistent temperature profile. This method, which exhibits resiliency to the intrinsic fluctuations arising in turbulence simulations, holds potential for integrating nonlocal gyrokinetic turbulence simulations into predictive, whole-device models.

Characterization with microturbulence simulations of the zero particle flux condition in case of a TCV discharge showing toroidal rotation reversal

In view of the stabilization effect of sheared plasma rotation on microturbulence, it is important to study the intrinsic rotation that develops in tokamaks that present negligible external toroidal torque, like ITER. Remarkable observations have been made on TCV, analysing discharges without NBI injection, as reported in [A. Bortolon et al. 2006 Phys. Rev. Lett. 97] and exhibiting a rotation inversion occurring in conjunction with a relatively small change in the plasma density. We focus in particular on a limited L-mode TCV shot published in [B. P. Duval et al. 2008 Phys. Plasmas 15], that shows a rotation reversal during a density ramp up. In view of performing a momentum transport analysis on this TCV shot, some constraints have to be considered to reduce the uncertainty on the experimental parameters. One useful constraint is the zero particle flux condition, resulting from the absence of direct particle fuelling to the plasma core. In this work, a preliminary study of the reconstruction of the zero particle flux hyper-surface in the physical parameters space is presented, taking into account the effect of the main impurity (carbon) and beginning to explore the effect of collisions, in order to find a subset of this hyper-surface within the experimental error bars. The analysis is done performing gyrokinetic simulations with the local (flux-tube) version of the Eulerian code GENE [Jenko et al 2000 Phys. Plasmas 7 1904], computing the fluxes with a Quasi-Linear model, according to [E. Fable et al. 2010 PPCF 52], and validating the QL results with Non-Linear simulations in a subset of cases.

Identifying microturbulence regimes in a TCV discharge making use of physical constraints on particle and heat fluxes

Reducing the uncertainty on physical input parameters derived from experimental measurements is essential towards improving the reliability of gyrokinetic turbulence simulations. This can be achieved by introducing physical constraints. Amongst them, the zero particle flux condition is considered here. A first attempt is also made to match as well the experimental ion/electron heat flux ratio. This procedure is applied to the analysis of a particular Tokamak a Configuration Variable discharge. A detailed reconstruction of the zero particle flux hyper-surface in the multi-dimensional physical parameter space at fixed time of the discharge is presented, including the effect of carbon as the main impurity. Both collisionless and collisional regimes are considered. Hyper-surface points within the experimental error bars are found. The analysis is done performing gyrokinetic simulations with the local version of the GENE code, computing the fluxes with a Quasi-Linear (QL) model and validating the QL results with non-linear simulations in a subset of cases.Reducing the uncertainty on physical input parameters derived from experimental measurements is essential towards improving the reliability of gyrokinetic turbulence simulations. This can be achieved by introducing physical constraints. Amongst them, the zero particle flux condition is considered here. A first attempt is also made to match as well the experimental ion/electron heat flux ratio. This procedure is applied to the analysis of a particular Tokamak a Configuration Variable discharge. A detailed reconstruction of the zero particle flux hyper-surface in the multi-dimensional physical parameter space at fixed time of the discharge is presented, including the effect of carbon as the main impurity. Both collisionless and collisional regimes are considered. Hyper-surface points within the experimental error bars are found. The analysis is done performing gyrokinetic simulations with the local version of the GENE code, computing the fluxes with a Quasi-Linear (QL) model and validating the QL results wi...