D. Zarzoso

H-mode pedestal scaling in DIII-D, ASDEX Upgrade, and JET

Multi device pedestal scaling experiments in the DIII-D, ASDEX Upgrade (AUG) and JET tokamaks are presented in order to test two plasma physics pedestal width models. The first model proposes a scaling of the pedestal width Δ/a ρ * 1/2 to ρ * based on the radial extent of the pedestal being set by the point where the linear turbulence growth rate exceeds the ExB velocity. In the multi device experiment where ρ * at the pedestal top was varied by a factor of four while other dimensionless parameters where kept fixed, it has been observed that the temperature pedestal width in real space coordinates scales with machine size, and that therefore the gyroradius scaling suggested by the model is not supported by the experiments. This density pedestal width is not invariant with ρ * which after comparison with a simple neutral fuelling model may be attributed to variations in the neutral fuelling patterns. The second model, EPED1, is based on kinetic ballooning modes setting the limit of the radial extent of the pedestal region and leads to Δ βp 1/2 . All three devices show a scaling of the pedestal width in normalised poloidal flux as Δ βp 1/2 , as described by the kinetic ballooning model, however on JET and AUG this could not be distinguished from an interpretation where the pedestal is fixed in real space. Pedestal data from all three devices have been compared with the predictive pedestal model EPED1 and the model produces pedestal height values that match the experimental data well.

Fully kinetic description of the linear excitation and nonlinear saturation of fast-ion-driven geodesic acoustic mode instability

We show in this paper that geodesic acoustic modes (GAMs) can be efficiently excited by a population of fast ions even when Landau damping on thermal ions is accounted for. We report in particular fully kinetic calculations of the GAM dispersion relation and its complete solution. Written under a variational form, the quasi-neutrality condition, together with the kinetic Vlasov equation, leads to the density of exchanged energy between particles and the mode. In particular, a linear threshold for the GAMs excitation is derived. Two examples of fast ion distribution have been discussed analytically. It turns out that particles with high perpendicular energy compared to the parallel resonance energy are most responsible for the excitation of the mode. Subsequent numerical simulations of circular plasmas using gysela code have been carried out. In particular, the linear kinetic threshold has been reproduced during the excitation phase, and a nonlinear saturation has been observed. Analysis in the phase space...

Comparison of hybrid and baseline ELMy H-mode confinement in JET with the carbon wall

The confinement in JET baseline type I ELMy H-mode plasmas is compared to that in so-called hybrid H-modes in a database study of 112 plasmas in JET with the carbon fibre composite (CFC) wall. The baseline plasmas typically have ?N???1.5?2, H98???1, whereas the hybrid plasmas have ?N???2.5?3, H98?<?1.5. The database study contains both low- (????0.2?0.25) and high-triangularity (????0.4) hybrid and baseline H-mode plasmas from the last JET operational campaigns in the CFC wall from the period 2008?2009. Based on a detailed confinement study of the global as well as the pedestal and core confinement, there is no evidence that the hybrid and baseline plasmas form separate confinement groups; it emerges that the transition between the two scenarios is of a gradual kind rather than demonstrating a bifurcation in the confinement. The elevated confinement enhancement factor H98 in the hybrid plasmas may possibly be explained by the density dependence in the ?98 scaling as n0.41 and the fact that the hybrid plasmas operate at low plasma density compared to the baseline ELMy H-mode plasmas. A separate regression on the confinement data in this study shows a reduction in the density dependence as n0.09?0.08. Furthermore, inclusion of the plasma toroidal rotation in the confinement regression provides a scaling with the toroidal Alfv?n Mach number as and again a reduced density dependence as n0.15?0.08. The differences in pedestal confinement can be explained on the basis of linear MHD stability through a coupling of the total and pedestal poloidal pressure and the pedestal performance can be improved through plasma shaping as well as high ? operation. This has been confirmed in a comparison with the EPED1 predictive pedestal code which shows a good agreement between the predicted and measured pedestal pressure within 20?30% for a wide range of ?N???1.5?3.5. The core profiles show a strong degree of pressure profile consistency. No beneficial effect of core density peaking on confinement could be identified for the majority of the plasmas presented here as the density peaking is compensated by a temperature de-peaking resulting in no or only a weak variation in the pressure peaking. The core confinement could only be optimized in case the ions and electrons are decoupled, in which case the ion temperature profile peaking can be enhanced, which benefits confinement. In this study, the latter has only been achieved in the low-triangularity hybrid plasmas, and can be attributed to low-density operation. Plasma rotation has been found to reduce core profile stiffness, and can explain an increase in profile peaking at small radius ?tor?=?0.3.

Pedestal width and ELM size identity studies in JET and DIII-D: implications for ITER

The dependence of the H-mode edge transport barrier width on normalized ion gyroradius (rho* = rho/a) in discharges with type I ELMs was examined in experiments combining data for the JET and DIII-D tokamaks. The plasma configuration as well as the local normalized pressure (beta), collisionality (nu*), Mach number and the ratio of ion and electron temperature at the pedestal top were kept constant, while rho* was varied by a factor of four. The width of the steep gradient region of the electron temperature (T-e) and density (n(e)) pedestals normalized to machine size showed no or only a weak trend with rho*. A rho(1/2) or rho(1) dependence of the pedestal width, given by some theoretical predictions, is not supported by the current experiments. This is encouraging for the pedestal scaling towards ITER as it operates at lower rho* than existing devices. Some differences in pedestal structure and ELM behaviour were, however, found between the devices; in the DIII-D discharges, the n(e) and T-e pedestal were aligned at high rho* but the ne pedestal shifted outwards in radius relative to T-e as rho* decreases, while on JET the profiles remained aligned while rho* was scanned by a factor of two. The energy loss at an ELM normalized to the pedestal energy increased from 10% to 40% as rho* increased by a factor of two in the DIII-D discharges but no such variation was observed in the case of JET. The measured pedestal pressures and widths were found to be consistent with the predictions from modelling based on peeling-ballooning stability theory, and are used to make projections towards ITER

