Graham McArdle

The MAST-U plasma control system platform had become increasingly difficult to maintain due to the presence of legacy devices lacking modern operating system support, resulting in a backlog of dependent software and execution framework upgrades. The original proposed approach to upgrading would have required multiple complex transitions to interim solutions that could not be easily reversed if unsuccessful, thus presenting a high risk to operations. To bypass these complexities, a complete “drop-in” platform replacement was developed, comprising a new control computer running the latest execution software framework on a new operating system with new hardware devices to interface with the same plant input and output connections. This “clean break” approach eliminated all the development bottlenecks caused by the chain of legacy dependencies and allowed us to skip over many years' worth of pending incremental upgrades. The strategy was fundamentally reliant on building a complete end-to-end replacement that could be installed and operated in parallel with and independently of the existing system, with pin-compatible pluggable input and output interfaces and retaining the control logic and user experience of the original system. By decoupling the installation from the commissioning phase, the team was able to perform exhaustive in-situ testing prior to the final switchover. This strategy transformed a high-risk system overhaul into a controlled, reversible event, resulting in a significantly shortened recommissioning window and a seamless transition to a modern, maintainable platform.

In this work, we present the experimental realization of a plasma with a steady negative triangularity (NT) in a spherical tokamak (MAST-U) for the first time. An average triangularity of − 0.1 was maintained for 100 ms in a double-null plasma with a plasma current of 600 kA, a neutral beam injected power of 3.2 MW, a normalized β of 2, and an elongation of 2 on the MAST-U spherical tokamak. Despite many constraints, this was made possible due to the development of a dedicated plasma shape control scheme and the validation and assessment of magnetic shape controllers in simulations. This effort lays the foundation for physics exploitation of the NT plasma in a spherical tokamak.

We present the FreeGSNKE pulse design tool (FPDT), an open-source, Python-based computational framework that enables in silico testing and predictive design of tokamak plasma scenarios and control strategies. The FPDT couples the FreeGSNKE evolutive equilibrium solver with a virtual plasma control system (PCS) containing modular and customisable controllers. Given a set of user-defined waveforms and control parameters, the virtual PCS uses feedback and feedforward control to modulate plasma current, position, and shape, while adhering to machine safety limits on poloidal field (PF) coil currents and voltages. The resulting framework allows simulation of the controlled dynamic evolution of plasma equilibria, along with the currents in both active PF coils and passive conducting structures, under the assumption of axisymmetry. The FPDT can be used to develop plasma scenarios, test control schemes, calibrate control parameters, and perform uncertainty quantification studies, thereby reducing iterative and expensive experimental testing on a physical tokamak. The FPDT is machine-agnostic and can be customised to implement different control algorithms tailored to the specific tokamak of interest. Here, we outline the overall framework and validate its performance on plasma discharges on the MAST Upgrade tokamak in the ‘flat-top’ phase. We demonstrate excellent quantitative agreement between the FPDT simulations, the desired control waveforms, and the experimental shot data. With this extension to the FreeGSNKE open-source suite of codes we aim to encourage more reproducible and collaborative research in plasma modeling and control.

The research programme performed on the Mega Amp Spherical Tokamak (MAST) Upgrade experiment has made significant advances in developing the physics understanding of low aspect ratio tokamaks in support of the operation of ITER and design of fusion powerplants. High performance plasma scenarios have been developed to facilitate a broad programme of experiments, in which confinement is constrained by the presence of m / n = 2/1 modes that cause substantial losses of fast ions. The onset of these modes coincides with the q = 2 surface residing in a local minimum in the toroidal current density profile. The maximum electron temperature at the pedestal top, T _e,ped is limited with gas fuelling to ∼350 eV to maintain regular ELMs; higher T _e,ped results in a transition to a non-stationary ELM-free regime. The operational space of spherical tokamaks has been expanded into small and ELM-free regimes. Strong shaping of the last closed flux surface can induce a transition from large to small ELMs, and ELM suppression with resonant magnetic perturbations has been observed for the first time in a low aspect ratio tokamak. Negative triangularity shaping has induced a transition from ELMy H-mode to a high-performance L-mode regime for the first time in a low aspect ratio tokamak. In studies of fast ion confinement, losses of fast particles due to Global Alfvén Eigenmodes have been identified. Interactions between fast ions generated by off-axis neutral beam injection and thermal neutrals can result in significant losses of fast ions. Experiments with on- and off-axis neutral beam injection exhibit a flux pumping mechanism, where the central safety factor is held to ∼1 in the absence of sawteeth. In studies of pedestal physics, it has been found that elevated main chamber neutral pressures result in an increase in the electron density and reduction in the temperature at the pedestal top. Advances in understanding plasma exhaust include the integration of a high-performance plasma core with detached outer divertors in the X -point target configuration. A newly commissioned lower divertor cryopump reduces the lower divertor neutral pressure by up to 50%, with minimal effect on the main chamber or upper divertor. New measurements and SOLPS-ITER simulations emphasise the importance of plasma–neutral interactions on divertor detachment in the conditions accessible in experiments. Real-time control of the ionisation front location in both divertor chambers independently has been demonstrated in double null experiments, enabled by the tightly baffled divertor chambers.

