A control-oriented model of the current profile in tokamak plasma

This paper proposes a control-oriented approach to the tokamak plasma current profile dynamics. It is established based on a consistent set of simplified relationships, in particular for the microwave current drive sources, rather than exact physical modelling. Assuming that a proper model for advanced control schemes can be established using the so-called cylindrical approximation and neglecting the diamagnetic effects, we propose a model that focuses on the flux diffusion (from which the current profile is inferred). Its inputs are some real-time measurements available on modern tokamaks and the effects of some major actuators, such as the magnetic coils, lower hybrid (LHCD), electron and ion cyclotron frequency (ECCD and ICRH) systems, are particularly taken into account. More precisely, the non-inductive current profile sources are modelled as 3-parameters functions of the control inputs derived either from approximate theoretical formulae for the ECCD and bootstrap terms or from experimental scaling laws specifically developed from hard x-ray Tore Supra data for the LHCD influence. The use of scaling laws in this model reflects the fact that the operation of future reactors will certainly depend upon a great number of scaling laws and specific engineering parameters. The discretization issues are also specifically addressed, to ensure robustness with respect to discretisation errors and the efficiency (in terms of computation time) of the associated algorithm. This model is compared with experimental results and the CRONOS solver for tore supra tokamak