Real-time magnetic control has been developed to deliver precise control of multiple plasma shape parameters for advanced divertor configurations, including double-null, Super-X, X-point target and X-divertor for the first time on the MAST Upgrade (MAST-U) spherical tokamak. Successful real-time magnetic equilibrium control of different plasma shape variables has been accomplished in the 2022–2023 MAST-U experimental campaign for the advanced MAST-U divertor configurations. Application of the MAST-U boundary reconstruction algorithm, LEMUR, is described and compared with off-line equilibrium reconstruction and diagnostic measurements. The process of design and verification of the axisymmetric magnetic control schemes using a suite of control analysis tools (known collectively as TokSys) is also described.
The DIII-D tokamak has elucidated crucial physics and developed projectable solutions for ITER and fusion power plants in the key areas of core performance, boundary heat and particle transport, and integrated scenario operation, with closing the core-edge integration knowledge gap being the overarching mission. New experimental validation of high-fidelity, multi-channel, non-linear gyrokinetic turbulent transport models for ITER provides strong confidence it will achieve Q 10 operation. Experiments identify options for easing H-mode access in hydrogen, and give new insight into the isotopic dependence of transport and confinement. Analysis of 2,1 islands in unoptimized low-torque IBS demonstration discharges suggests their onset time occurs randomly in the constant beta phase, most often triggered by non-linear 3-wave coupling, thus identifying an NTM seeding mechanism to avoid. Pure deuterium SPI for disruption mitigation is shown to provide favorable slow cooling, but poor core assimilation, suggesting paths for improved SPI on ITER. At the boundary, measured neutral density and ionization source fluxes are strongly poloidally asymmetric, implying a 2D treatment is needed to model pedestal fuelling. Detailed measurements of pedestal and SOL quantities and impurity charge state radiation in detached divertors has validated edge fluid modelling and new self-consistent 'pedestal-to-divertor' integrated modeling that can be used to optimize reactors. New feedback adaptive ELM control minimizes confinement reduction, and RMP ELM suppression with sustained high core performance was obtained for the first time with the outer strike point in a W-coated, compact and unpumped small-angle slot divertor. Advances have been made in integrated operational scenarios for ITER and power plants. Wide pedestal intrinsically ELM-free QH-modes are produced with more reactor-relevant conditions, Low torque IBS with W-equivalent radiators can exhibit predator-prey oscillations in T-e and radiation which need control. High-beta(P) scenarios with q(min) > 2, q(95)-7.9, beta(N) > 4, beta(T)-3.3% and H-98y2 > 1.5 are sustained with high density ((n) over bar = 7E19 m(-3), f(G)-1) for 6 tau(E), improving confidence in steady-state tokamak reactors. Diverted NT plasmas achieve high core performance with a non-ELMing edge, offering a possible highly attractive core-edge integration solution for reactors.
The paper details the process of developing the ITER Plasma Control System (PCS), that is, how to design and deploy it systematically, in the most efficient and effective manner. The integrated nature of the ITER PCS, with its multitude of coupled control functions, and its long-term development, calls for a different approach than the design and short-term deployment of individual controllers. It requires, in the first place, a flexible implementation strategy and system architecture that allows system re-configuration and optimization throughout its development. Secondly, a model-based system engineering approach is carried out, for the complete PCS development, i.e. both its design and deployment. It requires clear definitions for both the PCS role and its functionality, as well as definitions of the design and deployment process itself. The design and deployment process is shown to allow tracing the relationships of the many individual design and deployment aspects, such as system requirements, assumed operation use-cases and response models, and eventually verification and functional validation of the system design. The functional validation will make use of a dedicated PCS simulation platform that includes the description of the control function design as well as plant, actuator and sensor models that enable the simulation of these functions. By establishing a clear understanding of the interconnected steps involved in designing, implementing, commissioning, and operating the system, a more systematic approach is achieved. This ensures the completion of a comprehensive design that can be deployed efficiently, hence preventing the loss of precious operational time needed to debug and retune control functions and more importantly avoiding tokamak discharge disruptions.
Magnetic controller performance requirements and design solutions for the Spherical Tokamak for Energy Production (STEP) are driven by the need to produce equilibrium and scenario trajectories, maintaining steady plasma vertical stabilization and shape control for a period of thousands of seconds or longer while avoiding contact with the plasma-facing components (PFCs). Axisymmetric magnetic control schemes including vertical stabilization, plasma current control, plasma shape control and poloidal field coil current control for the STEP Prototype Reactor (SPR) concept are being developed using a suite of control analysis tools (known collectively as TokSys) supporting the integrated plasma control design process. The vertical growth rate based on a linear rigid plasma response model in TokSys is used for assessing the controllability of the vertical instability in SPR. TokSys closed-loop simulations with the axisymmetric non-linear, free-boundary evolution code, GSevolve, are performed for the assessment, identification, and verification of algorithm implementation and controller performance for the various axisymmetric control systems. Dynamic control performance and controllability are shown to be consistent with noise-affected scenario requirements under reasonable power supply and sensor performance assumptions.
HL-2M is a medium-sized tokamak constructed by the Southwestern Institute of Physics (SWIP) in China. A new plasma control system (PCS) based on the software framework of DIII-D PCS for HL-2M has been proposed to replace the original loose system. The new hardware system proposed in this paper consists of three computer nodes and peripheral equipment. Data acquisition is completed by a standing server and a reflective memory network that passes back data. The software architecture of DIII-D PCS is generally retained, but some interfaces are added to adapt to the existing system. In addition, the Linux operating system has also been upgraded and optimized so that real-time tasks can achieve better performance. The initial system involves some basic control algorithms satisfying early experiments. The new PCS has been preliminarily verified through simulation and can output control commands as expected in the integrated environment. This paper focuses on two aspects of DIII-D PCS on a new device, i.e., performance and interface adaptation. More efforts will be made to add more algorithms and use the system in real experiments.
An algorithm has been designed to simultaneously control the central safety factor (q0) and normalized beta (βN) while ensuring near-zero torque from neutral beam injection in DIII-D. Feedback control of q0 and βN in tokamaks can be beneficial due to the close relationship that these variables have with plasma performance and magneto-hydrodynamic stability. In addition, low neutral-beam-torque conditions are of special interest in present devices because future burning-plasma tokamaks such as ITER will most likely operate at very low plasma rotation. The control synthesis of the algorithm presented in this work is based on a linearized, one-dimensional (1D) model of the current-profile dynamics coupled with a zero-dimensional (0D) plasma-energy balance. The actuators considered are neutral beam injection and electron-cyclotron heating and current drive, and discrete logic determines the neutral-beam injection powers that deliver near-zero torque. The algorithm has been tested in nonlinear, 1D simulations using COTSIM (Control-Oriented Transport SIMulator) and in DIII-D experiments, demonstrating satisfactory performance.
For the problem of excessively high divertor heat flux, active impurity seeding is an effective method to radiate the plasma energy reaching the divertor and thus achieve the divertor detachment. Neon is a very effective radiation impurity on many current tokamaks, which is also a candidate species to be applied on ITER. In the EAST 2019 experimental campaign, a series of experiments were performed by seeding a mixture of neon and deuterium (Ne-D 2 ) for detachment and core-edge-divertor integration in H-mode plasmas. The divertor partial detachment with high-confinement core plasma has been achieved by using Ne-D 2 seeding in EAST with ITER-like tungsten divertor. Both the plasma stored energy and H 98, y 2 > 1.1 are maintained, with the divertor electron temperature, heat flux and the surface temperature near the strike point being all significantly reduced. The differences between Ne-D 2 seeding at the scrape-off layer (SOL) upstream and downstream have been experimentally investigated in detail. It is found that impurity seeding at SOL downstream is more beneficial to reducing the divertor electron temperature and peak heat flux. By comparison with experiments using divertor D 2 fueling, it is further demonstrated that gas seeding in the SOL downstream will enrich more particles near the strike point, while the seeding in the SOL upstream will influence the entire outer target more evenly. Furthermore, in most of the experiments, gas seeding does not cause obvious toroidal asymmetry in the divertor plasma. However, when D 2 is injected in an amount similar to that used to build the plasma, it causes the particle flux near the gas-puff to increase locally, i.e., much more than that at the toroidal location far from the gas-puff location. It is a competition between particle source and transport. When the particle source is stronger, it will naturally increase the local particles. In addition, dedicated experiments with different poloidal distances between impurity seeding and strike point on the radiation ability were carried out. Both experimental results and SOLPS simulation show that the seeding close to the strike point is more conducive to neon ionization and energy radiation.
The compatibility of efficient divertor detachment with high-performance core plasma is vital to the development of magnetically controlled fusion energy. The joint research on the EAST and DIII-D tokamaks demonstrates successful integration of divertor detachment with excellent core plasma confinement quality, a milestone towards solving the critical plasma–wall-interaction (PWI) issue and core-edge integration for ITER and future reactors. In EAST, actively controlled partial detachment with T et,div ∼ 5 eV around the strike point and H 98 > 1 in different H-mode scenarios including the high β P H-mode scenario have been achieved with ITER-like tungsten divertor, by optimizing the detachment access condition and performing detailed experiments for core-edge integration. For active long-pulse detachment feedback control, a 30 s H-mode operation with detachment-control duration being 25 s has been successfully achieved in EAST. DIII-D has achieved actively controlled fully detached divertor with low plasma electron temperature ( T et,div ⩽ 5 eV across the entire divertor target) and low particle flux (degree of detachment, DoD > 3), simultaneously with very high core performance ( β N ∼3, β P > 2 and H 98 ∼ 1.5) in the high β P scenario being developed for ITER and future reactors. The high- β P high confinement scenario is characterized by an internal transport barrier (ITB) at large radius and a weak edge transport barrier (ETB, or pedestal), which are synergistically self-organized. Both the high- β P scenario and impurity seeding facilitate divertor detachment. The detachment access leads to the reduction of ETB, which facilitates the development of an even stronger ITB at large radius in the high β P scenario. Thus, this strong large radius ITB enables the core confinement improvement during detachment. These significant joint DIII-D and EAST advances on the compatibility of high confinement core and detached divertor show a great potential for achieving a high-performance core plasma suitable for long-pulse operation of fusion reactors with controllable steady-state PWIs.
DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter- I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high- Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ∼ 8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co- I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.
Magnetic diagnostics in tokamaks are key to plasma equilibrium control (plasma current, plasma shape, and position) and amelioration of plasma instabilities. Thus, real-time identification of the anomalous sensor is mandatory. A new real-time system based on autoencoder (AE) neural networks (NNs) for anomaly detection in magnetics signals, including both flux loops and magnetic probes, has been successfully implemented on the plasma control system (PCS) of the DIII-D tokamak. The AE NN is trained on over 4000 plasma discharges, with an optimized latent space representation of the input signals while minimizing the reconstruction loss. An algorithm determining anomalous sensors based on excessive deviations from accurate reconstruction from the trained NNs is constructed in a MATLAB/Simulink environment and is deployed on the real-time PCS using the embedded MATLAB coder (EMC) environment. The reconstruction performance of the algorithm is quantified in hardware-in-the-loop (HIL) simulation by injecting artificially common fault types, including stuck-at-zero fault, saturated fault, noise fault, and drift fault. The AE NN algorithm has been successfully commissioned on the DIII-D PCS.
Reliable disruption avoidance techniques are critical for the development of safe and economically viable fusion reactors. This incipient need has driven the fusion community to spend substantial effort to develop machine learning (ML) models that aim to assist in this area. One area of concern is a sudden loss in global confinement associated with an unintended back-transition from H-mode confinement to L-mode confinement ("H–L back-transition"), which can sometimes lead to a disruption and increased risk of damage to a reactor-grade fusion device. A recent experiment on DIII-D demonstrated a real-time control method able to steer away from such unwanted back transitions. The present confinement mode of the plasma was inferred in real-time using a temporal convolutional neural network that outputs a continuous scalar estimate suitable for input to a closed-loop controller. This scalar estimate was then regulated to a value consistent with H-mode by adjusting power associated with the ratio of the plasma pressure to magnetic pressure ( $\beta _{N}$ ). The model was trained on hundreds of thousands of discrete time steps of expert-labeled data of periods of H-mode and L-mode confinement, taken from hundreds of DIII-D tokamak discharges. This work presents these experimental results, including the method of collecting expert-labeled H- and L-mode confinement data, and the ML model performance and efficacy on the mode prediction task. In addition, the first application of this tool for H–L back-transition prevention is demonstrated in this experiment. Details are presented on the training of the neural network model, including concerns such as hyperparameter tuning and the network architecture. In addition, the methodology for embedding the neural network into the control system for real-time inference and closed-loop control is discussed.
Next generation High Performance (HP) tokamaks risk damage from unmitigated disruptions at high current and power. Achieving reliable disruption prediction for a device’s HP operation based on its Low Performance (LP) data is a key to its success. In this letter, through explorative data analysis and dedicated numerical experiments on multiple existing tokamaks, we demonstrate how the operational regimes of tokamaks can affect the power of a trained disruption predictor. First, our results suggest data-driven disruption predictors trained on abundant LP discharges work poorly on the HP regime of the same tokamak, which is a consequence of the distinct distributions of the tightly correlated signals related to disruptions in these two regimes. Second, we find that matching operational parameters among tokamaks strongly improves cross-machine accuracy which implies our model learns from the underlying scalings of dimensionless physics parameters like q 95, β p and confirms the importance of these parameters in disruption physics and cross machine domain matching from the data-driven perspective. Finally, our results show in the absence of HP data from the target devices, the best predictivity of the HP regime for the target machine can be achieved by combining LP data from the target with HP data from other machines. These results provide a possible disruption predictor development strategy for next generation tokamaks, such as ITER and SPARC, and highlight the importance of developing baseline scenario discharges of future tokamaks on existing machines to collect more relevant disruptive data.
Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real-time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad s−1 with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I N < 0.3 before terminating, scaling to the 3 MA ITER requirement. This is in contrast to diverted shutdowns, the majority of which disrupt at I N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ∼100 ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real-time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the disruption free protocol, in DIII-D’s 2019–20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n = 1 or locked modes.
Integration of divertor detachment with a high-performance (βN ∼ 3, βp > 2, H98 ∼ 1.5) core plasma has been demonstrated in DIII-D high-βp (poloidal beta) plasmas associated with a sustained core internal transport barrier (ITB) and an H-mode edge transport barrier (ETB). Such good core-edge integration has been achieved for both neon and nitrogen seeding, for both favorable and unfavorable B-field directions, independently from the impurity puffing locations, though these variations play important roles on divertor characteristics. Compared to the standard H-mode plasmas, the high-βp plasma exhibits a much wider window of detachment compatible with high confinement core. Fully detached divertor plasmas with low plasma temperature (Te < 5 eV), low particle flux, and low heat flux across the entire divertor target plate were obtained by using nitrogen seeding. This detached high-βp plasma is compatible with a newly developed detachment control system which can help optimize the nitrogen gas flow rate. Several features, i.e., the high edge safety factor in the high-βp scenario, impurity injection, closed divertor and reduced heating power requirement due to the high confinement, facilitate the achievement of full divertor detachment at lower density. Instead of degrading global performance, the divertor detachment facilitates the access to an even stronger ITB at large radius with a relatively weak ETB through self-organized synergy between ITB and ETB, leading to sustained high confinement. The strengthening of the large-radius ITB compensates for the ETB degradation associated with divertor detachment. In addition, a weak ETB naturally has smaller edge localized modes (ELMs). In particular, with neon injection, a long-period no-ELM H-mode phase has been achieved simultaneously with high-performance core and partially detached divertor plasmas. These results demonstrate the possibility of integrating excellent core plasma performance with an effective divertor solution, an essential step toward steady-state operation of reactor-grade plasmas.
Future tokamaks will require robust technologies for the mitigation of heat exhaust onto the plasma-facing components. As a first step towards this development, a system has been developed at DIII-D that estimates and controls in real-time the heat flux to the PFCs. Real-time estimation of the peak power flux from this model -based approach is validated with off-line infra-red measurements for various DIII-D plasma discharges. A nonlinear free-boundary simulation code (GSevolve) is used for simulating the closed loop response and for the off-line determination of the control parameters. The implementation and first experimental results of the application during the DIII-D plasma campaign are reported.
Recent control advances at KSTAR enabled us to not only establish but also to sustain the ITER-similar shape (ISS) in a stationary manner, matching the ITER-like dimensionless shape parameters in lower single null configuration with vastly contrasting upper/lower triangularities. The kinetic parameters of the ISS are different from the typical KSTAR configurations seen in ELM-RMP suppression experiments. We experimentally demonstrate the ISS-compatible resonant magnetic perturbation (RMP) control in KSTAR for the first time, using the n = 2, +90 degrees phasing RMP.
Divertor detachment offers a promising solution to the challenge of plasma-wall interactions for steady-state operation of fusion reactors. Here, we demonstrate the excellent compatibility of actively controlled full divertor detachment with a high-performance ( β N ~ 3, H 98 ~ 1.5) core plasma, using high-β p (poloidal beta, β p > 2) scenario characterized by a sustained core internal transport barrier (ITB) and a modest edge transport barrier (ETB) in DIII-D tokamak. The high- β p high-confinement scenario facilitates divertor detachment which, in turn, promotes the development of an even stronger ITB at large radius with a weaker ETB. This self-organized synergy between ITB and ETB, leads to a net gain in energy confinement, in contrast to the net confinement loss caused by divertor detachment in standard H-modes. These results show the potential of integrating excellent core plasma performance with an efficient divertor solution, an essential step towards steady-state operation of reactor-grade plasmas.
Robust disruption avoidance techniques are critical for the development of reliable fusion reactor devices. A viable reactor will require non-disruptive, long pulse operation where simply shutting down a discharge is undesirable. To achieve such performance, the plasma must be controlled to continuously avoid hazardous regimes instead of asynchronously aborting. A recent experiment on DIII-D demonstrated for the first time real-time control of proximity to a disruptive instability boundary. In particular, the vertical growth rate, an eigenvalue that characterizes the degree of instability of the plasma's vertical position, was regulated so as not to exceed DIII-D's vertical controllability limit. The open-loop growth rate was estimated in real time on the DIII-D plasma control system using a neural network model trained with tens of thousands of DIII-D shots. The model was trained to replicate the results of RZRIG [1], a rigid displacement code for calculating the growth rate. Once trained, producing an estimate using the neural network is multiple orders of magnitude faster than RZRIG, thereby making the calculation suitable for real-time execution. The control system regulated the estimated growth rate by adjusting plasma elongation and distance to the inner wall of the vessel, and this regulation was shown to reliably avoid vertical displacement event disruptions (i.e. uncontrolled vertical oscillations) of the plasma. This work presents these experimental results, including the dynamic performance and the effectiveness of the control technique. Details are presented on the training of the neural network model, including concerns such as hyperparameter tuning and uncertainty quantification. Additionally, the methodology for embedding the neural network into the control system is discussed.