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

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

Development of a concept and basis for the DEMO diagnostic and control system

An initial concept for the plasma diagnostic and control (D&C) system has been developed as part of European studies towards the development of a demonstration tokamak fusion reactor (DEMO). The main objective is to develop a feasible, integrated concept design of the DEMO D&C system that can provide reliable plasma control and high performance (electricity output) over extended periods of operation. While the fusion power is maximized when operating near to the operational limits of the tokamak, the reliability of operation typically improves when choosing parameters significantly distant from these limits. In addition to these conflicting requirements, the D&C development has to cope with strong adverse effects acting on all in vessel components on DEMO (harsh neutron environment, particle fluxes, temperatures, electromagnetic forces, etc.). Moreover, space allocation and plasma access are constrained by the needs for first wall integrity and optimization of tritium breeding. Taking into account these boundary conditions, the main DEMO plasma control issues have been formulated, and a list of diagnostic systems and channels needed for plasma control has been developed, which were selected for their robustness and the required coverage of control issues. For a validation and refinement of this concept, simulation tools are being refined and applied for equilibrium, kinetic and mode control studies

Neural Networks for Modeling and Control of Particle Accelerators

We describe some of the challenges of particle accelerator control, highlight recent advances in neural network techniques, discuss some promising avenues for incorporating neural networks into particle accelerator control systems, and describe a neural network-based control system that is being developed for resonance control of an RF electron gun at the Fermilab Accelerator Science and Technology (FAST) facility, including initial experimental results from a benchmark controller.Comment: 21 p

Plasma position control in the STOR-M tokamak : a fuzzy logic approach

Adequate control of the position of the plasma column within the STOR-M tokamak is a chief requirement in order for experimental quality discharges to be obtained. Optimal control over tokamak discharge parameters, including the plasma position, is very difficult to achieve. This is due in large part to the difficulty in modelling the tokamak discharge parameters, as they are highly nonlinear and time varying in nature. The difficulty of modelling the tokamak discharge parameters suggests that a control system, such as a fuzzy logic based controller, which does not require a system model may be well suited to the control of fusion plasma. In order to improve the quality of control over the plasma position within the STOR-M tokamak, the existing analog PID controller was modified. These modifications facilitate the application of a digital controller by a personal computer via the Advantech PCL-711B data acquisition card. The performance of the modified plasma position controller and an Arbitrary Signal Generator developed by the author was evaluated. This modified plasma position controller was applied successfully to the STOR-M tokamak during both normal mode and A.C. mode operation. In both cases, the modified controller provided adequate control over the position of the plasma column within the discharge chamber. Furthermore, the modified controller was more convenient to optimize than the original, existing analog PID controller. By taking advantage of the modifications that were made to the plasma position controller, a fuzzy logic controller was developed by the author. The fuzzy logic based plasma position controller was also successfully applied to the STOR-M tokamak during both normal mode and A.C. operation. The fuzzy controller was demonstrated to reliably provide a higher degree of control over the position of the plasma column within the STOR-M tokamak than the modified PID controller

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

A Strict Control Lyapunov Function for a Diffusion Equation with Time-Varying Distributed Coefficients

International audienceIn this paper, a strict Lyapunov function is developed in order to show the exponential stability and input-to-state stability (ISS) properties of a diffusion equation for nonhomogeneous media. Such media can involve rapidly time-varying distributed diffusivity coefficients. Based on this Lyapunov function, a control law is derived to preserve the ISS properties of the system and improve its performance. A robustness analysis with respect to disturbances and estimation errors in the distributed parameters is performed on the system, precisely showing the impact of the controller on the rate of convergence and ISS gains. This is important in light of a possible implementation of the control since, in most cases, diffusion coefficient estimates involve a high degree of uncertainty. An application to the safety factor profile control for the Tore Supra tokamak illustrates and motivates the theoretical results. A constrained control law (incorporating nonlinear shape constraints in the actuation profiles) is designed to behave as closely as possible to the unconstrained version, albeit with the equivalent of a variable gain. Finally, the proposed control laws are tested under simulation, first in the nominal case and then using a model of Tore Supra dynamics, where they show adequate performance and robustness with respect to disturbances

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

Plasma boundary reconstruction and shape control of Tokamak discharges in RFX-mod

La tesi si occupa del controllo di forma del plasma in configurazione Tokamak a doppio nullo nell’esperimento RFX-mod. Vengono affrontate problematiche di ricostruzione della forma, con un apporto originale nell’utilizzo di una interpolante thin plate spline per la funzione flusso magnetico ricostruita. Per mezzo di un simulatore agli elementi finiti (Maxfea), si è studiata la configurazione di equilibrio e si è derivato un modello lineare con un metodo alle perturbazioni di flusso. Infine, si è sintetizzato il controllore con tecnica LQG. Gli algoritmi di ricostruzione e di controllo sono orientati a soddisfare requisiti di tempo reale e a rispettare i limiti fisici della macchina, in modo da rendere in prospettiva possibile l’effettiva implementazion