G. Saibene

ELM control strategies and tools: status and potential for ITER

Operating ITER in the reference inductive scenario at the design values of Ip = 15 MA and QDT = 10 requires the achievement of good H-mode confinement that relies on the presence of an edge transport barrier whose pedestal pressure height is key to plasma performance. Strong gradients occur at the edge in such conditions that can drive magnetohydrodynamic instabilities resulting in edge localized modes (ELMs), which produce a rapid energy loss from the pedestal region to the plasma facing components (PFC). Without appropriate control, the heat loads on PFCs during ELMs in ITER are expected to become significant for operation in H-mode at Ip = 6–9 MA; operation at higher plasma currents would result in a very reduced life time of the PFCs. Currently, several options are being considered for the achievement of the required level of ELM control in ITER; this includes operation in plasma regimes which naturally have no or very small ELMs, decreasing the ELM energy loss by increasing their frequency by a factor of up to 30 and avoidance of ELMs by actively controlling the edge with magnetic perturbations. Small/no ELM regimes obtained by influencing the edge stability (by plasma shaping, rotational shear control, etc) have shown in present experiments a significant reduction of the ELM heat fluxes compared to type-I ELMs. However, so far they have only been observed under a limited range of pedestal conditions depending on each specific device and their extrapolation to ITER remains uncertain. ELM control by increasing their frequency relies on the controlled triggering of the edge instability leading to the ELM. This has been presently demonstrated with the injection of pellets and with plasma vertical movements; pellets having provided the results more promising for application in ITER conditions. ELM avoidance/suppression takes advantage of the fact that relatively small changes in the pedestal plasma and magnetic field parameters seem to have a large stabilizing effect on large ELMs. Application of edge magnetic field perturbation with non-axisymmetric fields is found to affect transport at the plasma edge and thus prevent the uncontrolled rise of the plasma pressure gradients and the occurrence of type-I ELMs. This paper compiles a brief overview of various ELM control approaches, summarizes their present achievements and briefly discusses the open issues regarding their application in ITER.

Edge localized mode physics and operational aspects in tokamaks

Recent progress in experimental and theoretical studies of edge localized mode (ELM) physics is reviewed for the reactor relevant plasma regimes, namely the high confinement regimes, that is, H-modes and advanced scenarios.Theoretical approaches to ELM physics, from a linear ideal magnetohydrodynamic (MHD) stability analysis to non-linear transport models with ELMs are discussed with respect to experimental observations, in particular the fast collapse of pedestal pressure profiles, magnetic measurements and scrape-off layer transport during ELMs.High confinement regimes with different types of ELMs are addressed in this paper in the context of development of operational scenarios for ITER. The key parameters that have been identified at present to reduce the energy losses in Type I ELMs are operation at high density, high edge magnetic shear and high triangularity. However, according to the present experimental scaling for the energy losses in Type I ELMs, the extrapolation of such regimes for ITER leads to unacceptably large heat loads on the divertor target plates exceeding the material limits. High confinement H-mode scenarios at high triangularity and high density with small ELMs (Type II), mixed regimes (Type II and Type I) and combined advanced regimes at high βp are discussed for present-day tokamaks. The optimum combination of high confinement and small MHD activity at the edge in Type II ELM scenarios is of interest to ITER. However, to date, these regimes have been achieved in a rather narrow operational window and far from ITER parameters in terms of collisionality, edge safety factor and βp.The compatibility of the alternative internal transport barrier (ITB) scenario with edge pedestal formation and ELMs is also addressed. Edge physics issues related to the possible combination of small benign ELMs (Type III, Type II ELMs, quiescent double barrier) and high performance ITBs are discussed for present-day experiments (JET, JT-60U, DIII-D) in terms of their relevance for ITER. Successful plasma edge control, at high triangularity (~0.5) and high density (~0.7nGR), in ITB scenarios in JET is reported.Active control of ELMs by edge current, pellet injection, impurities and external magnetic perturbations creating an ergodic zone localized at the separatrix are discussed for present-day experiments and from the perspective of future reactors.

