Recent JET D-T campaigns opened the possibility of unique isotope studies to investigate the L-H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L-H transitions in the low and high-density branches of JET NBI-heated D-T plasmas. As discussed in the paper, L-H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux (Qi). We present the first power balance analysis of JET NBI-heated D-T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D-T Qi, with a lower Qi in D-T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L-H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Qi deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Qi and the power threshold minimum at JET. At low plasma density, NBI power dominates Qi, while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D-T NBI-heated plasmas at JET.
An important step on the way to future fusion power plants was the 2021 deuterium-tritium experimental campaign (DTE2) at the Joint European Torus (JET), in which crucial DT physics was investigated. In this study, we have reconstructed the fast-ion deuterium distribution function in JET discharge 99971 which broke the former fusion energy record. It is the first time that the fast-ion distribution has been reconstructed from experimental data in a DT discharge. The reconstruction shows that the fast-ion deuterium distribution is anisotropic, with a bias towards co-going ions (p > 0). The fast-ion deuterium distribution likely peaks in energy (E) at around E similar to 60-70 keV and has a marginal high-energy tail ( E greater than or similar to 180 keV). Furthermore, an orbit analysis shows that the fast-ion distribution is composed of mostly co-passing orbits ( 50%), trapped orbits ( 21%) and counter-passing orbits ( 27%), as well as a small population of potato orbits ( 1.7%) and counter-stagnation orbits ( 0.3%). The orbit-type constituents of the neutron measurements are distributed in similar fractions.
This paper presents a study on the dependence of the ion temperature stiffness on the plasma main ion isotope mass in JET ITER-like wall and C wall discharges. To this aim, a database of H, D and T shots is analyzed, including new dedicated shots, comparing experiments with lower and higher power injected by the NBI system. In order to characterize the turbulence dependence on the isotope mass, three of these discharges (two in T and one in D) with same external heating scheme are studied in detail and interpreted with gyrokinetic linear and nonlinear simulations. The analysis is performed at fixed radius rho tor=0.33, selected in order to maximize the electromagnetic stabilizing effects on turbulence, both from thermal and suprathermal particles. The experimental results show a clear ion temperature stiffness reduction when heavier isotopes are considered, thus moving from H to D to T, which is attributed to an increasing thermal electromagnetic stabilization with increasing main isotope mass.
Understanding the physics of the H-mode pedestal structure and dynamics and predicting its performance in future fusion reactors is key to reducing the uncertainties associated with the realization of burning plasma conditions. ITER and future fusion reactors will operate with mixed deuterium-tritium (D-T) plasmas. Perhaps the greatest challenge to date related to understanding the impact of hydrogen isotope mass on plasma transport lies in the processes governing the plasma edge region and H-mode pedestal. This paper, based on a review talk at the 19th H-mode workshop (Mito, Japan, September 2024), reviews recent progress in understanding the impact of isotope mass on pedestal structure and inter-ELM transport. Recent T and D-T experiments on JET have confirmed a strong mass dependence on the pedestal density and provided new and unique results. State of the art pedestal diagnostics were employed to measure the pedestal structure in D, D-T and T and the most advanced theoretical models were used to identify the physics mechanisms linked to the dependence of pedestal stability and inter-ELM transport on isotope mass. Different modelling and analysis workflows are consistent in finding that changes in pedestal stability and in inter-ELM transport with isotope mass are required to explain the isotope mass dependencies of pedestal density and temperature. The competing levels of turbulent and neoclassical transport with varying isotope mass has been quantified using the gyrokinetic code GENE for selected JET pedestals. However, progress in this area is hindered by the difficulty to obtain accurate measurements of the edge particle source. The isotope dependence of neutral fuelling in the H-mode pedestal region was examined quantitatively in DIII-D experiments using direct spectroscopic measurements of neutral hydrogen penetration in H vs D plasmas. Projecting to D-T H-modes in future reactors, which will be largely opaque to neutrals, the impact of isotope mass dependence of edge neutral fuelling is expected to be negligible, while the pedestal density structure will depend primarily on its transport properties, in turn affecting pedestal stability.
This work studies the influence of radio frequency (RF) waves in the ion cyclotron resonance heating (ICRH) range of frequencies on fusion alphas during the recent JET D-T campaign. Fusion alphas from D-T reactions are created with energies of about 3.5 MeV and therefore have significant Doppler shifts enabling synergistic interactions between them and RF waves at a broad range of frequencies, including the ones foreseen for future fusion machines in ITER (Schneider et al 2021 Nucl. Fusion 61 126058) and SPARC (Creely et al 2020 J. Plasma Phys. 86 865860502). Resonant interactions between RF waves and alphas, also called synergistic effects, will modify the alpha distribution and ultimately will have an impact on alpha orbit losses and heating. Data from JET 3.43 T/2.3 MA pulses based on the hybrid scenario (Hobirk et al 2023 Nucl. Fusion ; Hobirk et al 29th IAEA FEC23 Conf. ( 16–21 October 2023 ); Challis et al 48th EPS Conf. on Plasma Physics ( 27 June–1 July 2022 ) during the DTE2 campaign (Maggi et al 2023 Nucl. Fusion )) were used for the analysis in this study. The impact of synergistic effects on alpha orbit losses and alpha heating are assessed. The conclusions are based on the analysis of experimental data for fast alpha losses, i.e. measurements from neutral particle analyser (NPA), fast ion losses scintillator detector, Faraday cups (FCs), and TRANSP (Hawryluk et al 1980 Physics of Plasmas Close to Thermonuclear Conditions vol 1 (CEC) pp 19–46) simulations. Experimental data and TRANSP analysis indicates that there are indeed changes in the alpha distribution function (DF) due to interaction with RF waves. Data from the NPA show increased 4 He flux in the range from a few hundred keV up to 800 keV for pulses with RF power, while TRANSP clearly shows modifications in the fast alpha DF for these energies. Data from the scintillator detector and the FCs were compared for pulses with and without ICRH power and versus cases with enhanced alpha losses due to MHD activity. The trends from these diagnostics consistently show no additional alpha losses due to interaction with RF waves. TRANSP predictions for the impact of ynergistic effects on alpha heating show up to a 42% increase in alpha electron heating and up to a 25% increase in alpha ion heating. These effects, however, become negligibly small, less than 1%, when alpha heating is compared to the total auxiliary heating power in the investigated JET pulses.
A baseline scenario of deuterium-tritium (D-T) plasma with peripheral high-field-side fuelling pellets has been produced in JET in order to mimic the situation in ITER. The isotope mix ratio is controlled in order to target the value of 50%-50% by a combination of tritium gas puffing and deuterium pellet injection. Multiple factors controlling the fuelling efficiency of individual pellets are analysed, with the following findings: (1) prompt particle losses due to pellet-triggered edge-localised modes (ELMs) are detected, (2) the plasmoid drift velocity might be smaller than that predicted by simulation, (3) post-pellet particle loss is controlled by transient phases with ELMs.The overall pellet particle flux normalised to the heat flux is similar to that in previous pellet fuelling experiments in AUG and JET.
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).
Providing stable and clean energy sources is a necessity for the increasing demands of humanity. Energy produced by fusion reactions, in particular in tokamaks, is a promising path towards that goal. However, there is little experience with plasmas under conditions close to those expected in future fusion reactors, because it requires the fusion of Deuterium (D) and Tritium (T), while most of the experiments are currently performed in pure D. After more than 20 years, the Joint European Torus (JET) has carried out new D-T experiments with the aim of exploring the unique characteristics of burning D-T plasmas, such as the presence of highly energetic ions. A new stable, high confinement and impurity-free D-T regime, with strong reduction of energy losses with respect to D, has been found. Multiscale physics mechanisms critically determine the thermal confinement and the fusion power yield. These crucial achievements importantly contribute to the establishment of fusion energy generation as an alternative to fossil fuels.
The fusion reaction between deuterium and tritium, D(T,n)4 He is the main source of energy in future thermonuclear reactors. Alpha-particles (4 He-ions) born with an average energy of 3.5 MeV transferring energy to the thermal plasma during their slowing down, should provide the self-sustained D-T plasma burn. The adequate confinement of alpha-particles is essential to provide efficient heating of the bulk plasma and steady burning of a reactor plasma. That is why the fusion-born alpha-particle studies have been a priority task in the second D-T experiments (DTE2) on the Joint European Torus (JET) to understand the main mechanisms of their slowing down, redistribution and losses and to develop optimal plasma scenarios. JET with Be-wall and W-divertor, enhanced auxiliary heating systems and improved energetic-particle diagnostic capabilities, producing significant population of alpha-particles, provided the possibility for comprehensive studying of the alpha-particle behaviour. Selected results of the confined and lost alpha-particle measurements, evidence of alpha-particle self-heating and assessments of the fusion performance are presented in this paper giving an opportunity for further modelling and extrapolation to the International Thermonuclear Experimental Reactor and burning plasma reactors.
This work presents a study of the interaction between radio frequency (RF) waves used for ion cyclotron resonance heating and the fast deuterium (D) and tritium (T) neutral Beam injected (NBI) ions in DT plasma. The focus is on the effects of this interaction, also referred to as synergistic effects, on the fusion performance in the recent JET DTE2 campaign. Experimental data from dedicated pulses at 3.43 T/2.3 MA heated at (i) 51.4 MHz, giving the central minority H and n = 2 D, and at (ii) 32.2 MHz for the central minority 3 He and n = 2 T. Resonances are analysed and conclusions are drawn and supported by modelling of the synergistic effects. Modelling with transport code TRANSP runs with and without the RF kick operator predict a moderate increase, of about 10%, in DT rates for the case of the RF wave—fast D NBI ion interactions at the n = 2 harmonic of ion cyclotron resonance, and a negligible impact due to synergistic interaction between fast T NBI ions and RF waves. JETTO modelling gives a 29% enhancement in fusion rates due to the interction between RF waves and fast D NBI ions, and an 18% enhancement in fast T NBI ions. Analysis of experimental neutron rates compared to TRANSP predictions without synergistic effects and magnetic proton recoil neutron spectrometer indicate an enhancement of approximately 25%–28% in fusion rates due to RF interaction with fast D ions, and an enhancement of approximately 5%–8% when RF waves and fast T NBI ions are interacting. The contributions of various heating and fast ion sources are assessed and discussed.
Within the 9th European Framework programme, since 2021 EUROfusion is operating five tokamaks under the auspices of a single Task Force called ‘Tokamak Exploitation’. The goal is to benefit from the complementary capabilities of each machine in a coordinated way and help in developing a scientific output scalable to future largre machines. The programme of this Task Force ensures that ASDEX Upgrade, MAST-U, TCV, WEST and JET (since 2022) work together to achieve the objectives of Missions 1 and 2 of the EUROfusion Roadmap: i) demonstrate plasma scenarios that increase the success margin of ITER and satisfy the requirements of DEMO and, ii) demonstrate an integrated approach that can handle the large power leaving ITER and DEMO plasmas. The Tokamak Exploitation task force has therefore organized experiments on these two missions with the goal to strengthen the physics and operational basis for the ITER baseline scenario and for exploiting the recent plasma exhaust enhancements in all four devices (PEX: Plasma EXhaust) for exploring the solution for handling heat and particle exhaust in ITER and develop the conceptual solutions for DEMO. The ITER Baseline scenario has been developed in a similar way in ASDEX Upgrade, TCV and JET. Key risks for ITER such as disruptions and run-aways have been also investigated in TCV, ASDEX Upgrade and JET. Experiments have explored successfully different divertor configurations (standard, super-X, snowflakes) in MAST-U and TCV and studied tungsten melting in WEST and ASDEX Upgrade. The input from the smaller devices to JET has also been proven successful to set-up novel control schemes on disruption avoidance and detachment.
Analysis of comparable discharges fuelled by either deuterium or protium reveals a clear relationship between the isotope mass and the H-mode density limit. Notably, the density limit is significantly lower in protium, showing a reduction of up to 35 % compared to identical deuterium plasma conditions. Within mixed H-mode density limit (HDL) plasmas, the maximum achievable density, or H-mode density limit, decreases with increasing protium concentration, denoted as cH. For instance, the highest corresponding maximum Greenwald fraction (fGW) of about 1.02 was observed in the pulse with the lowest cH value of 4.4 %. This fGW decreases to 0.96 at cH = 48 %. The average atomic mass, A, of the plasma species decreases in these pulses from the value of 1.96 (cH = 4.4%) down to 1.52 (cH = 48 %). Interestingly, the maximum achievable density appears to be largely unaffected by the applied power value, regardless of whether deuterium or protium is used, as well as under mixed H/D fuelling conditions. Additionally, the measured Greenwald fractions are agreed with a heuristic model based on the SOL pressure threshold of an MHD instability, as proposed by Goldston. This comparison, especially concerning the model's dependence on isotopic mass, shows full consistency between the measured and predicted Greenwald fractions.
In the JET DTE2 deuterium-tritium campaign, neutron diagnostics were employed to measure 14 MeV neutrons originating from D(T,n)4He reactions. In discharge 99965, a diamond matrix detector (KM14) and a magnetic proton recoil (MPRu) detector with a vertical and an oblique line-of-sight were used, respectively. At the timepoints of interest, a significant decrease in the expected diagnostic signals can be observed as electromagnetic wave heating in the ion cyclotron range of frequencies (ICRF) is switched off. Utilizing only TRANSP simulation data, the fast-ion distribution is found to have been likely composed mostly of trapped orbits. In contrast, analysis performed using orbit weight functions revealed that the majority of neutrons in the KM14 Ed=9.3 MeV and MPRu Xcm=33 cm measurement bins are to have originated from fast deuterium ions on co-passing orbits. This work explains the perhaps surprising results and shows that the relative signal decrease as ICRF heating is switched off is largest for counter-passing orbits. Finally, for the magnetic equilibria of interest, it is shown how stagnation orbits, corresponding to similar to 1 % of the fast-ion distribution, were completely unobservable by the KM14 diagnostic.
In the JET DTE2 campaign a new method was successfully tested to detect the heating of bulk electrons by α-particles, using the dynamic response of the electron temperature T e to the modulation of ion cyclotron resonance heating (ICRH). A fundamental deuterium (D) ICRH scheme was applied to a tritium-rich hybrid plasma with D-neutral beam injection (NBI). The modulation of the ion temperature T i and of the ICRH accelerated deuterons leads to modulated α -heating with a large delay with respect to other modulated electron heating terms. A significant phase delay of ∼40° is measured between central T e and T i , which can only be explained by α -particle heating. Integrated modelling using different models for ICRH absorption and ICRH/NBI interaction reproduces the effect qualitatively. Best agreement with experiment is obtained with the European Transport Solver/Heating and Current Drive workflow.
Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power (P rad) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling ('guess and check') approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced chi 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98+0.03/ -0.29 and 1.01+0.02/-0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.
JET, the world’s largest operating tokamak with unique Be/W wall and tritium handling capability, completed a Deuterium-Tritium (D-T) campaign in 2021 (Maggi et al 29th Fusion Energy Conf. ) following a decade of preparatory experiments, dedicated enhancements, technical rehearsals and training (Horton et al 2016 Fusion Eng. Des. 109–111 925). Operation with tritium raises significant technical, safety and scientific challenges not encountered in standard protium or deuterium operation. This contribution describes the tritium operational requirements, pulses and technical preparations, new operating procedures, lessons learned and details on the achieved operational availability and performance. The preparation and execution of the recent JET tritium experiments benefitted from the previous experience in 1991 (Preliminary Tritium Experiment), 1997 (DTE1 campaign) and 2003 (Trace Tritium Campaigns) and consisted of the following five phases: technical rehearsals and scenario preparation, tritium commissioning, 100% tritium campaign, D-T campaign (DTE2), tritium clean-up. Following the clean-up JET resumed normal operation and is currently undertaking a further D-T campaign (DTE3).
To produce fusion reactions efficiently, thermonuclear plasmas have to reach extremely high temperatures, which is incompatible with their coming into contact with material surfaces. Confinement of plasmas using magnetic fields has progressed significantly in the last years, particularly in the tokamak configuration. Unfortunately, all tokamak devices, and particularly metallic ones, are plagued by catastrophic events called disruptions. Many disruptions are preceded by anomalies in the radiation patterns, particularly in ITER-relevant scenarios. These specific forms of radiation emission either directly cause or reveal the approaching collapse of the configuration. Detecting the localization of these radiation anomalies in real time requires an innovative and specific elaboration of bolometric measurements, confirmed by visible cameras and the inversion of sophisticated tomographic algorithms. The information derived from these measurements can be interpreted in terms of local power balances, which suggest a new quantity, the radiated power divided by the plasma internal energy, to determine the criticality of the plasma state. Combined with robust indicators of the temperature profile shape, the identified anomalous radiation patterns allow determination of the sequence of macroscopic events leading to disruptions. A systematic analysis of JET campaigns at high power in deuterium, full tritium, and DT, for a total of almost 2000 discharges, proves the effectiveness of the approach. The warning times are such that, depending on the radiation anomaly and the available actuators, the control system of future devices is expected to provide enough notice to enable deployment of effective prevention and avoidance strategies.
JET-ILW type I ELMy H-modes at 2.5 MA/2.8 T with constant NBI heating (23 MW) and gas fuelling rate were performed, utilising edge localised mode (ELM) pacing by vertical kicks and plasma shaping (triangularity, δ ) as tools to disentangle the effects of ELMs, inter-ELM transport and edge stability on the pedestal particle balance. In agreement with previous studies, the pedestal confinement improves with increasing δ , mostly due to a significant increase in pedestal density while the ELM frequency ( f E L M ) is decreased. Improved pedestal confinement with increasing δ was observed even when the pedestal MHD stability was degraded artificially by vertical kicks, implying that increased triangularity may favourably affect the inter-ELM pedestal recovery. The workflow developed to quantify the pedestal particle balance uses high time-resolution profile reflectometry to characterise the inter-ELM evolution of the plasma particle content ( d N / d t ), the NEO drift-kinetic solver to evaluate the neoclassical fluxes and interpretative EDGE2D-EIRENE simulations to estimate the edge particle source. The edge particle source is then constrained by deuterium Balmer- α line intensity measurements in the main chamber, which are, however, strongly affected by reflections from the metal walls. The reflections are accounted for by the CHERAB code taking the divertor emission (the brightest light source in the torus) distribution from imaging spectroscopy measurements as input. Our analysis shows that in the second half of the ELM cycle, the volume-integrated particle source is larger than d N / d t , indicating that transport plays a key role in the inter-ELM pedestal recovery.
This work applies the coupled JINTRAC and QuaLiKiz-neural-network (QLKNN) model on the ohmic current ramp-up phase of a JET D discharge. The chosen scenario exhibits a hollow T e profile attributed to core impurity accumulation, which is observed to worsen with the increasing fuel ion mass from D to T. A dynamic D simulation was validated, evolving j , n e , T e , T i , n Be , n Ni , and n W for 7.25 s along with self-consistent equilibrium calculations, and was consequently extended to simulate a pure T plasma in a predict-first exercise. The light impurity (Be) accounted for Z eff while the heavy impurities (Ni, W) accounted for P rad . This study reveals the role of transport on the T e hollowing, which originates from the isotope effect on the electron-ion energy exchange affecting T i . This exercise successfully affirmed isotopic trends from previous H experiments and provided engineering targets used to recreate the D q -profile in T experiments, demonstrating the potential of neural network surrogates for fast routine analysis and discharge design. However, discrepancies were found between the impurity transport behaviour of QuaLiKiz and QLKNN, which lead to notable T e hollowing differences. Further investigation into the turbulent component of heavy impurity transport is recommended.
The recent deuterium–tritium campaign in JET-ILW (DTE2) has provided a unique opportunity to study the isotope dependence of the L-H power threshold in an ITER-like wall environment (Be wall and W divertor). Here we present results from dedicated L-H transition experiments at JET-ILW, documenting the power threshold in tritium and deuterium–tritium plasmas, comparing them with the matching deuterium and hydrogen datasets. From earlier experiments in JET-ILW it is known that as plasma isotopic composition changes from deuterium, through varying deuterium/hydrogen concentrations, to pure hydrogen, the value of the line averaged density at which the threshold is minimum, nˉe,min , increases, leading us to expect that nˉe,min (T) < nˉe,min (DT) < nˉe,min (D) < nˉe,min (H). The new power threshold data confirms these expectations in most cases, with the corresponding ordering of the minimum power thresholds. We present a comparison of this data to power threshold scalings, used for extrapolation to future devices such as ITER and DEMO.