Impurity transport in temperature gradient driven turbulence

In the present paper the transport of impurities driven by trapped electron (TE) mode turbulence is studied. Non-linear (NL) gyrokinetic simulations using the code GENE are compared with results from quasilinear (QL) gyrokinetic simulations and a computationally efficient fluid model. The main focus is on model comparisons for electron temperature gra- dient driven turbulence regarding the sign of the convective impurity velocity (pinch) and the impurity density gradient R/LnZ (peaking factor) for zero impurity flux. In particular, the scaling of the impurity peaking factors with impurity charge Z and with driving temper- ature gradient is investigated and compared with the results for Ion Temperature Gradient (ITG) driven turbulence. In addition, the impurity peaking is compared to the main ion peaking obtained by a self-consistent fluid calculation of the density gradients corresponding to zero particle fluxes. For the scaling of the peaking factor with impurity charge Z, a weak dependence is obtained from NL GENE and fluid simulations. The QL GENE results show a stronger dependence for low Z impurities and overestimates the peaking factor by up to a factor of two in this region. As in the case of ITG dominated turbulence, the peaking factors saturate as Z increases, at a level much below neoclassical predictions. However, the scaling with Z is weak or reversed as compared to the ITG case. The scaling of impurity peaking with the background temperature gradients is found to be weak in the NL GENE and fluid simulations. The QL results are also here found to significantly overestimate the peaking factor for low Z values. For the parameters considered, the background density gradient for zero particle flux is found to be slightly larger than the corresponding impurity zero flux gradient.Comment: 23 pages, 13 figures. Submitted to AIP: Physics of Plasma

Constructing a regular histogram : a comparison of methods

Even for a well-trained statistician the construction of a histogram for a given real-valued set is a sifficult problem. It is even more difficult to construct a fully automatic procedure which specifies the number and widths of the binss in a satisfactory manner for a wide range of data sets. In this paper we compare several histogram construction methods by means of a simulation study. The study includes plug-in methods, cross-validation, penalized maximum likehood and the taut string procedure. Their performance on different test beds is measured by the Hellinger distance and the ability to identify the modes of the underlying density. --regular histogramm,model selection,penalized likehood,taut-string

Upgrade and benchmark of quasi-linear transport model EDWM

The verification of a new saturation rule applied to the quasi-linear fluid model EDWM (extended drift wave model) and the calibration of several other features are presented. As one of the computationally fastest first-principle-based core transport models, EDWM can include an arbitrary number of ions and charge states. This feature is especially important for experimental devices with plasma-facing components made of heavy elements, such as the upcoming ITER device. As a quasi-linear model, EDWM solves a linear dispersion relation to obtain the instabilities driving the turbulence and combines the linear description with an estimation of the saturation level of the electrostatic potential to determine the fluxes. A new saturation rule at the characteristic length combined with a spectral filter for the poloidal wavenumber dependency is developed. The shape of the filter has been fitted against the poloidal wavenumber dependency of the electrostatic potential from non-linear gyrokinetic simulations. Additionally, EDWM\u27s collision frequency and safety factor dependencies, as well as the electron heat flux level, have been calibrated against gyrokinetic and gyrofluid results. Finally, the saturation level has been normalized against non-linear gyrokinetic simulations and later validated against experimental measured fluxes from 12 discharges at JET

Constructing a regular histogram - a comparison of methods

Even for a well-trained statistician the construction of a histogram for a given real-valued data set is a difficult problem. It is even more difficult to construct a fully automatic procedure which specifies the number and widths of the bins in a satisfactory manner for a wide range of data sets. In this paper we compare several histogram construction methods by means of a simulation study. The study includes plug-in methods, cross-validation, penalized maximum likelihood and the taut string procedure. Their performance on different test beds is measured by the Hellinger distance and the ability to identify the modes of the underlying density

Zonal flow generation in collisionless trapped electron mode turbulence

In the present work the generation of zonal flows in collisionless trapped electron mode (TEM) turbulence is studied analytically. A reduced model for TEM turbulence is utilized based on an advanced fluid model for reactive drift waves. An analytical expression for the zonal flow growth rate is derived and compared with the linear TEM growth, and its scaling with plasma parameters is examined for typical tokamak parameter values.Comment: 20 pages, 4 figure

Impact of fast ions on density peaking in JET: fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position rho(t) = 0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile

Influence of the radio frequency ponderomotive force on anomalous impurity transport in tokamaks

Trace impurity transport in tokamaks is studied using an electrostatic, collisionless fluid model for ion-temperature-gradient and trapped-electron mode driven turbulence in the presence of radio frequency (rf) fields, and the results are compared with neoclassical predictions. It is shown that the inward impurity convective velocity (pinch) that is usually obtained can be reduced by the rf fields, in particular close to the wave resonance location where the rf ponderomotive force may be significant. However, the impurity diffusivity and convective velocity are usually similarly affected by the ponderomotive force, and hence the steady-state impurity density peaking factor -∇nz/nz is only moderately affected by the rf fields

Impact of fast ions on density peaking in JET: fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/ trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position ρt = 0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile.EURATOM 63305

Interpretative and predictive modelling of Joint European Torus collisionality scans

Transport modelling of Joint European Torus (JET) dimensionless collisionality scaling experiments in various operational scenarios is presented. Interpretative simulations at a fixed radial position are combined with predictive JETTO simulations of temperatures and densities, using the TGLF transport model. The model includes electromagnetic effects and collisions as well as (E)over-right-arrow x (b)over-right-arrow shear in Miller geometry. Focus is on particle transport and the role of the neutral beam injection (NBI) particle source for the density peaking. The experimental 3-point collisionality scans include L-mode, and H-mode (D and H and higher beta D plasma) plasmas in a total of 12 discharges. Experimental results presented in (Tala et al 2017 44th EPS Conf.) indicate that for the H-mode scans, the NBI particle source plays an important role for the density peaking, whereas for the L-mode scan, the influence of the particle source is small. In general, both the interpretative and predictive transport simulations support the experimental conclusions on the role of the NBI particle source for the 12 JET discharges.Peer reviewe