Density regulation with disruption avoidance in next-generation tokamaks using a safe reinforcement learning-based controller
Not Available
Engineering topics
Publications and source records attributed to Pajares, Andres.
Not Available
Simultaneous regulation of multiple properties in next-generation tokamaks like ITER and fusion pilot plant may require the integration of different plasma control algorithms. Such integration requires the conversion of individual controller commands into physical actuator requests while accounting for the coupling between different plasma properties. This work proposes a tokamak and scenario-agnostic actuator-sharing algorithm (ASA) to perform the above-mentioned command-request conversion and, hence, integrate multiple plasma controllers. The proposed algorithm implicitly solves a quadratic programming (QP) problem formulated to account for the saturation limits and the relation between the controller commands and physical actuator requests. Since the constraints arising in the QP program are linear, the proposed ASA is highly computationally efficient and can be implemented in the tokamak plasma control system in real time. Furthermore, the proposed algorithm is designed to handle real-time changes in the control objectives and actuators’ availability. Nonlinear simulations carried out using the Control Oriented Transport SIMulator illustrate the effectiveness of the proposed algorithm in achieving multiple control objectives simultaneously.
The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.
Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.