Search NASA⌕ Search

SEARCH · Search NASA

Results for “tokamak control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Plasma steering to avoid disruptions in ITER and tokamak power plants

Steering tokamak plasmas is commonly viewed as a way to avoid disruptions and runaway electrons. Plasma steering sounds as safe as driving to work but will be shown to more closely resemble driving at high speed through a dense fog on an icy road. The long time required to terminate an ITER discharge compared to time over which dangers can be foreseen is analogous to driving in a dense fog. Here, the difficulty of regaining plasma control if it is lost resembles driving on an icy road. Disruptions and runaways are associated with three issues—a solution to one tends to complicate the solution to the other two: loss of plasma position control, excessive heat deposition, and wall melting due to runaway electrons. All three risks must be addressed for ITER to achieve its mission and essentially eliminated before tokamak power plants can be deployed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First observation of RMP ELM mitigation on MAST Upgrade

Abstract The first experimental attempts at controlling edge localised modes (ELMs) via the application of resonant magnetic perturbations on the MAST Upgrade tokamak are reported. Using the linear MHD model MARS-F, the phase shift between the upper and lower coil rows ΔΦ was optimised for toroidal mode number n = 1 and n = 2 fields, to provide forward guidance to experiments. In low β N discharges, the application of n = 1 3D fields caused the ELM frequency f E L M to increase by over a factor 20 relative to the reference, and also induced a locked mode, which did not cause a plasma termination nor an H-L back transition. However when β N was raised, this induced locked mode caused plasma termination which precluded mitigation access. Initially, applying a numerically optimised n = 2 field had no effect. However applying a rigid toroidal shift to this field caused a locked mode disruption, demonstrating the presence of a substantial n = 2 error field. Coil current ramps were conducted with ΔΦ set at 6 different values, resulting in either locked mode disruptions or no effect, but mitigation with n = 2 fields was not established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stellarators with enhanced tritium confinement and edge radiation control

A stellarator design is described with the purpose of achieving three goals: (1) enhance the confinement time of tritium. (2) Have a sufficient density of high-Z impurities to radiate the thermal power escaping from the core while having an extremely low impurity density in the core. (3) Maintain a large fraction of the plasma in a burning plasma state with an optimal tritium fraction. Some features of this design could be used in tokamaks. Although having three confinement zones is natural for stellarators, it is not for tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantification of locked mode instability triggered by a change in confinement

Abstract This work presents the first analysis of the disruptive locked mode (LM) triggered by the dynamics of a confinement change. It shows that, under certain conditions, the LM threshold during the transient is significantly lower than expected from steady states. We investigate the sensitivity to a controlled n = 1 error field (EF) activated prior to the L-H transition in the COMPASS tokamak, at q 95 ∼ 3, β N ∼ 1, and using EF coils on the high-field side of the vessel. A threshold for EF penetration subsequent to the L-H transition is identified, which shows no significant trend with density or applied torque, and is an apparent consequence of the reduced intrinsic rotation of the 2/1 mode during this transient phase. This finding challenges the assumption made in theoretical and empirical works that natural mode rotation can be predicted by global plasma parameters and urges against using any parametric EF penetration scaling derived from steady-state experiments to define the EF correction strategy in the entire discharge. Furthermore, even at EFs below the identified penetration threshold, disruptive locking of sawtooth-seeded 2/1 tearing modes is observed after about 30% of L-H transitions without external torque.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toward fusion plasma scenario planning for NSTX-U using machine-learning-accelerated models

One of the most promising devices for realizing power production through nuclear fusion is the tokamak. To maximize performance, it is preferable that tokamak reactors achieve advanced operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic (MHD) stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high \emph{fusion gain} --- the ratio of produced fusion power to the external power provided through the plasma boundary. Precise and robust control of the evolution of the plasma boundary shape as well as the spatial distribution of the plasma current, density, temperature, and rotation will be essential to achieving and maintaining such scenarios. The complexity of the evolution of tokamak plasmas, arising due to nonlinearities and coupling between various parameters, motivates the use of model-based control algorithms that can account for the system dynamics. In this work, a learning-based accelerated model trained on data from the National Spherical Torus Experiment Upgrade (NSTX-U) is employed to develop planning and control strategies for regulating the density and temperature profile evolution around desired trajectories. The proposed model combines empirical scaling laws developed across multiple devices with neural networks trained on empirical data from NSTX-U and a database of first-principles-based computationally intensive simulations. The reduced execution time of the accelerated model will enable practical application of optimization algorithms and reinforcement learning approaches for scenario planning and control development. An initial demonstration of applying optimization approaches to the learning-based model is presented, including a strategy for mitigating the effect of leaving the finite validity range of the accelerated model. The approach shows promise for actuator planning between experiments and in real-time.

