Next Article in Journal
Numerical Study of the Aerodynamic Loads on Offshore Wind Turbines under Typhoon with Full Wind Direction
Next Article in Special Issue
TBM/MTM for HTS-FNSF: An Innovative Testing Strategy to Qualify/Validate Fusion Technologies for U.S. DEMO
Previous Article in Journal
Environment-Friendly Heterogeneous Alkaline-Based Mixed Metal Oxide Catalysts for Biodiesel Production
Previous Article in Special Issue
A Review of Dangerous Dust in Fusion Reactors: from Its Creation to Its Resuspension in Case of LOCA and LOVA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real Time Hybrid Model Predictive Control for the Current Profile of the Tokamak à Configuration Variable (TCV)

1
Faculty of Engineering, University of the Basque Country (UPV/EHU), Paseo Rafael Moreno 3, Bilbao 48013, Spain
2
Centre de Recherches en Physique des Plasmas, École Polytechnique Fédérale de Lausanne (CRPP-EPFL), CH-1015 Lausanne, Switzerland
*
Author to whom correspondence should be addressed.
Energies 2016, 9(8), 609; https://doi.org/10.3390/en9080609
Submission received: 27 February 2016 / Revised: 26 July 2016 / Accepted: 27 July 2016 / Published: 3 August 2016
(This article belongs to the Special Issue Fusion Power)

Abstract

:
Plasma stability is one of the obstacles in the path to the successful operation of fusion devices. Numerical control-oriented codes as it is the case of the widely accepted RZIp may be used within Tokamak simulations. The novelty of this article relies in the hierarchical development of a dynamic control loop. It is based on a current profile Model Predictive Control (MPC) algorithm within a multiloop structure, where a MPC is developed at each step so as to improve the Proportional Integral Derivative (PID) global scheme. The inner control loop is composed of a PID-based controller that acts over the Multiple Input Multiple Output (MIMO) system resulting from the RZIp plasma model of the Tokamak à Configuration Variable (TCV). The coefficients of this PID controller are initially tuned using an eigenmode reduction over the passive structure model. The control action corresponding to the state of interest is then optimized in the outer MPC loop. For the sake of comparison, both the traditionally used PID global controller as well as the multiloop enhanced MPC are applied to the same TCV shot. The results show that the proposed control algorithm presents a superior performance over the conventional PID algorithm in terms of convergence. Furthermore, this enhanced MPC algorithm contributes to extend the discharge length and to overcome the limited power availability restrictions that hinder the performance of advanced tokamaks.

Graphical Abstract

1. Introduction

