Fusion performance in tokamaks hinges critically on the efficacy of the edge transport barrier (ETB) in suppressing energy losses. The new concept of 'fingerprints' is introduced to identify the instabilities that cause transport losses in the ETBs of many of today's experiments, from among widely posited candidates. Analysis of the gyrokinetic-Maxwell equations and gyrokinetic simulations of experiments reveals that each mode type produces characteristic ratios of transport in the various channels: density, heat, and impurities. This, together with experimental observations of transport in some channel or of the relative size of the driving sources of channels, can identify or determine the dominant modes causing energy transport. In multiple H-mode cases with edge-localized modes that are examined, these fingerprints indicate that magnetohydrodynamic (MHD)-like modes are apparently not the dominant agent of energy transport; rather, this role is played by micro-tearing modes (MTMs) and electron temperature gradient (ETG) modes, and in addition, possibly by ion temperature gradient/ trapped electron modes (ITG/TEM) on JET (Joint European 'Torus). MHD-like modes may dominate the electron particle losses. Fluctuation frequency can also be an important means of identification, and is often closely related to the transport fingerprint. The analytical arguments unify and explain previously disparate experimental observations on multiple devices, including DIII-D, JET, and ASDEX-U. Detailed simulations of two DIII-D ETBs also demonstrate and corroborate this.

