The major increase in discharge duration and plasma energy in a next step DT fusion reactor will give rise to important plasma-material effects that will critically influence its operation, safety and performance. Erosion will increase to a scale of several centimetres from being barely measurable at a micron scale in today's tokamaks. Tritium co-deposited with carbon will strongly affect the operation of machines with carbon plasma facing components. Controlling plasma-wall interactions is critical to achieving high performance in present day tokamaks, and this is likely to continue to be the case in the approach to practical fusion reactors. Recognition of the important consequences of these phenomena stimulated an internationally co-ordinated effort in the field of plasma-surface interactions supporting the Engineering Design Activities of the International Thermonuclear Experimental Reactor project (ITER), and significant progress has been made in better understanding these issues. The paper reviews the underlying physical processes and the existing experimental database of plasma-material interactions both in tokamaks and laboratory simulation facilities for conditions of direct relevance to next step fusion reactors. Two main topical groups of interaction are considered: (i) erosion/redeposition from plasma sputtering and disruptions, including dust and flake generation and (ii) tritium retention and removal. The use of modelling tools to interpret the experimental results and make projections for conditions expected in future devices is explained. Outstanding technical issues and specific recommendations on potential R&D avenues for their resolution are presented.

/pdf/plasma-material-interactions-in-current-tokamaks-and-their-1a6iwv89x6.pdf

