Wenyu Shi

Toroidal current profile control during low confinement mode plasma discharges in DIII-D via first-principles-driven model-based robust control synthesis

In order for ITER to be capable of operating in advanced tokamak operating regimes, characterized by a high fusion gain, good plasma confinement, magnetohydrodynamic stability and a non-inductively driven plasma current, for extended periods of time, several challenging plasma control problems still need to be solved. Setting up a suitable toroidal current density profile in the tokamak is key for one possible advanced operating scenario characterized by non-inductive sustainment of the plasma current. At the DIII-D tokamak, the goal is to create the desired current profile during the ramp-up and early flat-top phases of the plasma discharge and then actively maintain this target profile for the remainder of the discharge. The evolution in time of the toroidal current profile in tokamaks is related to the evolution of the poloidal magnetic flux profile, which is modelled in normalized cylindrical coordinates using a first-principles, nonlinear, dynamic partial differential equation (PDE) referred to as the magnetic diffusion equation. The magnetic diffusion equation is combined with empirical correlations developed from physical observations and experimental data from DIII-D for the electron temperature, the plasma resistivity and the non-inductive current drive to develop a simplified, control-oriented, nonlinear, dynamic PDE model of the poloidal flux profile evolution valid for low confinement mode discharges. In this work, we synthesize a robust feedback controller to reject disturbances and track a desired reference trajectory of the poloidal magnetic flux gradient profile by employing the control-oriented model of the system. A singular value decomposition of the static gain matrix of the plant model is utilized to identify the most relevant control channels and is combined with the dynamic response of system around a given operating trajectory to design the feedback controller. A general framework for real-time feedforward + feedback control of magnetic and kinetic plasma profiles was implemented in the DIII-D Plasma Control System and was used to demonstrate the ability of the feedback controller to control the toroidal current profile evolution in the DIII-D tokamak. These experiments constitute the first time ever a first-principles-driven, model-based, closed-loop magnetic profile controller was successfully implemented and tested in a tokamak device.

Integrated magnetic and kinetic control of advanced tokamak plasmas on DIII-D based on data-driven models

The first real-time profile control experiments integrating magnetic and kinetic variables were performed on DIII-D in view of regulating and extrapolating advanced tokamak scenarios to steady-state devices and burning plasma experiments. Device-specific, control-oriented models were obtained from experimental data using a generic two-time-scale method that was validated on JET, JT-60U and DIII-D under the framework of the International Tokamak Physics Activity for Integrated Operation Scenarios (Moreau et al 2011 Nucl. Fusion 51 063009). On DIII-D, these data-driven models were used to synthesize integrated magnetic and kinetic profile controllers. The neutral beam injection (NBI), electron cyclotron current drive (ECCD) systems and ohmic coil provided the heating and current drive (H&CD) sources. The first control actuator was the plasma surface loop voltage (i.e. the ohmic coil), and the available beamlines and gyrotrons were grouped to form five additional H&CD actuators: co-current on-axis NBI, co-current off-axis NBI, counter-current NBI, balanced NBI and total ECCD power from all gyrotrons (with off-axis current deposition). Successful closed-loop experiments showing the control of (a) the poloidal flux profile, Ψ(x), (b) the poloidal flux profile together with the normalized pressure parameter, βN, and (c) the inverse of the safety factor profile, , are described.

Physics-based control-oriented modeling of the safety factor profile dynamics in high performance tokamak plasmas

The tokamak is a device that utilizes magnetic fields to confine a reactant gas to generate energy from nuclear fusion reactions. The next step towards the realization of a tokamak power plant is the ITER project, and extensive research has been conducted to find high performance operating scenarios characterized by a high fusion gain and plasma stability. A key property related to both the stability and performance of the plasma is the safety factor profile (q-profile). In this work, a general control-oriented physics-based modeling approach is developed, with emphasis on high performance scenarios, to convert the first-principles physics model that describes the q-profile evolution in the tokamak into a form suitable for control design, with the goal of developing closed-loop controllers to drive the q-profile to a desired target evolution. The DINA-CH& CRONOS and PTRANSP advanced tokamak simulation codes are used to tailor the first-principles-driven (FPD) model to the ITER and DIII-D tokamak geometries, respectively. The models prediction capabilities are illustrated by comparing the prediction to simulated data from DINA-CH&CRONOS for ITER and to experimental data for DIII-D.

