OpenStar Technologies is pursuing the levitated dipole (LD) as a highly modular, loosely-coupled system that leverages their expertise in high temperature superconductor (HTS) technology. The next generation experiment at OpenStar, Tahi (Ma¯ori for “first”), will demonstrate the generation and confinement of fast ions in a levitated dipole for the first time. Ion cyclotron range of frequency (ICRF) heating is a leading candidate for energetic ion formation in Tahi. A frequency in the 10 MHz range will be used for H minority heating or D majority heating with waves launched from an antenna located above the floating coil. Unlike a tokamak, where the targeted cyclotron resonance is typically a vertical path through the center of the plasma, in a dipole the resonance location follows a C-shaped path from the separatrix to the center of the plasma. The value of B also varies significantly within the confined plasma resulting in a large number of cyclotron harmonics present in the low field region. Furthermore, levitated dipoles contain a “first closed flux surface” surrounding the floating coil, in addition to the traditional separatrix /last closed flux surface. Simulation e fforts using full-wave ICRF codes show that ICRF heating of a levitated dipole reactor is feasible using a pair of toroidal current straps phased to launch the appropriate parallel ( i.e. poloidal) refractive index.

The tokamak offers a promising path to fusion energy, but disruptions pose a major economic risk, motivating solutions to manage their consequence. This work develops a reinforcement learning approach to this problem by training a policy to ramp-down the plasma current while avoiding limits on a number of quantities correlated with disruptions. The policy training environment is a hybrid physics and machine learning model trained on simulations of the SPARC primary reference discharge (PRD) ramp-down, an upcoming burning plasma scenario which we use as a testbed. To address physics uncertainty and model inaccuracies, the simulation is massively parallelized on GPU with randomized physics parameters during policy training. The trained policy is then run in feedback on a transport simulator as a demonstration. We also directly address the crucial issue of control validation by demonstrating that a constraint-conditioned policy can be a trajectory design assistant that designs a library of feed-forward trajectories to handle different physics conditions and user constraint settings, a promising approach for the sensitive context of burning plasma tokamaks. Finally, we demonstrate that the training environment can be a useful platform for feed-forward optimization approaches by optimizing feed-forward trajectories that are robust to physics uncertainty.

To achieve its performance goals, SPARC plans to operate in equilibrium configurations with a strong elongation of $\kappa_\mathrm{areal}\sim1.75$, destabilizing the $n=0$ vertical instability. However, SPARC also features a relatively thick conducting wall that is designed to withstand disruption forces, leading to lower vertical instability growth rates than usually encountered. In this work, we use the TokSyS framework to survey families of accessible shapes near the SPARC baseline configuration, finding maximum growth rates in the range of $\gamma\lesssim100\,$s$^{-1}$. The addition of steel vertical stability plates has only a modest ($\sim25\%$) effect on reducing the vertical growth rate and almost no effect on the plasma controllability when the full vertical stability system is taken into account, providing flexibility in the plate conductivity in the SPARC design. Analysis of the maximum controllable displacement on SPARC is used to inform the power supply voltage and current limit requirements needed to control an initial vertical displacement of $5\%$ of the minor radius. From the expected spectra of plasma disturbances and diagnostic noise, requirements for filter latency and vertical stability coil heating tolerances are also obtained. Small modifications to the outboard limiter location are suggested to allow for an unmitigated vertical disturbance as large as $5\%$ of the minor radius without allowing the plasma to become limited. Further, investigations with the 3D COMSOL code reveal that strategic inclusion of insulating structures within the VSC supports are needed to maintain sufficient magnetic response. The workflows presented here help to establish a model for the integrated predictive design for future devices by coupling engineering decisions with physics needs.

