Conceptual design study for heat exhaust management in the ARC fusion pilot plant

The ARC pilot plant conceptual design study has been extended beyond its initial scope [B. N. Sorbom et al., FED 100 (2015) 378] to explore options for managing ~525 MW of fusion power generated in a compact, high field (B_0 = 9.2 T) tokamak that is approximately the size of JET (R_0 = 3.3 m). Taking advantage of ARC's novel design - demountable high temperature superconductor toroidal field (TF) magnets, poloidal magnetic field coils located inside the TF, and vacuum vessel (VV) immersed in molten salt FLiBe blanket - this follow-on study has identified innovative and potentially robust power exhaust management solutions.Comment: Accepted by Fusion Engineering and Desig

Model predictive control of resistive wall mode for ITER

Active feedback stabilization of the dominant resistive wall mode (RWM) for an ITER H-mode scenario at high plasma pressure using infinite-horizon model predictive control (MPC) is presented. The MPC approach is closely-related to linear-quadratic-Gaussian (LQG) control, improving the performance in the vicinity of constraints. The control-oriented model for MPC is obtained with model reduction from a high-dimensional model produced by CarMa code. Due to the limited time for on-line optimization, a suitable MPC formulation considering only input (coil voltage) constraints is chosen, and the primal fast gradient method is used for solving the associated quadratic programming problem. The performance is evaluated in simulation in comparison to LQG control. Sensitivity to noise, robustness to changes of unstable RWM dynamics, and size of the domain of attraction of the initial conditions of the unstable modes are examined.Comment: Original manuscript as submitted to Fusion Engineering and Desig

Controllers for high-performance nuclear fusion plasmas

A succesful nuclear fusion reactor will confine plasma at hig temperatures and densities, with low thermal losses. The workhorse of the nuclear fusion community is the tokamak, a toroidal device in which plasmas are confined by poloidal and toroidal magnetic fields. Ideally, the confirming magnetic fields form a set of nested tori. A repetive magnetohydrodynamic (MHD) event in the plasma core (the sawteeth instability) perturbs the confirming magnetic field by producing seed islands. In low-pressure plasmas the seed islands will self-heal. In high-pressure plasmas the seed islands can grow and saturate. These neoclassical tearing models (NTMs) reduce the plasma performance or lead to plasma disruption. This sets the resistive pressure limit in tokamaks. High-performance operation in tokamaks therefore implies the control or amelioration of the NTMs. Controllers for the sawteeth and the NTMs will be discussed, with special emphasis on the development of dedicated sensors and models for MHD control

The modeling of a tokamak plasma discharge, from first principles to a flight simulator

A newly developed tool to simulate a tokamak full discharge is presented. The tokamak \u27flight simulator\u27 Fenix couples the tokamak control system with a fast and reduced plasma model, which is realistic enough to take into account several of the plasma non-linearities. A distinguishing feature of this modeling tool is that it only requires the pulse schedule (PS) as input to the simulator. The output is a virtual realization of the full discharge, whose time traces can then be used to judge if the PS satisfies control/physics goals or needs to be revised. This tool is envisioned for routine use in the control room before each pulse is performed, but can also be used off-line to correct PS in advance, or to develop and validate reduced models, control schemes for future machines like a commercial reactor, simulating realistic actuators and sensors behavior

Plasma Shape and Current Density Profile Control in Advanced Tokamak Operating Scenarios

The need for new sources of energy is expected to become a critical problem within the next few decades. Nuclear fusion has sufficient energy density to potentially supply the world population with its increasing energy demands. The tokamak is a magnetic confinement device used to achieve controlled fusion reactions. Experimental fusion technology has now reached a level where tokamaks are able to produce about as much energy as is expended in heating the fusion fuel. The next step towards the realization of a nuclear fusion tokamak power plant is ITER, which will be capable of exploring advanced tokamak (AT) modes, characterized by a high fusion gain and plasma stability. The extreme requirements of the advanced modes motivates researchers to improve the modeling of the plasma response as well as the design of feedback controllers. This dissertation focuses on several magnetic and kinetic control problems, including the plasma current, position and shape control, and data-driven and first-principles-driven modeling and control of plasma current density profile and the normalized plasma pressure ratio βN.The plasma is confined within the vacuum vessel by an external electromagnetic field, produced primarily by toroidal and poloidal field coils. The outermost closed plasma surface or plasma boundary is referred to as the shape of the plasma. A central characteristic of AT plasma regimes is an extreme elongated shape. The equilibrium among the electromagnetic forces acting on an elongated plasma is unstable. Moreover, the tokamak performance is improved if the plasma is located in close proximity to the torus wall, which guarantees an efficient use of available volume. As a consequence, feedback control of the plasma position and shape is necessary. In this dissertation, an H∞-based, multi-input-multi-output (MIMO) controller for the National Spherical Torus Experiment (NSTX) is developed, which is used to control the plasma position, shape, and X-point position.Setting up a suitable toroidal current profile is related to both the stability and performance of the plasma. The requirements of ITER motivate the research on plasma current profile control. Currently, physics-based control-oriented modeling techniques of the current profile evolution can be separated into two major classes: data-driven and first-principles-driven. In this dissertation, a two-timescale linear dynamic data-driven model of the rotational transform profile and βN is identified based on experimental data from the DIII-D tokamak. A mixed-sensitivity H∞ controller is developed and tested during DIII-D high-confinement (H-mode) experiments by using the heating and current drive (H&CD) systems to regulate the plasma rotational transform profile and βN around particular target values close to the reference state used for system identification. The preliminary experimental results show good progress towards routine current profile control in DIII-D. As an alternative, a nonlinear dynamic first-principles-driven model is obtained by converting the physics-based model that describes the current profile evolution in H-mode DIII-D discharges into a form suitable for control design. The obtained control-oriented model is validated by comparing the model prediction to experimental data. An H∞ control design problem is formulated to synthesize a stabilizing feedback controller, with the goal of developing a closed-loop controller to drive the current profile in DIII-D to a desirable target evolution. Simulations show that the controller is capable of regulating the system around the target rotational transform profile in the presence of disturbances. When compared to a previously designed data-driven model-based controller, the proposed first-principles-driven model-based controller shows potential for improving the control performance