Nuclear fusion has the potential capability to satisfy the world population’s increasing demand for energy. As is known, fusion is the process by which the nuclei of two light atoms are fused together to form a heavier nucleus, with energy produced as a by-product. In this sense, there exists a sufficient fuel supply for several thousand years since the necessary hydrogen isotopes can be generated from water, and lithium is available during the reaction cycle. Moreover, fusion would produce no air pollution or greenhouse gases during normal operation since the fusion reaction product is helium, and it constitutes an intrinsically safe process since a nuclear meltdown with a large, uncontrolled release of energy cannot occur. In spite of significant advances over the last twenty years, nuclear fusion suffers from plasma instabilities that provoke reaction decay. Therefore, the social and economic benefits of being able to achieve controlled fusion are invaluable.
A consortium of governments that includes the European Union, China, Japan, South Korea, Russia, and the United States has promoted the construction of the ITER tokamak (originally an acronym for International Thermonuclear Experimental Reactor), which will generate five to ten times as much energy as it consumes.
Tokamaks are toroidal fusion reactors that confine the hot plasma using a helical magnetic field. They are inherently pulsed devices because a limited power-supply induces a current that heats the plasma up to its operational temperature. Therefore, the aim is to maintain the reaction by self-heating from internal fusion reactions and minimal non-inductive current drive. Since tokamak plasmas are high-order, nonlinear distributed systems with a large number of instabilities, there are many challenging control problems whose solutions are required for ITER. In this sense, plasma control is one of the most basic problems in tokamaks for improving the plasma performance and the pulse duration by avoiding the premature reaction decay. Once the plasma has been created using the primary, electromagnetic fields generated by a set of poloidal field (PF) coils distributed around the vacuum vessel are used to control the plasma current, position and shape.
TCV is one of the existing tokamaks around the world that is being used as a test bench for developing control techniques in order to address the control problems that must be solved before fusion power production becomes viable. A model for plasma current, position and shape based on the RZIp model, and PID control schemes are used in the TCV tokamak. However, the control capabilities on some crucial parameters of the system, as is the case of the plasma current, could be improved. This set up provides the basis for implementing an additional control over the plasma current that, on the one hand should be fast enough to perform real time control of plasma instabilities which require a time response of milliseconds, and on the other hand minimizes the control action in order to adequately deal with the limited power available from actuators. Thus, a natural improvement of the existing magnetic control systems is to optimize these current drive actuators. Moreover, predicting the state changes that will be caused by the control actions, will allow for a smooth shape transition minimizing plasma instabilities. These are the reasons why Model Predictive Control (MPC) has been considered as an excellent control scheme candidate for tokamak operation due to its ability to optimize the future behavior while intrinsically minimizing the control action [1,2,3,4,5].
TCV is a tokamak that operates with a modular and scalable distributed digital control system able to study the hot-ion, high-β ITER-like plasmas (see Figure 1). Even when TCV plasmas would have Te/Ti in the range of 2–5, a new Neutral Beam Injection heating system (NBI 1 MW) [6], combined with an increment in the electron cyclotron plasma heating power at the third harmonic (2 MW) [7] will bring this ratio close to ITER and reactor relevant domain, Te/Ti~1. Therefore, TCV will bring the plasma at the β limit in H-mode (βN~2.8), which is of great interest for ITER and DEMO (Demonstration Power Station that is intended to build upon the ITER) [8], will provide variations in the momentum input to the plasma, and will generate fast ions that can be used to study wave-particle interaction phenomena of interest for burning plasmas [9]. The main parameters, are major radius 0.88 m, horizontal minor radius 0.25 m, toroidal field at major radius 1.5 T, and plasma current 1MA with elongation up to 2.8 and high triangularity from −0.7 to 1. This uniquely expansive shaping capability extends to limited and diverted plasmas, the latter featuring open divertor geometry and including single- and double-null as well as snowflake topologies [10]. Due to the flexible configuration and control, TCV is a tokamak that is especially active in control research, and all of this invaluable experience is expected to provide important input for safe and effective operation of ITER [11,12,13,14,15,16].
TCV has achieved stable plasma current and position control. In a typical plasma discharge, the plasma is feedback controlled using a sequence of PIDs with a multiloop structure [17,18]. PIDs are partly designed using decoupled models of the system to be controlled, and later, fine-tuned during plant operation. Decoupling of the controllers has often meant that the effect on a particular plasma parameter from coils other than the actuator has been neglected, (see Figure 2). Recent efforts seek to design more complex and robust controllers focusing the attention on these effects [19,20,21,22]. In this sense, MPC has been favored in this article as the controller choice due to the well-known advantages that this kind of controllers present and in particular because it is able to deal in an optimal sense with opposite behaviors in strongly coupled variables.
Common problems within real time optimization for distributed model predictive control when executing a control action with objectives are stability, constraint satisfaction, performance guarantee and real time execution guarantee. The key challenge is to obtain a fast, fixed time optimization. Interior point methods resolve within milliseconds, fast gradient methods in microseconds, while explicit methods solve in nanoseconds. Nowadays, CPLEX [23,24], is among the most popular optimization engines, although there are other very popular ones such as SeDuMi [25] or Multi-Parametric Toolbox [26].
The current interior point methods for solving real time MPC consider a warm start, exploiting the internal structure in MPC problems and terminating the online optimization early. Therefore, the solution obtained during online optimization steps is suboptimal, so that it may be infeasible, it can destabilize the system and it may cause steady-state offset. This is to say, the requirements to solve the system within a time window restricts the algorithm to maximum on-line optimization steps which, in turn, may destroy the stability properties of optimal MPC. Wang et al. developed a custom solver that yields a fast solution of the tracking problem with guaranteed stability for all suboptimal iterates for all time constraints [27]. They modified the matrix structure in the Newton step computation using the tracking formulation and Lyapunov constraint.
Fast gradient methods will require a number of iterations that are strongly dependent on the choice of initial iterate. Therefore, a warm starting is favored [28].
Explicit model predictive control uses offline computations to determine all operating regions in which the optimal control moves are determined by evaluating a linear function. Thus, these controllers require a maximum computational time per control interval equal to that needed to test each region, computing the violation term in each case, and then calculating the suboptimal control adjustment and are therefore useful for applications that require small sample times. However, explicit MPC may not satisfy given real time constraints and it is exponentially sensitive to model dimension and changes in system dynamics.
In the centralized MPC a single controller decides all actions to take. On the other hand, a distributed MPC approach facilitates the coordination of highly complex or independent system plants (see [5]). Reviews such as [29,30] present different compatible ways to label distributed control architectures. In this context, the main novelty of the paper relies on the application of the hybrid MPC control scheme to the particular case of the TCV.
The control method presented in this article is based on a multiloop concept where the traditional PID is combined with a fast MPC that only acts over a reduced system. The intrinsic minimization of the MPC further reduces the current drive source requirements [31]. Therefore, MPC composes an excellent control option since it not only considers the information related to the system at an operation instant as in other advanced controllers [32,33], but also takes into account the predicted future information so as to obtain an optimal control action. Furthermore, the Ohmic Heating (OH) coil may reduce the power requirements of the non-inductive current drive sources required for current profile control due to its ability to transiently manipulate the low order moments of the current profile. Until recently, this line of research [31,32,33] had been carried out with plant models of the type described in [34], which are two orders of magnitude smaller than the RZIp model. These models are excellent for real time control algorithms although the underlying plasma physics are reflected in a limited manner. The current choice of a more complex model such us RZIP is motivated from the need to faithfully represent the plasma dynamics. However, model size is paramount in real-time control on these systems, since it implies that there is a time window to gather data, process that data, and update the system. This is quite a strong limitation since the stability of the system degrades when the time window is missed. To achieve this restriction, all the MPC developments within this article will attempt to reduce the computational time as much as possible. Thus, only strictly needed computations will be performed and analytical expressions will be prioritized over numerical solvers. The main novelty of the paper relies on the application of novel hybrid MPC control scheme to the particular case of the TCV so that the performance improvements that this controller affords are of special relevance to ITER and DEMO, since these Advanced Tokamaks have to rely on having large bootstrap current fractions and/or a pulsed operation [35,36,37] with limited power available for non-inductive current drive actuators.
This article is organized as follows: The RZIp model for TCV plasma control is presented in Section 2. The plasma enhanced MPC for the TCV model and its block diagram are discussed in Section 3. Next, in Section 4 experimental results for the implementation of the novel enhanced controller that was developed in the previous section, are presented and compared with those obtained from a traditional PID control scheme. Comparisons of experimental data with simulations generated using different algorithms are also shown. Finally, Section 5 will present some concluding remarks to end the article.

