Integrated Control in Tokamaks using Nonlinear Robust Techniques and Actuator Sharing Strategies

Tokamaks are devices whose final purpose is obtaining energy by means of nuclear fusion reactions. To achieve such purpose, a gas is injected into the tokamak\u27s torus-shaped chamber and heated up to extremely high temperatures, giving birth to a plasma. When the necessary conditions of temperature, density, and confinement time are achieved, virtually inexhaustible energy can be produced in a tokamak. The main contribution of this dissertation is the development of integrated control strategies for tokamak plasmas. The development of integrated control architectures is necessary for tokamaks to become efficient and commercially competitive power plants. Because a tokamak plasma is a highly nonlinear, coupled dynamical system, the great diversity of complex control problems that coexist in a tokamak are indeed closely interrelated. However, this variety of control problems must be tackled by means of a limited number of actuators. A functional design for integrated tokamak-control architectures should employ multi-input multi-output controllers to simultaneously regulate as many plasma variables as possible with the available actuators. Supervisory and exception handling systems that monitor the plasma state arise as a necessity to ensure a safe tokamak-operation. Finally, actuator sharing and management capabilities should also exist in order to utilize the available actuators in an optimal way. Various control problems are tackled in this dissertation, including kinetic, magnetic, and instability control problems. Control-oriented, physics-based models that characterize some specific aspects of the plasma dynamics have been employed to develop new control-oriented simulation codes and integrated-control solutions that employ nonlinear, robust control techniques and optimization-based actuator-management strategies. Some of those control solutions have been experimentally tested in the DIII-D tokamak

Physics-model-based Optimization and Feedback Control of the Current Profile Dynamics in Fusion Tokamak Reactors

As the demand for energy continues to increase, the need to develop alternative energy sources to complement (and one day replace) conventional fossil fuels is becoming increasingly important. One such energy source is nuclear fusion, which has the potential to provide a clean source of energy and possesses an abundant fuel supply. However, due to the technological difficulty in creating the conditions necessary for controlled fusion to occur, nuclear fusion is not yet commercially viable. The tokamak is a device that utilizes magnetic fields to confine the reactants, which are in the plasma state, and it is one of the most promising devices capable of achieving controlled fusion. The ITER tokamak project is the next phase of tokamak development and will be the first tokamak reactor to explore the burning plasma (one with a significant amount of fusion reactions) operating regime.In order for ITER to meet its demanding goals, extensive research has been conducted to develop advanced tokamak operating scenarios characterized by a high fusion gain, good plasma confinement, magnetohydrodynamic stability, and a significant fraction of noninductively driven plasma current to maximize the plasma performance and potentially enable steady-state operation. As the dynamics of the tokamak plasma magnetic and kinetic states are highly coupled, distributed, nonlinear systems that exhibit many instabilities, it is extremely difficult to robustly achieve advanced operating scenarios. Therefore, active control of the plasma dynamics has significant potential to improve the ability to access advanced operating regimes. One of the key plasma properties investigated in the development of advanced scenarios is the plasma current profile because of its intimate relationship to plasma energy/particle transport and to plasma stability limits that are approached by increasing the plasma pressure. The plasma density and temperature profiles are also important parameters due to their close relationship to the amount of generated fusion power, to the total plasma stored energy, and to the amount of noninductive current drive. In tokamaks, the current and electron temperature profiles are coupled through resistive diffusion, noninductive current drive, and plasma energy/particle transport. As a result, integrated algorithms for current profile and electron temperature profile control will be necessary to maintain plasma stability, optimize plasma performance, and respond to changing power demand in ITER, and eventually a commercial, power producing tokamak reactor.In this work, model-based feedforward and feedback algorithms are developed to control the plasma current profile and thermal state dynamics with the goal of improving the ability to achieve robust tokamak operation. A first-principles-driven (FPD), physics-based approach is employed to develop models of the plasma response to the available actuators, which provides the freedom to handle the trade-off between the physics accuracy and the tractability for control design of the models. A numerical optimization algorithm to synthesize feedforward trajectories for the tokamak actuators that steer the plasma through the tokamak operating space to achieve a predefined target scenario (characterized by a desired current profile and total stored energy), subject to the plasma dynamics (described by the developed physics-based model), actuator constraints, and plasma state constraints, is developed. Additionally, robust feedback control algorithms for current profile, combined current profile + total stored energy, and simultaneous current profile + electron temperature profile control are synthesized for various tokamaks by embedding a FPD model into the control design process.Examples of the performance of the controllers in simulations (DIII-D, ITER, and TCV tokamaks) and DIII-D experiments are presented to illustrate the potential and versatility of the employed control methodology. The DIII-D experimental tests demonstrate the potential physics-model-based profile control has to provide a systematic approach for the development and robust sustainment of advanced scenarios. The ITER simulations demonstrate the ability to drive the current profile to a stationary target while simultaneously modulating the amount of fusion power that is generated. Finally, the TCV simulations demonstrate the ability to drive the current and electron temperature profiles to a self consistent target, as well as to maintain the current profile in a stationary condition while simultaneously modulating the electron temperature profile between equilibrium points

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