Plasma Edge Kinetic-MHD Modeling in Tokamaks Using Kepler Workflow for Code Coupling, Data Management and Visualization

A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localized modes (ELMs) is presented in this report. This tool brings together, in a coordinated and effective manner, several first-principles physics simulation codes, stability analysis packages, and data processing and visualization tools. A Kepler workflow is used in order to carry out an edge plasma simulation that loosely couples the kinetic code, XGC0, with an ideal MHD linear stability analysis code, ELITE, and an extended MHD initial value code such as M3D or NIMROD. XGC0 includes the neoclassical ion-electron-neutral dynamics needed to simulate pedestal growth near the separatrix. The Kepler workflow processes the XGC0 simulation results into simple images that can be selected and displayed via the Dashboard, a monitoring tool implemented in AJAX allowing the scientist to track computational resources, examine running and archived jobs, and view key physics data, all within a standard Web browser. The XGC0 simulation is monitored for the conditions needed to trigger an ELM crash by periodically assessing the edge plasma pressure and current density profiles using the ELITE code. If an ELM crash is triggered, the Kepler workflow launches the M3D code on a moderate-size Opteron cluster to simulate the nonlinear ELM crash and to compute the relaxation of plasma profiles after the crash. This process is monitored through periodic outputs of plasma fluid quantities that are automatically visualized with AVS/Express and may be displayed on the Dashboard. Finally, the Kepler workflow archives all data outputs and processed images using HPSS, as well as provenance information about the software and hardware used to create the simulation. The complete process of preparing, executing and monitoring a coupled-code simulation of the edge pressure pedestal buildup and the ELM cycle using the Kepler scientific workflow system is described in this paper

Nonlinear Burn Condition and Kinetic Profile Control in Tokamak Fusion Reactors

One of the most promising devices for realizing power production through nuclear fusion is the tokamak. In order to maximize performance, it is preferable that tokamaks achieve operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high fusion gain, the ratio of fusion power produced to the external heating power needed to sustain reactions. There are many experimental tokamaks around the world, each exploring different facets of plasma physics and fusion technology. These experiments have reached the point where the power released from fusion is nearly equal to the power input required to heat the plasma. The next experimental step is ITER, which aims to reach a fusion gain exceeding ten for short pulses, and to sustain a gain of five for longer pulses (around 1000 s). In order for ITER to be a success, several challenging control engineering problems must be addressed.Among these challenges is to precisely regulate the plasma density and temperature, or burn condition. Due to the nonlinear and coupled dynamics of the system, modulation of the burn condition (either during ramp-up/shut-down or in response to changing power demands) without a well designed control scheme could result in undesirable transient performance. Feedback control will also be necessary for responding to unexpected changes in plasma confinement, impurity content, or other parameters, which could significantly alter the burn condition during operation. Furthermore, although stable operating points exist for most confinement scalings, certain conditions can lead to thermal instabilities. Such instabilities can either lead to quenching or a thermal excursion in which the system moves to a higher temperature equilibrium point. In any of these situations, disruptive plasma instabilities could be triggered, stopping operation and potentially causing damage to the confinement vessel.In this work, the problem of burn condition control is addressed through the design of a nonlinear control law guaranteeing stability of desired equilibria. Multiple actuation methods, including auxiliary heating, isotopic fueling, and impurity injection, are used to ensure the burn condition is regulated even when actuators saturate. An adaptive control scheme is used to handle model uncertainty, and an online optimization scheme is proposed to ensure that the plasma is driven to an operating point that minimizes an arbitrary cost function. Due to the possible limited availability of diagnostic systems in ITER and future reactors, an output feedback control scheme is also proposed that combines the nonlinear controller with an observer that estimates the states of the burning plasma system based on available measurements. Finally, the control scheme is tested using the integrated modeling code METIS.The control of spatial profiles of parameters, including current, density, and temperature, is also an important challenge in fusion research, due to their effect on MHD stability, non-inductive current drive, and fusion power. In this work, the problem of kinetic profile control in burning plasmas is addressed through a nonlinear boundary feedback control law designed using a technique called backstepping. A novel implementation of the backstepping technique is used that enables the use of both boundary and interior actuation. The backstepping technique is then applied to the problem of current profile control in both low-confinement and high-confinement mode discharges in the DIII-D tokamak based on a first-principles-driven model of the current profile evolution. Both designs are demonstrated in simulations and experimental tests

MHD stability and disruptions in the SPARC tokamak

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) approximate to 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 similar to 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