Edge Localized Modes Suppression Research Articles

Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak

The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψN) or large (>4% ψN) island width depending on the magnitude of the edge bootstrap current, where ψN is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χǁ/χ⊥) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.

Non-linear free boundary simulations of the plasma response to resonant magnetic perturbations in ASDEX Upgrade plasmas

A promising method for the control of Edge Localized Modes (ELMs) in H-Mode tokamak plasmas is the application of Resonant Magnetic Perturbations (RMPs), where small helical field perturbations are introduced into the plasma via a set of external coils. While RMPs are used for suppression of ELMs in many present-day tokamaks, the mechanisms that lead to RMP-ELM control are still subject of debate.Here, we use the non-linear MHD code JOREK to investigate the penetration of the magnetic perturbation fields into ASDEX Upgrade (AUG) plasmas. We present an extension of the coupled JOREK-STARWALL code, that replaces the commonly used fixed boundary treatment with a free boundary treatment. Instead of prescribing the magnetic field at the boundary according to the vacuum field using Dirichlet boundary conditions, natural boundary conditions are applied, so that the magnetic field and plasma current density are evolving freely at the boundary. This allows a fully self-consistent development of the plasma response and the magnetic perturbation in the whole computational domain. The direct comparison of both approaches demonstrates that the artificial suppression of the plasma response with the fixed boundary treatment reduces the excitation of marginally stable modes. An overall larger perturbation is observed using the free boundary approach.The presented simulations are performed in realistic geometry with fully realistic plasma parameters and plasma flows based on reconstructions of experimental AUG equilibria. While the use of realistic plasma parameters makes the simulations particularly challenging, it also allows for quantitative comparisons to the experiment. When the RMP induced corrugation of the boundary is compared to electron density measurements from the lithium beam emission spectroscopy, only the free boundary approach shows excellent agreement with the experiment.

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

We present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

Recent EAST Experimental Results and Systems Upgrade in Support of Long-Pulse Steady-State Plasma Operation

In recent years, EAST has achieved a 60-s scale steady-state plasma operation (SSO), which is electron dominant heated by the electron cyclotron wave and low hybrid wave. Substantial progresses have been made in the development of long-pulse operation, such as active control of radiative divertor, edge localized mode (ELM) suppression, and He plasma experiment. To further develop the steady-state operation scenario and study the related critical physical issues, the upgrade of EAST has been completed in 2020. The auxiliary heating systems were rearranged and upgraded with a total heating power capability of 33 MW. The 2.45-GHz lower hybrid current drive (LHCD) launcher has been updated from full active multijunction (FAM) to passive active multijunction (PAM) to avoid damage caused by plasma–material interaction. The injection power capability of electron cyclotron resonance heating (ECRH) system is upgraded to 1.4 MW. The upgraded ion cyclotron resonance frequency (ICRF) antenna injects low parallel wavenumber ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$k_{\vert \vert }$ </tex-math></inline-formula> ) spectra with good single pass absorption. The neutral beam injection (NBI) system is moved from F-port to D-port and is reoriented from counterclockwise to clockwise injection. The armor material of lower divertor has been upgraded from graphite to tungsten improving the steady-state heat handling capability from 2 to 10 MW/m2. Based on the upgraded pumping system and improved conductance, the effective pumping speed for <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$D_{2}$ </tex-math></inline-formula> of the EAST lower divertor increased twice. With all these upgrades, 100-s scale steady-state operation ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$T_{e0} &gt;10$ </tex-math></inline-formula> keV) and 1056-s high-performance operation with 1.73-GJ injection energy were achieved in 2021 campaigns. The experience in high-performance long-pulse operation obtained on EAST will lay a foundation for the ITER and Chinese Fusion Engineering Testing Reactor (CFETR) operation.

Nonlinear MHD modeling of n = 1 RMP-induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

Fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

Study on divertor heat flux under n = 3 and n = 4 resonant magnetic perturbations using infrared thermography diagnostic in EAST

Resonant magnetic perturbations (RMPs) with high toroidal mode number n are considered for controlling edge-localized modes (ELMs) and divertor heat flux in future ITER H-mode operations. In this paper, characteristics of divertor heat flux under high-n RMPs (n = 3 and 4) in H-mode plasma are investigated using newly upgraded infrared thermography diagnostic in EAST. Additional splitting strike point (SSP) accompanying with ELM suppression is observed under both RMPs with n = 3 and n = 4, the SSP in heat flux profile agrees qualitatively with the modeled magnetic footprint. Although RMPs suppress ELMs, they increase the stationary heat flux during ELM suppression. The dependence of heat flux on during ELM suppression is preliminarily investigated, and further splitting in the original strike point is observed at during ELM suppression. In terms of ELM pulses, the presence of RMPs shows little influence on transient heat flux distribution.