In this paper, we analyze heat transport in the JET tokamak using data from its high resolution ECE diagnostic and analyses based on the transfer entropy (TE). The analysis reveals that heat transport is not smooth and continuous, but is characterized by 'trapping regions' separated by `minor transport barriers'. Meat may 'jump over' these barriers and when the heating power is raised, this 'jumping' behavior becomes more prominent. To check that our results are relevant for global heat transport, we deduced an effective diffusion coefficient from the TE results. Both its value and overall radial variation are consistent with heat diffusivities reported in literature. The detailed radial structure of the effective diffusion coefficient was shown to be linked to the mentioned minor transport barriers.

Results from two different sets of JET experiments are presented.In experiments in which toroidicity-induced Alfvén eigenmodes (TAEs) localized at different radial locations had the same frequencies and toroidal mode numbers, the occurrence of enhanced losses after the excitation of TAEs in the core of the plasma was observed.On the contrary, enhanced losses were not observed if the TAEs localized at different radial locations had different frequencies and toroidal mode numbers.Numerical modeling indicates that, in the first set of experiments, the enhanced losses were caused by a combined effect of the TAEs localized at different radial locations.The TAEs localized in the plasma core transported energetic ions from the core to outer regions of the plasma.Then, the TAEs localized in outer regions of the plasma interacted with these ions just transported by the core-localized TAEs causing a further radial displacement of the ions to the plasma edge.This process eventually ends up causing the loss of the resonant ions.In the second set of experiments, it was found that TAEs localized in the plasma core and in outer regions did not interact with the same ions and so no enhanced losses were measured.Sheared profiles of the safety factor combined with flat mass density profiles lead to larger differences on the frequencies of the TAEs localized at different radial locations, eventually avoiding loss of energetic ions through the described mechanism.

In this study, beryllium tiles from Joint European Torus (JET) vacuum vessel wall were analysed and compared regarding their position in the vacuum vessel and differences in the exploitation conditions during two campaigns of ITER-Like-Wall (ILW) in 2011-2012 (ILW1) and 2013-2014 (ILW2) Tritium content in beryllium samples were assessed. Two methods were used to measure tritium content in the samples - dissolution under controlled conditions and tritium thermal desorption. Prior to desorption and dissolution experiments, scanning electron microscopy and energy dispersive x-ray spectroscopy were used to study structure and chemical composition of plasma-facing-surfaces of the beryllium samples. Experimental results revealed that tritium content in the samples is in range of 2.10(11)-2.10(13) tritium atoms per square centimetre of the surface area with its highest content in the samples from the outer wall of the vacuum vessel (up to 1.9.10(13) atoms/cm(2) in ILW1 campaign and 2.4.10(13) atoms/cm(2) in ILW2). The lowest content of tritium was found in the upper part of the vacuum vessel (2.0.10(12) atoms/cm(2) and 2.0.10(11) atoms/cm(2) in ILW1 and ILW2, respectively). Results obtained from scanning electron microscopy has shown that surface morphology is different within single tile, however if to compare two campaigns main tendencies remains similar.

The first divertor was installed in the JET machine between 1992 and 1994 and was operated with carbon tiles and then beryllium tiles in 1994-5. Post-mortem studies after these first experiments demonstrated that most of the impurities deposited in the divertor originate in the main chamber, and that asymmetric deposition patterns generally favouring the inner divertor region result from drift in the scrape-off layer. A new monolithic divertor structure was installed in 1996 which produced heavy deposition at shadowed areas in the inner divertor corner, which is where the majority of the tritium was trapped by co-deposition during the deuterium-tritium experiment in 1997. Different divertor geometries have been tested since such as the Gas-Box and High-Delta divertors; a principle objective has been to predict plasma behaviour, transport and tritium retention in ITER. Transport modelling experiments were carried out at the end of four campaigns by puffing C-13-labelled methane, and a range of diagnostics such as quartz-microbalance and rotating collectors have been installed to add time resolution to the post-mortem analyses. The study of material migration after D-D and D-T campaigns clearly revealed important consequences of fuel retention in the presence of carbon walls. They gave a strong impulse to make a fundamental change of wall materials. In 2010 the carbon divertor and wall tiles were removed and replaced with tiles with Be or W surfaces for the ITER-Like Wall Project.

Data on erosion and melting of beryllium upper limiter tiles, so-called dump plates (DP), are presented for all three campaigns in the JET tokamak with the ITER-like wall. High-resolution images of the upper wall of JET show clear signs of flash melting on the ridge of the roofshaped tiles. The melt layers move in the poloidal direction from the inboard to the outboard tile, ending on the last DP tile with an upward going waterfall-like melt structure. Melting was caused mainly by unmitigated plasma disruptions. During three ILW campaigns, around 15% of all 12376 plasma pulses were catalogued as disruptions. Thermocouple data from the upper dump plates tiles showed a reduction in energy delivered by disruptions with fewer extreme events in the third campaign, ILW-3, in comparison to ILW-1 and ILW-2. The total Be erosion assessed via precision weighing of tiles retrieved from JET during shutdowns indicated the increasing mass loss across campaigns of up to 0.6 g from a single tile. The mass of splashed melted Be on the upper walls was also estimated using the high-resolution images of wall components taken after each campaign. The results agree with the total material loss estimated by tile weighing (similar to 130 g). Morphological and structural analysis performed on Be melt layers revealed a multilayer structure of re-solidified material composed mainly of Be and BeO with some heavy metal impurities Ni, Fe, W. IBA analysis performed across the affected tile ridge in both poloidal and toroidal direction revealed a low D concentration, in the range 1-4 x 10(17) D atoms cm(-2).

Beryllium limiter marker tiles were exposed to plasma in the Joint European Torus to diagnose the erosion of main chamber wall materials. A limiter marker tile consists of a beryllium coating layer (7-9 mu m) on the top of bulk beryllium, with a nickel interlayer (2-3 mu m) between them. The thickness variation of the beryllium coating layer, after exposure to plasma, could indicate the erosion measured by ion beam analysis with backscattering spectrometry. However, interpretations from broad beam backscattering spectra were limited by the non-uniform surface structures. Therefore, micro-ion beam analysis (mu-IBA) with 3 MeV proton beam for Elastic back scattering spectrometry (EBS) and PIXE was used to scan samples. The spot size was in the range of 3-10 mu m. Scanned areas were analysed with scanning electron microscopy (SEM) as well. Combining results from mu-IBA and SEM, we obtained local spectra from carefully chosen areas on which the surface structures were relatively uniform. Local spectra suggested that the scanned area (approximate to 600 mu m x 1200 mu m) contained regions with serious erosion with only 2-3 mu m coating beryllium left, regions with intact marker tile, and droplets with 90% beryllium. The nonuniform erosion, droplets mainly formed by beryllium, and the possible mixture of beryllium and nickel were the major reasons that confused interpretation from broad beam EBS.

EDGE2D-EIRENE (the ‘code’) simulations show that radial electric field, Er, in the near scrape-off layer (SOL) of tokamaks can have large variations leading to a strong local E × B shear greatly exceeding that in the core region. This was pointed out in simulations of JET plasmas with varying divertor geometry, where the magnetic configuration with larger predicted near SOL Er was found to have lower H-mode power threshold, suggesting that turbulence suppression in the SOL by local E × B shear can be a player in the L–H transition physics (Delabie et al 2015 42nd EPS Conf. on Plasma Physics (Lisbon, Portugal, 22–26 June 2015) paper O3.113 (http://ocs.ciemat.es/EPS2015PAP/pdf/O3.113.pdf), Chankin et al 2017 Nucl. Mater. Energy 12 273). Further code modeling of JET plasmas by changing hydrogen isotopes (H–D–T) showed that the magnitude of the near SOL Er is lower in H cases in which the H-mode threshold power is higher (Chankin et al 2017 Plasma Phys. Control. Fusion 59 045012). From the experiment it is also known that hydrogen plasmas have poorer particle and energy confinement than deuterium plasmas, consistent with the code simulation results showing larger particle diffusion coefficients at the plasma edge, including SOL, in hydrogen plasmas (Maggi et al 2018 Plasma Phys. Control. Fusion 60 014045). All these experimental observations and code results support the hypothesis that the near SOL E × B shear can have an impact on the plasma confinement. The present work analyzes neutral ionization patterns of JET plasmas with different hydrogen isotopes in L-mode cases with fixed input power and gas puffing rate, and its impact on target electron temperature, Te, and SOL Er. The possibility of a self-feeding mechanism for the increase in the SOL Er via the interplay between poloidal E × B drift and target Te is discussed. It is also shown that reducing anomalous turbulent transport coefficients, particle diffusion and electron and ion heat conductivities, leads to higher peak target Te and larger Er, suggesting the possibility of a positive feedback loop, under an implicitly made assumption that the E × B shear in the SOL is capable of suppressing turbulence.

Hot spots on the divertor tiles at JET result in overestimation of the tile surface temperature which causes unnecessary termination of pulses. However, the appearance of hot spots can also indicate the condition of the divertor tile surfaces. To analyse the behaviour of the hot spots in the outer divertor tiles of JET, a simple image processing algorithm is developed. The algorithm isolates areas of bright pixels in the camera image and compares them to previously identified hot spots. The activity of the hot spots is then linked to values of other signals and parameters in the same time intervals. The operation of the detection algorithm was studied in a limited pulse range with high hot spot activity on the divertor tiles 5, 6 and 7. This allowed us to optimise the values of the controlling parameters. Then, the wider applicability of the method has been demonstrated by the analysis of the hot spot behaviour in a whole experimental campaign.

The impurity concentration in the tokamak divertor plasma is a necessary input for predictive scaling of divertor detachment, however direct measurements from existing tokamaks in different divertor plasma conditions are limited. To address this, we have applied a recently developed spectroscopic N II line ratio technique for measuring the N concentration in the divertor to a range of H-mode and L-mode plasma from the ASDEX Upgrade and JET tokamaks, respectively. The results from both devices show that as the power crossing the separatrix, Psep, is increased under otherwise similar core conditions (e.g. density), a higher N concentration is required to achieve the same detachment state. For example, the N concentrations at the start of detachment increase from ≈ 2% to ≈ 9% as Psep is increased from ≈ 2.5 MW to ≈ 7 MW. These results tentatively agree with scaling law predictions (e.g. Goldston et al.) motivating a further study examining the parameters which affect the N concentration required to reach detachment. Finally, the N concentrations from spectroscopy and the ratio of D and N gas valve fluxes agree within experimental uncertainty only when the vessel surfaces are fully-loaded with N.

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.

Recent achievements in studies of the effects of ICRF (Ion Cyclotron Range of Frequencies) power on the SOL (Scrape-Off Layer) and PWI (Plasma Wall Interactions) in ASDEX Upgrade (AUG), Alcator C-Mod, and JET-ILW are reviewed. Capabilities to diagnose and model the effect of DC biasing and associated impurity production at active antennas and on magnetic field connections to antennas are described. The experiments show that ICRF near-fields can lead not only to E x B convection, but also to modifications of the SOL density, which for Alcator C-Mod are limited to a narrow region near antenna. On the other hand, the SOL density distribution along with impurity sources can be tailored using local gas injection in AUG and JET-ILW with a positive effect on reduction of impurity sources. The technique of RF image current cancellation at antenna limiters was successfully applied in AUG using the 3-strap AUG antenna and extended to the 4-strap Alcator C-Mod field-aligned antenna. Multiple observations confirmed the reduction of the impact of ICRF on the SOL and on total impurity production when the ratio of the power of the central straps to the total antenna power is in the range 0.6 < P-cen / P-total < 0.8. Near-field calculations indicate that this fairly robust technique can be applied to the ITER ICRF antenna, enabling the mode of operation with reduced PWI. On the contrary, for the A2 antenna in JET-ILW the technique is hindered by RF sheaths excited at the antenna septum. Thus, in order to reduce the effect of ICRF power on PWI in a future fusion reactor, the antenna design has to be optimized along with design of plasmafacing components.

Post-disruption runaway electron (RE) beams in tokamaks with large current can cause deep melting of the vessel and are one of the major concerns for ITER operations. Consequently, a considerable effort is provided by the scientific community in order to test RE mitigation strategies. We present an overview of the results obtained at FTU and TCV controlling the current and position of RE beams to improve safety and repeatability of mitigation studies such as massive gas (MGI) and shattered pellet injections (SPI). We show that the proposed RE beam controller (REB-C) implemented at FTU and TCV is effective and that current reduction of the beam can be performed via the central solenoid reducing the energy of REs, providing an alternative/parallel mitigation strategy to MGI/SPI. Experimental results show that, meanwhile deuterium pellets injected on a fully formed RE beam are ablated but do not improve RE energy dissipation rate, heavy metals injected by a laser blow off system on low-density flat-top discharges with a high level of RE seeding seem to induce disruptions expelling REs. Instabilities during the RE beam plateau phase have shown to enhance losses of REs, expelled from the beam core. Then, with the aim of triggering instabilities to increase RE losses, an oscillating loop voltage has been tested on RE beam plateau phase at TCV revealing, for the first time, what seems to be a full conversion from runaway to ohmic current. We finally report progresses in the design of control strategies at JET in view of the incoming SPI mitigation experiments.

Geodesic acoustic modes (GAMs) may generate strong oscillations in the radial electric field and therefore are considered as a possible trigger mechanism for the L-H transition.This contribution focuses on the characterization of GAMs in JET plasmas when approaching the L-H transition aiming at understanding their possible role in triggering the transition.GAM and turbulence characteristics are measured at the plasma edge using Doppler backscattering for different plasma current and line-averaged densities.The radial location of the GAM often moves further inside when neutral beam injection is applied possibly as a response to changes in the turbulence drive.GAMs are found to have modest amplitude at the transition except for high density discharges where GAMs are stronger, suggesting that the GAM is not responsible for facilitating the transition as the L-H power threshold also increases with density in the high density branch of the L-H transition.Our results suggest that the GAM alone does not play a leading role for causing the L-H transition at JET.

The need for predictive capabilities greater than 95% with very limited false alarms are demanding requirements for reliable disruption prediction systems in tokamaks such as JET or, in the near future, ITER.The prediction of an upcoming disruption must be provided sufficiently in advance in order to apply effective disruption avoidance or mitigation actions to prevent the machine from being damaged.In this paper, following the typical machine learning workflow, a generative topographic mapping (GTM) of the operational space of JET has been built using a set of disrupted and regularly terminated discharges.In order to build the predictive model, a suitable set of dimensionless, machine-independent, physics-based features have been synthesized, which make use of 1D plasma profile information, rather than simple zero-D time series.The use of such predicting features, together with the power of the GTM in fitting the model to the data, obtains, in an unsupervised way, a 2D map of the multi-dimensional parameter space of JET, where it is possible to identify a boundary separating the region free from disruption from the disruption region.In addition to helping in operational boundaries studies, the GTM map can also be used for disruption prediction exploiting the potential of the developed GTM toolbox to monitor the discharge dynamics.Following the trajectory of a discharge on the map throughout the different regions, an alarm is triggered depending on the disruption risk of these regions.The proposed approach to predict disruptions has been evaluated on a training and an independent test set and achieves very good performance with only one tardive detection and a limited number of false detections.The warning times are suitable for avoidance purposes and, more important, the detections are consistent with physical causes and mechanisms that destabilize the plasma leading to disruptions.

In plasmas heated with deuterium beams a deficit of the expected fusion neutron rate is an indicator of the deterioration of the fast-ion confinement, caused, for instance, by magnetohydrodynamic instabilities. The capability of predicting this deficit during the discharge relies on the availability of real-time estimates of the neutron rate from NBI codes which must be fast and accurate at the same time. Therefore, the recently developed real-time RABBIT code for neutral beam injection (NBI) simulations has been extended to output the distribution function and calculate the neutron emission. After the description of this newly installed diagnostics in RABBIT, benchmarks with NUBEAM, a massively used and validated Monte Carlo NBI solver, are discussed on ASDEX-Upgrade and JET cases. A first application for control-room intershot analysis on DIII-D is presented, and the results are compared on a large database with a slower NUBEAM analysis. Further application possibilities, e.g. for real-time control of Alfvén eigenmodes, are outlined.

Preliminary results investigating the microstructure, bonding and effect of beryllium oxide formation on retention in the JET ITER-like wall beryllium tiles, are presented. The tiles have been investigated by several techniques: Scanning Electron Microscopy (SEM) equipped with Energy Dispersive X-ray (EDX), Transmission Electron microscopy (TEM) equipped with EDX and Electron Energy Loss Spectroscopy (EELS), Raman Spectroscopy and Thermal Desorption Spectroscopy (TDS). This paper focuses on results from melted materials of the dump plate tiles in JET. From our results and the literature, it is concluded, beryllium can form micron deep oxide islands contrary to the nanometric oxides predicted under vacuum conditions. The deepest oxides analyzed were up to 2-micron thicknesses. The beryllium Deuteroxide (BeOxDy) bond was found with Raman Spectroscopy. Application of EELS confirmed the oxide presence and stoichiometry. Literature suggests these oxides form at temperatures greater than 700 °C where self-diffusion of beryllium ions through the surface oxide layer can occur. Further oxidation is made possible between oxygen plasma impurities and the beryllium ions now present at the wall surface. Under Ultra High Vacuum (UHV) nanometric Beryllium oxide layers are formed and passivate at room temperature. After continual cyclic heating (to the point of melt formation) in the presence of oxygen impurities from the plasma, oxide growth to the levels seen experimentally (approximately two microns) is proposed. This retention mechanism is not considered to contribute dramatically to overall retention in JET, due to low levels of melt formation. However, this mechanism, thought the result of operation environment and melt formation, could be of wider concern to ITER, dependent on wall temperatures.

An effect due to fusion born triton production has been observed in JET high-performance deuterium plasma discharges with neutral beam injection (NBI) and H-minority ion cyclotron range of frequencies (ICRF) heating, using DD and deuterium tritium (DT) neutron spectrometry as well as fusion product loss measurements. The observations show that a decrease of the second harmonic omega = 2 omega(cD) enhancement of the DD neutron rate correlates with an increase of the triton burnup rate. An acceleration of tritons due to absorbing ICRH power at the third harmonic omega = 3 omega(cT) has been observed. This effect could indicate a redistribution of ICRH power absorption at omega approximate to omega(cH) = 2 omega(cD) = 3 omega(cT) with increasing triton concentration at the ion cyclotron resonance layer. Also, the reduction of the second harmonic enhancement of the DD neutron rate can be caused by burning of the accelerated deuterium as the tritium concentration grows. This is an extremely non-linear process as both mechanisms intensify with triton concentration. It determines the necessity to consider the ICRH power absorption omega = 3 omega(cT) in modelling of high-performance deuterium discharges with simultaneous NBI and H-minority ICRF heating as well as the assessment of enhanced burnup of ICRF accelerated deuterium for the development of high-performance plasma scenarios and DT fusion rate predictions.

The measurement of the energy spectra and densities of α-particles and other fast ions are part of the ITER measurement requirements, highlighting the importance of energy-resolved energetic-particle measurements for the mission of ITER. However, it has been found in recent years that the velocity-space interrogation regions of the foreseen energetic-particle diagnostics do not allow these measurements directly. We will demonstrate this for γ-ray spectroscopy (GRS), collective Thomson scattering (CTS), neutron emission spectroscopy and fast-ion Dα spectroscopy by invoking energy and momentum conservation in each case, highlighting analogies and differences between the different diagnostic velocity-space sensitivities. Nevertheless, energy spectra and densities can be inferred by velocity-space tomography which we demonstrate using measurements at JET and ASDEX Upgrade. The measured energy spectra agree well with corresponding simulations. At ITER, α-particle energy spectra and densities can be inferred for energies larger than 1.7 MeV by velocity-space tomography based on GRS and CTS. Further, assuming isotropy of the α-particles in velocity space, their energy spectra and densities can be inferred by 1D inversion of spectral single-detector measurements down to about 300 keV by CTS. The α-particle density can also be found by fitting a model to the CTS measurements assuming the α-particle distribution to be an isotropic slowing-down distribution.

The gyrokinetic toroidal code (GTC) has been used to study toroidal Alfvén eigenmodes (TAEs) in high-performance plasmas. Experiments performed at the Joint European Torus (JET), where TAEs were driven by energetic particles arising from neutral beams, ion cyclotron resonant heating, and resonantly excited by dedicated external antennas, have been simulated. Modes driven by populations of energetic particles are observed, matching the TAE frequency seen with magnetic probes in JET experiments. A synthetic antenna, composed of one toroidal and two neighboring poloidal harmonics has been used to probe the modes' damping rates and quantify mechanisms for this damping in GTC simulations. This method was also applied to frequency and damping rate measurements of stable TAEs made by the Alfvén eigenmode active diagnostic in these discharges.