machine learning↗
DOE OSTI2023

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Multi-device study of temporal characteristics of magnetohydrodynamic modes initiating disruptions

Disruptions in tokamaks are often preceded by magnetohydrodynamic (MHD) instabilities that can rotate or become locked to the wall. Measurements from magnetic diagnostics in the presence of MHD mode precursors to disruptions can yield potentially valuable input to the plasma control system, with a view to disruption avoidance, prediction and mitigation. This report presents an exploratory analysis of the growth of MHD modes and corresponding time scales on the basis of magnetic measurements in multiple tokamaks. To this end, a database was compiled using disruptive discharges from COMPASS, ASDEX Upgrade, DIII-D and JET, manually classified according to disruption root cause, and characterized by a great diversity of operational conditions and mode dynamics. The typical time during which a mode can be detected using saddle coils and the duration of the locked mode phase in the database both extend over several orders of magnitude, but generally the time scales increase with plasma size. Several additional factors are discussed that can influence these durations, including the disruption root cause. A scaling law for the locked phase duration was estimated, yielding predictions toward ITER of the order of hundreds of milliseconds or even seconds. In addition, a scaling law for the mode amplitude at the disruption onset, proposed earlier by de Vries et al. (2016), is applied to the database, and its predictive capabilities are assessed. Despite significant uncertainty on the predictions from both scaling laws, encouraging trends are observed of the fraction of disruptions that may be detected with sufficient warning time to allow mitigation or even avoidance, based solely on observations of MHD mode dynamics. When combined with similar analysis of measurements from diagnostics that are sensitive to other disruption precursors, our analysis methods and results may contribute to the reliability, robustness and generalization of disruption warning schemes for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Control of Alfvén Eigenmodes by Externally Applied 3D Magnetic Perturbations

The suppression and excitation of Alfvén eigenmodes have been experimentally obtained, for the first time, by means of externally applied 3D perturbative fields with different spatial spectra in a tokamak plasma. The applied perturbation causes an internal fast-ion redistribution that modifies the phase-space gradients responsible for driving the modes, determining, ultimately their existence. Hybrid kinetic-magnetohydrodynamic simulations reveal an edge resonant transport layer activated by the 3D perturbative field as the responsible mechanism for the fast-ion redistribution. Finally, the results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multi-machine validation of plasma initiation modelling and prospects for future devices: Predicting plasma initiation using only hardware design and control room input data

This paper reports on the generic prediction capability of full electromagnetic plasma initiation modelling with DYON, which was carried out for the first time in fusion research by the joint modelling of the International Tokamak Physics Activity—Integrating Operation Scenario group. The following devices were included in the experiment database: VEST (spherical torus, copper coils, Stainless steel wall, R/a = 0.3 m/0.2 m, V v = 3.7 m 3 ), MAST-U (spherical torus, copper coils, C wall, R/a = 0.7 m/0.5 m, V v = 55 m 3 ), EAST (conventional tokamak, superconducting coils, metallic wall, R/a = 1.85 m/0.5 m, V v = 38 m 3 ), DIII-D (conventional tokamak, copper coils, C wall, R/a = 1.67 m/0.65 m, V v = 35 m 3 ), and KSTAR (conventional tokamak, superconducting coils, C wall, R/a = 1.8 m/0.5 m, V v = 55 m 3 ). Despite the different hardware features of the devices, the required operating spaces of the loop voltage induction and prefill gas pressure for inductive plasma initiation in each device were successfully reproduced by the predictive simulations with DYON using only the individual hardware design and the control room input data for each discharge. This successful validation across multiple machines demonstrates that the full electromagnetic DYON modelling can capture the essential physics of inductive plasma initiation. The simulation settings commonly employed for all modelling and the modifications necessary to account for the discrepancies between individual devices are reported. Predictions for ITER based on the multi-machine validation indicate that a wide range of prefill gas pressures exists for the Townsend breakdown and the plasma burn-through (0.01–1.5 mPa).

