Edge Localized Modes Energy Loss Research Articles

The simulation of ELM control by the advanced divertor configuration in EAST

Abstract Edge localized modes (ELMs) are effectively suppressed in the ‘quasi-snowflake’ (QSF) divertor discharges, which has been observed in the Experimental Advanced Superconducting Tokamak (EAST). To obtain the physical mechanism of ELM suppression, the numerical simulations are carried out using the BOUT++ turbulence model. The simulations reveal that the large local magnetic shear near the outer mid-plane (OMP) induced by QSF divertor plays a key role in the ELM suppression. Using the EFIT code, a series of plasma equilibria with different 2nd X-points and nearly fixed last closed flux surfaces (LCFSs) are generated to analyze the effects of the different magnetic configurations on ELMs. Here we mainly discuss the standard single-null (SN), snowflake plus (SF+), and snowflake minus (SF-) divertors. The simulation results indicate that: (1) for linear instability, compared to SN, SF+ is more unstable, while SF- is more stable. Essentially, the local magnetic shear formed by different divertor geometries can alter the growth rate of the peeling-ballooning (P-B) mode. Through statistical analysis, there is an inverse correlation between the strength of local magnetic shear and the growth rate of P-B mode; (2) for ELM energy loss, SN is 4.60%, SF+ is 7.50%, and SF- is 0.35%. The SF+ divertor triggers a larger ELM, which is consistent with the TCV experiments; while the SF- divertor reduces the ELM amplitude, which is similar to the QSF experiments in EAST. Further analysis shows that the Reynolds stress determines the ELM size under different divertor configurations. The Reynolds stress can redistribute energy to fluctuations and cause the growth of low-n modes. What’s more, the SF- divertor not only suppresses the radial transport, but also has large magnetic flux expansion and connection length, which can reduce the target heat flux effectively. The conclusion of this paper shows that the advanced divertor configurations are promising for the future fusion.

Effect of negative triangularity on peeling-ballooning instability

Experiments have achieved high confinement discharges in tokamaks with a negative triangularity (NT) plasma shape accompanied by a lower pedestal and smaller and more frequent edge localized modes (ELMs) compared with positive triangularity (PT). Some existing theories emphasize the linear instability variations result from the change of pedestal. However, NT can directly bring significant changes on magnetic field structures which may also influence the instability of ELMs. Based on a series of equilibria constructed with different triangularities and pressure profiles, the influence of NT on peeling-ballooning mode (P–B mode) is investigated. It is found that NT can increase the growth rates of low to intermediate n (toroidal mode number) modes in the linear stage and lead to a larger pedestal collapse in the nonlinear stage if its pressure profile is the same with the PT shape. Further analyses demonstrate that NT enlarges the unfavorable curvature area, which provides stronger driving source and larger unstable region for the instability. Meanwhile, the diamagnetic effect and local magnetic shear helps to stabilize high n modes in the linear phase, and the E × B shearing rate at the top of the pedestal contributes to suppress the transport of turbulence into the plasma core in the nonlinear phase for the NT shape. What’s more, further simulations with different pedestal heights demonstrate that there exists a threshold value of pressure ratio, below which the ELM energy loss in NT shapes can be smaller than that in PT shapes, suggesting that the smaller energy loss with NT in experiment mainly results from the lower pedestal heigh. The results reveal behaviors of P–B modes and provide possible mechanisms for the phenomenon of lower pedestal height with negative triangularities in experiments.

Effect of pedestal density on the formation of small edge localized modes

Recent experiments have demonstrated that a high separatrix density and a large ratio of separatrix density to pedestal top density are two crucial conditions for achieving high confinement operation with small edge localized modes (ELMs). In order to identify the underlying physics of this phenomenon, a series of equilibria with different separatrix and pedestal top densities are constructed, and their peeling–ballooning (P–B) instabilities are analyzed through simulation. It is found that there is a threshold value of pedestal top density which comes from competition between ion inertia and diamagnetic effect, and ELM energy loss can be minimized at the threshold value for a fixed separatrix density. When the pedestal top density is smaller than the threshold value, the ion inertial effect induced by the density profile has a significant influence on the growth of ELMs, resulting in an increased linear growth rate and more ELM energy loss by trigging low-n modes (n being the toroidal mode number) in the nonlinear phase. When the pedestal top density is larger than the threshold value, the diamagnetic effect is the main factor determining the mode spectrum, which moves to the high-n region with a larger growth rate and the nonlinear ELM energy loss increases. However, for a fixed pedestal top density, a larger separatrix density leads to a wider mode spectrum with a smaller growth rate; thus ELM energy loss is reduced. The results of this research provide a new mechanism, namely that the P–B mode is possibly transferred to a resistive ballooning mode, to interpret the experimental findings during high pedestal density operation.