Multivariable robust control of the plasma rotational transform profile for advanced tokamak scenarios in DIII-D

The tokamak is a high order, distributed parameter, nonlinear system with a large number of instabilities. Therefore, accurate theoretical plasma models are difficult to develop. However, linear plasma response models around a particular equilibrium can be developed by using data-driven modeling techniques. This paper introduces a linear model of the rotational transform ι profile evolution based on experimental data from the DIII-D tokamak. The model represents the response of the ι profile to the electric field due to induction as well as to heating and current drive (H&CD) systems. The control goal is to use both induction and H&CD systems to regulate the plasma ι profile around a particular target profile. A singular value decomposition (SVD) of the plasma model at steady state is carried out to decouple the system and identify the most relevant control channels. A mixed sensitivity H∞ control design problem is formulated to synthesize a stabilizing feedback controller without input constraint that minimizes the reference tracking error and rejects external disturbances with minimal control energy. The feedback controller is then augmented with an anti-windup compensator, which keeps the given profile controller well-behaved in the presence of magnitude constraints in the actuators and leaves the nominal closed-loop unmodified when no saturation is present. Finally, computer simulations and experimental results illustrate the performance of the model-based profile controller.

A two-time-scale model-based combined magnetic and kinetic control system for advanced tokamak scenarios on DIII-D

System identification techniques have been successfully used to obtain linear dynamic plasma response models around a particular equilibrium in different tokamaks. This paper identifies a two-time-scale dynamic model of the rotational transform ι profile and βN in response to the electric field due to induction as well as to heating and current drive (H&CD) systems based on experimental data from DIII-D. The control goal is to regulate the plasma ι profile and βN around a particular target value. A singular value decomposition (SVD) of the plasma model at steady state is carried out to decouple the system and identify the most relevant control channels. A mixed sensitivity H∞ control design problem is solved to determine a stabilizing feedback controller that minimizes the reference tracking error and rejects external disturbances with minimal control energy. The feedback controller is augmented with an anti-windup compensator, which keeps the given controller well-behaved in the presence of magnitude constraints in the actuators and leaves the nominal closed-loop unmodified when no saturation is present. Experimental results illustrate the performance of the proposed controller, which is one of the first profile controllers integrating magnetic and kinetic variables ever implemented in DIII-D.

Optimal feedback control of the poloidal magnetic flux profile in the DIII-D tokamak based on identified plasma response models

First-principles predictive models based on flux-averaged transport equations often yield complex expressions not suitable for real-time control implementations. It is however always possible to reduce these models to forms suitable for control design while preserving the dominant physics of the system. If further model simplification is desired at the expense of less model accuracy and controller capability, data-driven modeling emerges as an alternative to first-principles modeling. System identification techniques have the potential of producing low-complexity, linear models that can capture the system dynamics around an equilibrium point. This paper focuses on the control of the poloidal magnetic flux profile evolution in response to the heating and current drive (H&CD) systems and the total plasma current. Open-loop data for model identification is collected during the plasma current flattop in a high-confinement scenario (H-mode). Using this data a linear state-space plasma response model for the poloidal magnetic flux profile dynamics around a reference profile is identified. The control goal is to use the H&CD systems and the plasma current to regulate the magnetic profile around a desired target profile in the presence of disturbances. The target profile is defined close enough to the reference profile used for system identification in order to stay within the range of validity of the identified model. An optimal state feedback controller with integral action is designed for this purpose. Experimental results showing the performance of the proposed controller implemented in the DIII-D tokamak are presented.

Physics-model-based nonlinear actuator trajectory optimization and safety factor profile feedback control for advanced scenario development in DIII-D