DYON↗
DOE OSTI2021

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. Here we have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Compelling Need for a Mid-Scale Stellarator Facility

In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. A new mid-size stellarator is needed to retire risks and innovate towards a high performance, economically attractive, stellarator Fusion Pilot Plant. In this presentation we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new and modern mid-size stellarator research facility. A new mid-scale stellarator is needed to realize the potential predicted by a solid body of theory and simulation along with advances in computational tools for optimization and non-linear turbulence modeling. Notably, it is possible to combine the advantages of the stellarator (steady state, no current drive, no disruptions) with the good confinement regularly achieved in tokamaks. The top priorities for experimental work, and the motivation for a mid-scale stellarator experiment are: turbulence control, non-resonant divertor, MHD stability at large beta, confinement of fast particles, and coil simplification. A new mid-size quasi-symmetric stellarator, built as a user facility, would complement existing research at Wendelstein 7-X and LHD and strongly augment private industry. It would provide a program of innovative research, concept validation, theoretical advancement, and workforce development. Growing support and interest for stellarators by the fusion community and private industry affirms this rationale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

In-Field Testing of Components for Feedback and Control of the ITER Disruption Mitigation System

Here, the shattered pellet injection (SPI) method has been chosen as the disruption mitigation system (DMS) for ITER. To protect the device from plasma disruptions that cause damaging heat and electromagnetic loads, SPI is used to inject high-Z material into the plasma. The process of SPI utilizes cryogenic cooling to form solid pellets. Pellets are accelerated down a barrel and into an angled surface, causing the pellet to shatter prior to entering the tokamak chamber. For the DMS to function reliably, the 27 separate shattered pellet injectors planned for ITER must rely on many components to provide accurate feedback data and for control functions. Each component in the DMS is exposed to an elevated background magnetic field depending on its placement and proximity to the plasma chamber. A Helmholtz coil test stand that is operated at Oak Ridge National Laboratory was utilized to test the components in relevant background field levels to assess component performance. This paper details the test design and results for in-field component operation for a variety of components. This list includes the following components: two different network switches for camera connectivity, a VAT fast shutter valve intended to reduce the flow of SPI propellant gas into the torus, a solenoid control valve intended for use in the pellet formation process, pressure/vacuum switches to be used for feedback and control, a printed circuit board piezo pressure sensor to be used to measure breech pressure, and various relays for the high-voltage pulsed power supply used to drive the SPI propellant valve.

Disruption mitigation↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2020

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of edge localized mode suppression by static resonant magnetic perturbations in the DIII-D tokamak

According to a recent paper [Hu et al ., Phys. Plasmas 26, 120702 (2019)], mode penetration at the top of the pedestal is a necessary and sufficient condition for the suppression of edge localized modes (ELMs) by a resonant magnetic perturbation (RMPs) in an H-mode tokamak discharge. This paper employs asymptotic matching theory to model a particular DIII-D discharge in which ELMs were suppressed by an externally generated, static, n ¼ 2, RMP whose amplitude was modulated at a frequency of 1 Hz. It is demonstrated that the response of the plasma to the applied RMP, in the immediate vicinities of the rational (i.e., resonant) surfaces, is governed by nonlinear, rather than by linear, physics. This is the case because the magnetic island widths associated with driven reconnection exceed the linear layer widths, even in cases where driven reconnection is strongly suppressed by plasma rotation. The natural frequency at a given rational surface (i.e., the helical frequency at which the locally resonant component of the RMP would need to propagate in order to maximize driven reconnection) is found to be offset from the local E B frame in the ion diamagnetic direction. The size of the offset is mostly determined by neoclassical poloidal rotation. Finally, the predictions of a fully nonlinear plasma response model are found to be broadly consistent with the DIII-D experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