Model-based Optimization and Feedback Control of the Current Density Profile Evolution in NSTX-U

Unlike nuclear fission in present nuclear power plants, where energy is generated by splitting heavy atoms like uranium, nuclear fusion generates energy by fusing light nuclei like hydrogen isotopes under high-temperature and high-pressure conditions, at which the reactants (hydrogen isotopes) separate from their electrons and form an ionized gas called plasma, which is considered as the fourth state of matter. Contrary to fission, fusion provides more energy density, poses almost no risk of a catastrophic nuclear accident, and produces mostly short-term, low-level radioactive waste.The main difficulty in maintaining fusion reactions is the development of a device that can confine the hot plasma for sufficiently long time while preventing it from hitting the walls of the confining device. Among several techniques, magnetic confinement appears as the most promising approach. In particular, the tokamak device is a toroidal device surrounded by large magnetic coils responsible for the magnetic fields that confine the plasma. A spherical tokamak, or a spherical torus (ST), is a variation of the conventional tokamak concept. Compared to a standard tokamak, the ST device extrapolates to a more compact, potentially lower-cost reactor with higher efficiency of confinement. Nuclear fusion research is a highly challenging, multidisciplinary field seeking contributions from both plasma physics and multiple engineering areas. As an application of plasma control engineering, this dissertation mainly explores methods to control the current density profile evolution within the National Spherical Torus eXperiment-Upgrade (NSTX-U), which is a substantial upgrade based on the NSTX device, which is located in Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ. Active control of the toroidal current density profile is among those plasma control milestones that the NSTX-U program must achieve to realize its next-step operational goals, which are characterized by high-performance, long-pulse, MHD-stable plasma operation with neutral beam heating. Therefore, the aim of this work is to develop model-based, feedforward and feedback controllers that can enable time regulation of the current density profile in NSTX-U by actuating the total plasma current, electron density, and the powers of the individual neutral beam injectors.Motivated by the coupled, nonlinear, multivariable, distributed-parameter plasma dynamics, the first step towards control design is the development of a physics-based, control-oriented model for the current profile evolution in NSTX-U in response to non-inductive current drives and heating systems. Numerical simulations of the proposed control-oriented model show qualitative agreement with the high-fidelity physics code TRANSP. The next step is to utilize the proposed control-oriented model to design an open-loop actuator trajectory optimizer. Given a desired operating state, the optimizer produces the actuator trajectories that can steer the plasma to such state. The objective of the feedforward control design is to provide a more systematic approach to advanced scenario planning in NSTX-U since the development of such scenarios is conventionally carried out experimentally by modifying the tokamak’s actuator trajectories and analyzing the resulting plasma evolution.Finally, the proposed control-oriented model is embedded in feedback control schemes based on optimal control and Model Predictive Control (MPC) approaches. Integrators are added to the standard Linear Quadratic Gaussian (LQG) and MPC formulations to provide robustness against various modeling uncertainties and external disturbances. The effectiveness of the proposed feedback controllers in regulating the current density profile in NSTX-U is demonstrated in closed-loop nonlinear simulations. Moreover, the optimal feedback control algorithm has been implemented successfully in closed-loop control simulations within TRANSP through the recently developed Expert routine

Full tokamak discharge simulation and kinetic plasma profile control for ITER

