Large ELMs Research Articles

Nonlinear dynamics and effects of fast-ion driven instabilities in HL-2A NBI heating high-βN H-mode plasmas

High-performance plasmas had been obtained with hybrid scenarios on HL-2A in past several years, and the quality factor G=H89βN/q952 is around 0.3-0.4. In these high-βN H-mode plasmas with βN>2.5 and H98>1.2 Chen et al. (2022) [24], there are multiple energetic-ion driven instabilities, such as toroidal Alfvén eigenmode (TAE) and beta-induced Alfvén eigenmode (BAE), and there are complex magnetohydrodynamic (MHD) dynamics, e.g. nonlinear mode-mode and mode-particle interactions. The TAE amplitude increases explosively, and the frequency is fast chirping in the range of f=90−150 kHz. The TAE toroidal mode-number is n=3−4, and they localize in outer plasma region. The BAE amplitude also increases explosively, and the BAE frequency is fast chirping in the range of f=50−90 kHz. The toroidal mode-number of dominant BAEs is n=2, and these modes localize at q=2 rational surface where there is a large radius internal transport barrier (ITB). Generally, the energetic-ion fraction is low in outer plasma region. It indicates that thermal ions at large radius ITB location possibly provide additional drive and destabilize the BAE in company with energetic-ions. These TAEs/BAEs occur regularly in high-βN region and before large type-I ELM burst. To further illuminate the characteristics of these Alfvénic modes, the 2D and radial mode structures had been calculated by using a well-benchmarked general Alfvénic mode eigenvalue code (MAS). The experimental results suggest that these Alfvénic modes may have important effects on the pedestal transport and ELM-onset.