The simulation of ELM suppression by ion cyclotron resonance heating in EAST using BOUT++

Edge localized modes (ELMs) are completely suppressed by ion cyclotron resonant heating (ICRH) during H-mode discharge, which is first observed in the Experimental Advanced Superconducting Tokamak (EAST). The electromagnetic turbulence model in the BOUT++ framework is applied to the simulations, with the aim of obtaining the key physics mechanisms of ELM suppression by ICRH. The simulations, for the EAST experiment with no ELMs after ICRH, reveal that the radio frequency (RF) sheath plays a key role in the suppression of ELMs. The simulation results show that: (a) using the flow-balanced electric field, the ELM sizes are about 3.4% before ICRH, and 2.1% after ICRH. Although the smaller ELM is obtained with ICRH, it still belongs to relatively large ELMs; (b) instead of the flow-balance condition, the radial electric field E r is calculated using the RF sheath model, and the ELM size is about 0.36% with ICRH, much closer to the small/no ELM regime. After considering the effect of ICRH, the larger E × B shear rate in the pedestal and scrape-off layer induced by the RF sheath reduces the linear growth rate, and enhances the nonlinear mode coupling. Furthermore, we find that there is an effective range of the sheath potential, in which the ELM can be well suppressed, but beyond this range the ELM may be only mitigated rather than suppressed. This can reasonably explain why the probability of ELM suppression by ICRH appearing in the experiments is low.

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. The main aim of this work was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ∼20 cm in inner divertor and ∼40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ∼5 MW m−2 to ∼1 MW m−2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

Design and experimental demonstration of feedback adaptive RMP ELM controller toward complete long pulse ELM suppression on KSTAR

Operation of a fusion power plant requires robust edge localized mode (ELM) suppression simultaneously with high plasma performance. In this paper, we describe a novel feedback adaptive resonant magnetic perturbation (RMP) ELM controller designed to address this problem by achieving optimized ELM suppression through the advanced application of 3D RMPs. From real-time Dα data, the controller is able to achieve robust ELM suppression while simultaneously minimizing the applied RMP in order to enhance plasma performance. In real-time, the instantaneous ELM-frequency is analyzed with an adaptive feedback algorithm to determine amplitudes and phases of RMP coil currents that will maximize plasma performance while maintaining ELM suppression. When applied through the KSTAR plasma control system in several experiments using n = 1 RMPs, robust ELM suppression is achieved and sustained in feedback while reducing the RMP strength to ∼65% of its initial value. Minimization of the RMP strength in this manner not only allows for operation of longer discharges due to a decrease in flux consumption but also allows for a strong recovery of up to ∼60% of βN throughout the ELM-free period.

Dependence of the L–H power threshold on the alignment of external non-axisymmetric magnetic perturbations in ASDEX Upgrade

Experiments at ASDEX Upgrade have been conducted to study the impact of the alignment of external magnetic perturbations (MPs) with n = 2 toroidal mode symmetry on the power threshold of the transition from L- to H-mode (PLH). This is interesting for ITER, since its operation will rely on an MP field configuration that securely suppresses edge localized modes (ELMs) while avoiding an increase in the L–H power threshold. PLH can increase up to ≈80% when the external MP field is aligned to maximize the associated plasma response at the edge. This alignment deviates from the equilibrium field pitch and is similar to the one that is required to suppress ELMs by MPs in H-mode. The relative radial magnetic field perturbation δBr/BT has to exceed a critical value, roughly 1.7−2×10−4 (at the q = 5 surface), to cause an increase in PLH. This value is above the one that is required to sustain ELM suppression at ASDEX Upgrade. Combining plasma response calculations with a critical value of δBr/BT captures well the measured dependencies of PLH on the alignment of the applied MP field. The increase in PLH is accompanied by a change in the E×B flow velocity (vE×B) profile. The application of MPs in L-mode flattens and reverses the edge vE×B profile when at least δBr/BT≈1.9×10−4 is applied. This is the same value above which PLH increases. More heating power is needed to roughly achieve the same shear in the vE×B profile at the transition to H-mode as without MPs but with an vE×B profile that is shifted into the ion diamagnetic direction.

Optimization of 3D controlled ELM-free state with recovered global confinement for KSTAR with n = 1 resonant magnetic field perturbation

