I. Voitsekhovitch

Physics of transport in tokamaks

This paper is an overview of recent results relating to turbulent particle and heat transport, and to the triggering of internal transport barriers (ITBs). The dependence of the turbulent particle pinch velocity on plasma parameters has been clarified and compared with experiment. Magnetic shear and collisionality are found to play a central role. Analysis of heat transport has made progress along two directions: dimensionless scaling laws, which are found to agree with the prediction for electrostatic turbulence, and analysis of modulation experiments, which provide a stringent test of transport models. Finally the formation of ITBs has been addressed by analysing electron transport barriers. It is confirmed that negative magnetic shear, combined with the Shafranov shift, is a robust stabilizing mechanism. However, some well established features of internal barriers are not explained by theory.

Integration of a radiative divertor for heat load control into JET high triangularity ELMy H-mode plasmas

Experiments on JET with a carbon-fibre composite wall have explored the reduction of steady-state power load in an ELMy H-mode scenario at high Greenwald fraction similar to 0.8, constant power and close to the L to H transition. This paper reports a systematic study of power load reduction due to the effect of fuelling in combination with seeding over a wide range of pedestal density ((4-8) x 10(19) m(-3)) with detailed documentation of divertor, pedestal and main plasma conditions, as well as a comparative study of two extrinsic impurity nitrogen and neon. It also reports the impact of steady-state power load reduction on the overall plasma behaviour, as well as possible control parameters to increase fuel purity. Conditions from attached to fully detached divertor were obtained during this study. These experiments provide reference plasmas for comparison with a future JET Be first wall and an all W divertor where the power load reduction is mandatory for operation.

Global and pedestal confinement in JET with a Be/W metallic wall

Type I ELMy H-mode operation in JET with the ITER-like Be/W wall (JET-ILW) generally occurs at lower pedestal pressures compared to those with the full carbon wall (JET-C). The pedestal density is similar but the pedestal temperature where type I ELMs occur is reduced and below to the so-called critical type I–type III transition temperature reported in JET-C experiments. Furthermore, the confinement factor H98(y,2) in type I ELMy H-mode baseline plasmas is generally lower in JET-ILW compared to JET-C at low power fractions Ploss/Pthr,08 2, the confinement in JET-ILW hybrid plasmas is similar to that in JET-C. A reduction in pedestal pressure is the main reason for the reduced confinement in JET-ILW baseline ELMy H-mode plasmas where typically H98(y,2) = 0.8 is obtained, compared to H98(y,2) = 1.0 in JET-C. In JET-ILW hybrid plasmas a similarly reduced pedestal pressure is compensated by an increased peaking of the core pressure profile resulting in H98(y,2) ≤ 1.25. The pedestal stability has significantly changed in high triangularity baseline plasmas where the confinement loss is also most apparent. Applying the same stability analysis for JET-C and JET-ILW, the measured pedestal in JET-ILW is stable with respect to the calculated peeling–ballooning stability limit and the ELM collapse time has increased to 2 ms from typically 200 µs in JET-C. This indicates that changes in the pedestal stability may have contributed to the reduced pedestal confinement in JET-ILW plasmas. A comparison of EPED1 pedestal pressure prediction with JET-ILW experimental data in over 500 JET-C and JET-ILW baseline and hybrid plasmas shows a good agreement with 0.8 < (measured pped)/(predicted pped,EPED) < 1.2, but that the role of triangularity is generally weaker in the JET-ILW experimental data than in the model predictions.

First scenario development with the JET new ITER-like wall

In the recent JET experimental campaigns with the new ITER-like wall (JET-ILW), major progress has been achieved in the characterization and operation of the H-mode regime in metallic environments: (i) plasma breakdown has been achieved at the first attempt and X-point L-mode operation recovered in a few days of operation; (ii) stationary and stable type-I ELMy H-modes with beta(N) similar to 1.4 have been achieved in low and high triangularity ITER-like shape plasmas and are showing that their operational domain at H = 1 is significantly reduced with the JET-ILW mainly because of the need to inject a large amount of gas (above 10(22) Ds(-1)) to control core radiation; (iii) in contrast, the hybrid H-mode scenario has reached an H factor of 1.2-1.3 at beta(N) of 3 for 2-3 s; and, (iv) in comparison to carbon equivalent discharges, total radiation is similar but the edge radiation is lower and Z(eff) of the order of 1.3-1.4. Strong core radiation peaking is observed in H-mode discharges at a low gas fuelling rate (i. e. below 0.5 x 10(22) Ds(-1)) and low ELM frequency (typically less than 10 Hz), even when the tungsten influx from the diverter is constant. High-Z impurity transport from the plasma edge to the core appears to be the dominant factor to explain these observations. This paper reviews the major physics and operational achievements and challenges that an ITER-like wall configuration has to face to produce stable plasma scenarios with maximized performance.

Current ramps in tokamaks: from present experiments to ITER scenarios

