Modelling of turbulent fluctuations measured in the TCV tokamak with gyrokinetic simulations and a synthetic phase contrast imaging diagnostic

Proper comparison of simulated turbulent fluctuations with experimental observations is neces- sary for code validation, and for establishing a link between experimental fluctuation measure- ments and the underlying physics. In this contribution we report on the progress of modelling localised measurements of turbulent electron density fluctuations, obtained with the tangential phase contrast imaging (TPCI) diagnostic installed on the TCV tokamak. The modelling is done in two steps. First, nonlinear flux-tube gyrokinetic simulations are performed with the Eulerian (grid-based) GENE code, taking into account realistic TCV geometry and profiles. A synthetic diagnostic that has been developed in MATLAB is then applied to the simulated density fluctua- tions to model the experimental measurement procedure, allowing detailed comparison between simulation and experiment. Using these tools, we extended a previous computational analysis of the role of plasma triangu- larity in improved confinement in TCV. In previous such studies, the focus was on comparing heat flux levels between simulation and experiments. Large discrepancies were found, which were attributed to the large experimental error bars on the input profile gradients, to the sim- plification of not including electron scales in the simulations, to possible global effects, or to a combination of all three. With the synthetic diagnostic a more in depth comparison is performed, investigating the importance of such effects also in the observed wavenumber and frequency spectrum. After motivating and verifying the described procedure for modelling density fluctuations at TCV, we make initial predictions on the role of electron scale turbulence on transport levels in highly shaped plasmas. This is performed in preparation of first electron scale measurements that are envisioned to be performed with TPCI, as a result of an upgrade currently underway.

Neutral atom dynamics and plasma turbulence in the tokamak periphery

Christoph Wersal

Understanding the physical mechanisms at play in the interaction between turbulent plasma and neutral particles is a crucial issue that we approach in this Thesis by using a first-principles self-consistent model of the tokamak periphery implemented in the GBS code. While the plasma is modeled by the drift-reduced two-fluid Braginskii equations, a kinetic model for the neutrals is developed, valid in short and in long mean free path scenarios. The model includes ionization, charge-exchange, recombination, and elastic collisional processes. The neutral kinetic equation is solved by using the method of characteristics. We identify the key elements determining the interaction between neutrals and the turbulent plasma focusing on a tokamak with a toroidal rail limiter on the high-field side equatorial midplane. For this purpose, we simulate the dynamics of the plasma and the neutrals in a domain that includes both the confined edge region and the scrape-off layer (SOL). It turns out that, in the considered plasma conditions, neither the fluctuations of the neutral moments, nor the friction between neutrals and the plasma impact the time-averaged plasma profiles significantly. Thanks to this study, we derive a simple model for the neutral-plasma interaction, which is helpful to identify and understand the principal physical processes at play in the tokamak periphery. By studying the dynamics of the neutral-plasma interplay along the magnetic field lines in the SOL, we derive a refined two-point model from the drift-reduced Braginskii equations that balances the parallel and perpendicular transport of plasma and heat, and takes into account the plasma-neutral interaction. The model estimates the electron temperature drop along a field line, from a region far from the limiter to the limiter plates. The refined two-point model is shown to be in very good agreement with the simulation results. Finally, we self-consistently simulate a diagnostic neutral gas puff, which is often used experimentally as a tool to learn about the turbulence properties in the tokamak periphery. In particular, we investigate the impact of neutral density fluctuations on the D-Î± light emission, finding that at a radial distance from the gas puff smaller than the neutral mean free path, neutral density fluctuations are anti-correlated with plasma density, electron temperature, and D-Î± fluctuations, while at distances from the gas puff larger than the neutral mean free path, a non-local shadowing effect influences the neutrals, and the D-Î± fluctuations are correlated with the neutral density fluctuations.

EPFL2017

Turbulence studies in TCV using the Correlation ECE diagnostic

Matteo Fontana