Mitigation of deleterious heat flux from edge-localized modes (ELMs) on fusion reactors is often attempted with 3D perturbations of the confining magnetic fields. However, the established technique of resonant magnetic perturbations (RMPs) also degrades plasma performance, complicating implementation on future fusion reactors. In this paper, we introduce an adaptive real-time control scheme on the KSTAR tokamak as a viable approach to achieve an ELM-free state and simultaneously recover high-confinement (β N ∼ 1.91, β p ∼ 1.53, and H 98 ∼ 0.9), demonstrating successful handling of a volatile complex system through adaptive measures. We show that, by exploiting a salient hysteresis process to adaptively minimize the RMP strength, stable ELM suppression can be achieved while actively encouraging confinement recovery. This is made possible by a self-organized transport response in the plasma edge which reinforces the confinement improvement through a widening of the ion temperature pedestal and promotes control stability, in contrast to the deteriorating effect on performance observed in standard RMP experiments. These results establish the real-time approach as an up-and-coming solution toward an optimized ELM-free state, which is an important step for the operation of ITER and reactor-grade tokamak plasmas.

Edge localized modes suppression via edge E × B velocity shear induced by RF sheath of ion cyclotron resonance heating in EAST

Edge localized modes suppression via edge E × B velocity shear induced by RF sheath of ion cyclotron resonance heating in EAST

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

Simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

First demonstration of full ELM suppression in low input torque plasmas to support ITER research plan using n = 4 RMP in EAST

Full suppression of type-I edge localized modes (ELMs) using n = 4 resonant magnetic perturbations (RMPs) as planned for ITER has been demonstrated for the first time (n is the toroidal mode number of the applied RMP). This is achieved in EAST plasmas with low input torque and tungsten divertor, and the target plasma for these experiments in EAST is chosen to be relevant to the ITER Q = 10 operational scenario, thus also addressing significant scenario issues for ITER. In these experiments the lowest neutral beam injection (NBI) input torque is around T NBI ∼ 0.44 Nm, which extrapolates to around 14 Nm in ITER (compared to a total torque input of 35 Nm when 33 MW of NBI are used for heating). The q 95 is around 3.6 and normalized plasma beta β N ∼ 1.5–1.8, similar to that in the ITER Q = 10 scenario. Suppression windows in both q 95 and plasma density are observed; in addition, lower plasma rotation is found to be favourabe to access ELM suppression. ELM suppression is maintained with line averaged density up to 60%n GW (Greenwald density limit) by feedforward gas fuelling after suppression is achieved. It is interesting to note that in addition to an upper density, a low density threshold for ELM suppression of 40%n GW is also observed. In these conditions energy confinement does not significantly drop (<10%) during ELM suppression when compared to the ELMy H-mode conditions, which is much better than previous results using low n (n = 1 and 2) RMPs in higher q 95 regimes. In addition, the core plasma tungsten concentration is clearly reduced during ELM suppression demonstrating an effective impurity exhaust. MHD response modelling using the MARS-F code shows that edge magnetic field stochasticity has a peak at q 95 ∼ 3.65 for the odd parity configuration, which is consistent to the observed suppression window around 3.6–3.75. These results expand the physical understanding of ELM suppression and demonstrate the effectiveness of n = 4 RMPs for reliable control ELMs in future ITER high Q plasma scenarios with minimum detrimental effects on plasma confinement.

Diverted negative triangularity plasmas on DIII-D: the benefit of high confinement without the liability of an edge pedestal

Diverted discharges at negative triangularity on the DIII-D tokamak sustain normalized confinement and pressure levels typical of standard H-mode scenarios (H 98y2 ≃ 1, β N ≃ 3) without developing an edge pressure pedestal, despite the auxiliary power far exceeding the L → H power threshold expected from conventional scaling laws. The power degradation of confinement is substantially weaker than the ITER-89P scaling, resulting in a confinement factor that improves with increasing auxiliary power. The absence of the edge pedestal is beneficial in several ways, such as eliminating the need for active mitigation or suppression of edge localized modes, low impurity retention and a reconstructed scrape-off layer heat flux width at the mid-plane that exceeds the ITPA multi-machine scaling law by up to 50%. Together with technological advantages granted by placing the divertor at larger radii, plasmas at negative triangularity without an edge pedestal feature both core confinement and power handling characteristics that are potentially suitable for operation in future fusion reactors.

Integrated ELM and divertor power flux control using RMPs with low input torque in EAST in support of the ITER research plan