Plasma{material interactions in current tokamaks and their implications for next step fusion reactors

Progress in the area of MHD stability and disruptions, since the publication of the 1999 ITER Physics Basis document (1999 Nucl. Fusion 39 2137-2664), is reviewed. Recent theoretical and experimental research has made important advances in both understanding and control of MHD stability in tokamak plasmas. Sawteeth are anticipated in the ITER baseline ELMy H-mode scenario, but the tools exist to avoid or control them through localized current drive or fast ion generation. Active control of other MHD instabilities will most likely be also required in ITER. Extrapolation from existing experiments indicates that stabilization of neoclassical tearing modes by highly localized feedback-controlled current drive should be possible in ITER. Resistive wall modes are a key issue for advanced scenarios, but again, existing experiments indicate that these modes can be stabilized by a combination of plasma rotation and direct feedback control with non-axisymmetric coils. Reduction of error fields is a requirement for avoiding non-rotating magnetic island formation and for maintaining plasma rotation to help stabilize resistive wall modes. Recent experiments have shown the feasibility of reducing error fields to an acceptable level by means of non-axisymmetric coils, possibly controlled by feedback. The MHD stability limits associated with advanced scenarios are becoming well understood theoretically, and can be extended by tailoring of the pressure and current density profiles as well as by other techniques mentioned here. There have been significant advances also in the control of disruptions, most notably by injection of massive quantities of gas, leading to reduced halo current fractions and a larger fraction of the total thermal and magnetic energy dissipated by radiation. These advances in disruption control are supported by the development of means to predict impending disruption, most notably using neural networks. In addition to these advances in means to control or ameliorate the consequences of MHD instabilities, there has been significant progress in improving physics understanding and modelling. This progress has been in areas including the mechanisms governing NTM growth and seeding, in understanding the damping controlling RWM stability and in modelling RWM feedback schemes. For disruptions there has been continued progress on the instability mechanisms that underlie various classes of disruption, on the detailed modelling of halo currents and forces and in refining predictions of quench rates and disruption power loads. Overall the studies reviewed in this chapter demonstrate that MHD instabilities can be controlled, avoided or ameliorated to the extent that they should not compromise ITER operation, though they will necessarily impose a range of constraints.

https://iopscience.iop.org/article/10.1088/0029-5515/47/6/S03/pdf

Chapter 3: MHD stability, operational limits and disruptions

The ITER Physics Basis presents and evaluates the physics rules and methodologies for plasma performance projections, which provide the basis for the design of a tokamak burning plasma device whose goal is to demonstrate the scientific and technological feasibility of fusion energy for peaceful purposes. This Chapter summarizes the physics basis for burning plasma projections, which is developed in detail by the ITER Physics Expert Groups in subsequent chapters. To set context, the design guidelines and requirements established in the report of ITER Special Working Group 1 are presented, as are the specifics of the tokamak design developed in the Final Design Report of the ITER Engineering Design Activities, which exemplifies burning tokamak plasma experiments. The behaviour of a tokamak plasma is determined by the interaction of many diverse physics processes, all of which bear on projections for both a burning plasma experiment and an eventual tokamak reactor. Key processes summarized here are energy and particle confinement and the H-mode power threshold; MHD stability, including pressure and density limits, neoclassical islands, error fields, disruptions, sawteeth, and ELMs; power and particle exhaust, involving divertor power dispersal, helium exhaust, fuelling and density control, H-mode edge transition region, erosion of plasma facing components, tritium retention; energetic particle physics; auxiliary power physics; and the physics of plasma diagnostics. Summaries of projection methodologies, together with estimates of their attendant uncertainties, are presented in each of these areas. Since each physics element has its own scaling properties, an integrated experimental demonstration of the balance between the combined processes which obtains in a reactor plasma is inaccessible to contemporary experimental facilities: it requires a reactor scale device. It is argued, moreover, that a burning plasma experiment can be sufficiently flexible to permit operation in a steady state mode, with non-inductive plasma current drive, as well as in a pulsed mode where current is inductively driven. Overall, the ITER Physics Basis can support a range of candidate designs for a tokamak burning plasma facility. For each design, there will remain a significant uncertainty in the projected performance, but the projection methodologies outlined here do suffice to specify the major parameters of such a facility and form the basis for assuring that its phased operation will return sufficient information to design a prototype commercial fusion power reactor, thus fulfilling the goal of the ITER project.

https://iopscience.iop.org/article/10.1088/0029-5515/39/12/301/pdf

Chapter 1: Overview and summary

Nonlinear gyrokinetic equations play a fundamental role in our understanding of the long-time behavior of strongly magnetized plasmas. The foundations of modern nonlinear gyrokinetic the- ory are based on three important pillars: (1) a gyrokinetic Vlasov equation written in terms of a gyrocenter Hamiltonian with quadratic low-frequency ponderomotive-like terms; (2) a set of gyrokinetic Maxwell (Poisson-Ampere) equations written in terms of the gyrocenter Vlasov dis- tribution that contain low-frequency polarization (Poisson) and magnetization (Ampere) terms derived from the quadratic nonlinearities in the gyrocenter Hamiltonian; and (3) an exact energy conservationlaw for the gyrokineticVlasov-Maxwell equations that includes all the relevant linear and nonlinear coupling terms. The foundations of nonlinear gyrokinetic theory are reviewed with an emphasis on the rigorous applications of Lagrangian and Hamiltonian Lie-transform perturba- tion methods used in the variationalderivationof nonlineargyrokineticVlasov-Maxwell equations. The physical motivations and applications of the nonlinear gyrokinetic equations, which describe the turbulent evolution of low-frequency electromagnetic fluctuations in a nonuniform magnetized plasmas with arbitrary magnetic geometry, are also discussed.

/pdf/foundations-of-nonlinear-gyrokinetic-theory-3l78lnlpai.pdf

Foundations of nonlinear gyrokinetic theory

Progress, since the ITER Physics Basis publication (ITER Physics Basis Editors et al 1999 Nucl. Fusion 39 2137–2664), in understanding the processes that will determine the properties of the plasma edge and its interaction with material elements in ITER is described. Experimental areas where significant progress has taken place are energy transport in the scrape-off layer (SOL) in particular of the anomalous transport scaling, particle transport in the SOL that plays a major role in the interaction of diverted plasmas with the main-chamber material elements, edge localized mode (ELM) energy deposition on material elements and the transport mechanism for the ELM energy from the main plasma to the plasma facing components, the physics of plasma detachment and neutral dynamics including the edge density profile structure and the control of plasma particle content and He removal, the erosion of low- and high-Z materials in fusion devices, their transport to the core plasma and their migration at the plasma edge including the formation of mixed materials, the processes determining the size and location of the retention of tritium in fusion devices and methods to remove it and the processes determining the efficiency of the various fuelling methods as well as their development towards the ITER requirements. This experimental progress has been accompanied by the development of modelling tools for the physical processes at the edge plasma and plasma–materials interaction and the further validation of these models by comparing their predictions with the new experimental results. Progress in the modelling development and validation has been mostly concentrated in the following areas: refinement in the predictions for ITER with plasma edge modelling codes by inclusion of detailed geometrical features of the divertor and the introduction of physical effects, which can play a major role in determining the divertor parameters at the divertor for ITER conditions such as hydrogen radiation transport and neutral–neutral collisions, modelling of the ion orbits at the plasma edge, which can play a role in determining power deposition at the divertor target, models for plasma–materials and plasma dynamics interaction during ELMs and disruptions, models for the transport of impurities at the plasma edge to describe the core contamination by impurities and the migration of eroded materials at the edge plasma and its associated tritium retention and models for the turbulent processes that determine the anomalous transport of energy and particles across the SOL. The implications for the expected performance of the reference regimes in ITER, the operation of the ITER device and the lifetime of the plasma facing materials are discussed.

/pdf/chapter-4-power-and-particle-control-3l68lj0fxw.pdf

Chapter 4: Power and particle control

The goal of this work is to generate large statistically representative data sets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student $t$ process regression. We apply Student $t$ process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via colouring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics and classic machine learning clustering algorithms.

Data augmentation for disruption prediction via robust surrogate models

Recent Alcator C-Mod experimental campaigns have focused upon the study of the Advanced Tokamak regimes, which includes characterization of the RF heating, the formation and dynamics of internal barriers, H-mode edge pedestal, and divertor and scrape-off physics. The ICRF system has been recently upgraded with the improved performance of the 4-strap antenna. Total ICRF power in excess of 5 MW has been launched successfully into the plasma during this campaign. Due to the compact nature of C-Mod, the power feeds for the antenna are vacuum strip lines. Their orientation, to the tokamak B-field, is governed by maintaining E<15 kV/cm in locations where the RF E-field is parallel to tokamak B-field. Other modifications included improved protection tile grounding and installation of protective shields for Faraday screen ceramic isolators. The antennas also make use of BN protection tiles to eliminate high Z impurities from the antennas. The present empirical power limit results from arcing in a region of the antenna strap where E/spl sim/15 kV/cm and parallel to B and injections from the metallic fasteners used to attach the BN tiles to the antenna.

https://ieeexplore.ieee.org/iel5/7981/22072/01027697.pdf

Results and status of the Alcator C-Mod tokamak

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c, that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c, due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

https://iopscience.iop.org/article/10.1088/1741-4326/acdd11/pdf

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

Asymmetric halo currents (HCs) can exert large net forces on the vacuum vessel and other components during disruptions on tokamaks. The displacements caused by these forces can then be amplified if these asymmetric forces rotate at frequencies resonant with the vessel. This paper reports on the investigation of a recently proposed scaling law for the disruption HC rotation frequency [Saperstein et al., “Halo current rotation scaling in post-disruption plasmas,” Nucl. Fusion 62, 026044 (2022)] that combines measurements on Alcator C-Mod with those on HBT-EP. We find that a new non-circular version of the scaling law [[Formula: see text]] takes into consideration the dependence of frot on the poloidal structure of the MHD instability ( m) driving the asymmetry and describes the disruption-averaged rotation frequency on C-Mod. Disruption rotation is also found to be insensitive to the vertical position and impurity content of the plasma at the onset of the disruption. However, a stagnation in the time evolution of frot is occasionally observed. Observations are consistent with the dominance of poloidal rotation during the disruption, which is motivated by the poloidal drift nature of the scaling law.

Disruption halo current rotation scaling on Alcator C-Mod and HBT-EP

Halo currents generated during disruptions on Alcator C-Mod have been measured with Langmuir "rail" probes These rail probes are embedded in a lower outboard divertor module in a closely-spaced vertical (poloidal) array The dense array provides detailed resolution of the spatial dependence (~1 cm spacing) of the halo current distribution in the plasma scrape-off region with high time resolution (400 kHz digitization rate) As the plasma limits on the outboard divertor plate, the contact point is clearly discernible in the halo current data (as an inversion of current) and moves vertically down the divertor plate on many disruptions These data are consistent with filament reconstructions of the plasma boundary, from which the edge safety factor of the disrupting plasma can be calculated Additionally, the halo current "footprint" on the divertor plate is obtained and related to the halo flux width The voltage driving halo current and the effective resistance of the plasma region through which the halo current flows to reach the probes are also investigated Estimations of the sheath resistance and halo region resistivity and temperature are given This information could prove useful for modeling halo current dynamics

Robert Granetz

Papers

Data augmentation for disruption prediction via robust surrogate models

Results and status of the Alcator C-Mod tokamak

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

Disruption halo current rotation scaling on Alcator C-Mod and HBT-EP

High-resolution disruption halo current measurements using Langmuir probes in Alcator C-Mod