Numerical study on the peeling–ballooning modes with electron cyclotron wave injection

The influence of electron cyclotron wave (ECW) injection with different deposition positions and injection powers on the evolution of peeling–ballooning (P–B) modes is studied with the BOUT++ code, in which the energy deposition and current drive of the ECW are calculated using a ray tracing code. It is shown that the changes in the profiles of plasma pressure, current density, and resistivity induced by ECW injection can significantly influence the linear property and the nonlinear evolution of P–B modes. For the linear simulation, the ECW deposited at the top of the pedestal makes the high toroidal mode number (n) P–B modes more unstable; however, it stabilizes the medium-n to high-n P–B modes when the ECW is deposited at the middle of the pedestal, and the ECW deposited at the bottom of the pedestal decreases the growth rate of P–B modes with medium-n. Further investigation shows that the injected ECW influences the characteristic of linear P–B modes by changing the diamagnetic effect, magnetic shear, pressure gradient, current density, resistivity, and so on. It is known from the nonlinear simulation that the energy loss caused by the edge localized mode (ELM) with ECW injection deposited at the top of the pedestal is nearly the same as that in the case without ECW injection, while the ECW deposited at the middle and bottom of the pedestal is helpful to decrease ELM energy loss. According to the analyses of the time evolution of the P–B mode toroidal spectrum, the physical mechanism of the decrease in ELM energy loss in the simulation is that ECW injection suppresses the most unstable toroidal harmonic of the P–B mode. On the other hand, the influence of ECW injection on P–B modes becomes more obvious when the power of the injected ECW increases. Moreover, the influence of current driven by the ECW on P–B modes is studied separately in this paper, which plays a different role from the bootstrap current.

The simulation of ELMs mitigation by pedestal coherent mode in EAST using BOUT++

A general phenomenon that the edge localized modes (ELMs) can be effectively mitigated with the enhanced coherent modes (CMs) has been observed on EAST. For this phenomenon, the experimental statistical analysis and electromagnetic (EM) simulations have been performed. There is a threshold value of the CM intensity in the experiments, which plays a key role in ELMs mitigation. Through the ELITE and conventional BOUT++ analysis, we found that when the insignificant ELM and enhanced CM co-exist, the pedestal is located in unstable P–B region and the ELM is relatively large. The simulation results only using the experimental profiles without considering other factors cannot reproduce the no significant ELM experiment. The CM enhances the edge turbulence, which can control ELMs. Therefore, the effects of CM are considered to explain the ELM mitigation. Modifying the three-field reduced model in BOUT++, an imposed perturbation is added as the CM. The simulation results indicate that: without the CM, the ELM size belongs to the relative large ELM region; after considering the CM, the ELM is mitigated and the energy loss is reduced by about 44.5%. Analysis shows that the CM enhances the three-wave nonlinear interactions in the pedestal and reduces the phase coherence time (PCT) between the pressure and potential, which lead the perturbation to tend to be ‘multiple-mode’ coupling. The competition of free energy between the multiple modes leads to the lack of obvious filament structures and the decreased energy loss. The above reveals that there is a competitive relationship between turbulence and ELMs, and the CM-enhanced turbulence can effectively reduce ELM energy loss. In addition, through the parameter scanning, there is a threshold of the amplitude A, which is consistent with the statistical results in the experiments.

The simulation of ELMs mitigation by pedestal coherent mode in EAST using BOUT++

