In nuclear fusion reactors, plasmas are heated to very high temperatures of more than 100 million kelvin and, in so-called tokamaks, they are confined by magnetic fields in the shape of a torus. Light nuclei, such as deuterium and tritium, undergo a fusion reaction that releases energy, making fusion a promising option for a sustainable and clean energy source. Tokamak plasmas, however, are prone to disruptions as a result of a sudden collapse of the system terminating the fusion reactions. As disruptions lead to an abrupt loss of confinement, they can cause irreversible damage to present-day fusion devices and are expected to have a more devastating effect in future devices. Disruptions expected in the next-generation tokamak, ITER, for example, could cause electromagnetic forces larger than the weight of an Airbus A380. Furthermore, the thermal loads in such an event could exceed the melting threshold of the most resistant state-of-the-art materials by more than an order of magnitude. To prevent disruptions or at least mitigate their detrimental effects, empirical models obtained with artificial intelligence methods, of which an overview is given here, are commonly employed to predict their occurrence—and ideally give enough time to introduce counteracting measures.

Alpha particles with energies on the order of megaelectronvolts will be the main source of plasma heating in future magnetic confinement fusion reactors. Instead of heating fuel ions, most of the energy of alpha particles is transferred to electrons in the plasma. Furthermore, alpha particles can also excite Alfvénic instabilities, which were previously considered to be detrimental to the performance of the fusion device. Here we report improved thermal ion confinement in the presence of megaelectronvolts ions and strong fast ion-driven Alfvénic instabilities in recent experiments on the Joint European Torus. Detailed transport analysis of these experiments reveals turbulence suppression through a complex multi-scale mechanism that generates large-scale zonal flows. This holds promise for more economical operation of fusion reactors with dominant alpha particle heating and ultimately cheaper fusion electricity.

The JET 2019-2020 scientific and technological programme exploited the results of years of concerted scientific and engineering work, including the ITER-like wall (ILW: Be wall and W divertor) installed in 2010, improved diagnostic capabilities now fully available, a major neutral beam injection upgrade providing record power in 2019-2020, and tested the technical and procedural preparation for safe operation with tritium. Research along three complementary axes yielded a wealth of new results. Firstly, the JET plasma programme delivered scenarios suitable for high fusion power and alpha particle (alpha) physics in the coming D-T campaign (DTE2), with record sustained neutron rates, as well as plasmas for clarifying the impact of isotope mass on plasma core, edge and plasma-wall interactions, and for ITER pre-fusion power operation. The efficacy of the newly installed shattered pellet injector for mitigating disruption forces and runaway electrons was demonstrated. Secondly, research on the consequences of long-term exposure to JET-ILW plasma was completed, with emphasis on wall damage and fuel retention, and with analyses of wall materials and dust particles that will help validate assumptions and codes for design and operation of ITER and DEMO. Thirdly, the nuclear technology programme aiming to deliver maximum technological return from operations in D, T and D-T benefited from the highest D-D neutron yield in years, securing results for validating radiation transport and activation codes, and nuclear data for ITER.

A global heat flux model based on a fractional derivative of plasma pressure is proposed for the heat transport in fusion plasmas. The degree of the fractional derivative of the heat flux, α, is defined through the power balance analysis of the steady state. The model was used to obtain the experimental values of α for a large database of the Joint European Torus (JET) carbon-wall as well as ITER like-wall plasmas. The fractional degrees of the electron heat flux are found to be α<2, for all the selected pulses in the database, suggesting a deviation from the diffusive paradigm. Moreover, the results show that as the volume integrated input power is increased, the fractional degree of the electron heat flux converges to α∼0.8, indicating a global scaling between the net heating and the pressure profile in the high-power JET plasmas. The model is expected to provide insight into the proper kinetic description for the fusion plasmas and improve the accuracy of the heat transport predictions.