2. RZIP Model for the TCV Tokamak

RZIp considers a rigid plasma displacement model, where the plasma is assumed to be a unit that is allowed to move either vertically or radially, and the vessel eddy currents will use an eigenmode representation. The model is derived from the plasma and circuit equations, perturbing the vertical and radial force balance equations. It was first developed by Lister and Sharma [27,28,38] with the aim of extending the stability region of the resulting closed-loop plasma dynamical system via the implementation of adequate control techniques.

2.1. RZIp System Formulation

The model is represented in state space form where state variables are the plasma current, the structure currents and the position variation of the plasma current:
q ˙ = [ Ι e Ι s r ˙ ] ,   where   r = [ R z ]
The entries of the input vector are the effective voltages applied to each element and a vector of externally applied poloidal field coil voltages:
U = [ V e V s 0 ]
The model addresses flux and energy conservation issues that may be extended into four vector equations, the Kirchoff voltage law for the plasma elements, the Kirchoff voltage law for the structural and poloidal circuits and the force balance in the R and z directions, respectively
d ( L e I e + M e s I s + R E I e ) d t + Ω e I e = V e
d ( L s I s + M s e I e ) d t + Ω s I s = V s
d ( m e R ˙ ) d t = 1 2 I e L e R I e + I s M s e R I e + E I e 2 2
d ( m e z ˙ ) d t = 1 2 I e L e z I e + I s M s e z I e

2.2. Lumping, Linearization and Eigenmode Reduction

A lumped model may be derived from Equations (3)–(6) by defining averaged plasma quantities. The total plasma current will be denoted by Ip, while the density distribution j(R,z) is calculated from the Grad-Shafranov equation by an inverse equilibrium reconstruction code. The Grad-Shafranov equilibrium equation describes the shape and current profile of a toroidally symmetric plasma in an externally applied magnetic field. The equilibrium equations are found by considering the plasma almost charge neutral so that the condition for balance is that the thermodynamic pressure and magnetic force satisfy:
j × B = p
subject to toroidal symmetry, where j is the current density, B is the magnetic field and p is the pressure. Thus, the magnetic surfaces are surfaces of constant pressure, and the lines of current lie on the magnetic surfaces. For computational simplicity the variables (R,z) are replaced with ( R I p 0 , z I p 0 ) , where I p 0 denotes the constant equilibrium plasma current. The system is solved in incremental form, so that the state vector is rewritten as:
x = [ I s I s 0 ( z z 0 ) I p 0 ( R R 0 ) I p 0 I p I p 0 ] = [ δ I s δ z I p 0 δ R I p 0 δ I p ]
The resulting equations may be linearized around the chosen operation point. The plasma circuit equation, structure circuit equation, and both vertical and radial force balance equations yield a state space model whose values may be found in Supplementary Materials:
x ˙ = A x + B u y = C x + D u
where A = M–1·R, M being the mutual inductance matrix and R a diagonal matrix of structure resistances. However, since the passive structure model is very large, an eigenmode reduction is performed by considering only those associated with the largest eigenvalues. In particular, the passive structure will be expressed in terms of the current eigenvalues instead of the passive structure currents themselves. A very positive side effect of this procedure is the increment of the stability of the system by removing modes closer to the imaginary axis.
Some relevant values of the RZIp-based TCV plant model that have been used are shown in Table 1. In particular, values for the actuators Poloidal Field coils may be found, where the inner ones are responsible for the fastest control action. The control action of these coils is explicitly stated in the model through Equations (3) and (4). Besides, an Ohmic Heating current peak value and its variation are also given since these both are relevant actuators for the case of the plasma current. Their effect is reflected in the internal energy and inductances of Equations (5) and (6). Finally, the plasma current average values are given since this is the output variable which value has been improved with the MPC controller.

3. Enhanced MPC

The proposed hybrid control scheme uses combined proportional integral derivative (PID) control of the full state space system coupled with a predictive control model for the plasma current. The excellent performance of the scheme relies on the matrices that decouple the PID control scheme. Decoupling of the control involves a process of modifying the dynamic system matrix of the system so as to reduce some of the off-diagonal elements to zero. For the axisymmetric system (9) characterized by the state vector (8) decoupling can also be regarded as a process of equalizing the eigenmodes of inductances and Green functions (between the structure and the plasma). The corresponding plasma current decoupling matrix can be viewed as the best fit port voltages under the condition that produces a zero field when a specific mode is excited. Thus, PIDs are partly designed using decoupled models of the system to be controlled, and later, fine-tuned during plant operation. The PID scheme gain values are provided jointly with the data of the TCV pulse 49,626 in Supplementary Materials. Decoupling of the controllers means that the effect on a particular plasma parameter from coils other than the actuator has been neglected. The aim is now to combine these system couplings by adding an MPC over a reduced system in a multiloop structure.
Advanced Tokamak control suffers from high power requirements, and, most importantly, in future tokamaks both the total and peak actuator power is strongly limited. In this context, MPC provides a systematic method of dealing with constraints on inputs and states. In MPC these limitations are explicitly accounted for by solving a constrained optimization problem in real-time to determine the optimal predicted control. In addition to this, nonlinear plant dynamics can be similarly incorporated in the prediction mode.