Effect of toroidal field ripple on plasma rotation in JET

Dedicated experiments on TF ripple effects on the performance of tokamak plasmas have been carried out at JET. The TF ripple was found to have a profound effect on the plasma rotation. The central Mach number, M, defined as the ratio of the rotation velocity and the thermal velocity, was found to drop as a function of TF ripple amplitude (3) from an average value of M = 0.40-0.55 for operations at the standard JET ripple of 6 = 0.08% to M = 0.25-0.40 for 6 = 0.5% and M = 0.1-0.3 for delta = 1%. TF ripple effects should be considered when estimating the plasma rotation in ITER. With standard co-current injection of neutral beam injection (NBI), plasmas were found to rotate in the co-current direction. However, for higher TF ripple amplitudes (delta similar to 1%) an area of counter rotation developed at the edge of the plasma, while the core kept its co-rotation. The edge counter rotation was found to depend, besides on the TF ripple amplitude, on the edge temperature. The observed reduction of toroidal plasma rotation with increasing TF ripple could partly be explained by TF ripple induced losses of energetic ions, injected by NBI. However, the calculated torque due to these losses was insufficient to explain the observed counter rotation and its scaling with edge parameters. It is suggested that additional TF ripple induced losses of thermal ions contribute to this effect.

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

JET intrinsic rotation studies in plasmas with a high normalized beta and varying toroidal field ripple

Understanding the origin of rotation in ion cyclotron resonance frequency (ICRF) heated plasmas is important for predictions for burning plasmas sustained by alpha particles, being characterized by a large population of fast ions and no external momentum input. The angular velocity of the plasma column has been measured in JET H-mode plasmas with pure ICRF heating both for the standard low toroidal magnetic ripple configuration, of about similar to 0.08% and, for increased ripple values up to 1.5% (Nave et al 2010 Phys. Rev. Lett. 105 105005). These new JET rotation data were compared with the multi-machine scaling of Rice et al (2007 Nucl. Fusion 47 1618) for the Alfven-Mach number and with the scaling for the velocity change from L-mode into H-mode. The JET data do not fit well any of these scalings that were derived for plasmas that are co-rotating with respect to the plasma current. With the standard low ripple configuration, JET plasmas with large ICRF heating power and normalized beta, beta(N) approximate to 1.3, have a very small co-current rotation, with Alfven-Mach numbers significantly below those given by the rotation scaling of Rice et al (2007 Nucl. Fusion 47 1618). In some cases the plasmas are actually counter-rotating. No significant difference between the H-mode and L-mode rotation is observed. Typically the H-mode velocities near the edge are lower than those in L-modes. With ripple values larger than the standard JET value, between 1% and 1.5%, H-mode plasmas were obtained where both the edge and the core counter-rotated.

Comparison between dominant NB and dominant IC heated ELMy H-mode discharges in JET

The experiment described in this paper is aimed at characterization of ELMy H-mode discharges with varying momentum input, rotation, power deposition profiles and ion to electron heating ratio obtained by varying the proportion between ion cyclotron (IC) and neutral beam (NB) heating. The motivation for the experiment was to verify if the basic confinement and transport properties of the baseline ITER H-mode are robust to these changes, and similar to those derived mostly from dominant NB heated H-modes. No significant difference in the density and temperature profiles or in the global confinement were found. Although ion temperature profiles were seen to be globally stiff, some variation of stiffness was obtained in the experiment by varying the deposition profiles, but not one that could significantly affect the profiles in terms of global confinement. This analysis shows the thermal plasma energy confinement enhancement factor to be independent of the heating mix, for the range of conditions explored. Moreover, the response of the global confinement to changes in density and power were also independent of heating mix, reflecting the changes in the pedestal, which is in agreement with globally stiff profiles. Consistently, the pedestal characteristics (pressure and width) and their dependences on global parameters such as density and power were the same during NB only or with predominant IC heating.