Physics research on the TCV tokamak facility: from conventional to alternative scenarios and beyond

The research program of the TCV tokamak ranges from conventional to advanced-tokamak scenarios and alternative divertor configurations, to exploratory plasmas driven by theoretical insight, exploiting the device’s unique shaping capabilities. Disruption avoidance by real-time locked mode prevention or unlocking with electron-cyclotron resonance heating (ECRH) was thoroughly documented, using magnetic and radiation triggers. Runaway generation with high-Z noble-gas injection and runaway dissipation by subsequent Ne or Ar injection were studied for model validation. The new 1 MW neutral beam injector has expanded the parameter range, now encompassing ELMy H-modes in an ITER-like shape and nearly non-inductive H-mode discharges sustained by electron cyclotron and neutral beam current drive. In the H-mode, the pedestal pressure increases modestly with nitrogen seeding while fueling moves the density pedestal outwards, but the plasma stored energy is largely uncorrelated to either seeding or fueling. High fueling at high triangularity is key to accessing the attractive small edge-localized mode (type-II) regime. Turbulence is reduced in the core at negative triangularity, consistent with increased confinement and in accord with global gyrokinetic simulations. The geodesic acoustic mode, possibly coupled with avalanche events, has been linked with particle flow to the wall in diverted plasmas. Detachment, scrape-off layer transport, and turbulence were studied in L- and H-modes in both standard and alternative configurations (snowflake, super-X, and beyond). The detachment process is caused by power ‘starvation’ reducing the ionization source, with volume recombination playing only a minor role. Partial detachment in the H-mode is obtained with impurity seeding and has shown little dependence on flux expansion in standard single-null geometry. In the attached L-mode phase, increasing the outer connection length reduces the in–out heat-flow asymmetry. A doublet plasma, featuring an internal X-point, was achieved successfully, and a transport barrier was observed in the mantle just outside the internal separatrix. In the near future variable-configuration baffles and possibly divertor pumping will be introduced to investigate the effect of divertor closure on exhaust and performance, and 3.5 MW ECRH and 1 MW neutral beam injection heating will be added

Real-time-capable prediction of temperature and density profiles in a tokamak using RAPTOR and a first-principle-based transport model

The RAPTOR code is a control-oriented core plasma profile simulator with various applications in control design and verification, discharge optimization and real-time plasma simulation. To date, RAPTOR was capable of simulating the evolution of poloidal flux and electron temperature using empirical transport models, and required the user to input assumptions on the other profiles and plasma parameters. We present an extension of the code to simulate the temperature evolution of both ions and electrons, as well as the particle density transport. A proof-of-principle neural-network emulation of the quasilinear gyrokinetic QuaLiKiz transport model is coupled to RAPTOR for the calculation of first-principle-based heat and particle turbulent transport. These extended capabilities are demonstrated in a simulation of a JET discharge. The multi-channel simulation requires ∼0.2 s to simulate 1 second of a JET plasma, corresponding to ∼20 energy confinement times, while predicting experimental profiles within the limits of the transport model. The transport model requires no external inputs except for the boundary condition at the top of the H-mode pedestal. This marks the first time that simultaneous, accurate predictions of T e, T i and n e have been obtained using a first-principle-based transport code that can run in faster-than-real-time for present-day tokamaks

Profile control simulations and experiments on TCV: A controller test environment and results using a model-based predictive controller

The successful performance of a model predictive profile controller is demonstrated in simulations and experiments on the TCV tokamak, employing a profile controller test environment. Stable high-performance tokamak operation in hybrid and advanced plasma scenarios requires control over the safety factor profile (q-profile) and kinetic plasma parameters such as the plasma beta. This demands to establish reliable profile control routines in presently operational tokamaks. We present a model predictive profile controller that controls the q-profile and plasma beta using power requests to two clusters of gyrotrons and the plasma current request. The performance of the controller is analyzed in both simulation and TCV L-mode discharges where successful tracking of the estimated inverse q-profile as well as plasma beta is demonstrated under uncertain plasma conditions and the presence of disturbances. The controller exploits the knowledge of the time-varying actuator limits in the actuator input calculation itself such that fast transitions between targets are achieved without overshoot. A software environment is employed to prepare and test this and three other profile controllers in parallel in simulations and experiments on TCV. This set of tools includes the rapid plasma transport simulator RAPTOR and various algorithms to reconstruct the plasma equilibrium and plasma profiles by merging the available measurements with model-based predictions. In this work the estimated q-profile is merely based on RAPTOR model predictions due to the absence of internal current density measurements in TCV. These results encourage to further exploit model predictive profile control in experiments on TCV and other (future) tokamaks

The science program of the TCV tokamak: Exploring fusion reactor and power plant concepts

TCV is acquiring a new 1 MW neutral beam and 2 MW additional third-harmonic electron cyclotron resonance heating (ECRH) to expand its operational range. Its existing shaping and ECRH launching versatility was amply exploited in an eclectic 2013 campaign. A new sub-ms real-time equilibrium reconstruction code was used in ECRH control of NTMs and in a prototype shape controller. The detection of visible light from the plasma boundary was also successfully used in a position-control algorithm. A new bang-bang controller improved stability against vertical displacements. The RAPTOR real-time transport simulator was employed to control the current density profile using electron cyclotron current drive. Shot-by-shot internal inductance optimization was demonstrated by iterative learning control of the current reference trace. Systematic studies of suprathermal electrons and ions in the presence of ECRH were performed. The L-H threshold power was measured to be ∼50-75% higher in both H and He than D, to increase with the length of the outer separatrix, and to be independent of the current ramp rate. Core turbulence was found to decrease from positive to negative edge triangularity deep into the core. The geodesic acoustic mode was studied with multiple diagnostics, and its axisymmetry was confirmed by a full toroidal mapping of its magnetic component. A new theory predicting a toroidal rotation component at the plasma edge, driven by inhomogeneous transport and geodesic curvature, was tested successfully. A new high-confinement mode (IN-mode) was found with an edge barrier in density but not in temperature. The edge gradients were found to govern the scaling of confinement with current, power, density and triangularity. The dynamical interplay of confinement and magnetohydrodynamic modes leading to the density limit in TCV was documented. The heat flux profile decay lengths and heat load profile on the wall were documented in limited plasmas. In the snowflake (SF) divertor configuration the heat flux profiles were documented on all four strike points. SF simulations with the EMC3-EIRENE code, including the physics of the secondary separatrix, underestimate the flux to the secondary strike points, possibly resulting from steady-state E × B drifts. With neon injection, radiation in a SF was 15% higher than in a conventional divertor. The novel triple-null and X-divertor configurations were also achieved in TCV