A power-balance model, with radiation losses from impurities and neutrals, gives a unified description of the density limit (DL) of the stellarator, the L-mode tokamak, and the reversed field pinch (RFP). The model predicts a Sudo-like scaling for the stellarator, a Greenwald- like scaling, alpha I-p(8/9), for the RFP and the ohmic tokamak, a mixed scaling, alpha (PIp4/9)-I-4/9, for the additionally heated L-mode tokamak. In a previous paper (Zanca et al 2017 Nucl. Fusion 57 056010) the model was compared with ohmic tokamak, RFP and stellarator experiments. Here, we address the issue of the DL dependence on heating power in the L-mode tokamak. Experimental data from high-density disrupted L-mode discharges performed at JET, as well as in other machines, arc taken as a term of comparison. The model fits the observed maximum densities better than the pure Greenwald limit.

Block-structured mesh generation techniques have been well addressed in the CFD community for automobile and aerospace studies, and their applicability to magnetic fusion is highly relevant, due to the complexity of the plasma-facing wall structures inside a tokamak device. Typically applied to non-linear simulations of MHD instabilities relevant to magnetically confined fusion, the JOREK code was originally developed with a 2D grid composed of isoparametric bi-cubic Bezier finite elements, that are aligned to the magnetic equilibrium of tokamak plasmas (the third dimension being represented by Fourier harmonics). To improve the applicability of these simulations, the grid-generator has been generalised to provide a robust extension method, using a block-structured mesh approach, which allows the simulations of arbitrary domains of tokamak vacuum vessels. Such boundary-aligned grids require the adaptation of boundary conditions along the edge of the new domain. Demonstrative non-linear simulations of plasma edge instabilities are presented to validate the robustness of the new grid, and future potential physics applications for tokamak plasmas are discussed. The methods presented here may be of interest to the wider community, beyond tokamak physics, wherever imposing arbitrary boundaries to quadrilateral finite elements is required.

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, P-sep, 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 approximate to 2% to approximate to 9% as P-sep, is increased from approximate to 2.5 MW to approximate 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.

Inconel-600 blocks and stainless steel covers for quartz microbalance crystals from remote corners in the JET-ILW divertor were studied with time-of-flight elastic recoil detection analysis and nuclear reaction analysis to obtain information about the areal densities and depth profiles of elements present in deposited material layers. Surface morphology and the composition of dust particles were examined with scanning electron microscopy and energy-dispersive X-ray spectroscopy. The analyzed components were present in JET during three ITER-like wall campaigns between 2010 and 2017. Deposited layers had a stratified structure, primarily made up of beryllium, carbon and oxygen with varying atomic fractions of deuterium, up to more than 20%. The range of carbon transport from the ribs of the divertor carrier was limited to a few centimeters, and carbon/deuterium co-deposition was indicated on the Inconel blocks. High atomic fractions of deuterium were also found in almost carbon-free layers on the quartz microbalance covers. Layer thicknesses up to more than 1 micrometer were indicated, but typical values were on the order of a few hundred nanometers. Chromium, iron and nickel fractions were less than or around 1% at layer surfaces while increasing close to the layer-substrate interface. The tungsten fraction depended on the proximity of the plasma strike point to the divertor corners. Particles of tungsten, molybdenum and copper with sizes less than or around 1 micrometer were found. Nitrogen, argon and neon were present after plasma edge cooling and disruption mitigation. Oxygen-18 was found on component surfaces after injection, indicating in-vessel oxidation. Compensation of elastic recoil detection data for detection efficiency and ion-induced release of deuterium during the measurement gave quantitative agreement with nuclear reaction analysis, which strengthens the validity of the results.

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 analysis of the polarimetry measurements has the aim of validating models (De Marco and Segre 1972 Plasma Phys. 14 245), with a careful attention to the clarification of their limits of application. In this paper a new approximation method is introduced, the so-called special constant Omega direction (SCOD), which gives an analytical solution to the polarimetry exact Stokes model equations. The available approximate solutions (including SCOD) of the polarimetry propagation equations are presented, compared and their application limits determined, using a reference tokamak configuration, which is a simplified equilibrium for a circular tokamak. The SCOD approximation is compared successfully to the Stokes model in the context also of equilibria evaluated for two JET discharges. The approximation methods are analytical or very simple mathematical expressions which can also be used in equilibrium codes for their optimization.