Abstract A general phenomenon that the edge localized modes (ELMs) can be effectively mitigated with the enhanced coherent modes (CMs) has been observed on EAST. For this phenomenon, the experimental statistical analysis and electromagnetic (EM) simulations have been performed. There is a threshold value of the CM intensity in the experiments, which plays a key role in ELMs mitigation. Through the ELITE and conventional BOUT++ analysis, we found that when the insignificant ELM and enhanced CM co-exist, the pedestal is located in unstable P-B region and the ELM is relatively large. The simulation results only using the experimental profiles without considering other factors cannot reproduce the no significant ELM experiment. The CM enhances the edge turbulence, which can control ELMs. Therefore, the effects of CM are considered to explain the ELM mitigation. Modifying the three-field reduced model in BOUT++, an imposed perturbation is added as the CM. The simulation results indicate that: without the CM, the ELM size belongs to the relative large ELM region; after considering the CM, the ELM is mitigated and the energy loss is reduced by about 44.5%. Analysis shows that the CM enhances the three-wave nonlinear interactions in the pedestal and reduces the phase coherence time (PCT) between the pressure and potential, which lead the perturbation to tend to be ‘multiple-mode’ coupling. The competition of free energy between the multiple modes leads to the lack of obvious filament structures and the decreased energy loss. The above reveals that there is a competitive relationship between turbulence and ELMs, and the CM-enhanced turbulence can effectively reduce ELM energy loss. In addition, through the parameter scanning, there is a threshold of the amplitude A, which is consistent with the statistical results in the experiments.

Type-I ELM power loads on the closed outer divertor targets in the HL-2A tokamak

The HL-2A tokamak has a very closed divertor geometry, and a new infrared camera has been installed for high resolution studies of edge-localized mode (ELM) heat load onto the outer divertor targets. The characteristics of power deposition patterns on the lower outer divertor target plates during ELMs are systematically analysed with infrared thermography. The ELM energy loss is in the range of 3%–8% of the total plasma stored energy. The peak heat flux on the outer divertor targets during ELMs currently achieved in HL-2A is about 1.5–3.2 MW m−2, the wetted area is about 0.5–0.7 m2, and the corresponding integrated power decay length at the midplane is about 25–40 mm. The rise time of the ELM power deposition is in the range of about 100 μs to 400 μs, and the decay time is typically 1.5 to 4 times longer than the corresponding rise time. Convective transport along open field lines during the ELM rise phase from the midplane towards the divertor targets is implied due to the correlation of parallel transport time in the scrape-off layer (SOL) and ELM power rise time. The peak ELM energy fluence is compared with those predicted by models and with experimental data from JET, ASDEX Upgrade, MAST, and COMPASS. The results, as a whole, show a good agreement.

Kinetic-fluid coupling simulations of ITER Type I ELM

Edge localized modes (ELMs)-induced transient heat loads on the divertor targets represent an important threat to target lifetime and can lead to the need to replace them with a frequency that has a major impact in the execution of the ITER Research Plan. Predicting the impact of such large transient heat loads through modeling is especially challenging and is often attempted through the use of fluid plasma boundary modeling codes, such as SOLPS-ITER, in which the ELM is crudely approximated as a fixed large, but limited in time, increase in anomalous cross-field transport coefficients for particles and heat to mimic a specified total ELM energy loss. However, one problem with this approach is that the boundary conditions at the target sheath entrance are expected to vary strongly in time through the ELM transient, whilst fixed kinetic target sheath heat transmission factors (SHTF), and more generally, constant heat flux limiters, are typically applied in the fluid codes.This contribution describes the first results of efforts to address ELMs issues for ITER simulations under high performance conditions using the 1D3V electrostatic parallel Particle-in-Cell (PIC) code BIT1, to study the kinetic effects and to provide time-dependent kinetic target sheath heat transmission factors. In a second step of the work, these are used in the formulation of fluid boundary conditions for calculations of ELM target heat loads using the SOLPS-ITER code.The BIT1-SOLPS-ITER coupling allows us to investigate the kinetic effects on the targets, by comparing power and particle fluxes from time-dependent simulations of ITER Type I ELMs.

Edge localized mode characteristics and divertor heat flux during stationary and transient phase for CFETR hybrid scenario