DIII-D experimental results are reported to demonstrate the potential of physics-model-based safety factor profile control for robust and reproducible sustainment of advanced scenarios. In the absence of feedback control, variability in wall conditions and plasma impurities, as well as drifts due to external disturbances, can limit the reproducibility of discharges with simple pre-programmed scenario trajectories. The control architecture utilized is a feedforward + feedback scheme where the feedforward commands are computed off-line and the feedback commands are computed on-line. In this work, a first-principles-driven (FPD), physics-based model of the q profile and normalized beta () dynamics is first embedded into a numerical optimization algorithm to design feedforward actuator trajectories that steer the plasma through the tokamak operating space to reach a desired stationary target state that is characterized by the achieved q profile and . Good agreement between experimental results and simulations demonstrates the accuracy of the models employed for physics-model-based control design. Second, a feedback algorithm for q profile control is designed following an FPD approach, and the ability of the controller to achieve and maintain a target q profile evolution is tested in DIII-D high confinement (H-mode) experiments. The controller is shown to be able to effectively control the q profile when is relatively close to the target, indicating the need for integrated q profile and control to further enhance the ability to achieve robust scenario execution. The ability of an integrated q profile + feedback controller to track a desired target is demonstrated through simulation.

Experimental and Simulation Testing of Physics-model-based Safety Factor Profile and Internal Energy Feedback Controllers in DIII-D Advanced Tokamak Scenarios

Abstract Active closed-loop control of the plasma safety factor profile ( q -profile) and internal energy dynamics in nuclear fusion tokamak devices has the potential to significantly impact the success of the ITER project. These plasma properties are related to both the stability and performance of a given plasma operating scenario. In this work, we develop integrated feedback control algorithms to control the q -profile and internal energy dynamics in DIII-D advanced tokamak (high performance) scenarios. The feedback controllers are synthesized by embedding a nonlinear, physics-based, control-oriented partial differential equation model of the plasma dynamics into the control design and to be robust to uncertainties in the plasma electron density, electron temperature, and plasma resistivity profiles. The auxiliary heating and current-drive system and the total plasma current are the actuators utilized by the feedback controllers to control the plasma dynamics. Finally, the feedback controllers are tested both through simulations based on the physics-based model and experimentally in the DIII-D tokamak.

Nonlinear Physics-Model-Based Actuator Trajectory Optimization for Advanced Scenario Planning in the DIII-D Tokamak

Abstract Extensive research has been conducted to find operating scenarios that optimize the plasma performance in nuclear fusion tokamak devices with the goal of enabling the success of the ITER project. The development, or planning, of these advanced scenarios is traditionally investigated experimentally by modifying the tokamaks actuator trajectories, such as the auxiliary heating/current-drive (H&CD) scheme, and analyzing the resulting plasma evolution. In this work, a numerical optimization algorithm is developed to complement the experimental effort of advanced scenario planning in the DIII-D tokamak. Two properties related to the plasma stability and performance are the safety factor profile ( q -profile) and the normalized plasma beta ( β N ). The optimization algorithm goal is to design actuator trajectories that steer the plasma to a target q -profile and plasma β N , such that the achieved state is stationary in time, subject to the plasma dynamics (described by a physics-based, nonlinear, control-oriented partial differential equation model) and practical plasma state and actuator constraints, such as the maximum available amount of H&CD power. This defines a nonlinear, constrained optimization problem that we solve by employing sequential quadratic programming. The optimized trajectories are then tested through simulation with the physics-based model and experimentally in DIII-D.

Simultaneous Boundary and Distributed Feedback Control of the Current Profile in H-mode Discharges on DIII-D

Abstract Control of the current profile in tokamak plasmas has been shown to play an important role in achieving advanced scenarios that could enable steady-state operation. The nonlinearity and spatially distributed nature of the current profile dynamics motivate the use of model-based control designs. In this work, we consider a control-oriented model of the current profile evolution in DIII-D high-confinement (H-mode) discharges, and the problem of regulating the current profile around a desired trajectory. The PDE model is discretized in space with a finite difference method and a backstepping design is applied to obtain a transformation from the original system into a particular target system with desirable properties. The resulting boundary condition control law is complemented with control laws for the available distributed actuators. The combined control strategy uses nonlinear combinations of the total plasma current, total power, and line averaged density as actuators. Simulation and experimental results show the ability of the controller to track desired targets and to reject input disturbances.