In 2021 JET exploited its unique capabilities to operate with T and D-T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D-T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements-new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors-as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D-T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D-T plasmas to date and expanding our understanding of isotopes and D-T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D-T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D-T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D-T operations is presented in (King et al Nucl. Fusion submitted).

This paper presents results of JET ITER-like wall L-mode experiments in hydrogen and deuterium (D) plasmas, dedicated to the study of the isotope dependence of ion heat transport by determination of the ion critical gradient and stiffness by varying the ion cyclotron resonance heating power deposition. When no strong role of fast ions in the plasma core is expected, the main difference between the two isotope plasmas is determined by the plasma edge and the core behavior is consistent with a gyro-Bohm scaling. When the heating power (and the fast ion pressure) is increased, in addition to the difference in the edge region, also the plasma core shows substantial changes. The stabilization of ion heat transport by fast ions, clearly visible in D plasmas, appears to be weaker in H plasmas, resulting in a higher ion heat flux in H with apparent anti-gyro-Bohm mass scaling. The difference is found to be caused by the different fast ion pressure between H and D plasmas, related to the heating power settings and to the different fast ion slowing down time, and is completely accounted for in non-linear gyrokinetic simulations. The application of the TGLF quasi-linear model to this set of data is also discussed.

The present paper offers an overview of the potential of ion cyclotron resonance heating (ICRH) or radio frequency heating for the DEMO machine. It is found that various suitable heating schemes are available. Similar to ITER and in view of the limited bandwidth of about 10 MHz that can be achieved to ensure optimal functioning of the launcher, it is proposed to make core second harmonic tritium heating the key ion heating scheme, assisted by fundamental cyclotron heating He-3 in the early phase of the discharge; for the present design of DEMO-with a static magnetic field strength of B-o = 5.855 T-that places the T and 3He layers in the core for f = 60 MHz and suggests centering the bandwidth around that main operating frequency. In line with earlier studies for hot, dense plasmas in large-size magnetic confinement machines, it is shown that good single pass absorption is achieved but that the size as well as the operating density and temperature of the machine cause the electrons to absorb a non-negligible fraction of the power away from the core when core ion heating is aimed at. Current drive and alternative heating options are briefly discussed and a dedicated computation is done for the traveling wave antenna, proposed for DEMO in view of its compatibility with substantial antenna-plasma distances. The various tasks that ICRH can fulfill are briefly listed. Finally, the impact of transport and the sensitivity of the obtained results to changes in the machine parameters is commented on.

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.

The First Mirror Test for ITER has been carried out in JET with mirrors exposed during: (i) the third ILW campaign (ILW-3, 2015-2016, 23.33 h plasma) and (ii) all three campaigns, i.e. ILW-1 to ILW-3: 2011-2016, 63,52 h in total. All mirrors from main chamber wall show no significant changes of the total reflectivity from the initial value and the diffuse reflectivity does not exceed 3% in the spectral range above 500 nm. The modified layer on surface has very small amount of impurities such as D, Be, C, N, O and Ni. All mirrors from the divertor (inner, outer, base under the bulk W tile) lost reflectivity by 20-80% due to the beryllium-rich deposition also containing D, C, N, O, Ni and W. In the inner divertor N reaches 5 x 10(17) cm(-2), W is up to 4.3 x 10(17) cm(-2), while the content of Ni is the greatest in the outer divertor: 3.8 x 10(17) cm(-2). Oxygen-18 used as the tracer in experiments at the end of ILW-3 has been detected at the level of 1.1 x 10(16) cm(-2). The thickness of deposited layer is in the range of 90 nm to 900 nm. The layer growth rate in the base (2.7 pm s(-1)) and inner divertor is proportional to the exposure time when a single campaign and all three are compared. In a few cases, on mirrors located at the cassette mouth, flaking of deposits and erosion occurred.

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.

