The aim of this paper is to present a signal processing algorithm that, applied to the raw Locked Mode signal, allows us to obtain a disruption indicator in principle exploitable on different tokamaks. A common definition of such an indicator for different machines would facilitate the development of portable systems for disruption prediction, which is becoming of increasingly importance for the next tokamak generations. Moreover, the indicator allows us to overcome some intrinsic problems in the diagnostic system such as drift and offset. The behavior of the proposed indicator as disruption predictor, based on crossing optimized thresholds of the signal amplitude, has been analyzed using data of both JET and ASDEX Upgrade experiments. A thorough analysis of the disruption prediction performance shows how the indicator is able to recover some missed and tardy detections of the raw signal. Moreover, it intervenes and corrects premature or even wrong alarms due to, e.g., drifts and/or offsets.

Developing successful tokamak operation scenarios, as well as confident extrapolation of present-day knowledge requires a rigorous understanding of plasma turbulence, which largely determines the quality of the confinement. In particular, accurate particle transport predictions are essential due to the strong dependence of fusion power or bootstrap current on the particle density details. Here, gyrokinetic turbulence simulations are performed with physics inputs taken from a JET power scan, for which a relatively weak degradation of energy confinement and a significant density peaking is obtained with increasing input power. This way physics parameters that lead to such increase in the density peaking shall be elucidated. While well-known candidates, such as the collisionality, previously found in other studies are also recovered in this study, it is furthermore found that edge E x B shearing may adopt a crucial role by enhancing the inward pinch. These results may indicate that a plasma with rotational shear could develop a stronger density peaking as compared to a non-rotating one, because its inward convection is increased compared to the outward diffusive particle flux as long as this rotation has a significant on E x B flow shear stabilization. The possibly significant implications for future devices, which will exhibit much less torque compared to present day experiments, are discussed.

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.

For many years, machine learning tools have proved to be very powerful disruption predictors in tokamaks. On the other hand, the vast majority of the techniques deployed assume that the input data is independent and is sampled from exactly the same probability distribution for the training set, the test set and the final real time deployment. This hypothesis is certainly not verified in practice, since the experimental programmes evolve quite rapidly, resulting typically in ageing of the predictors and consequent suboptimal performance. This paper describes various adaptive training strategies that have been tested to maintain the performance of disruption predictors in non-stationary conditions. The proposed approaches have been implemented using new ensembles of classifiers, explicitly developed for the present application. The improvements in performance are unquestionable and, given the difficulties encountered so far in translating predictors from one device to another, the proposed adaptive methods from scratch can therefore be considered a useful option in the arsenal of alternatives envisaged for the next generation of devices, particularly at the very beginning of their operation.

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 analysed 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 mu m were indicated, but typical values were on the order of a few hundred nm. 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 mu m 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.

This paper outlines an approach towards improved rigour in tokamak turbulence transport model validation within integrated modelling. Gaussian process regression (GPR) techniques were applied for profile fitting during the preparation of integrated modelling simulations allowing for rigourous sensitivity tests of prescribed initial and boundary conditions as both fit and derivative uncertainties are provided. This was demonstrated by a JETTO integrated modelling simulation of the JET ITER-like-wall H-mode baseline discharge #92436 with the QuaLiKiz quasilinear turbulent transport model, which is the subject of extrapolation towards a deuterium-tritium plasma. The simulation simultaneously evaluates the time evolution of heat, particle, and momentum fluxes over similar to 10 confinement times, with a simulation boundary condition at rho(tor) = 0.85. Routine inclusion of momentum transport prediction in multi-channel flux-driven transport modelling is not standard and is facilitated here by recent developments within the QuaLiKiz model. Excellent agreement was achieved between the fitted and simulated profiles for n(e), T-e, T-i, and Omega(tor) within 2 sigma, but the simulation underpredicts the mid-radius Ti and overpredicts the core n(e) and T-e profiles for this discharge. Despite this, it was shown that this approach is capable of deriving reasonable inputs, including derivative quantities, to tokamak models from experimental data. Furthermore, multiple figures-of-merit were defined to quantitatively assess the agreement of integrated modelling predictions to experimental data within the GPR profile fitting framework.

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.