3.1. Model Predictive Plasma Current Control

This section is concerned mainly with the case of plasma current model predictive control for the discrete-time linear system. The real system has a time constant in the order of milliseconds so that the discretization time step has to be of that order of magnitude to obtain an acceptable difference between the continuous and the discretized optimal control profiles. However, such a small discretization time step creates a need to significantly reduce the computational time at each optimization step. The formulations discussed in this paper adopt a discretization time step of 50 msec, which is found to be sufficiently small to approximate the continuous optimal control profile. For this purpose, the plant dynamics are linearized at each time step with the last computed state space variable values other than plasma current. The resulting optimization problem is solved with an analytical solution equation approach using an exact Hessian approximation. Thus, the future response of the controlled tokamak plasma current is predicted using a reduced dynamic model:
x d ( k + 1 ) = A d x d ( k ) + B d u d ( k ) y ( k ) = C d x d ( k )  
where y ( k ) is the plasma current, while x d ( k ) and u d ( k ) are the states and input vectors at the kth sampling instant. In particular, the discrete matrix C d is a reduced matrix from the original system output matrix that corresponds to the row vector for the plasma current.
There is a tradeoff between the plant model size and the accuracy of the solution. In this sense, on the one side, most of the transport models are not useful for real-time control due to the fact that they are computationally expensive. On the other side, the reduced state-space models do not fully reflect the plant behavior. In the particular case of the RZIp model considered it would not be feasible to implement a real-time MPC involving the complete set of variables, so that it is necessary to adequately choose the variables to be considered.
Given a predicted input sequence, the corresponding sequence of plasma current prediction values is generated by simulating the model forward over the prediction horizon. For notational convenience, these predicted sequences are often stacked into vectors that may be denoted by:
x ( k ) = [ x ( k + 1 | k ) x ( k + 2 | k ) x ( k + N p | k ) ] ,  u ( k ) = [ u ( k | k ) u ( k + 1 | k ) u ( k + N c 1 | k ) ]
where x ( k + i | k ) and u ( k + i | k ) denote plasma current value and ohmic heating voltage input at time k + i that are predicted at time k . Therefore, x d ( k + i | k ) evolves according to the prediction model:
x d ( k + i + 1 | k ) = A d x d ( k + i | k ) + B d u d ( k + i | k )
with an initial condition at the beginning of the prediction horizon defined x d ( k | k ) = x d ( k ) k = 0 , 1 , and a control horizon, N c , whose value is smaller than or equal to that of the prediction horizon, N p .
The predictive control feedback law is computed by minimizing a predicted performance cost, which is defined in terms of the predicted sequences x ( k ) and u ( k ) . In particular, rewriting the output current in compact form as:
y ^ = Fx ( k ) + G u ^ ( k )
u ^ ( k ) denotes the ohmic heating voltage for a given reference signal r ( k ) at a given time. Within a prediction horizon the control system aims to bring the plasma current as close as possible to the reference:
r ^ = [ r ( k ) r ( k + N p 1 ) ] T
The corresponding quadratic cost function:
J = ( r ^ y ^ ) T ( r ^ y ^ ) + u ^ T R u ^
is clearly a function of the predicted OH voltages and a positive definite weighing matrix R . The optimal voltage sequence for the problem of minimizing this cost function takes into account the input restrictions as equivalent constraints on u ^ in the optimization.
Only the first element of the N c optimal predicted voltage sequence is input back into the plant. This procedure, known as the receding horizon strategy, provides a degree of robustness to modeling errors and uncertainty since the input at each step takes into account the current state values to compute future events. Additionally, by continually shifting the horizon to where future inputs are optimized, it attempts to compensate for the fact that this horizon is finite. Therefore, a receding horizon strategy aims to ensure that the performance of the closed-loop system is at least as good as that of the optimal prediction.

3.2. Control System Multiloop Structure

This novel hybrid control system for the TCV has been developed with the aim of achieving closed-loop performance objectives that cannot be accomplished using a PID or a MPC on its own.
A block diagram with a brief summary of the proposed scheme may be seen in Figure 3. The error between all the references and the estimated values serves as an input to the block that contains the PID control. The output is composed by a set of commands for different coils that are given in Dall. In parallel, the reference error for the plasma current is the input to the MPC. This block performs the MPC of the plant computed as shown in Equation (10). The computational load within this block is reduced as much as possible, since the controller involves the analytical solution of the minimization for the functional given in Equation (15). The corresponding control command output is provided in DIp. The block denoted by “Fast Control” is a third block that commands an independent set of coils located inside the vessel to perform vertical stabilization, which requires a very fast control response of about 5 μs. These three blocks determine a control action over the different coils. Moreover, the block “Substitute FPS” represents the data transmission from the Fast Power Supply (FPS) channels to the internal coil system. Finally, the estimated values obtained from the diagnostics close the control loop.
Considering different time–scales, makes it possible to compete with homogeneous control schemes. On the other hand, the benefits derived from the model predictive control horizon, such as stability and robustness, do automatically apply. Additionally, the hybrid scheme allows for computational cost efficiency by restricting the MPC usual overhead to the plasma current computation.
The PID and receding horizon control algorithms are implemented in two separate loops which run in parallel. The MPC is implemented over a reduced plant model and the overall execution time remains unchanged with respect to the implementation of the PID without MPC. Assuming that both loops start at time tk, the Pseudo code is as follows:
  • PID loop
    1:
    Solve PID for j = k and obtain vk, where two of the entries will be corrected with uk.
    2:
    Wait until uk.
    3:
    Update the entries of vk.
    4:
    Apply inputs to the model.
    5:
    k = k + 1.
    6:
    post process.
    7:
    Jump to 1.
  • MPC loop
    1:
    Wait until k + 1.
    2:
    Reduce the plant model to compute Ip.
    3:
    Solve (15) for j = k; … ; k + Np and obtain uk; … ; uk+Nc.
    4:
    Apply uk to the inputs.
    5:
    Jump to 1.