The erosion of tungsten (W), induced by the bombardment of plasma and impurity particles, determines the lifetime of plasma-facing components as well as impacting on plasma performance by the influx of W into the confined region. The screening of W by the divertor and the transport of W in the plasma determines largely the W content in the plasma core, but the W source strength itself has a vital impact on this process. The JET tokamak experiment provides access to a large set of W erosion-determining parameters and permits a detailed description of the W source in the divertor closest to the ITER one: (i) effective sputtering yields and fluxes as function of impact energy of intrinsic (Be, C) and extrinsic (Ne, N) impurities as well as hydrogenic isotopes (H, D) are determined and predictions for the tritium (T) isotope are made. This includes the quantification of intra- and inter-edge localised mode (ELM) contributions to the total W source in H-mode plasmas which vary owing to the complex flux compositions and energy distributions in the corresponding phases. The sputtering threshold behaviour and the spectroscopic composition analysis provides an insight in the dominating species and plasma phases causing W erosion. (ii) The interplay between the net and gross W erosion source is discussed considering (prompt) re-deposition, thus, the immediate return of W ions back to the surface due to their large Larmor radius, and surface roughness, thus, the difference between smooth bulk-W and rough W-coating components used in the JET divertor. Both effects impact on the balance equation of local W erosion and deposition. (iii) Post-mortem analysis reveals the net erosion/deposition pattern and the W migration paths over long periods of plasma operation identifying the net W transport to remote areas. This W transport is related to the divertor plasma regime, e.g. attached operation with high impact energies of impinging particles or detached operation, as well as to the applied magnetic configuration in the divertor, e.g. close divertor with good geometrical screening of W or open divertor configuration with poor screening. JET equipped with the ITER-like wall (ILW) provided unique access to the net W erosion rate within a series of 151 subsequent H-mode discharges (magnetic field: T, plasma current: MA, auxiliary power MW) in one magnetic configuration accumulating 900 s of plasma with particle fluences in the range of 5– in the semi-detached inner and attached outer divertor leg. The comparison of W spectroscopy in the intra-ELM and inter-ELM phases with post-mortem analysis of W marker tiles provides a set of gross and net W erosion values at the outer target plate. ERO code simulations are applied to both divertor leg conditions and reproduce the erosion/deposition pattern as well as confirm the high experimentally observed prompt W re-deposition factors of more than 95% in the intra- and inter-ELM phase of the unseeded deuterium H-mode plasma. Conclusions to the expected divertor conditions in ITER as well as to the JET operation in the DT plasma mixture are drawn on basis of this unique benchmark experiment.

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.

ERO is a 3D Monte-Carlo impurity transport and plasma-surface interaction code. In 2011 it was applied for the ITER first wall (FW) life time predictions [1] (critical blanket module BM11). After that the same code was significantly improved during its application to existing fusion-relevant plasma devices: the tokamak JET equipped with an ITER-like wall and linear plasma device PISCES-B. This has allowed testing the sputtering data for beryllium (Be) and showing that the "ERO-min" fit based on the large (50%) deuterium (D) surface content is well suitable for plasma-wetted areas (D plasma). The improved procedure for calculating of the effective sputtering yields for each location along the plasma-facing surface using the recently developed semi-analytical sheath approach was validated. The re-evaluation of the effective yields for BM11 following the similar revisit of the JET data has indicated significant increase of erosion and motivated the current re-visit of ERO simulations.

In this contribution we present OES measurements in the JET tokamak of the deuterated NH (ND) radical and the correlation between results of those experiments and measurement of ammonia production. The observation region covers most of the divertor and its outer throat. Measurements are performed in different magnetic configurations. The results include temporal and spatial dependence of the molecular emission intensity and study of the emission band shape (vibrational and rotational temperatures) during different JET pulses, with or without nitrogen seeding. Results are a step towards the understanding of nitrogen-containing molecule creation and destruction in the divertor plasma.

Measurements of neutron streaming through penetrations in biological shields are being carried out at JET fusion device by means of thermoluminescence (TL) detectors with the objective to validate the neutronics codes and nuclear data widely applied in ITER nuclear analyses in a real fusion environment. TLDs response due to the neutron component of the radiation field is related to the neutron fluence in a well-defined neutron energy spectrum. Therefore, a TLDs calibration in real fusion radiation fields is necessary to allow neutron fluence from TL measurements at JET to be more precisely calculated. Hence, a MCP-N and MCP-7 TLDs produced at the IFJ PAN in Kraków were calibrated at the ENEA facilities of Frascati and Casaccia laboratories. The obtained results have been analysed and new calibration factors are proposed. The detection system based on TLDs developed and calibrated for JET experiments can then be generally applied not only to fusion neutron fields, but also to ITER to monitor the neutron fluence outside the biological shield.

Experimental observations in JET tokamak plasmas and gyrokinetic simulations point to an important role, for electron heat transport, of electron-scale instabilities and of their interaction with ion-scale instabilities. Since these effects are maximized for strong electron heating and ion-scale modes close to marginal stability, these findings are of high relevance for ITER plasmas, featuring both conditions. Gyrokinetic and quasi-linear transport models accounting for multi-scale effects are assessed against JET experimental results.

This paper analyses the final phase of hybrid scenario discharges at JET, the reduction of auxiliary heating towards finally the Ohmic phase.The here considered Ohmic phase is mostly still in the current flattop but may also be in the current ramp down.For this purpose a database is created of 54 parameters in 7 phases distributed in time of the discharge.It is found that the occurrence of a locked mode is in most cases preceded by a radiation peaking after the main heating phase either in a low power phase and/or in the Ohmic phase.To gain insight on the importance of different parameters in this process a correlation analysis to the radiation peaking in the Ohmic phase is done.The first finding is that the further away in time the analysed phases are the less the correlation is.This means in the end that a good termination scenario might also be able to terminate unhealthy plasmas safely.The second finding is that remaining impurities in the plasma after reducing the heating power in the termination phase are the most important reason for generating a locked mode which can lead to a disruption.

