No Arabic abstract
Gyrokinetic simulations of ion temperature gradient mode and trapped electron mode driven impurity transport in a realistic tokamak geometry are presented and compared with results using simplified geometries. The gyrokinetic results, obtained with the GENE code in both linear and non-linear modes are compared with data and analysis for a dedicated impurity injection discharge at JET. The impact of several factors on heat and particle transport is discussed, lending special focus to tokamak geometry and rotational shear. To this end, results using s-alpha and concentric circular equilibria are compared with results with magnetic geometry from a JET experiment. To further approach experimental conditions, non-linear gyrokinetic simulations are performed with collisions and a carbon background included. The impurity peaking factors, computed by finding local density gradients corresponding to zero particle flux, are discussed. The impurity peaking factors are seen to be reduced by a factor of ~2 in realistic geometry compared with the simplified geometries, due to a reduction of the convective pinch. It is also seen that collisions reduce the peaking factor for low-Z impurities, while increasing it for high charge numbers, which is attributed to a shift in the transport spectra towards higher wavenumbers with the addition of collisions. With the addition of roto-diffusion, an overall reduction of the peaking factors is observed, but this decrease is not sufficient to explain the flat carbon profiles seen at JET.
The turbulence-induced quasi-linear particle flux of a highly-charged, collisional impurity species is calculated from the electrostatic gyrokinetic equation including collisions with the bulk ions and the impurities themselves. The equation is solved by an expansion in powers of the impurity charge number $Z$. In this formalism, the collision operator only affects the impurity flux through the dynamics of the impurities in the direction parallel to the magnetic field. At reactor-relevant collisionality, the parallel dynamics is dominated by the parallel electric field, and collisions have a minor effect on the turbulent particle flux of highly-charged, collisional impurities.
Barriers have been identified in magnetically confined plasmas reducing the particle transport and improving the confinement. One of them, the primary shearless barriers are associated to extrema of non-monotonic plasma profiles. Previously, we identified these barriers in a model described by a map that allows the integration of charged particles motion in drift waves for a long time scale. In this work, we show how the existence of these robust barriers depends on the fluctuation amplitude and on the electric shear. Moreover, we also find control parameter intervals for which these primary barriers onset and break-up are recurrent. Another noticeable feature, in these transitions, is the appearance of a layer of particle trajectory stickiness after the shearless barrier break-up or before its onset. Besides the mentioned primary barriers, we also observe sequences of secondary shearless barriers, not reported before, created and destroyed by a sequence of bifurcations as the main control parameters, the fluctuation amplitude and electric shear, are varied. Furthermore, in these bifurcations, we also find hitherto unknown double and triple secondary shearless barriers which constitute a noticeable obstacle to the chaotic transport.
Particle transport due to Ion Temperature Gradient/Trapped Electron (ITG/TE) mode turbulence is investigated using the gyrokinetic code GENE. Both a reduced quasilinear (QL) treatment and nonlinear (NL) simulations are performed for typical tokamak parameters corresponding to ITG dominated turbulence. A selfconsistent treatment is used, where the stationary local profiles are calculated corresponding to zero particle flux simultaneously for electrons and trace impurities. The scaling of the stationary profiles with magnetic shear, safety factor, electron-to-ion temperature ratio, collisionality, toroidal sheared rotation, triangularity, and elongation is investigated. In addition, the effect of different main ion mass on the zero flux condition is discussed. The electron density gradient can significantly affect the stationary impurity profile scaling. It is therefore expected, that a selfconsistent treatment will yield results more comparable to experimental results for parameter scans where the stationary background density profile is sensitive. This is shown to be the case in scans over magnetic shear, collisionality, elongation, and temperature ratio, for which the simultaneous zero flux electron and impurity profiles are calculated. A slight asymmetry between hydrogen, deuterium and tritium with respect to profile peaking is obtained, in particular for scans in collisionality and temperature ratio.
In tokamaks, internal transport barriers, produced by modifications of the plasma current profile, reduce particle transport and improve plasma confinement. The triggering of the internal transport barriers and their dependence on the plasma profiles is a key nonlinear dynamics problem still under investigation. We consider the onset of shearless invariant curves inside the plasma which create internal transport barriers. A non-integrable drift-kinetic model is used to describe particle transport driven by drift waves and to investigate these shearless barriers onset in tokamaks. We show that for some currently observed plasma profiles shearless particle transport barriers can be triggered by properly modifying the electric field profile and the influence of non-resonant modes in the barriers onset. In particular, we show that a broken barrier can be restored by enhancing non-resonant modes.
Neutral atoms can strongly influence the intrinsic rotation and radial electric field at the tokamak edge. Here, we present a framework to investigate these effects when the neutrals dominate the momentum transport. We explore the parameter space numerically, using highly flexible model geometries and a state of the art kinetic solver. We find that the most important parameters controlling the toroidal rotation and electric field are the major radius where the neutrals are localized and the plasma collisionality. This offers a means to influence the rotation and electric field by, for example, varying the radial position of the X-point to change the major radius of the neutral peak.