Overview of recent experimental results on the EAST Tokamak

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (∼1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high β P has been significantly extended (H 98y2 > 1.3, β P ∼ 4.0, β N ∼ 2.4 and n e/n GW ∼ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q 95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m−1) with EC pre-ionization is developed. A beneficial role on the lower hybrid wave injection to control the tungsten concentration in the NBI discharge is observed for the first time in EAST suggesting a potential way toward steady-state H-mode NBI operation.

High-frequency fluctuation and EHO-like mode in the H-mode pedestal on the EAST tokamak

In the pedestal region of the Experimental Advanced Superconducting Tokamak (EAST) during high confinement mode plasma operations with radio-frequency heating, two distinct fluctuations are observed: high-frequency fluctuations (HFFs) and edge harmonic oscillation-like (EHO-like) modes. The HFFs are characterized by intermittent fluctuations with a broadband frequency range of 1−3 MHz and a poloidal wave number ( kθ ) greater than 0.9 cm−1 . On the other hand, the EHO-like mode exhibits characteristics similar to magnetohydrodynamics (MHD)-like modes with n = 1−5 and lower poloidal wave numbers ( kθ⩽0.12 cm−1 ). During the pedestal establishing phase following the L–H transition, a significant concurrent presence of HFF and EHO-like modes in high-density pedestal regions has been noted. In this phase, the EHO-like mode not only modulates the amplitude of the HFF but also engages in nonlinear interactions. The occurrence of EHO-like mode and HFF is associated with particle transport toward the divertor, though it is notably less than that caused by edge coherent modes. During the inter-edge localized mode (ELM) period, a significant decrease in the Dα baseline is observed whenever the low frequency fluctuation (LFF) weakens and the HFF grows, prior to each large ELM. One possible explanation is that the rapid increase of E×B shear stabilizes the LFF and destabilizes the HFF, which lowers the pedestal transport and enables the further growth of the pedestal until the onset of the ELM.

Progress in the development of the ITER baseline scenario in TCV

Under the auspices of EUROfusion, the ITER baseline (IBL) scenario has been jointly investigated on AUG and TCV in the past years and this paper reports on the developments on TCV. Three ITER shapes, namely the JET, AUG and ITER IBL have been reproduced in TCV, illustrating that the higher the triangularity the larger the ELM perturbation and the more difficult it is to reach stationary states with q 95< 3.6. It is found that the performance of TCV IBL is mainly limited by (neoclassical) tearing modes, in particular 2/1 modes which are triggered after a large ELM. It is demonstrated that the shorter the ELM period the larger β N at the NTM onset. We show that these modes can be avoided with central X3 EC heating at relatively high q 95 and moderate β N . However, the lack of significant ECH at the high central densities obtained in TCV IBL scenario limits the duration of low q 95 cases to about four confinement times. During this time, density usually keeps peaking until (neoclassical) tearing modes are triggered. Nevertheless, the TCV IBL database covers the ITER target values ( H98y2∼ 1, βN∼ 1.8 at q95∼ 3) and a slightly better confinement than requested for ITER is reported. Integrated modelling results show that ITG modes are the dominant instabilities, and show that, in TCV, fuelling also plays a role to sustain peaked density profiles. The role of profiles, sawteeth and ELMs regarding MHD stability are also discussed.

Compatibility of large ELM control and stable partial detachment with neon/argon seeding in EAST

It is necessary to achieve simultaneous exhaust of excessive transient and steady-state heat fluxes on the divertor target for the divertor protection in the future fusion reactors. The sustained large ELM control and stable partial detachment have been achieved concurrently with argon (Ar) or neon (Ne) seeding in EAST. With Ne seeding, the large ELMs with frequency f ELM ∼ 100 Hz disappear and a stable ELM-free state with H 98,y2 > 1 is maintained. Meanwhile, the electron temperature T et around the lower outer strike point decreases from more than 70 eV during the large ELM burst to less than 5 eV in the stable ELM-free phase. In addition, a slight improvement of plasma confinement is observed in the partially detached state, mainly attributed to the increased electron density n e and ion temperature T i in the core region. In the pedestal region, the density gradient and the electron temperature show subtle variation. The effective charge number Z eff increases significantly after Ne seeding, leading to a decrease in the edge bootstrap current and the pedestal pressure gradient, and thus the stabilization of ELMs. With Ar seeding, the large ELMs are also suppressed at first, but soon transit to type-III ELMs with a high f ELM ∼ 1 kHz, highly correlated with the energy confinement degradation. The steady-state and transient heat fluxes on the divertor can be both well reduced with Ar/Ne seeding in EAST.

Physical mechanisms for the transition from type-III to large ELMs induced by impurity injection on EAST

Transition from type-III to large-amplitude ELMs induced by neon injection has been observed in the EAST tokamak at overlapping q95 space between large and small ELMs. With neon injection, pedestal density gradient shows a remarkable increase accompanied by some decrease of pedestal electron temperature, and consequently the pressure gradient increases moderately and edge bootstrap current has minimal change. Further experiment demonstrates that the occurrence of large ELMs after neon injection is highly correlated with the change in edge density. Linear peeling-ballooning stability analysis indicates that the large ELM case is more unstable than the type-III ELM case during the ELM transition. A scan of pedestal density gradient in linear stability analysis shows that the direct destabilizing effect of steep pedestal density gradient on peeling-ballooning instabilities via two-fluid effects could also facilitate the transition to large ELMs. These results could provide more insight into the role of pedestal density gradient on pedestal stability and ELM behavior.

Study on pedestal fluctuations in H-modes without large ELMs during the transition to a detached tungsten divertor in EAST

Simultaneous control of the damaging erosion induced by the transient and steady-state heat/particle fluxes on the divertor target material is one of the critical issues for next-step magnetic fusion devices. H-mode operation without large edge-localized modes has been achieved in EAST with an ITER-like tungsten divertor, while being compatible with the partial and pronounced detachment in divertor, via either ramping-up of bulk density or injection of low/high-Z impurities. The pedestal characteristics during the transition from the attached to the detached divertor and the reversed transition (detached to attached) under different detachment methods are studied in detail, where the evolution of multi fluctuating structures commonly residing in the H-mode pedestal of EAST (edge coherent mode (ECM), magnetic coherent mode (MCM) and high frequency mode (HFM)) is highlighted. In addition, the possible mechanisms that affect the behavior of these modes, such as the pedestal pressure gradient and the collisionality, have also been discussed. The radial structures of ECM, MCM and HFM are detected, for the first time, in one discharge. Relevant research may provide contribution to obtaining an integrated small/no ELM and radiative divertor scenario in the next step.

Towards understanding reactor relevant tokamak pedestals

The physics of the tokamak pedestal is still not fully understood, for example there is no fully predictive model for the pedestal height and width. However, the pedestal is key in determining the fusion power for a given scenario. If we can improve our understanding of reactor relevant pedestals we will improve our confidence in designing potential fusion power plants. Work has been carried out as part of a collaboration on reactor relevant pedestal physics. We report some of the results in detail here and review some of the wider work which will be reported in full elsewhere. First, we attempt to use a gyrokinetic-based calculation to eliminate the pedestal top density as a model input for Europed/EPED pedestal predictions. We assume power balance at the top of the pedestal, that is, the heat flux crossing the separatrix must be equal to the heat source at the top of the pedestal and investigate the consequences of this assumption. Unfortunately, the transport assumptions of the EPED model mean that this method does not discriminate between different pairs of density and temperature profiles for a given pressure profile. Second, we investigate the effects of non flux surface density on the bootstrap current. Third, type I ELMs will not be tolerable for a reactor relevant regime due to the damage that they are expected to cause to plasma facing components. In recent years various methods of running tokamak plasmas without large ELMs have been developed. These include small and no ELM regimes, the use of resonant magnetic perturbations and the use of vertical kicks. We discuss the quiescent H-mode here. Finally we give a summary and directions for future work.

Рост и состояние смешанных древостоев лесных насаждений Каменной Степи

By now, a vast amount of experimental material has been accumulated in the field of protective afforestation of the Central Chernozem Region of Russia, which allows estimating the stability and longevity of various trees and shrubs in different environmental conditions. Good growth and high vitality of woody vegetation largely depend on growing conditions and are determined primarily by the relationship between species (individuals within a species) when they are placed in the forest area and the type of terrain. The researchdifferent combinations in tree and shrub mixing schemes and unequal location by types of terrain. The study was carried out in mature (115–118-year-old) forest strips laid out by the staff of the Kamennaya Steppe experimental forestry G.F. Morozov and N.A. Mikhailov on the territory of the Kamennaya Steppe (Voronezh region, Talovsky district) using archival materials, scientific works of the authors of this article and the employees of the Department of Agroforestry. A comparative analysis of the materials of forest survey work carried out in protective forest plantations created according to the tree-shrub type of mixing on different types of terrain revealed an excess of biometric indicators of tree species growing on the upland type of terrain over those on the slope. It is shown that the initial percentage of participation of English oak is a significant, but not always decisive factor in the creation of oak plantations. They can also be grown with a smaller proportion of this species in the culture, but in this case, timely thinning will be of paramount importance. Over time, there have been significant changes in the composition of plantations, the number of trees and their valuation indicators. There are no viable ash and elm specimens left in the upper tier of the stands. They are severely damaged by stem pests. Certain tree species are damaged by the following pests: common ash by large ash bark beetle (Hylesinus crenatus) and ash bark beetle (Hylesinus fraxini); elm species by cambium-feeding beetles (large elm bark beetle (Scolytus scolytus), European elm bark beetle (Scolytus multistriatus), and pygmy elm bark beetle (Scolytus pygmaeus); English oak by gold pit oak splendour beetle (Chrysobothris affinis), oak borer (Agrilus angustulus), longhorn beetle (Plagionotus detritus), and European oak bark beetle (Scolytus intricatus). The predominant number of trees of all tree species belongs to the categories of limited viable and inviable.

Evolution of ELMs, pedestal profiles and fluctuations in the inter-ELM period in NBI- and ECH-dominated discharges in DIII-D

In DIII-D, it has been observed that ELM frequency decreases by 40% and ELM spacing becomes more regular in time when heating is changed from pure neutral beam injection (NBI) to predominantly electron cyclotron heating (ECH) in ITER similar shape plasmas. In comparison with the pure NBI discharges, pedestal fluctuations in magnetics and density increase in the ECH-dominated discharges. Recovery of the pedestal profiles like electron density (n e), temperature (T e) and pressure (p e) shows marked differences for these two heating schemes. Average profiles in the last 30% of the ELM cycle show higher T e, lower n e, and similar p e at the pedestal top for the ECH discharge when compared to the NBI discharge. The gradient of T e (∇T e) is also steeper at the pedestal in the ECH discharge. Magnetic fluctuations show three distinct modes at 13–116 kHz in the ECH discharges only. n e fluctuations show two modes evolving in the inter-ELM period of the ECH discharge, a low-frequency (400 kHz) quasi-coherent mode (LFQC) and high-frequency (∼2 MHz) broadband (HFB) fluctuations. Evolution of these modes has marked correspondence with the inter-ELM ∇T e recovery. A sharp decrease in the D α baseline is observed whenever the LFQC weakens and the HFB grows, prior to each large ELM. Transport coefficients obtained from TRANSP show that MTM and/or TEM are plausible candidates for the observed fluctuations. Linear gyrofluid simulation (TGLF) corroborates this characterization. TGLF shows that the linear growth rate of the most dominant mode peaks at ion-scale (k θ ρ s ∼ 0.4) at the pedestal steep gradient and the frequency is in the electron diamagnetic direction. It is proposed that increased fluctuations in the ECH-dominated case, due to increased ∇T e, caused an increase in fluctuation-driven transport in the pedestal and slowed the pedestal recovery between ELMs, leading to a reduction in the ELM frequency.

Impact of plasma-wall interaction and exhaust on the EU-DEMO design

Abstract In the present work, the role of plasma facing components protection in driving the EU-DEMO design will be reviewed, focusing on steady-state and, especially, on transients. This work encompasses both the first wall (FW) as well as the divertor. In fact, while the ITER divertor heat removal technology has been adopted, the ITER FW concept has been shown in the past years to be inadequate for EU-DEMO. This is due to the higher foreseen irradiation damage level, which requires structural materials (like Eurofer) able to withstand more than 5 dpa of neutron damage. This solution, however, limits the tolerable steady-state heat flux to ~1 MW/m2, i.e. a factor 3–4 below the ITER specifications. For this reason, poloidally and toroidally discontinuous protection limiters are implemented in EU-DEMO. Their role consists in reducing the heat load on the FW due to charged particles, during steady state and, more importantly, during planned and off-normal plasma transients. Concerning the divertor configuration, EU-DEMO currently assumes an ITER-like, lower single null (LSN) divertor, with seeded impurities for the dissipation of the power. However, this concept has been shown by numerous simulations in the past years to be marginal during steady-state (where a detached divertor is necessary to maintain the heat flux below the technological limit and to avoid excessive erosion) and unable to withstand some relevant transients, such as large ELMs and accidental loss of detachment. Various concepts, deviating from the ITER design, are currently under investigation to mitigate such risks, for example in-vessel coils for strike point sweeping in case of reattachment, as well as alternative divertor configurations. Finally, a broader discussion on the impact of divertor protection on the overall machine design is presented.

Transport at high and development of candidate steady state scenarios for ITER

On DIII-D, the high scenario has an internal transport barrier (ITB), , , and very high normalized confinement . Recently, plasmas starting with these conditions have been dynamically driven to and , where we find the ITB and high performance persist for five energy confinement times. These conditions are projected to meet the ITER steady-state goal of Q = 5. The ITB is maintained at lower with a strong reverse shear, consistent with predictions that negative central shear can lower the threshold for the ITB. There are two observed confinement states in the high scenario: H-mode confinement state with a high edge pedestal, and an enhanced confinement state with a low pedestal and an ITB. It has been observed in a scan of external resonant magnetic perturbation amplitude that when there are no large type-I ELMs, there is no transition to enhanced confinement. This is consistent with the proposed mechanism for ITB formation being a type-I ELM. Quasilinear gyro-Landau fluid predictive modeling of ITER suggests that only a modest reverse shear is required to achieve the ITB formation necessary for Q = 5 when electromagnetic physics including the kinetic ballooning mode (KBM) is incorporated.

Design of a PAM launcher and comparative analysis with the old FAM launcher for 2.45 GHz LHCD system on EAST

In order to support EAST to realize the main mission of long-pulse operation lasting 1000s, the existing full active multijunction (FAM) launcher will be upgraded to a passive active multijunction (PAM) for the 2.45 GHz lower hybrid current drive (LHCD) system. The RF components were designed and optimized with consideration of the port size, power density inside the waveguides to avoid RF breakdown and power spectrum for efficient current drive. The designed PAM has a power handling capability of 4 MW CW at 2.45 GHz with matched load. The coupling properties and power spectrum were calculated using the ALOHA code and a comparative analysis with the old FAM was carried out. The main feature of this new launcher is that the density for optimum coupling and power spectrum directivity (ne ˜ 0.3 × 1017/m3) lies below the cut-off value (ne_co ˜ 0.74 × 1017/m3), which is favored for long-pulse steady-state operation on EAST. It has a capability to deliver 3.6 MW of power during real operation (ne ˜ 2 × 1017/m3) and 2 MW during large ELMs (ne = 10 ne_co). For the nominal phasing of 180 °s between modules, the power spectrum is peaked at N// = 1.8 and a maximum power directivity of 71.2% and reflection coefficient less than 1% can be obtained via simulations. Weaker hot spot phenomenon is expected with the designed PAM than with the FAM for their respective optimum densities. Finally, a thermal analysis was performed and the maximum temperature at the front of the antenna is expected to be ˜ 260 °C which is safe enough for operation. The PAM is scheduled to be installed on EAST in 2021.

Observation of a fully non-inductive H-mode regime dominated by the sporadic-small edge-localized modes in EAST with a tungsten divertor

Here we report a fully non-inductive H-mode operation with a sporadic-small-edge-localized mode(ELM)-dominant pedestal realized in EAST with a tungsten divertor, where the normalized electron collisionality at the pedestal top is estimated below unity. This is regarded as the expansion and/or supplement of our previous research having been preliminarily reported (section 4.2 in (B N Wan et al 2017 Nucl. Fusion 57 102019)) in the EAST overview article contributed to the special issue on FEC 2016 summaries and overviews. In this operation, spontaneous avoidance of large ELMs is realized reproducibly, for the first time on EAST, at a lower collisionality (below unity) at the pedestal top, profiting from high heating power. Besides, this regime shows the well controlled bulk plasma density and light impurity concentration, good energy confinement with H98, (y2) ≃ 1, together with a loop voltage very close to zero. All of these features are considered critical to a steady-state operation with high performance, which EAST is on target to demonstrate in the near future in support to the International Thermonuclear Experimental Reactor (ITER). However, this operation regime suffers a continuous accumulation of medium-/high-Z impurities, which may seriously limit its sustainability. The resonant magnetic perturbation (RMP) with a dominant toroidal mode number n = 1 has been employed for assisting exhaust of heavy impurities, which is, unfortunately, not successful since the RMP with an effective poloidal spectrum always brings back large ELMs. This scenario transits to an ELM H-mode regime as plasma current is decreased, and elimination of large ELMs is well reproduced by ELITE simulation. Finally, we present a discussion on the potential mechanism eliminating large ELMs, which seems to highlight the LHW-induced profound change in edge magnetic topology.

Impact of impurity mixture gas seeded by supersonic molecular beam injection on edge-localized modes in the HL-2A tokamak

Impurity seeding experiments in ELMy H-mode plasmas by using supersonic molecular beam injection (SMBI) have been performed in the HL-2A tokamak. Besides pure impurity gas (neon, argon), impurities mixed with the main ion fuelling gas by different ratios were seeded into the plasma. It has been observed that the ELM behavior and plasma confinement are very sensitive to the impurity ratio. For 30% Ne impurity seeding, large ELMs are replaced by high frequency bursts (HFBs), which have higher frequency and smaller amplitude compared to the ELMs and the peak heat load on the divertor plate was considerably reduced. Analyses indicate that HFBs originate from the pedestal region and enhance the pedestal particle transport, resulting in a softened pedestal gradient after SMBI. However, for the pure Ne and Ar SMBI seeding, the ELM frequency is decreased and the plasma confinement is moderately improved. Experimental observations suggest that the impurity mixture gas plays a role in changing the pedestal dynamics and there should be an optimal impurity ratio for efficient heat load control in ELMy H-mode plasmas.

Bursaphelenchus michalskii sp. n. (Nematoda: Aphelenchoididae), a nematode associate of the large elm bark beetle, Scolytus scolytus Fabr. (Coleoptera: Curculionidae), in Dutch elm disease-affected elm, Ulmus laevis Pall.

Summary Bursaphelenchus michalskii sp. n. is described from the bark of the European white elm, Ulmus laevis. All propagative stages of the nematode were found in larval galleries of the large elm bark beetle, Scolytus scolytus, and in overlapping gallery systems of this species and the small European elm bark beetle, S. multistriatus. Dauer juveniles of the new nematode are transmitted to new breeding trees under elytra of adult S. scolytus. Bursaphelenchus michalskii sp. n. is characterised by the female body length of 953 (838-1108) μm and male body length of 893 (811-971) μm, very slender body (a = 53.9 (46.1-58.5) and 60.9 (52.2-72.0) in female and male, respectively), lateral fields with three incisures (two bands), excretory pore usually located anterior to the median bulb, lack of vulval flap, long post-uterine sac, relatively small spicules 12.3 (10.8-13.3) μm long with no cucullus and with distinct, somewhat thorn-like, dorsally bent or reflexed condylus and a conical or digitate rostrum, and the arrangement of the seven male caudal papillae (i.e., a single precloacal ventromedian papilla (P1), one pair of adcloacal ventrosublateral papillae (P2) at or just anterior to cloacal slit, one ventrosublateral, postcloacal pair (P3) located at ca 60% of the tail length, posterior to cloacal slit, and one pair (P4) of ventrosublateral papillae located near the base of the bursa). The newly described species shares most of the key morphological characters with members of the eremus-group (sensu Braasch et al., 2009). However, B. michalskii sp. n. is unique amongst Bursaphelenchus species by a combination of female tail and spicule shape, excretory pore position, and other morphometric characters. These findings were confirmed by DNA sequencing and phylogenetic analysis of the 18S and 28S rDNA regions and by the unique molecular profile of the ITS region (ITS-RFLP).

Experimental study on the magnetic coherent mode in the H-mode pedestal of EAST

Here we report the experimental analysis on the low-n (mostly n = 1, sometimes n = 2) magnetic coherent mode (MCM) at a characteristic frequency 20–, which has been frequently observed in various H-mode discharges on EAST. This mode can be easily identified in the magnetic fluctuations measured by the fast Mirnov coils mounted on the vacuum vessel wall, but is detected by the local measurements of edge electrostatic fluctuations only when the mode is sufficiently strong. The apperance of the MCM is summarized covering broad ranges of discharge parameters, in particular, the different heating schemes including pure neutral beams injected in either co- or counter-current direction as well as pure ratio-frequency waves. This may rule out the possibility of fast particle driven modes. Radial distribution and poloidal propagation of the MCM are investigated using the Doppler backscattering system and Langmuir probes inserted at the outer midplane, respectively. Temporal evolution of MCM amplitude during large ELM crashes is evaluated in detail, may suggesting the mode is closely correlated with pedestal buildup. Dedicated experiments reveal the possible correlations of MCM’s frequencies with edge line-averaged density and edge safety factor q95. We also present the observation of multi MCMs at relatively high q95, which are speculated locating at different rational surfaces in the pedestal via analyzing their mode structures and nonlinear interactions. Finally, effect of the MCM on edge particle transport is explored via surveying the correlation between the intermittent events of the mode and the particle fluxes deposited on the divertor target plates, utilizing the conditional analysis method. Corresponding results suggest that the MCM seems to primarily result in a notable poloidal redistribution of the divertor particle flux, rather than a considerable net increase of the total flux.

Intrinsic ELMing in ASDEX Upgrade and global control system-plasma self-entrainment

It is well established that edge localized modes can be entrained to the frequency of applied global magnetic perturbations. These perturbations are delivered to the plasma using the vertical control system field coil currents. These field coils are part of an active control system that is required to maintain the plasma in a steady state. We perform time domain timeseries analysis of natural ELMing when there are no applied perturbations in the ASDEX Upgrade tokamak. We find that the plasma can transition into a state in which the control system field coil currents continually oscillate and are synchronized with oscillations in characteristic plasma parameters such as plasma edge position and total MHD energy. These synchronous oscillations have a one-to-one correlation with the naturally occurring ELMs; the ELMs all occur when the control system coil current is around a specific temporal phase. Large and small ELMs may be distinguished by the amplitude of inward movement of the edge following an ELM. Large ELMs are then found to occur preferentially around a specific temporal phase of the vertical position control coil current. Small ELMs are most likely in antiphase to this. The large and small natural ELMs occur at the opposite extrema of the oscillations in the control system vertical position control coil current. The control system coil current phase may thus provide a useful parameter to order the observed ELM dynamics. We have identified a class of natural ELMing which is a self-entrained state, in which there is a continual non-linear feedback between the global plasma dynamics and the active control system that is intrinsic to the cyclic dynamics of naturally occurring ELMs. Control system-plasma feedback thus becomes an essential component for integration into future models of natural ELM dynamics.

Grassy-ELM regime with edge resonant magnetic perturbations in fully noninductive plasmas in the DIII-D tokamakaaEdge-localized-mode.

Resonant magnetic perturbations (n = 3 RMPs) are used to suppress large amplitude ELMs and mitigate naturally occurring ‘grassy’-ELMs in DIII-D plasmas relevant to the ITER steady-state mission. Fully non-inductive discharges in the ITER shape and pedestal collisionality ( ≈ 0.05–0.15) are routinely achieved in DIII-D with RMP suppression of the Type-I ELMs. The residual grassy-ELMs deliver a low peak heat flux to the divertor as low as 1.2× the inter-ELM heat flux in plasmas with sustained high H-factor (H98y2 ≈ 1.2). The operating window for the RMP grassy-ELM regime is q95 = 5.3–7.1 and external torque in the range 9–0.7 Nm in the co-Ip direction, which is in the range required for a steady-state tokamak reactor. The RMP grassy-ELM regime is associated with a two-step pedestal, with strong flattening of the density around the zero crossing in the E × B shear. The edge magnetic response of the plasma to the n = 3 RMP is found to be ≈2–3× larger than for comparable ITER baseline plasmas (βN ≈ 1.8, q95 ≈ 3.1). The amplification of the RMP is consistent with the weak magnetic perturbation level (δB/B ≈ 1 × 10−4) required for effective Type-I ELM suppression. Cyclic variations in the pedestal pressure, width, and toroidal rotation are observed in these plasmas, correlated with cyclic variations in the strength and frequency of the grassy-ELMs. Extended MHD analysis and magnetic measurements indicate that these pedestal pulsations are driven by cyclic variations in the resonant field strength at the top of the pedestal. These pedestal pulsations reveal that the grassy-ELMs is correlated with the proximity of the pedestal to the low-n peeling-ballooning mode stability boundary. The use of low amplitude magnetic fields to access grassy-ELM conditions free of Type-I ELMs in high beta poloidal plasmas (βP ≈ 1.5–2.0) opens the possibility for the further optimization of the steady-state tokamak by use of edge resonant magnetic perturbations.

Investigation of the role of pedestal pressure and collisionality on type-I ELM divertor heat loads in DIII-D

A non-dimensional collisionality scan conducted on DIII-D confirms a model for ELM energy densities recently put forward by Eich (2017 Nucl. Mater. Energy 12 84–90), but also reveals key effects that may explain the large scatter typically observed about the scaling. Electron cyclotron heating (ECH) close to the plasma edge was used to raise electron temperatures at the pedestal top and lower collisionality to ITER level, while the power of neutral beam injection (NBI) was decreased during discharges to operate closer to the L–H transition threshold. The scan reveals no explicit pedestal pressure dependence of the ELM energy densities. While collisionality does not play a decisive role, the ratio of heating power to the H-mode threshold power is identified as parameter determining the agreement with the model, with discharges marginally above the threshold showing the largest scatter in the database and exceeding the predicted ELM energy up to twofold. Operation close to the L–H-threshold is accompanied by low ELM frequency and large ELM heat loads. Using linear stability calculations, ELM energy densities are shown to scale inversely with to the most unstable linear mode number before the ELM crash. There are indications that the scatter in the data when compared with the Eich model prediction is caused by including only a limited set out of all quantities considered by linear stability analysis. While further ELM studies near the LH threshold are of great priority, the overall agreement of DIII-D with the Eich model recommends its use in extrapolations towards ITER.