The aim of this paper is to present a signal processing algorithm that, applied to the raw Locked Mode signal, allows us to obtain a disruption indicator in principle exploitable on different tokamaks. A common definition of such an indicator for different machines would facilitate the development of portable systems for disruption prediction, which is becoming of increasingly importance for the next tokamak generations. Moreover, the indicator allows us to overcome some intrinsic problems in the diagnostic system such as drift and offset. The behavior of the proposed indicator as disruption predictor, based on crossing optimized thresholds of the signal amplitude, has been analyzed using data of both JET and ASDEX Upgrade experiments. A thorough analysis of the disruption prediction performance shows how the indicator is able to recover some missed and tardy detections of the raw signal. Moreover, it intervenes and corrects premature or even wrong alarms due to, e.g., drifts and/or offsets.
In ITER and DEMO, various control objectives related to plasma control must be simultaneously achieved by the plasma control system (PCS), in both normal operation as well as off-normal conditions. The PCS must act on off-normal events and deviations from the target scenario, since certain sequences (chains) of events can precede disruptions. It is important that these decisions are made while maintaining a coherent prioritization between the real-time control tasks to ensure high-performance operation. In this paper, a generic architecture for task-based integrated plasma control is proposed. The architecture is characterized by the separation of state estimation, event detection, decisions and task execution among different algorithms, with standardized signal interfaces. Central to the architecture are a plasma state monitor and supervisory controller. In the plasma state monitor, discrete events in the continuous-valued plasma state are modeled using finite state machines. This provides a high-level representation of the plasma state. The supervisory controller coordinates the execution of multiple plasma control tasks by assigning task priorities, based on the finite states of the plasma and the pulse schedule. These algorithms were implemented on the TCV digital control system and integrated with actuator resource management and existing state estimation algorithms and controllers. The plasma state monitor on TCV can track a multitude of plasma events, related to plasma current, rotating and locked neoclassical tearing modes, and position displacements. In TCV experiments on simultaneous control of plasma pressure, safety factor profile and NTMs using electron cyclotron heating (ECH) and current drive (ECCD), the supervisory controller assigns priorities to the relevant control tasks. The tasks are then executed by feedback controllers and actuator allocation management. This work forms a significant step forward in the ongoing integration of control capabilities in experiments on TCV, in support of tokamak reactor operation.
In this paper we validate the finite element code RAPLICASOL, which models radiofrequency wave propagation in edge plasmas near ICRF antennas, against calculations with the TOPICA code. We compare the output of both codes for the ASDEX Upgrade 2-strap antenna, and for a 4-strap WEST-like antenna. Although RAPLICASOL requires considerably fewer computational resources than TOPICA, we find that the predicted quantities of experimental interest (including reflection coefficients, coupling resistances, - and -matrix entries, optimal matching settings, and even radiofrequency electric fields) are in good agreement provided we are careful to use the same geometry in both codes.
Edge localized modes (ELMs) are repetitive instabilities driven by the large pressure gradients and current densities in the edge of H‐mode plasmas. Type‐I ELMs lead to a fast collapse of the H‐mode pedestal within several hundred microseconds to a few milliseconds. Localized transient heat fluxes to divertor targets are expected to exceed tolerable limits for ITER, requiring advanced insights into ELM physics and applicable mitigation methods. This paper describes how non‐linear magneto‐hydrodynamic (MHD) simulations can contribute to this effort. The JOREK code is introduced, which allows the study of large‐scale plasma instabilities in tokamak X‐point plasmas covering the main plasma, the scrape‐off layer, and the divertor region with its finite element grid. We review key physics relevant for type‐I ELMs and show to what extent JOREK simulations agree with experiments and help reveal the underlying mechanisms. Simulations and experimental findings are compared in many respects for type‐I ELMs in ASDEX Upgrade. The role of plasma flows and non‐linear mode coupling for the spatial and temporal structure of ELMs is emphasized, and the loss mechanisms are discussed. An overview of recent ELM‐related research using JOREK is given, including ELM crashes, ELM‐free regimes, ELM pacing by pellets and magnetic kicks, and mitigation or suppression by resonant magnetic perturbation coils (RMPs). Simulations of ELMs and ELM control methods agree in many respects with experimental observations from various tokamak experiments. On this basis, predictive simulations become more and more feasible. A brief outlook is given, showing the main priorities for further research in the field of ELM physics and further developments necessary.
Electron cyclotron emission imaging (ECEI) provides measurements of electron temperature (Te) and its fluctuations (δTe). However, when measuring at the plasma edge, in the steep gradient region, radiation transport effects must be taken into account. It is shown that due to these effects, the scrape-off layer region is not accessible to the ECEI measurements in steady state conditions and that the signal is dominated by the shine-through emission. Transient effects, such as filaments, can change the radiation transport locally, but cannot be distinguished from the shine-through. Local density measurements are essential for the correct interpretation of the electron cyclotron emission, since the density fluctuations influence the temperature measurements at the plasma edge. As an example, a low frequency 8 kHz mode, which causes 10%–15% fluctuations in the signal level of the ECEI, is analysed. The same mode has been measured with the lithium beam emission spectroscopy density diagnostic, and is very well correlated in time with high frequency magnetic fluctuations. With radiation transport modelling of the electron cyclotron radiation in the ECEI geometry, it is shown that the density contributes significantly to the radiation temperature (Trad) and the experimental observations have shown the amplitude modulation in both density and temperature measurements. The poloidal velocity of the low frequency mode measured by the ECEI is 3 km s–1. The calculated velocity of the high frequency mode measured with the magnetic pick-up coils is about 25 km s–1. Velocities are compared with the E × B background flow velocity and possible explanations for the origin of the low frequency mode are discussed.
Edge localized modes (ELMs) are magnetohydrodynamic (MHD) instabilities that cause fast periodic relaxations of the strong edge pressure gradient in tokamak fusion plasmas. A novel diagnostic method allows the extraction of toroidal mode numbers, rotation velocities and spatial information during the ELM cycle including the crash. While mode number branches n=3–6 and n=8–10 are dominant just before the ELM crash, during the ELM crash n=2–5 are observed in typical discharges with type-I ELMs in the tokamak experiment. These findings are compared to results from nonlinear MHD simulations. Although n=6 is linearly dominant, nonlinear coupling in which n=1 is particularly important leads to the dominance of n=3–5 during the ELM crash, in excellent agreement with experimental observations. The simultaneous occurrence of these modes over a wide radial region leads to high stochasticity and thus increased transport.
In ASDEX Upgrade, the propagation of cold pulses induced by type-I edge localized modes (ELMs) is studied using electron cyclotron emission measurements, in a dataset of plasmas with moderate triangularity. It is found that the edge safety factor or the plasma current are the main determining parameters for the inward penetration of the Te perturbations. With increasing plasma current the ELM penetration is more shallow in spite of the stronger ELMs. Estimates of the heat pulse diffusivity show that the corresponding transport is too large to be representative of the inter-ELM phase. Ergodization of the plasma edge during ELMs is a possible explanation for the observed properties of the cold pulse propagation, which is qualitatively consistent with non-linear magneto-hydro-dynamic simulations.
The explanation of the existence of the rotating MHD modes in the pedestal region before Type I edge localized mode (ELM) crash and in the inter-ELM periods (ELM precursors) observed in KSTAR is provided for the first time in the present paper. The dynamics of ELMs, observed using electron cyclotron emission imaging (ECEI) in KSTAR tokamak, is compared to the modelling results of the non-linear reduced resistive MHD code JOREK. The realistic KSTAR pulse parameters and geometry including X-point and scrape off layer (SOL) were used. The full ELM crash modelling was performed using JOREK code for single and multi-harmonic representation and in multi-cycles ELMy regimes including relevant flows. The most unstable toroidal modes numbers (n = 5–8), velocity (~5 km s−1 for n = 8 mode) and the direction of the mode rotation were reproduced in modelling. The two fluid diamagnetic effects and toroidal rotations included in the model were found to be the most important factors in explaining the experimentally observed rotation of the ballooning modes before the ELM crash and in the inter-ELM phase. In multi-harmonic multi-cycle simulations the spectrum of temperature fluctuations is similar to the experimental one in the inter-ELM phase, where several rotating modes with medium n numbers were detected in 5–30 kHz frequency range. The rotating modes can contain single or several harmonics which last from 0.2 ms to few ms in time, and can appear and disappear in the inter ELM period or persist until a new ELM crash.
The plasma response to resonant magnetic perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which edge localized modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces q = m / n , the coupling between the kink component ( m > nq ) and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.
Integrating the plasma core performance with an edge and scrape-off layer (SOL) that leads to tolerable heat and particle loads on the wall is a major challenge. The new European medium size tokamak task force (EU-MST) coordinates research on ASDEX Upgrade (AUG), MAST and TCV. This multi-machine approach within EU-MST, covering a wide parameter range, is instrumental to progress in the field, as ITER and DEMO core/pedestal and SOL parameters are not achievable simultaneously in present day devices. A two prong approach is adopted. On the one hand, scenarios with tolerable transient heat and particle loads, including active edge localised mode (ELM) control are developed. On the other hand, divertor solutions including advanced magnetic configurations are studied. Considerable progress has been made on both approaches, in particular in the fields of: ELM control with resonant magnetic perturbations (RMP), small ELM regimes, detachment onset and control, as well as filamentary scrape-off-layer transport. For example full ELM suppression has now been achieved on AUG at low collisionality with n = 2 RMP maintaining good confinement H-H(98,H-y2) approximate to 0.95. Advances have been made with respect to detachment onset and control. Studies in advanced divertor configurations (Snowflake, Super-X and X-point target divertor) shed new light on SOL physics. Cross field filamentary transport has been characterised in a wide parameter regime on AUG, MAST and TCV progressing the theoretical and experimental understanding crucial for predicting first wall loads in ITER and DEMO. Conditions in the SOL also play a crucial role for ELM stability and access to small ELM regimes.
The plasma response to Resonant Magnetic Perturbations (RMPs) and the interplay between Edge Localized Modes (ELMs) and RMPs is modeled with the non-linear resistive MHD code JOREK, using data from ASDEX Upgrade experimental discharges. The RMP spectrum leading to a maximum amplification of edge kink modes also induce the largest resonant response due to poloidal coupling between modes. In this case, the toroidal coupling between ELMs and RMPs prevents the higher n modes to grow into a large ELM. However when the plasma response to RMPs is minimum, the RMP drive is not strong enough to affect the growth of the ELMs.
Future devices like JT-60SA, ITER and DEMO require quantitative predictions of pedestal density and temperature levels, as well as (inter-ELM and ELM) divertor heat fluxes, in order to improve global confinement capabilities while preventing divertor erosion/melting in the planning of future experiments. Such predictions can be obtained from dedicated pedestal models like EPED, and from non-linear MHD codes like JOREK, for which systematic validation against current experiments is necessary. In this paper, we show progress in the validation of the JOREK code using MAST, KSTAR, ASDEX Upgrade, DIII-D, JT-60U and JET simulations, with both qualitative and quantitative comparisons to experiments, and we present the latest achievements of EUROPED as an extension of the EPED model, to clarify the pedestal width description based on kinetic ballooning modes and turbulence. In addition, we describe how JOREK and EUROPED can interact to improve pedestal predictions in cases where ideal MHD fails to describe experimental observations, which is the case for many type-I ELMs in JET-ILW.
diamagnetic drift effects A. Lessig1, M. Holzl1, I. Krebs4, S. Gunter1, E. Franck1, F. Orain2, J. Morales2, M. Becoulet2, G. Huysmans3, ASDEX Upgrade Team 1 Max-Planck-Institut fur Plasmaphysik, Boltzmannstrasse 2, 85748 Garching, Germany 2 CEA-IRFM, Cadarache, 13108 Saint-Paul-Lez-Durance, France 3 ITER Organization, 13067 Saint-Paul-Lez-Durance, France 4 Max-Planck/Princeton Center for Plasma Physics