Observation of Confined Current Ribbon in JET Plasmas

We report the identification of a localized current structure inside the JET plasma. It is a field-aligned closed helical ribbon, carrying current in the same direction as the background current profile (cocurrent), rotating toroidally with the ion velocity (corotating). It appears to be located at a flat spot in the plasma pressure profile, at the top of the pedestal. The structure appears spontaneously in low density, high rotation plasmas, and can last up to 1.4 s, a time comparable to a local resistive time. It considerably delays the appearance of the first edge localized mode.

A 5D gyrokinetic full-f global semi-lagrangian code for flux-driven ion turbulence simulations

This paper addresses non-linear gyrokinetic simulations of ion temperature gradient (ITG) turbulence in tokamak plasmas. The electrostatic Gysela code is one of the few international 5D gyrokinetic codes able to perform global, full-f and flux-driven simulations. Its has also the numerical originality of being based on a semi-Lagrangian (SL) method. This reference paper for the Gysela code presents a complete description of its multi-ion species version including: (i) numerical scheme, (ii) high level of parallelism up to 500k cores and (iii) conservation law properties.

The linear tearing instability in three dimensional, toroidal gyro-kinetic simulations

Linear gyro-kinetic simulations of the classical tearing mode in three-dimensional toroidal geometry were performed using the global gyro-kinetic turbulence code, GKW. The results were benchmarked against a cylindrical ideal MHD and analytical theory calculations. The stability, growth rate, and frequency of the mode were investigated by varying the current profile, collisionality, and the pressure gradients. Both collisionless and semi-collisional tearing modes were found with a smooth transition between the two. A residual, finite, rotation frequency of the mode even in the absence of a pressure gradient is observed, which is attributed to toroidal finite Larmor-radius effects. When a pressure gradient is present at low collisionality, the mode rotates at the expected electron diamagnetic frequency. However, the island rotation reverses direction at high collisionality. The growth rate is found to follow a η1∕7 scaling with collisional resistivity in the semi-collisional regime, closely following the se...

Relation between Energetic and Standard Geodesic Acoustic Modes

Geodesic Acoustic Modes (GAMs) are electrostatic, axisymmetric modes which are non-linearly excited by turbulence. They can also be excited linearly by fast-particles; they are then called Energetic-particle-driven GAMs (EGAMs). Do GAMs and EGAMs belong to the same mode branch? Through a linear, analytical model, in which the fast particles are represented by a Maxwellian bump-on-tail distribution function, we find that the answer depends on several parameters. For low values of the safety factor q and for high values of the fast ion energy, the EGAM originates from the GAM. On the contrary, for high values of q and for low values of the fast ion energy, the GAM is not the mode which becomes unstable when fast particles are added: the EGAM then originates from a distinct mode, which is strongly damped in the absence of fast particles. The impact of other parameters is further explored: ratio of the ion temperature to the electron temperature, width of the fast particle distribution, mass and charge of the fast ions. The ratio between the EGAM and the GAM frequencies was found in experiments (DIII-D) and in non-linear numerical simulations (code GYSELA) to be close to 1/2: the present analytical study allows one to recover this ratio.

Turbulent momentum transport in core tokamak plasmas and penetration of scrape-off layer flows

The turbulent transport of toroidal angular momentum in the core of a tokamak plasma is investigated in global, full-f gyrokinetic simulations, performed with the GYSELA code in the flux-driven regime. During the initial turbulent phase, a front of positive Reynolds stress propagates radially, generating intrinsic toroidal rotation from a vanishing initial profile. This is also accompanied by a propagating front of turbulent heat flux. In the statistical steady-state regime, turbulent transport exhibits large-scale avalanche-like events which are found to transport both heat and momentum, and similar statistical properties are obtained for both fluxes. The impact of scrape-off layer flows is also investigated by modifying the boundary conditions in the simulations. The observed impact is radially localized for L-mode like poloidal profiles of parallel velocity at the edge, while a constant velocity at the edge can modify the core toroidal rotation profile in a large fraction of the radial domain.

Entropy production and collisionless fluid closure

A novel method is proposed to construct collisionless fluid closures accounting for some kinetic properties. The idea consists in optimizing the agreement between the fluid and kinetic quasi-linear entropy production rates, so as to constrain the closure coefficients. This procedure is applied to the slab branch of the ion temperature gradient driven instability. Focusing on the kinetic regime characterized by slow waves, the closure proposed by Hammett and Perkins (Hammett and Perkins 1990 Phys. Rev. Lett. 64 3019) naturally emerges from the systematic identification of the kinetic and fluid entropy production rates. This closure is revealed to be extremely powerful well beyond the kinetic regime. Besides, it reconciles the fluid and kinetic linear stability diagrams in the two-dimensional space of the density and temperature gradient lengths. Such a method is systematic and generic. As such, it is applicable to other models and classes of instabilities.