Experiments have been carried out at the EAST tokamak to study ITER-relevant scenario integration issues, related to edge localized mode (ELM) control in H-mode plasmas by the application of three-dimensional (3D) resonant magnetic perturbations (RMPs), which have a large impact on the execution of the ITER research plan. The EAST experiments have successfully demonstrated ELM suppression at normalized torque inputs similar to ITER. The application of RMP fields with high toroidal mode number (n = 4) reduces the impact of ELM control on energy and particle confinement compared to those use lower n (n = 1, 2) RMPs. Injection of successive pellets is found to be effective in increasing the plasma density in ELM-suppressed H-modes and reducing the divertor power without triggering large ELMs at EAST. Access to high recycling and radiative divertor conditions while maintaining ELM suppression has been demonstrated in EAST by the use of gas fuelling and neon impurity seeding. Both approaches have been found to be effective in reducing power fluxes to the divertor strike points in near-separatrix lobes for both n = 2 and n = 4 RMPs. However, reduction of power fluxes in off-separatrix lobes is only effective for n = 4 RMP application, which is consistent with magnetic topology modelling (including plasma response) results showing a shallow penetration into the confined plasma region of field lines connected to these lobes compared to n = 2. The EAST results support the use of high n 3D fields for ELM suppression in ITER high Q DT scenarios since they provide optimum integration features regarding energy and particle confinement, pellet fuelling, radiative divertor operation while eliminating ELM transient power loads and being compatible with low torque input.

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q 95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95. The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

Plasma response to resonant magnetic perturbations near rotation zero-crossing in low torque plasmas

The plasma response to resonant magnetic perturbations (RMPs) near the pedestal top is crucial for accessing edge localized modes (ELMs) suppression in tokamaks. Since the radial location of the rotation zero-crossing plays a key role in determining the threshold for field penetration of RMPs, the plasma response may be different in low input torque plasmas. In this work, the linear MHD code MARS-F is applied to reveal the dependence of the plasma response to RMPs on the location of the rotation zero-crossing by a scan of E × B rotation profiles based on an EAST equilibrium. It is shown that the plasma response is enhanced when the E × B rotation zero-crossing occurs near the rational surfaces. The dependence of the plasma response on the location of the E × B rotation zero-crossing is well fitted by a double Gaussian, indicating two effects in this enhancement. One is induced by a rotation screening effect shown as a wide base (with a width around 10–20 krad/s), and the other is related to a resistive singular layer effect characterized by a localized peak (with a width around 3–4 krad/s). The peak of each resonant harmonic in the plasma response appears always at the rotation zero-crossing. The width of the peak scales with the resistive singular layer width. The plasma displacement suggests that the response is tearing-like when the zero-crossing is within the singular layer, while it is kink-like when the zero-crossing is far from the layer. The reduction in the magnetic islands width is only by a factor of two, when the absolute value of the local E × B rotation increases from 0 to 10–20 krad/s. This is further confirmed in modeling of the plasma response in an EAST ELM suppression discharge. Although there is a zero-crossing in the E × B rotation but none in the electron perpendicular rotation, no significant difference in the plasma response is obtained using these two rotation profiles. This suggests that the rotation near the pedestal top should not be far away from zero but it may not be necessary to have a zero-crossing for accessing ELM suppression.

Effects of resonant magnetic perturbations on radial electric fields in DIII-D tokamak

Gyrokinetic simulations of DIII-D tokamak equilibrium find that resonant magnetic perturbation (RMP) drives a neoclassical non-ambipolar electron particle flux, which causes a rapid change of equilibrium radial electric fields consistent with experimental observations during the suppression of the edge localized mode (ELM). The simulation results provide a support for the conjecture that RMP-induced changes of radial electric fields lead to the enhanced turbulent transport at the pedestal top during the ELM suppression (Taimourzadeh et al 2019 Nucl. Fusion 59 046005). Furthermore, gyrokinetic simulations of collisionless damping of zonal flows show that resonant responses to the RMP decrease the residual level of the zonal flows and damp the geodesic acoustic mode.

Modeling of q95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the KSTAR tokamak

The EPEC code is employed to model the q95 windows for n = 2 and n = 1 resonant magnetic perturbation (RMP)-induced edge localized mode (ELM) suppression in typical KSTAR H-mode discharges. The plasma equilibria used in the study are derived by rescaling the experimental plasma equilibrium in KSTAR discharge #18594 measured at time t = 6450 ms. The n = 2 and n = 1 ELM suppression windows predicted by EPEC are comparatively narrow in q95 (i.e., Δq95≃0.1) and are such that (assuming a relative pedestal pressure reduction of 15% is required to trigger ELM suppression), even in the middle of the windows, there is a threshold RMP coil current of about 1–3 kA/turn that must be exceeded before ELM suppression occurs. The n = 2 and n = 1 ELM suppression windows calculated by the EPEC code are consistent with the experimental observations and are also broadly similar to those recently calculated by the TM1 code [Hu et al., Phys. Plasmas 28, 052505 (2021)].