A new three-dimensional electromagnetic modeling tool ThinCurr has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vessel and other conducting structural elements with realistic material resistivities. This model was leveraged to support the design of a passive runaway electron mitigation coil (REMC), studying the effect of various design parameters, including coil resistivity, current quench duration, and plasma vertical position, on the effectiveness of the coil. The REMC is a non-axisymmetric coil designed to passively drive large non-axisymmetric fields during the plasma disruption thereby destroying flux surfaces and deconfining RE seed populations. These studies indicate that current designs should apply substantial 3D fields at the plasma surface during future plasma current disruptions as well as highlight the importance of having the REMC conductors away from the machine midplane in order to ensure they are robust to off-normal disruption scenarios.

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 tokamak offers a promising path to fusion energy, but plasma disruptions pose a major economic risk, motivating considerable advances in disruption avoidance. This work develops a reinforcement learning approach to this problem by training a policy to safely ramp-down the plasma current while avoiding limits on a number of quantities correlated with disruptions. The policy training environment is a hybrid physics and machine learning model trained on simulations of the SPARC primary reference discharge (PRD) ramp-down, an upcoming burning plasma scenario which we use as a testbed. To address physics uncertainty and model inaccuracies, the simulation environment is massively parallelized on GPU with randomized physics parameters during policy training. The trained policy is then successfully transferred to a higher fidelity simulator where it successfully ramps down the plasma while avoiding user-specified disruptive limits. We also address the crucial issue of safety criticality by demonstrating that a constraint-conditioned policy can be used as a trajectory design assistant to design a library of feed-forward trajectories to handle different physics conditions and user settings. As a library of trajectories is more interpretable and verifiable offline, we argue such an approach is a promising path for leveraging the capabilities of reinforcement learning in the safety-critical context of burning plasma tokamaks. Finally, we demonstrate how the training environment can be a useful platform for other feed-forward optimization approaches by using an evolutionary algorithm to perform optimization of feed-forward trajectories that are robust to physics uncertainty

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ∼1MA compared to ∼8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a<0.3 and with kinetic energies ∼0.2 – 15MeV , contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ∼18m2s−1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field’s effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

Accurate reconstruction of the plasma equilibrium is imperative for successful operation of the SPARC tokamak. In order to assess the expected reconstruction accuracy throughout the duration of design-point discharges, the EFIT equilibrium reconstruction code was deployed for SPARC. Reconstructions from SPARC baseline scenarios were compared with free-boundary equilibria generated by FreeGS, Toksys, and the Tokamak Simulation Code. The key geometric areas of interest, where design constraints are imposed, included: the inner and outer midplane gaps, the X-point locations, as well as the strike point locations. Successful reconstructions of various reference discharges, using deviations in these key geometric quantities as metrics, were calculated from synthetic signals considering an optimized equilibrium magnetic sensor set. The optimization process for this sensor set combined a scan of randomized sensor placement with a linear perturbation analysis to determine critical sensor locations, while simultaneously conforming to design constraints on the sensor placement. This optimized set was also successful in performing equilibrium reconstructions with the addition of error to synthetic measurements of magnetic flux and magnetic field, as well as contributions from eddy currents in conducting structures. These methods represent a workflow of optimization and validation that balances the engineering constraints of sensor placement with achieving sufficient reconstruction fidelity for science and operations missions for SPARC.

While fusion reactors known as tokamaks hold promise as a firm energy source, advances in plasma control, and handling of events where control of plasmas is lost, are needed for them to be economical. A significant bottleneck towards applying more advanced control algorithms is the need for better plasma simulation, where both physics-based and data-driven approaches currently fall short. The former is bottle-necked by both computational cost and the difficulty of modelling plasmas, and the latter is bottle-necked by the relative paucity of data. To address this issue, this work applies the neural ordinary differential equations (ODE) framework to the problem of predicting a subset of plasma dynamics, namely the coupled plasma current and internal inductance dynamics. As the neural ODE framework allows for the natural inclusion of physics-based inductive biases, we train both physics-based and neural network models on data from the Alcator C-Mod fusion reactor and find that a model that combines physics-based equations with a neural ODE performs better than both existing physics-motivated ODEs and a pure neural ODE model.

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.

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q -profile evolution is seen to be critical to predicting the later time effects of the coil.

The SPARC tokamak project, currently in engineering design, aims to achieve breakeven and burning plasma conditions in a compact device, thanks to new developments in high-temperature superconductor technology. With a magnetic field of 12.2 T on axis and 8.7 MA of plasma current, SPARC is predicted to produce 140 MW of fusion power with a plasma gain of Q ≈ 11, providing ample margin with respect to its mission of Q > 2. All tokamak systems are being designed to produce this landmark plasma discharge, thus enabling the study of burning plasma physics and tokamak operations in reactor relevant conditions to pave the way for the design and construction of a compact, high-field fusion power plant. Construction of SPARC is planned to begin by mid-2021.

The potential formation of multi-mega-ampere beams of relativistic 'runaway' electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium-tritium plasmas.

SPARC is being designed to operate with a normalized beta of $\beta _N=1.0$, a normalized density of $n_G=0.37$ and a safety factor of $q_{95}\approx 3.4$, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $\beta _p=0.19$ at the safety factor $q=2$ surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ${\sim }80$ %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order $10^{-2}$ that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.