The study of edge localized mode (ELM) behavior, stationary heat flux and transient heat flux during the ELM crash phase is performed for a China Fusion Engineering Test Reactor (CFETR) 1 GW hybrid mode operation scenario ( = 7.2 m, = 6.5 T, = 13.78 MA). Modeling and simulation start with a scenario obtained by multi-code integrated modeling on the one modeling framework for integrated tasks framework. Linear stability and nonlinear simulations of ELM dynamics are carried out using the BOUT++ six-field reduced magnetohydrodynamic module, which show a much smaller ELM energy loss (ΔELM ~ 0.13%) compared to that of a Type-I ELM. Parametric analysis of the weak linear growth rate and small ELM energy loss characteristics shed light on physics corresponding to a grassy ELM regime for CFETR 1 GW hybrid scenario. The transient heat flux on the divertor target during this small ELM phase is investigated using BOUT++. We found that upstream radial transport in the scrape-off-layer (SOL) induced by small ELMs is weak, which keeps it in the drift-dominated region. However, the heat flux width is still broadened to = 4.64 mm by the increase of separatrix temperature during ELM nonlinear evolution. The impact of transient peak heat load and ELM energy fluence on tungsten melting and net erosion rate of divertor target is evaluated for the first-time using physics-based transport that connects the pedestal with the SOL. Energy fluence caused by a single ELM pulse is below the tungsten melting limit, while tungsten erosion would exceed the material requirements. We conclude that external mitigation methods, such as divertor detachment and advanced divertor geometry are likely needed for the steady state operation of CFETR.

Nonlinear dependence of plasma potential and ion impinging energy on energy loss of edge localized modes

Particle-in-cell (PIC) modelling has been performed to investigate the impact of energy loss during edge localized modes (ELMs) on the plasma potential and ion impinging energy on the divertor target. A double-peak structure of the ion impinging energy has been identified under JET-relevant ELM conditions. The ELM burst leads to a strong increase in the potential drop in front of the target plate, which accelerates the cold ions from the downstream divertor and accordingly causes a peak value of ion impinging energy. Moreover, the great potential drop helps confine the fast electrons and leads to a reduction in the potential drop and ion impinging energy. The arrival of the upstream hot ions results in the second peak value of ion impinging energy. The maximum potential drop and ion impact energy show a linear dependence on the pedestal temperature. Further, a nonlinear dependence of the peak potential drop and ion impact energy on the ELM energy loss can be ascertained based on the PIC simulations.

Toward integrated multi-scale pedestal simulations including edge-localized-mode dynamics, evolution of edge-localized-mode cycles, and continuous fluctuations

The high-fidelity BOUT++ two-fluid code suite has demonstrated significant recent progress toward integrated multi-scale simulations of tokamak pedestal, including Edge-Localized-Mode (ELM) dynamics, evolution of ELM cycles, and continuous fluctuations, as observed in experiments. Nonlinear ELM simulations show three stages of an ELM event: (1) a linear growing phase; (2) a fast crash phase; and (3) a slow inward turbulence spreading phase lasting until the core heating flux balances the ELM energy loss and the ELM is terminated. A new coupling/splitting model has been developed to perform simulations of multi-scale ELM dynamics. Simulation tracks five ELM cycles for 10 000 Alfvén times for small ELMs. The temporal evolution of the pedestal pressure is similar to that of experimental measurements for the pedestal pressure profile collapses and recovers to a steep gradient during ELM cycles. To validate BOUT++ simulations against experimental data and develop physics understanding of the fluctuation characteristics for different tokamak operation regimes, both quasi-coherent fluctuations (QCFs) in ELMy H-modes and Weakly Coherent Modes in I-modes have been simulated using three dimensional 6-field 2-fluid electromagnetic model. The H-mode simulation results show that (1) QCFs are localized in the pedestal region having a predominant frequency at f≃300−400 kHz and poloidal wavenumber at kθ≃0.7 cm−1, and propagate in the electron diamagnetic direction in the laboratory frame. The overall signatures of simulation results for QCFs show good agreement with C-Mod and DIII-D measurements. (2) The pedestal profiles giving rise to QCFs are near the marginal instability threshold for ideal peeling-ballooning modes for both C-Mod and DIII-D, while the collisional electromagnetic drift-Alfvén wave appears to be dominant for DIII-D. (3) Particle diffusivity is either smaller than the heat diffusivity for DIII-D or similar to the heat diffusivity for C-Mod. Key I-mode simulation results are that (1) a strong instability exists at n≥20 for resistive ballooning mode and drift-Alfvén wave; (2) the frequency spectrum of nonlinear BOUT++ simulation features a peak around 300 kHz for the mode number n = 20, consistent with a reflectometer measurement at nearby position; (3) the calculated particle diffusivity is larger than the calculated heat diffusivity, which is consistent with a key feature of the I-mode pedestal with no particle barrier.

