CRONOS is a suite of numerical codes for the predictive/interpretative simulation of a full tokamak discharge. It integrates, in a modular structure, a 1D transport solver with general 2D magnetic equilibria, several heat, particle and impurities transport models, as well as heat, particle and momentum sources. This paper gives a first comprehensive description of the CRONOS suite: overall structure of the code, main available models, details on the simulation workflow and numerical implementation. Some examples of applications to the analysis of experimental discharges and the predictions of ITER scenarios are also given.

Disruptions are a major operational concern for next generation tokamaks, including ITER. They may generate excessive heat loads on plasma facing components, large electromagnetic forces in the machine structures and several MA of multi-MeV runaway electrons. A more complete understanding of the runaway generation processes and methods to suppress them is necessary to ensure safe and reliable operation of future tokamaks. Runaway electrons were studied at JET-ILW showing that their generation dependencies (accelerating electric field, avalanche critical field, toroidal field, MHD fluctuations) are in agreement with current theories.

The role played by electron density fluctuations near the plasma edge on rf current drive in tokamaks is assessed quantitatively. For this purpose, a general framework for incorporating density fluctuations in existing modelling tools has been developed. It is valid when rf power absorption takes place far from the fluctuating region of the plasma. The ray-tracing formalism is modified in order to take into account time-dependent perturbations of the density, while the Fokker-Planck solver remains unchanged. The evolution of the electron distribution function in time and space under the competing effects of collisions and quasilinear diffusion by rf waves is determined consistently with the time scale of fluctuations described as a statistical process.Using the ray-tracing code C3PO and the 3D linearized relativistic bounce-averaged Fokker-Planck solver LUKE, the effect of electron density fluctuations on the current driven by the lower hybrid (LH) and the electron cyclotron (EC) waves is estimated quantitatively. A thin fluctuating layer characterized by electron drift wave turbulence at the plasma edge is considered. The effect of fluctuations on the LH wave propagation is equivalent to a random scattering process with a broadening of the poloidal mode spectrum proportional to the level of the perturbation. However, in the multipass regime, the LH current density profile remains sensitive to the ray chaotic behaviour, which is not averaged by fluctuations. The effect of large amplitude fluctuations on the EC driven current is found to be similar to an anomalous radial transport of the fast electrons. The resulting lower current drive efficiency and broader current profile are in better agreement with experimental observations. Finally, applied to the ITER ELMy H-mode regime, the model predicts a significant broadening of the EC driven current density profile with the fluctuation level, which can make the stabilization of neoclassical tearing mode potentially more challenging.