In order to prepare adequate current ramp-up and ramp-down scenarios for ITER, present experiments from various tokamaks have been analysed by means of integrated modelling in view of determining relevant heat transport models for these operation phases. A set of empirical heat transport models for L-mode (namely, the Bohm–gyroBohm model and scaling based models with a specific fixed radial shape and energy confinement time factors of H96−L = 0.6 or HIPB98 = 0.4) has been validated on a multi-machine experimental dataset for predicting the li dynamics within ±0.15 accuracy during current ramp-up and ramp-down phases. Simulations using the Coppi–Tang or GLF23 models (applied up to the LCFS) overestimate or underestimate the internal inductance beyond this accuracy (more than ±0.2 discrepancy in some cases). The most accurate heat transport models are then applied to projections to ITER current ramp-up, focusing on the baseline inductive scenario (main heating plateau current of Ip = 15 MA). These projections include a sensitivity study to various assumptions of the simulation. While the heat transport model is at the heart of such simulations (because of the intrinsic dependence of the plasma resistivity on electron temperature, among other parameters), more comprehensive simulations are required to test all operational aspects of the current ramp-up and ramp-down phases of ITER scenarios. Recent examples of such simulations, involving coupled core transport codes, free-boundary equilibrium solvers and a poloidal field (PF) systems controller are also described, focusing on ITER current ramp-down.

Neutron spectroscopy measurements and modeling of neutral beam heating fast ion dynamics

The energy spectrum of the neutron emission from beam-target reactions in fusion plasmas at the Joint European Torus (JET) has been investigated. Different beam energies as well as injection angles were used. Both measurements and simulations of the energy spectrum were done. The measurements were made with the time-of-flight spectrometer TOFOR. Simulations of the neutron spectrum were based on first-principle calculations of neutral beam deposition profiles and the fast ion slowing down in the plasma using the code NUBEAM, which is a module of the TRANSP package. The shape of the neutron energy spectrum was seen to vary significantly depending on the energy of the beams as well as the injection angle and the deposition profile in the plasma. Cross validations of the measured and modeled neutron energy spectra were made, showing a good agreement for all investigated scenarios.

The European Integrated Tokamak Modelling (ITM) effort: achievements and first physics results

A selection of achievements and first physics results are presented of the European Integrated Tokamak Modelling Task Force (EFDA ITM-TF) simulation framework, which aims to provide a standardized platform and an integrated modelling suite of validated numerical codes for the simulation and prediction of a complete plasma discharge of an arbitrary tokamak. The framework developed by the ITM-TF, based on a generic data structure including both simulated and experimental data, allows for the development of sophisticated integrated simulations (workflows) for physics application.The equilibrium reconstruction and linear magnetohydrodynamic (MHD) stability simulation chain was applied, in particular, to the analysis of the edgeMHDstability of ASDEX Upgrade type-I ELMy H-mode discharges and ITER hybrid scenario, demonstrating the stabilizing effect of an increased Shafranov shift on edge modes. Interpretive simulations of a JET hybrid discharge were performed with two electromagnetic turbulence codes within ITM infrastructure showing the signature of trapped-electron assisted ITG turbulence. A successful benchmark among five EC beam/ray-tracing codes was performed in the ITM framework for an ITER inductive scenario for different launching conditions from the equatorial and upper launcher, showing good agreement of the computed absorbed power and driven current. Selected achievements and scientific workflow applications targeting key modelling topics and physics problems are also presented, showing the current status of the ITM-TF modelling suite.

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.

The beta scaling of energy confinement in ELMy H-modes in JET

The disagreement between the weak dependence of the energy confinement time on normalized pressure, β, observed in dedicated scans and the strongly negative dependence in the confinement scaling laws used for the design of next step tokamaks and future reactors, remains an outstanding problem. As such, scans of β have been undertaken in single null, low triangularity (δ ≈ 0.2) ELMy H-mode plasmas in JET with the MarkIIGB-SRP divertor. The scans varied β by a factor of 2.8 (normalized β from 0.72 to 2.04) and covered a range of magnetic fields (1.5–2.3 T), plasma currents (1.5–2.75 MA) and safety factors (q95 = 2.8 and 3.3). A weak β dependence was observed both globally (B0τE varied less than 9% across any one scan) and locally. A scan within Type I ELMy H-modes suggests that this weaker dependence is not due to ELM regimes. A statistical analysis indicates that these results are consistent with log–linear regressions performed on a wide JET database of ELMy H-modes, if correlations in this database are considered.

Internal transport barrier dynamics with plasma rotation in JET

At JET the dynamics of internal transport barriers (ITBs) has been explored by trying to decouple the effects of heating on the one hand and torque on the other with the ultimate objective of identifying the minimum torque required for the formation of transport barriers. The experiments shed light on the physics behind the initial trigger for ITBs, which often shows to be linked to the shape of the q profile and magnetic shear, while the further development was influenced by the strength of the rotational shear. In discharges with a small amount of rotational shear ITBs were triggered, which suggest that the overall rotational shear is not the dominant factor in the triggering process. However, the subsequent growth of the barrier was limited if the rotational shear was too low at the time of triggering. This growth phase may be highly non-linear, with several possible positive feedback loops, such as the increases in the toroidal and poloidal component of the rotational shear caused by the ITB itself.