Gyrokinetic Large Eddy Simulations

The Large Eddy Simulation (LES) approach is adapted to the study of plasma microturbulence in a fully three-dimensional gyrokinetic system. Ion temperature gradient driven turbulence is studied with the {\sc GENE} code for both a standard resolution and a reduced resolution with a model for the sub-grid scale turbulence. A simple dissipative model for representing the effect of the sub-grid scales on the resolved scales is proposed and tested. Once calibrated, the model appears to be able to reproduce most of the features of the free energy spectra for various values of the ion temperature gradient

Free energy cascade in gyrokinetic turbulence

In gyrokinetic theory, the quadratic nonlinearity is known to play an important role in the dynamics by redistributing (in a conservative fashion) the free energy between the various active scales. In the present study, the free energy transfer is analyzed for the case of ion temperature gradient driven turbulence. It is shown that it shares many properties with the energy transfer in fluid turbulence. In particular, one finds a forward (from large to small scales), extremely local, and self-similar cascade of free energy in the plane perpendicular to the background magnetic field. These findings shed light on some fundamental properties of plasma turbulence, and encourage the development of large eddy simulation techniques for gyrokinetics.Comment: 4 pages, 2 Postscript figure

Corrigendum: Electromagnetic turbulence suppression by energetic particle driven modes (2019 Nucl. Fusion 59 124001), Nuclear Fusion 60, 089501 (2020)

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

Benchmark of quasi-linear models against gyrokinetic single scale simulations in deuterium and tritium plasmas for a JET high beta hybrid discharge

A benchmark of the reduced quasi-linear models QuaLiKiz and TGLF with GENE gyrokinetic simulations has been performed for parameters corresponding to a JET high performance hybrid pulse in deuterium. Given the importance of the study of such advanced scenarios in view of ITER and DEMO operations, the dependence of the transport on the ion isotope mass has also been assessed, by repeating the benchmark changing the ion isotope to tritium. TGLF agrees better with GENE on the linear spectra and the flux levels. However, concerning the isotope dependence, only QuaLiKiz reproduces the GENE radial trend of a basically gyro-Bohm (gB) scaling at inner radii and instead anti-gB at outer radii. The physics effects which are responsible of the antigB effect in GENE simulations have been singled out

A flux-matched gyrokinetic analysis of DIII-D L-mode turbulence

Previous nonlinear gyrokinetic simulations of specific DIII-D L-mode cases have been found to significantly underpredict the ion heat transport and associated density and temperature fluctuation levels by up to almost one of order of magnitude in the outer-core domain, i.e., roughly in the last third of the minor radius. Since then, this so-called shortfall issue has been subject to various speculations on possible reasons and furthermore motivation for a number of dedicated comparisons for L-mode plasmas in comparable machines. However, only a rather limited number of simulations and gyrokinetic codes has been applied to the original scenario, thus calling for further dedicated investigations in order to broaden the scientific basis. The present work contributes along these lines by employing another well-established gyrokinetic code in a numerically and physically comprehensive manner. Contrary to the previous studies, only a mild underprediction is observed at the outer radial positions which can furthermore be overcome by varying the ion temperature gradient within the error bars associated with the experimental measurement. The significance and reliability of these simulations are demonstrated by benchmarks, numerical convergence tests, and furthermore by extensive validation studies. The latter involve cross-phase and cross-power spectra analyses of various fluctuating quantities and confirm a high degree of realism. The code discrepancies come as a surprise since the involved software packages had been benchmarked repeatedly and very successfully in the past. Further collaborative effort in identifying the underlying difference is hence required.European Union. Horizon 2020 Research and Innovation Programme (Grant 633053