Plasma tomography consists of reconstructing a two-dimensional radiation profile of a poloidal cross section of a fusion device based on line-integrated measurements along several lines of sight. The reconstruction process is computationally intensive, and in practice, only a few reconstructions are usually computed per pulse. In this work, we trained a deep neural network based on a large collection of sample tomograms that have been produced at JET over several years. Once trained, the network is able to reproduce those results with high accuracy. More importantly, it can compute all the tomographic reconstructions for a given pulse in just a few seconds. This makes it possible to visualize several phenomena—such as plasma heating, disruptions, and impurity transport—over the course of the entire pulse.

This paper presents the first application to real JET data of the new equilibrium code NICE which enables the consistent resolution of the inverse equilibrium reconstruction problem in the framework of non-linear free-boundary equilibrium coupled to the Stokes model equation for polarimetry. The conducted numerical experiments enable first of all to validate NICE by comparing it to the well-established EFIT code on 4 selected high performance shots. Secondly the results indicate that the fit to polarimetry measurements clearly benefits from the use of Stokes vector measurements compared to the classical case of Faraday measurements, and that the reconstructed p' and ff' profiles are better constrained with smaller error bars and are closer to the profiles reconstructed by EFTM, the EFIT JET code using internal MSE constraints.

Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machine-specific minimum current quench time, tau(CQ); (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, < f(h)>. These data reveal that the normalized rotation duration, t(rot)/tau(CQ), and the average rotation velocity, < v(h)>, are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above < f(h)> = 20 Hz is to be expected. More importantly, depending on the projected value of tau(CQ) in ITER, substantial rotation could also occur in the resonant frequency range of 6-20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out.

Progress in understanding lower hybrid current drive (LHCD) at high density has been made through experiments and modeling, which is encouraging given the need for an efficient off-axis current profile control technique in burning plasma. By reducing the wall recycling of neutrals, the edge temperature is increased and the effect of parametric instability (PI) and collisional absorption (CA) is reduced, which is beneficial for increasing the current drive efficiency. Strong single pass absorption is preferred to prevent CA and high LH operating frequency is essential for wave propagation to the core region at high density, presumably to mitigate the effect of PI. The dimensionless parameter that characterizes LH wave accessibility and wave refraction for the experiments in this joint study is shown to bracket the region in parameter space where ITER LHCD experiments will operate in the steady state scenario phase. Further joint experiments and cross modeling are necessary to understand the LHCD physics in weak damping regimes which would increase confidence in predictions for ITER where the absorption is expected to be strong.

In preparation for the upcoming deuterium-tritium campaign at the JET tokamak, the ex-vessel fission chamber neutron diagnostics and the neutron activation system will be calibrated in absolute terms at 14 MeV neutron energy, to a required accuracy of less than 10%. Two nominally identical DT neutron generators were chosen as the calibration sources, both of which were fully calibrated and characterized at the UK’s National Physical Laboratory. The neutron activation method was adopted as a complementary method for the purpose of determining the absolute value of the neutron yield from the neutron generators and to provide a means of cross check for the active detection methods being employed. The work being presented here shows the derivation of the neutron emission rate from the neutron generators based upon experimental activation foil measurements.

Several important ELM control techniques are in large part motivated by the empirically observed inverse relationship between average ELM energy loss and ELM frequency in a plasma. However, to ensure a reliable effect on the energy released by the ELMs, it is important that this relation is verified for individual ELM events. Therefore, in this work the relation between ELM energy loss (WELM) and waiting time (ΔtELM) is investigated for individual ELMs in a set of ITER-like wall plasmas in JET. A comparison is made with the results from a set of carbon-wall and nitrogen-seeded ITER-like wall JET plasmas. It is found that the correlation between W ELM and ΔtELM for individual ELMs varies from strongly positive to zero. Furthermore, the effect of the extended collapse phase often accompanying ELMs from unseeded JET ILW plasmas and referred to as the slow transport event (STE) is studied on the distribution of ELM durations, and on the correlation between W ELM and ΔtELM . A high correlation between W ELM and ΔtELM , comparable to CW plasmas is only found in nitrogen-seeded ILW plasmas. Finally, a regression analysis is performed using plasma engineering parameters as predictors for determining the region of the plasma operational space with a high correlation between W ELM and ΔtELM .

In metallic machines ICRH heating is playing an increasingly important role. One of its most recent applications on the Joint Europena Torus (JET) is sawtooth control by ICRH modulation, for avoiding triggering dangerous neo-classical tearing modes (NTMs) and counteracting impurity accumulation. Some of the main difficulties of these experiments are the assessment of the synchronization efficiency and the understanding of the main physical mechanisms at play. In this paper, three independent classes of statistical indicators are introduced to address these issues: Recurrence Plots, Convergent Cross Mapping and Transfer Entropy. The application to JET experiments with the ILW shows that the proposed indicators agree quite well among themselves and provide sound estimates of the efficiency of the synchronisation scheme investigated. They also support, with a shot to shot basis analysis and an estimate of the uncertainties, the interpretation that the fast ions play a fundamental role in the stabilization of the sawteeth, in both L and H mode. Proposals for experiments to be carried out in the future to consolidate the interpretation of the results are discussed.