4. Numerical vs. Experimental Results

A study of the MPC control action on a real shot for the already described TCV is performed. As previously mentioned, a correct modeling of the tokamak has a large influence over the control cost-effectiveness because the minimization of the cost function is based on the system state space matrices. Therefore, it can be emphasized that obtaining an accurate and computationally efficient system model is one of the main challenges in order to successfully implement a MPC. In this article the control oriented RZIp model has been considered, because it is widely recognized as a good control oriented model by the fusion research community.
Once the prediction and control horizons have been defined in the algorithm with N p = 5 and N c = 3 , the TCV system response is obtained, (see left hand side of Figure 4). It can be seen how the plasma current follows the given reference. The values of the parameters of the MPC have been set in order to get an adequate adjustment with the set point. The results for the TCV pulse 49,626 are shown in the right hand side of Figure 4, where a fine-tuned PID controller was used, instead of the hybrid MPC.
It can be clearly seen from these two figures that using an enhanced MPC (red line in the left hand side figure), it gives a better fit than the PID without the enhancement (red line in the right hand side figure) to the real measure (blue line). This effect is associated with the fact that the MPC has the plant model embedded, which enables the controlled TCV to yield a better fit of the tracking variable. Other plant model outputs are not shown because there are no significant differences between them when using either a PID or the PID enhanced with the MPC, as should be expected since the plasma current was the variable to improve.

Prediction Horizon and Control Horizon

In this section the effect of different horizons on the TCV pulse 49626 will be studied. There are two kinds of horizons: the prediction horizon N p , and the control horizon, N c . The value N p denotes the number of time-steps ahead in the horizon where the expanded system solution will predict the change in the plasma current that will be caused by changes in the independent variables. These chosen dependent variables are determined by the value of N c , that denotes when the control law stops trying, and which might be smaller than the prediction horizon. High values of horizons are supposed to provide a better system response while small ones may present difficulties for achieving the control objectives; also, both horizon values have to be balanced. Furthermore, there is a tradeoff between the computational effort and the desired accuracy.
To analyze the effect of different horizons on the controlled horizon, let us recall the satisfactory results achieved for the TCV system by an enhanced MPC with N p = 5 and N c = 3 (see Figure 4).
First at all, the horizon values will be increased up to N p = 9 and N c = 5 . It may be seen that even when the left of Figure 5 shows an improvement of the plasma current response with respect to that displayed in the left of Figure 4, the added control complexity might not fully justify the improved accuracy in the controlled plasma current.
A very small number for both the control horizon and the predictive horizon has also been considered, with only two steps. The result corresponds to a system that has not been adequately controlled. As might be expected, the response of the system when using a short prediction horizon and a short control horizon is not stable. The plasma current does not even reach the set point. The control is unable to anticipate its response quickly enough so that the plasma current slides out of bounds while the TCV system clearly moves out of its stability region. In the other cases, the system stays within the stable region.
Finally, a long predictive horizon together with a small control horizon, shown in the right hand side of Figure 5, provides a poor response. This is due to the propagation of the last calculated manipulated variable value over the rest of the prediction horizon, which becomes too long. It is important to maintain an adequate relation between the two horizons.

5. Conclusions

An enhanced MPC control system has been designed for the TCV that takes advantage of a global PID control algorithm. The RZIp has been chosen to implement the control because it is widely accepted within the fusion research community as a reputed control oriented modeling tool. RZIp considers both active coils and the passive structure, which are strongly coupled in the physical system, so that it requires an eigenmode reduction in the passive structure. This procedure not only reduces the system size but also contributes to system stabilization by removing those modes corresponding to eigenvalues close to the imaginary axis. Although the PID scheme has been designed to work over the whole system, the system reduction gives the opportunity to implement a multiloop control system, such as is the case of the one presented in this article considering an MPC on the plasma current.
This novel enhanced MPC system meets the plasma current control target. Experimental evidence has also been provided, showing that the closed loop system yields improved results with respect to the global PID control. Furthermore, these satisfactory results are obtained for relatively short horizons. It is, therefore, concluded that the enhanced MPC implementation in RZIp for the TCV tokamak provides an excellent complementarity with respect to a more traditional PID control system. The benefits are twofold because it not only provides a smoother fit to the state of interest but also minimizes the control action, saving much needed energy from the non-inductive limited resources.
There are two distinct future lines of research. The obvious way forward is to extend the MPC to other variables, by considering different time control actions. Also, it would be desirable to validate other models based on RZIp simplifications, with the aim to reduce the computational load while still representing the particular plasma dynamics at hand. Finally, it would also be interesting to consider parareal type algorithm solvers to deal with the optimization problem. These kinds of time domain predictor-corrector solvers allow high performance computations so that the prediction part of the MPC algorithm could be implicitly treated within the solver itself, allowing the solution of more complex models.