In this thesis work, the flexibility of the Tokamak à Configuration Variable (TCV) is exploited to study the influence of plasma parameters on turbulent fluctuations. The correlation electron cyclotron emission (CECE) diagnostic is used to measure low amplitude, large bandwidth radiative temperature fluctuations, associated with the anomalous transport terms that constitute the largest contribution to heat and particle loss in tokamak plasmas. In a series of L-mode limited plasmas, fluctuations are characterized over a large range of plasma parameters, obtained varying collisionality, Te/Ti and plasma triangularity. The influence of negative triangularity on confinement and fluctuations is also investigated. Temperature fluctuations profiles have been measured in low density, ohmic, positive and negative triangularity discharges, with matched current and density profiles. Strong suppression of fluctuations is observed in negative triangularity discharges from the plasma edge to the plasma mid radius. This is despite triangularity quickly decreasing when moving away from the edge. The existence of a region of low stiffness at the plasma edge is postulated to be responsible for this. Negative triangularity is confirmed having beneficial effects on confinement and reducing fluctuations amplitude in regimes with high Te/Ti. The neutral beam injector (NBI) has been used on positive and negative triangularity discharges in order to study the effect of negative triangularity on plasmas where Te/Ti~1. This condition is of particular interest since future reactor-like tokamaks are expected to work with thermalized electrons and ions. A pair of discharges with positive and negative triangularity has been realized, where matched temperature profiles are attained with ~385 kW and ~145 kW of NBI power respectively for the two shapes. This is evidence of negative triangularity improving plasma confinement also in conditions of low Te/Ti. In these discharges, CECE measurements find suppressed fluctuations in the negative triangularity discharges. Linear, flux tube, gyrokinetic simulations show that, in the same plasmas, the dominant turbulent regime is a mix of ion temperature gradient (ITG) instabilities and trapped electron modes (TEM). This is different with respect to all previous works in TCV, in which no ion heating had been available and only pure TEM regimes were observed. The results of the linear simulations indicate that, also in this mixed turbolence regime, negative triangularity has a stabilizing effect on instabilities, more visible for radial positions closer to the plasma edge, where the magnitude of triangularity is higher. A database of fluctuations measurements has been constructed from a series of discharges with positive and negative triangularity covering a large range of plasma conditions, particularly focusing on the effects of different combinations of collisionality and the Te/Ti ratio. Collisionality is found to be the main parameter to control the relative suppression of fluctutations in negative triangularity plasmas. No direct effect of Te/Ti on the fluctuation amplitudes is observed in the explored parameter range. Measurements taken in positive and negative triangularity plasmas, with comparable conditions and matched normalized temperature scale lengths still show reduced fluctuation levels for negative triangularity, suggesting influence on the threshold gradients for the onset of fluctuations.

EPFL2018

Gyrokinetic particle-in-cell global simulations of ion-temperature-gradient and collisionless-trapped-electron-mode turbulence in tokamaks

Sébastien Jolliet

The goal of thermonuclear fusion research is to provide power plants, that will be able to produce one gigawatt of electricity. Among the different ways to achieve fusion, the tokamak, based on magnetic confinement, is the most promising one. A gas is heated up to hundreds of millions of degrees and becomes a plasma, which is maintained – or confined – in a toroidal vessel by helical magnetic field lines. Then, deuterium and tritium are injected and fuse to create an α particle and an energetic neutron. In order to have a favorable power balance, the power produced by fusion reactions must exceed the power needed to heat the plasma and the power losses. This can be cast in a very simple expression which stipulates that the product of the density, the temperature and the energy confinement time must exceed some given value. Unfortunately, present-days tokamaks are not able to reach this condition, mostly due to plasma turbulence. The latter phenomenon enhances the heat losses and degrades the energy confinement time, which cannot be predicted by analytical theories such as the so-called neoclassical theory in which the heat losses are caused by Coulomb collisions. Therefore, numerical simulations are being developed to model plasma turbulence, mainly caused by the Ion and Electron Temperature-Gradient and the Trapped-Electron-Mode instabilities. The plasma is described by a distribution function which evolves according to the Vlasov equation. The electromagnetic fields created by the particles are self-consistently obtained through Maxwell's equations. The resulting Vlasov-Maxwell system is greatly simplified by using the gyrokinetic theory, which consists, through an appropriate ordering, of eliminating the fast gyromotion (compared to the typical frequency of instabilities). Nevertheless, it is still extremely difficult to solve this system numerically due to the large range of time and spatial scales to be resolved. In this thesis, the Vlasov-Maxwell system is solved in the electrostatic and collisionless limit with the Particle-In-Cell (PIC) ORB5 code in global tokamak geometry. This Monte-Carlo approach suffers from statistical noise which unavoidably degrades the quality of the simulation. Consequently, the first part of this work has been devoted to the optimization of the code with a view to reduce the numerical noise. The code has been rewritten in a new coordinate system which takes advantage of the anisotropy of turbulence, which is mostly aligned with the magnetic field lines. The overall result of the optimization is that for a given accuracy, the CPU time has been decreased by a factor two thousand, the total memory has been decreased by a factor ten and the numerical noise has been reduced by a factor two hundred. In addition, the scaling of the code with respect to plasma size is presently optimal, suggesting that ORB5 could compute heat transport for future fusion devices such as ITER. The second part of this thesis presents the validation of the code with numerical convergence tests, linear (including dispersion relations) and nonlinear benchmarks. Furthermore, the code has been applied to important issues in gyrokinetic theory. It is shown for the first time that a 5D global delta-f PIC code can achieve a thermodynamic steady state on the condition that some dissipation is present. This is a fundamental result as the main criticism against delta-f PIC codes is their inability to deal with long time simulations. Next, the role of the parallel nonlinearity is studied and it is demonstrated in this work that this term has no real influence on turbulence, provided the numerical noise is sufficiently low. This result should put an end to the controversy that recently occurred, in which gyrokinetic simulations using different numerical approaches yielded contradictory results. Finally, thanks to the optimization of the code, the gyrokinetic model has been extended to include the kinetic response of trapped-electrons, in place to the usual adiabatic (Boltzmann) approximation. For the first time, global TEM nonlinear simulations are presented, and the role of the zonal flow on heat transport is analyzed. This study will help in acquiring some knowledge on the less-known TEM turbulence (as compared to ITG). In conclusion, this thesis is one of the main steps of the development of ORB5, which is now a state-of-the-art gyrokinetic code for collisionless ITG and TEM turbulence, and has brought several contributions to the understanding of these phenomena.

EPFL2009