Nuclear fusion could offer clean, abundant energy. However, managing the power exhausted from the core fusion plasma towards the reactor wall remains a major challenge. This is compounded in emerging compact reactor designs promising more cost-effective pathways towards commercial fusion energy. Alternative Divertor Configurations (ADCs) are a potential solution. In this work, we demonstrate exhaust control in ADCs, employing a novel method to diagnose the neutral gas buffer, which shields the target. Our work on the Mega Ampere Spherical Tokamak Upgrade shows that ADCs tackle key risks and uncertainties for fusion energy. Their highly reduced sensitivity to perturbations enables active exhaust control in otherwise unfeasible situations and facilitates an increased passive absorption of transients, which would otherwise damage the divertor. We observe a strong decoupling of each divertor from other reactor regions, enabling near-independent control of the divertors and core plasma. Our work showcases the real-world benefits of ADCs for effective heat load management in fusion power reactors.

In this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

A key aspect in the modelling of magnetohydrodynamic (MHD) equilibria in tokamak devices is having access to fast, accurate, and stable numerical simulation methods. There is an increasing demand for reliable methods that can be used to develop traditional or machine learning-based shape control feedback systems, optimise scenario designs, and integrate with other plasma edge or transport modelling codes. To handle such applications, these codes need to be flexible and, more importantly, they need to have been validated against both analytically known and real-world tokamak equilibria to ensure they are consistent and credible. In this paper, we are interested in solving the static forward Grad–Shafranov (GS) problem for free–boundary MHD equilibria. Our focus is on the validation of the static forward solver in the Python-based equilibrium code FreeGSNKE by solving equilibria from magnetics-only EFIT ++ reconstructions of MAST-U shots. In addition, we also validate FreeGSNKE against equilibria simulated using the well-established MATLAB-based equilibrium code Fiesta. To do this, we develop a computational pipeline that allows one to load the same (a)symmetric MAST-U machine description into each solver, specify the required inputs (active/passive conductor currents, plasma profiles and coefficients, etc.) from EFIT++, and solve the GS equation for all available time slices across a shot. For a number of different MAST-U shots, we demonstrate that both FreeGSNKE and Fiesta can successfully reproduce various poloidal flux quantities and shape targets (e.g. midplane radii, magnetic axes, separatrices, X-points, and strikepoints) in agreement with EFIT++ calculations to a very high degree of accuracy. We also provide public access to the code/data required to load the MAST-U machine description in FreeGSNKE/Fiesta and reproduce the equilibria in the shots shown.

A plasma scenario with negative triangularity (NT) shaping is achieved on MAST-U for the first time. While edge localized modes (ELMs) are eventually suppressed as the triangularity is decreased below δ ≲ − 0.06 , an extended period of H-mode operation with Type-III ELMs is sustained at less negative δ even through access to the second stability region for ideal ballooning modes is closed. This documents a qualitative difference from the ELM-free access conditions documented in NT scenarios on conventional aspect ratio machines. The electron temperature at the pedestal top drops across the transition to ELM-free operation, but a steady rise in core temperature as δ is decreased allows for similar normalized β in the ELM-free NT and H-mode positive triangularity shapes.

Recent results from MAST Upgrade are presented, emphasising understanding the capabilities of this new device and deepening understanding of key physics issues for the operation of ITER and the design of future fusion power plants. The impact of MHD instabilities on fast ion confinement have been studied, including the first observation of fast ion losses correlated with Compressional and Global Alfvén Eigenmodes. High-performance plasma scenarios have been developed by tailoring the early plasma current ramp phase to avoid internal reconnection events, resulting in a more monotonic q profile with low central shear. The impact of m/n = 3/2, 2/1 and 1/1 modes on thermal plasma confinement and rotation profiles has been quantified, and scenarios optimised to avoid them have transiently reached values of normalised beta approaching 4.2. In pedestal and ELM physics, a maximum pedestal top temperature of ∼350 eV has been achieved, exceeding the value achieved on MAST at similar heating power. Mitigation of type-I ELMs with n = 1 RMPs has been observed. Studies of plasma exhaust have concentrated on comparing conventional and Super-X divertor configurations, while X-point target, X-divertor and snowflake configurations have been developed and studied in parallel. In L-mode discharges, the separatrix density required to detach the outer divertors is approximately a factor 2 lower in the Super-X than the conventional configuration, in agreement with simulations. Detailed analysis of spectroscopy data from studies of the Super-X configuration reveal the importance of including plasma-molecule interactions and D₂ Fulcher band emission to properly quantify the rates of ionisation, plasma-molecule interactions and volumetric recombination processes governing divertor detachment. In H-mode with conventional and Super-X configurations, the outer divertors are attached in the former and detached in the latter with no impact on core or pedestal confinement.

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.