Supplementary Materials

The following are available online at www.mdpi.com/1996-1073/9/8/609/s1. TCV Model—Hybrid MPC Data.

Acknowledgments

The authors would like to acknowledge the collaboration of the Basque Energy Board (EVE) through Agreement UPV/EHUEVE23/6/2011 and the Spanish National Fusion Laboratory (CIEMAT) UPV/EHUCIEMAT08/190. Authors are also very grateful for the partial financial contribution to the UPV/EHU and EJ/GV through projects GIU14/07 and IT987-16, as well as to the FEDER and MINECO through project DPI2015-70075-R.

Author Contributions

Izaskun Garrido and Aitor Garrido at UPV/EHU conceived and designed the MPC control loops and the experimental work, while Stefano Coda, Hoang B. Le and Jean Marc Moret at EPFL designed de experiments and carried out the plant model identification with the help of the TCV team.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ΒPlasma pressure/magnetic pressure ratio
Te, TiElectron and ion temperatures
βNβ coefficient
Ι e Plasma elements currents
Ι s Structure currents
R Radial current position
z Vertical current position
V e Plasma elements effective voltages
V s Poloidal field coil voltages
L e , M e s Self and mutual inductance matrices
Ω e Plasma elements resistance matrix
R E I e 2 Plasma internal energy
E Plasma internal energy constant matrix
m e Plasma elements mass
j Current density
B Magnetic field
p Thermodynamic pressure
I p 0 Constant equilibrium plasma current
I p Total plasma current
M Mutual inductance matrix
R Diagonal matrix of structure resistances
N p Prediction horizon
N c Control horizon
E1–E8; F1–F8Poloidal field shaping coils
κ Plasma elongation parameter
OHOhmnic heating coils
DOHOH1-OH2
ZIPdslow Slow vertical observer plasma current
ZIPpFast vertical observer plasma current
RZIP(R, z, Ip) plasma model
MPCModel Predictive Control
PIDProportional Integral Derivative global scheme
MIMOMultiple Input Multiple output
TCVTokamak à Configuration Variable
ITERInternational Thermonuclear Experimental Reactor
PFPoloidal field
βRatio of the plasma pressure to the magnetic pressure
NBINeutral beam injection
H-modeHigh confinement mode
DEMODEMOnstration Power Station

