E. Maljaars

Control of the tokamak safety factor profile with time-varying constraints using MPC

A controller is designed for the tokamak safety factor profile that takes real-time-varying operational and physics limits into account. This so-called model predictive controller (MPC) employs a prediction model in order to compute optimal control inputs that satisfy the given limits. The use of linearized models around a reference trajectory results in a quadratic programming problem that can easily be solved online. The performance of the controller is analysed in a set of ITER L-mode scenarios simulated with the non-linear plasma transport code RAPTOR. It is shown that the controller can reduce the tracking error due to an overestimation or underestimation of the modelled transport, while making a trade-off between residual error and amount of controller action. It is also shown that the controller can account for a sudden decrease in the available actuator power, while providing warnings ahead of time about expected violations of operational and physics limits. This controller can be extended and implemented in existing tokamaks in the near future.

Numerical optimization of actuator trajectories for ITER hybrid scenario profile evolution

Optimal actuator trajectories for an ITER hybrid scenario ramp-up are computed using a numerical optimization method. For both L-mode and H-mode scenarios, the time trajectory of plasma current, EC heating and current drive distribution is determined that minimizes a chosen cost function, while satisfying constraints. The cost function is formulated to reflect two desired properties of the plasma q profile at the end of the ramp-up. The first objective is to maximize the ITG turbulence threshold by maximizing the volume-averaged s/q ratio. The second objective is to achieve a stationary q profile by having a flat loop voltage profile. Actuator and physics-derived constraints are included, imposing limits on plasma current, ramp rates, internal inductance and q profile. This numerical method uses the fast control-oriented plasma profile evolution code RAPTOR, which is successfully benchmarked against more complete CRONOS simulations for L-mode and H-mode mode ITER hybrid scenarios. It is shown that the optimized trajectories computed using RAPTOR also result in an improved ramp-up scenario for CRONOS simulations using the same input trajectories. Furthermore, the optimal trajectories are shown to vary depending on the precise timing of the L–H transition.