Understanding non-linearly coupled physics between plasma transport and free-boundary equilibrium evolution is essential to operating future tokamak devices, such as ITER and DEMO, in the advanced tokamak operation regimes. To study the non-linearly coupled physics, we need a simulation tool which can self-consistently calculate all the main plasma physics, taking the operational constraints into account. As the main part of this thesis work, we have developed a full tokamak discharge simulator by combining a non-linear free-boundary plasma equilibrium evolution code, DINA-CH, and an advanced transport modelling code, CRONOS. This tokamak discharge simulator has been used to study the feasibility of ITER operation scenarios and several specific issues related to ITER operation. In parallel, DINA-CH has been used to study free-boundary physics questions, such as the magnetic triggering of edge localized modes (ELMs) and plasma dynamic response to disturbances. One of the very challenging tasks in ITER, the active control of kinetic plasma profiles, has also been studied. In the part devoted to free-boundary tokamak discharge simulations, we have studied dynamic responses of the free-boundary plasma equilibrium to either external voltage perturbations or internal plasma disturbances using DINA-CH. Firstly, the opposite plasma behaviour observed in the magnetic triggering of ELMs between TCV and ASDEX Upgrade has been investigated. Both plasmas experience similar local flux surface expansions near the upper G-coil set and passive stabilization loop (PSL) when the ELMs are triggered, due to the presence of the PSLs located inside the vacuum vessel of ASDEX Upgrade. Secondly, plasma dynamic responses to strong disturbances anticipated in ITER are examined to study the capability of the feedback control system in rejecting the disturbances. Specified uncontrolled ELMs were controllable with the feedback control systems. However, the specifications for fast H-L mode transitions were not fully achievable due to a vertical displacement event (VDE) caused by a strong inward plasma movement. In the part dedicated to full tokamak discharge simulations, firstly, we have introduced the combined DINA-CH/CRONOS tokamak discharge simulator. DINA-CH self-consistently calculates the non-linear evolution of the free-boundary plasma equilibrium with the plasma current diffusion, in response to both controlled poloidal field (PF) coil currents and inductively driven currents in the surrounding conducting system. CRONOS provides the evolution of the plasma profiles by self-consistently solving heat and particle transport with source profiles. Secondly, we have successfully simulated ITER operation scenario 2 as a demonstration of the capabilities of the combined simulator, as well as being a design study in itself. The fusion power ratio to the total auxiliary power Q was about 10 with the application of 53MW of auxiliary heating and current drive (H&CD) power. We have investigated several specific issues related to the tokamak operation, such as the vertical instability, PF coil current limits and poloidal flux consumption during the current ramp-up. Lower hybrid (LH) applied from the initial phase of the plasma current ramp-up increased the safety margins in operating the superconducting PF coils both by reducing resistive ohmic flux consumption and by providing non-inductively driven plasma current. Lastly, we have studied ITER hybrid mode operation, focusing on the operational capability of obtaining a stationary at safety factor (q/) profile at the start of at-top (SOF) phase and sustaining it as long as possible by combining various non-inductively driven current sources. Application of a near on-axis electron cyclotron current drive (ECCD) appears to be effective compared to the far off-axis lower hybrid current drive (LHCD), at least on short time scales. In the active plasma profile control part, we have developed a robust control technique that simplifies the active real-time control of several kinetic plasma profiles in ITER. The response of the plasma profiles to power changes of auxiliary H&CD systems is modelled. To allow real-time update of the plasma profile response model, the related physics are simplified with several assumptions. The electron temperature profile response is modelled by simplifying the electron heat transport equation. The q profile response is modelled by directly relating it to the changes of source current density profiles. The required actuator power changes are calculated using the singular value decomposition (SVD) technique, taking the saturation of the actuator powers into account. The potential of this control technique has been shown by applying it to simulations of the ITER hybrid mode operation

Free-Boundary Simulations of ITER Advanced Scenarios

The successful operation of ITER advanced scenarios is likely to be a major step forward in the development of controlled fusion as a power production source. ITER advanced scenarios raise specific challenges that are not encountered in presently-operated tokamaks. In this thesis, it is argued that ITER advanced operation may benefit from optimal control techniques. Optimal control ensures high performance operation while guaranteeing tokamak integrity. The application of optimal control techniques for ITER operation is assessed and it is concluded that robust optimisation is appropriate for ITER operation of advanced sce- narios. Real-time optimisation schemes are discussed and it is concluded that the necessary conditions of optimality tracking approach may potentially be appropriate for ITER operation, thus offering a viable closed-loop optimal control approach. Simulations of ITER advanced operation are necessary in order to assess the present ITER design and uncover the main difficulties that may be encountered during advanced operation. The DINA-CH&CRONOS full tokamak simulator is used to simulate the operation of the ITER hybrid and steady-state scenarios. It is concluded that the present ITER design is appropriate for performing a hybrid scenario pulse lasting more than 1000s, with a flat-top plasma current of 12MA, and a fusion gain of Q ∼= 8. Similarly, a steady-state scenario without internal transport barrier, with a flat-top plasma current of 10MA, and with a fusion gain of Q ∼= 5 can be realised using the present ITER design. The sensitivity of the advanced scenarios with respect to transport models and physical assumption is assessed using CRONOS. It is concluded that the hybrid scenario and the steady- state scenario are highly sensitive to the L-H transition timing, to the value of the confinement enhancement factor, to the heating and current drive scenario during ramp-up, and, to a lesser extent, to the density peaking and pedestal pressure