References

  1. Xiong, M.; Gao, F.; Liu, K.; Chen, S.; Dong, J. Optimal Real-Time Scheduling for Hybrid Energy Storage Systems and Wind Farms Based on Model Predictive Control. Energies 2015, 8, 8020–8051. [Google Scholar] [CrossRef]
  2. Nguyen, T.-T.; Yoo, H.-J.; Kim, H.-M. Application of Model Predictive Control to BESS for Microgrid Control. Energies 2015, 8, 8798–8813. [Google Scholar] [CrossRef]
  3. Vu, T.-V.; Chen, C.-K.; Hung, C.-W. A Model Predictive Control Approach for Fuel Economy Improvement of a Series Hydraulic Hybrid Vehicle. Energies 2014, 7, 7017–7040. [Google Scholar] [CrossRef]
  4. Moradzadeh, M.; Boel, R.; Vandevelde, L. Anticipating and Coordinating Voltage Control for Interconnected Power Systems. Energies 2014, 7, 7017–7040. [Google Scholar] [CrossRef] [Green Version]
  5. Rowe, M.; Yunusov, T.; Haben, S.; Potter, W.H.B. The Real-Time Optimisation of DNO Owned Storage Devices on the LV Network for Peak Reduction. Energies 2014, 7, 1027–1047. [Google Scholar] [CrossRef]
  6. Karpushov, A.N.; Duval, B.P.; Chavan, R.; Fable, E.; Mayor, J.M.; Sauter, O.; Weisen, H. A scoping study of the application of neutral beam heating on the TCV tokamak. Fusion Eng. Des. 2011, 86, 868–871. [Google Scholar] [CrossRef]
  7. Fasoli, A.; the TCV Team. Overview of physics research on the TCV tokamak. Nucl. Fusion 2009, 49, 104005. [Google Scholar] [CrossRef]
  8. Zohm, H.; Angioni, C.; Fable, E.; Federici, G.; Gantenbein, G.; Hartmann, T.; Lackner, K.; Poli, E.; Porte, L.; Sauter, O. The physics guidelines for a tokamak DEMO. Nucl. Fusion 2013, 53, 073019. [Google Scholar] [CrossRef]
  9. Fasoli, A.; Alberto, S.; Chavan, R.; Duval, B.P.; Karpushov, A.; Martin, Y.; Sauter, O.; Toussaint, M.; Weisen, H.; the TCV Team. An upgraded TCV for tokamak physics in view of ITER and DEMO. In Proceedings of the 40th EPS Conference on Plasma Physics, Espoo, Finland, 1–5 July 2013; Volume 37D.
  10. Coda, S.; the TCV Team. The science program of the TCV tokamak: Exploring fusion reactor and power plant concepts. Nucl. Fusion 2015, 55, 104004. [Google Scholar] [CrossRef]
  11. Le, H.B.; Felic, F.; Paley, J.I.; Duval, B.P.; Moret, J.-M.; Coda, S.; Sauter, O.; Fasel, D.; Marmillod, P.; the TCV team. Distributed digital real-time control system for TCV tokamak. Fusion Eng. Des. 2014, 89, 155–164. [Google Scholar] [CrossRef]
  12. Moret, J.-M.; Duval, B.P.; Le, H.B.; Coda, S.; Felici, F.; Reimerdes, H. Tokamak equilibrium reconstruction code LIUQE and its real time implementation. Fusion Eng. Des. 2015, 91, 1–15. [Google Scholar] [CrossRef]
  13. Paley, J.I.; Felici, F.; Berrino, J.; Coda, S.; Cruz, N.; Duval, B.P.; Goodman, T.P.; Martin, Y.; Moret, J.M.; Piras, F. Real time control of plasmas and ECRH systems on TCV. Nucl. Fusion 2009, 49, 085017. [Google Scholar] [CrossRef]
  14. Kim, D.; Goodman, T.P.; Sauter, O. Real-time sawtooth control and neoclassical tearing mode preemption in ITER. Phys. Plasmas 2014, 21, 061503. [Google Scholar] [CrossRef]
  15. Hommen, G.; De Baar, M.; Duval, B.P.; Andrebe, Y.; Le, H.B.; Klop, M.A.; Doelman, N.J.; Witvoet, G.; Steinbuch, M.; the TCV Team. Real-time optical plasma boundary reconstruction for plasma position control at the TCV Tokamak. Nucl. Fusion 2014, 54, 073018. [Google Scholar] [CrossRef]
  16. Felici, F.; Sauter, O.; Coda, S.; Duval, B.P.; Goodman, T.P.; Moret, J.-M.; Paley, J.I.; the TCV Team. Real-time physics-model-based simulation of current density profile in tokamak plasmas. Nucl. Fusion 2011, 51, 083052. [Google Scholar] [CrossRef]
  17. Sharma, A.S.; Limebeer, D.J.N.; Jaimoukha, I.; Lister, J.B. Modelling and Control of TCV. IEEE Trans. Control Syst. 2005, 13, 356–369. [Google Scholar] [CrossRef] [Green Version]
  18. Lister, J.B.; Sharma, A.; Limebeer, D.J.N.; Nakamura, Y.; Wainwright, J.P.; Yoshino, R. Plasma equilibrium response modeling and validation on JT-60U. Nucl. Fusion 2002, 42. [Google Scholar] [CrossRef]
  19. Sevillano, M.G.; Garrido, I.; Garrido, A.J. ASTRA-matlab integration for the control of tokamaks. Informatica 2011, 22, 241–257. [Google Scholar]
  20. García, C.E.; Prett, D.M.; Morari, M. Model predictive control: Theory and practice—A survey. Automatica 1989, 25, 335–348. [Google Scholar] [CrossRef]
  21. Sevillano, M.G.; Garrido, I.; Garrido, A.J. Sliding-mode loop voltage control using ASTRA-Matlab integration in Tokamak reactors. Int. J. Innov. Comput. Inf. Control 2012, 8, 6473–6489. [Google Scholar]
  22. Sevillano, M.G.; Garrido, I.; Garrido, A.J. Control-oriented automatic system for transport analysis (ASTRA)-matlab integration for Tokamaks. Energy 2011, 36, 2812–2819. [Google Scholar] [CrossRef]
  23. Carrascal, E.; Garrido, I.; Garrido, A.J.; Sala, J.M. Optimization of the Heating System Use in Aged Public Buildings via Model Predictive Control. Energies 2016, 9, 251. [Google Scholar] [CrossRef]
  24. CPLEX Optimization Software Package. IBM. Available online: http://www-01.ibm.com/software/commerce/optimization/cplex-optimizer/ (accessed on 6 April 2016).
  25. Sturm, J.F. Using SeDuMi 1.02, a MATLAB Toolbox for Optimization over Symmetric Cones. Optim. Methods Softw. 1999, 11, 625–653. [Google Scholar] [CrossRef]
  26. Kvasnica, M.; Grieder, P.; Baotic, M.; Morari, M. Multi-Parametric Toolbox (MPT). In Hybrid Systems: Computation and Control; Springer: Berlin/Heidelberg, Geramny, 2004; pp. 448–462. [Google Scholar]
  27. Wang, Y.; Boyd, S. Fast model predictive control using online optimization. IEEE Trans. Control Syst. Technol. 2010, 18, 267–278. [Google Scholar] [CrossRef]
  28. Richter, S.; Jones, C.N.; Morari, M. Real-time input-constrained MPC using fast gradient methods. In Proceedings of the 48th IEEE Conference on Decision and Control, 2009 Held Jointly with the 28th Chinese Control Conference, CDC/CCC 2009, Shanghai, China, 15–18 December 2009.
  29. Christofides, P.D.; Scattolini, R.; Muñoz, D.; Liue, J. Distributed model predictive control: A tutorial review and future research directions. Comput. Chem. Eng. 2013, 51, 21–41. [Google Scholar] [CrossRef]
  30. Negenborn, R.R.; Maestre, J.M. Distributed Model Predictive Control: An overview and roadmap of future research opportunities. IEEE Control Syst. 2014, 34, 87–97. [Google Scholar] [CrossRef]
  31. Garrido, I.; Garrido, A.J.; Romero, J.; Carrascal, E.; Sevillano, M.G. Low Effort Li Nuclear Fusion Plasma Control Using Model Predictive Control Laws. Math. Prob. Eng. 2015, 2015. [Google Scholar] [CrossRef]
  32. Romero, J.A.; Coda, S.; Felici, F.; Moret, J.-M.; Paley, J.; Sevillano, G.; Garrido, I.; Le, H.B. Sliding mode control of a tokamak transformer. In Proceedings of the 2012 IEEE 51th Conference on Decision and Control (CDC), Maui, HI, USA, 10–13 December 2012.
  33. Garrido, I.; Garrido, A.J.; Sevillano, M.G.; Romero, J.A. Robust sliding mode control for tokamaks. Math. Probl. Eng. 2012, 2012. [Google Scholar] [CrossRef]
  34. Romero, J.A.; Moret, J.-M.; Coda, S.; Felici, F.; Garrido, I. Development and validation of a tokamak skin effect transformer model. Nucl Fusion 2012, 52, 023019. [Google Scholar] [CrossRef]
  35. Luce, T.C. Realizing steady-state tokamak operation for fusion energy. Phys. Plasmas 2011, 18, 030501. [Google Scholar] [CrossRef]
  36. Cunningham, G. High Performance Plasma Vertical Position Control System for Upgraded MAST. Fusion Eng. Des. 2013, 88, 3238–3247. [Google Scholar] [CrossRef]
  37. Yuan, Q.P.; Xiao, B.J.; Luo, Z.P.; Walker, M.L.; Welander, A.S.; Hyatt, A.; Qian, J.P.; Zhang, R.R.; Humphreys, D.A.; Leuer, J.A.; et al. Plasma current, position and shape feedback control on EAST. Nucl. Fusion 2013, 53. [Google Scholar] [CrossRef]
  38. Sharma, A.S. Tokamak Modelling and Control. Ph.D. Thesis, Imperial College of Science, Technology and Medicine, London, UK, 2002. [Google Scholar]