Effect of resonant magnetic perturbations on low collisionality discharges in MAST and a comparison with ASDEX Upgrade

Sustained edge localized mode (ELM) mitigation has been achieved on MAST and AUG using resonant magnetic perturbations (RMPs) with various toroidal mode numbers over a wide range of low to medium collisionality discharges. The ELM energy loss and peak heat loads at the divertor targets have been reduced. The ELM mitigation phase is typically associated with a drop in plasma density and overall stored energy. In one particular scenario on MAST, by carefully adjusting the fuelling it has been possible to counteract the drop in density and to produce plasmas with mitigated ELMs, reduced peak divertor heat flux and with minimal degradation in pedestal height and confined energy. While the applied resonant magnetic perturbation field can be a good indicator for the onset of ELM mitigation on MAST and AUG there are some cases where this is not the case and which clearly emphasize the need to take into account the plasma response to the applied perturbations. The plasma response calculations show that the increase in ELM frequency is correlated with the size of the edge peeling-tearing like response of the plasma and the distortions of the plasma boundary in the X-point region. In many cases the RMPs act to increase the frequency of type I ELMs, however, there are examples where the type I ELMs are suppressed and there is a transition to a small or type IV ELM-ing regime.

Impact of the pedestal plasma density on dynamics of edge localized mode crashes and energy loss scaling

The latest BOUT++ studies show an emerging understanding of dynamics of edge localized mode (ELM) crashes and the consistent collisionality scaling of ELM energy losses with the world multi-tokamak database. A series of BOUT++ simulations are conducted to investigate the scaling characteristics of the ELM energy losses vs collisionality via a density scan. Linear results demonstrate that as the pedestal collisionality decreases, the growth rate of the peeling-ballooning modes decreases for high n but increases for low n (1 < n < 5), therefore the width of the growth rate spectrum γ(n) becomes narrower and the peak growth shifts to lower n. Nonlinear BOUT++ simulations show a two-stage process of ELM crash evolution of (i) initial bursts of pressure blob and void creation and (ii) inward void propagation. The inward void propagation stirs the top of pedestal plasma and yields an increasing ELM size with decreasing collisionality after a series of micro-bursts. The pedestal plasma density plays a major role in determining the ELM energy loss through its effect on the edge bootstrap current and ion diamagnetic stabilization. The critical trend emerges as a transition (1) linearly from ballooning-dominated states at high collisionality to peeling-dominated states at low collisionality with decreasing density and (2) nonlinearly from turbulence spreading dynamics at high collisionality into avalanche-like dynamics at low collisionality.

Simulation of peeling-ballooning modes with pellet injection

The influence of pellet ablation on the evolution of peeling-ballooning (P-B) modes is studied with BOUT++ code. The atoms coming from pellet ablation can significantly reshape the plasma pressure profile, so the behaviors of P-B modes and edge localized mode (ELM) are modified dramatically. This paper shows that the energy loss associated with an ELM increases substantially over that without the pellet, if the pellet is deposited at the top of the pedestal. On the contrary, for pellet deposition in the middle of the pedestal region the ELM energy loss can be less.

Progress on the application of ELM control schemes to ITER scenarios from the non-active phase to DT operation

Progress in the definition of the requirements for edge localized mode (ELM) control and the application of ELM control methods both for high fusion performance DT operation and non-active low-current operation in ITER is described. Evaluation of the power fluxes for low plasma current H-modes in ITER shows that uncontrolled ELMs will not lead to damage to the tungsten (W) divertor target, unlike for high-current H-modes in which divertor damage by uncontrolled ELMs is expected. Despite the lack of divertor damage at lower currents, ELM control is found to be required in ITER under these conditions to prevent an excessive contamination of the plasma by W, which could eventually lead to an increased disruptivity. Modelling with the non-linear MHD code JOREK of the physics processes determining the flow of energy from the confined plasma onto the plasma-facing components during ELMs at the ITER scale shows that the relative contribution of conductive and convective losses is intrinsically linked to the magnitude of the ELM energy loss. Modelling of the triggering of ELMs by pellet injection for DIII-D and ITER has identified the minimum pellet size required to trigger ELMs and, from this, the required fuel throughput for the application of this technique to ITER is evaluated and shown to be compatible with the installed fuelling and tritium re-processing capabilities in ITER. The evaluation of the capabilities of the ELM control coil system in ITER for ELM suppression is carried out (in the vacuum approximation) and found to have a factor of ∼2 margin in terms of coil current to achieve its design criterion, although such a margin could be substantially reduced when plasma shielding effects are taken into account. The consequences for the spatial distribution of the power fluxes at the divertor of ELM control by three-dimensional (3D) fields are evaluated and found to lead to substantial toroidal asymmetries in zones of the divertor target away from the separatrix. Therefore, specifications for the rotation of the 3D perturbation applied for ELM control in order to avoid excessive localized erosion of the ITER divertor target are derived. It is shown that a rotation frequency in excess of 1 Hz for the whole toroidally asymmetric divertor power flux pattern is required (corresponding to n Hz frequency in the variation of currents in the coils, where n is the toroidal symmetry of the perturbation applied) in order to avoid unacceptable thermal cycling of the divertor target for the highest power fluxes and worst toroidal power flux asymmetries expected. The possible use of the in-vessel vertical stability coils for ELM control as a back-up to the main ELM control systems in ITER is described and the feasibility of its application to control ELMs in low plasma current H-modes, foreseen for initial ITER operation, is evaluated and found to be viable for plasma currents up to 5–10 MA depending on modelling assumptions.

Non-linear MHD simulation of ELM energy deposition

The mechanisms of edge-localized mode (ELM) energy deposition are studied by means of non-linear magnetohydrodynamic (MHD) simulation of ELMs. The footprint of the ELM heat flux at the divertor is found to increase approximately linearly with the total ELM energy loss for JET-scale plasmas, which is similar to the experimentally observed broadening of the ELM energy deposition with ELM energy loss. For these relatively large ELMs, in which conductive losses dominate, the divertor footprint broadening is due to an increase in the magnetic perturbation of the ballooning mode with increasing ELM energy loss, which results in a widening of the homoclinic tangles intersecting the target. The first results from ELM simulations in the ITER Q = 10 scenario indicate that on the ITER scale the broadening is similar for conductive and convective ELMs at least up to an ELM energy loss of 4 MJ. For the larger conductive-type ELMs the magnetic perturbation and its homoclinic tangles determine the pattern of the ELM heat flux at the divertor target similar to the JET-scale results. For the smaller convective ELMs, the ELM footprint is determined by the radial distance travelled by plasma filaments expelled by the ELM and the loss of the plasma energy in the filaments along the magnetic field lines.

Reduction of ELM energy loss by pellet injection for ELM pacing

The energy loss caused by the edge-localized mode (ELM) needs to be reduced for ITER operations with ELMy H-mode plasmas. The reduction in ELM energy loss by pellet injection for ELM pacing is studied by an integrated core/scrape-off layer/divertor transport code TOPICS-IB with a magnetohydrodynamic stability code and a pellet model taking account of the E × B drift of the pellet plasma cloud. It is found that the energy loss can be significantly reduced by a pellet injected to the pedestal plasma equivalent to that at the middle timing in the natural ELM cycle, whose pressure height is only about 5% lower than that of the natural ELM onset. In this case, pellet injection from the low-field side enables a small pellet, with about 1–2% of pedestal particle content and a speed high enough to approach the pedestal top, to reduce the energy loss significantly. With the above suitable conditions for ELM pacing, a pellet penetrates deep into the pedestal and triggers high-n ballooning modes with localized eigenfunctions near the pedestal top, where n is the toroidal mode number. Under suitable conditions, ELM pacing with reduced energy loss is successfully demonstrated in simulations, in which the gas puff reduction and the enhancement of divertor pumping can compensate for the core density increase due to additional particle fuelling by the pacing pellet.

Overview of physics results from MAST towards ITER/DEMO and the MAST Upgrade

New diagnostic, modelling and plant capability on the Mega Ampère Spherical Tokamak (MAST) have delivered important results in key areas for ITER/DEMO and the upcoming MAST Upgrade, a step towards future ST devices on the path to fusion currently under procurement. Micro-stability analysis of the pedestal highlights the potential roles of micro-tearing modes and kinetic ballooning modes for the pedestal formation. Mitigation of edge localized modes (ELM) using resonant magnetic perturbation has been demonstrated for toroidal mode numbers n = 3, 4, 6 with an ELM frequency increase by up to a factor of 9, compatible with pellet fuelling. The peak heat flux of mitigated and natural ELMs follows the same linear trend with ELM energy loss and the first ELM-resolved Ti measurements in the divertor region are shown. Measurements of flow shear and turbulence dynamics during L–H transitions show filaments erupting from the plasma edge whilst the full flow shear is still present. Off-axis neutral beam injection helps to strongly reduce the redistribution of fast-ions due to fishbone modes when compared to on-axis injection. Low-k ion-scale turbulence has been measured in L-mode and compared to global gyro-kinetic simulations. A statistical analysis of principal turbulence time scales shows them to be of comparable magnitude and reasonably correlated with turbulence decorrelation time. Te inside the island of a neoclassical tearing mode allow the analysis of the island evolution without assuming specific models for the heat flux. Other results include the discrepancy of the current profile evolution during the current ramp-up with solutions of the poloidal field diffusion equation, studies of the anomalous Doppler resonance compressional Alfvén eigenmodes, disruption mitigation studies and modelling of the new divertor design for MAST Upgrade. The novel 3D electron Bernstein synthetic imaging shows promising first data sensitive to the edge current profile and flows.

DIII-D research towards resolving key issues for ITER and steady-state tokamaks

The DIII-D research program is addressing key ITER research needs and developing the physics basis for future steady-state tokamaks. Pellet pacing edge-localized mode (ELM) control in the ITER configuration reduces ELM energy loss in proportion to 1/fpellet by inducing ELMs at up to 12× the natural ELM rate. Complete suppression of ELMs with resonant magnetic perturbations has been extended to the q95 expected for ITER baseline scenario discharges, and long-duration ELM-free QH-mode discharges have been produced with ITER-relevant co-current neutral-beam injection (NBI) using external n = 3 coils to generate sufficient counter-Ip torque. ITER baseline discharges at βN ∼ 2 and scaled NBI torque have been maintained in stationary conditions for more than four resistive times using electron cyclotron current drive (ECCD) for tearing mode suppression and disruption avoidance; active tracking with steerable launchers and feedback control catch these modes at small amplitude, reducing the ECCD power required to suppress them. Massive high-Z gas injection into disruption-induced 300–600 kA 20 MeV runaway electron (RE) beams yield dissipation rates ∼10× faster than expected from e–e collisions and demonstrate the possibility of benign dissipation of such REs should they occur in ITER. Other ITER-related experiments show measured intrinsic plasma torque in good agreement with a physics-based model over a wide range of conditions, while first-time main-ion rotation measurements show it to be lower than expected from neoclassical theory. Core turbulence measurements show increased temperature fluctuations correlated with sharply enhanced electron transport when exceeds a critical-gradient scale length. In H-mode, data show the pedestal height and width growing between ELMs with ∇P at the computed kinetic-ballooning limit, in agreement with the EPED model. Successful modification of a neutral-beam line to provide 5 MW of adjustable off-axis injection has enabled sustained operation at βN ∼ 3 with broader current and pressure profiles at higher qmin than previously possible, though energy confinement is lower than expected. Initial experiments aimed at developing integrated core and boundary solutions demonstrated heat flux reduction using enhanced edge radiation from neon injection and innovative divertor geometries (e.g. snowflake configuration).

Understanding edge-localized mode mitigation by resonant magnetic perturbations on MAST

Sustained edge-localized mode (ELM) mitigation has been achieved using resonant magnetic perturbations (RMPs) with a toroidal mode number of n = 4 and n = 6 in lower single null and with n = 3 in connected double null plasmas on MAST. The ELM frequency increases by up to a factor of eight with a similar reduction in ELM energy loss. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. A comparison of the filament structures observed during the ELMs in the natural and mitigated stages shows that the mitigated ELMs have the characteristics of type I ELMs even though their frequency is higher, their energy loss is reduced and the pedestal pressure gradient is decreased. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. The size of these lobes is correlated with the increase in ELM frequency observed. The RMPs produce a clear 3D distortion to the plasma and it is likely that these distortions explain why ELMs are destabilized and hence why ELM mitigation occurs.