Figure 1. Partial view of the installations of Tokamak à Configuration Variable, TCV, at CRPP-EPFL.
Figure 1. Partial view of the installations of Tokamak à Configuration Variable, TCV, at CRPP-EPFL.
Energies 09 00609 g001
Figure 2. A poloidal cross-section of TCV showing the single (A and D) and multi (B and C) turn loops used for diamagnetic measurements. The ohmic transformer coils A1, B1, B2, C1, C2, D1 and D2 and the poloidal field shaping coils E1–E8 and F1–F8.
Figure 2. A poloidal cross-section of TCV showing the single (A and D) and multi (B and C) turn loops used for diamagnetic measurements. The ohmic transformer coils A1, B1, B2, C1, C2, D1 and D2 and the poloidal field shaping coils E1–E8 and F1–F8.
Energies 09 00609 g002
Figure 3. Block diagram of the enhanced MPC for the TCV.
Figure 3. Block diagram of the enhanced MPC for the TCV.
Energies 09 00609 g003
Figure 4. Plasma current: RZIp measured values (PID), blue (-) lines. Hybrid MPC RZIp computed values with N p = 5 and N c = 3 , red (- -) (A); PID RZIp computed values, red (- -) (B).
Figure 4. Plasma current: RZIp measured values (PID), blue (-) lines. Hybrid MPC RZIp computed values with N p = 5 and N c = 3 , red (- -) (A); PID RZIp computed values, red (- -) (B).
Energies 09 00609 g004
Figure 5. Plasma current: RZIp measured values (PID), blue (-) lines. Hybrid MPC RZIp computed values with N p = 9 and N c = 5 , red (- -) (A) and using a hybrid MPC with N p = 9 and N c = 2 (B).
Figure 5. Plasma current: RZIp measured values (PID), blue (-) lines. Hybrid MPC RZIp computed values with N p = 9 and N c = 5 , red (- -) (A) and using a hybrid MPC with N p = 9 and N c = 2 (B).
Energies 09 00609 g005
Table 1. TCV relevant values.
Table 1. TCV relevant values.
Inner PF Coils (V)Outer PF Coils (V)Plasma Current & OH Current (A)
E2−715F2c−4914OH1.7 × 104
E4−3817F3c−1930DOH87 (max. 150)
E5−1443F5c−1662IP2.5 × 105
E7−204F6c−316ZIPp0.17 × 104
E8−1571F8c−578ZIPslow4.5 × 104

Share and Cite

MDPI and ACS Style

Garrido, I.; Garrido, A.J.; Coda, S.; Le, H.B.; Moret, J.M. Real Time Hybrid Model Predictive Control for the Current Profile of the Tokamak à Configuration Variable (TCV). Energies 2016, 9, 609. https://doi.org/10.3390/en9080609

AMA Style

Garrido I, Garrido AJ, Coda S, Le HB, Moret JM. Real Time Hybrid Model Predictive Control for the Current Profile of the Tokamak à Configuration Variable (TCV). Energies. 2016; 9(8):609. https://doi.org/10.3390/en9080609

Chicago/Turabian Style

Garrido, Izaskun, Aitor J. Garrido, Stefano Coda, Hoang B. Le, and Jean Marc Moret. 2016. "Real Time Hybrid Model Predictive Control for the Current Profile of the Tokamak à Configuration Variable (TCV)" Energies 9, no. 8: 609. https://doi.org/10.3390/en9080609

APA Style

Garrido, I., Garrido, A. J., Coda, S., Le, H. B., & Moret, J. M. (2016). Real Time Hybrid Model Predictive Control for the Current Profile of the Tokamak à Configuration Variable (TCV). Energies, 9(8), 609. https://doi.org/10.3390/en9080609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop