Advances in Global MHD Mode Stabilization Research on NSTX

Similar documents
Resistive Wall Mode Stabilization and Plasma Rotation Damping Considerations for Maintaining High Beta Plasma Discharges in NSTX

Global Mode Control and Stabilization for Disruption Avoidance in High-β NSTX Plasmas *

Characterization and Forecasting of Unstable Resistive Wall Modes in NSTX and NSTX-U *

KSTAR Equilibrium Operating Space and Projected Stabilization at High Normalized Beta

Resistive Wall Mode Control in DIII-D

RESISTIVE WALL MODE STABILIZATION RESEARCH ON DIII D STATUS AND RECENT RESULTS

Effect of Resonant and Non-resonant Magnetic Braking on Error Field Tolerance in High Beta Plasmas

DIII D. by M. Okabayashi. Presented at 20th IAEA Fusion Energy Conference Vilamoura, Portugal November 1st - 6th, 2004.

GA A26247 EFFECT OF RESONANT AND NONRESONANT MAGNETIC BRAKING ON ERROR FIELD TOLERANCE IN HIGH BETA PLASMAS

Analysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission

RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations

GA A27444 PROBING RESISTIVE WALL MODE STABILITY USING OFF-AXIS NBI

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations a

Resistive Wall Mode Observation and Control in ITER-Relevant Plasmas

The Effects of Noise and Time Delay on RWM Feedback System Performance

Plasma Stability in Tokamaks and Stellarators

STABILIZATION OF THE RESISTIVE WALL MODE IN DIII D BY PLASMA ROTATION AND MAGNETIC FEEDBACK

Supported by. Role of plasma edge in global stability and control*

Requirements for Active Resistive Wall Mode (RWM) Feedback Control

Extended Lumped Parameter Model of Resistive Wall Mode and The Effective Self-Inductance

Effects of Noise in Time Dependent RWM Feedback Simulations

Effect of an error field on the stability of the resistive wall mode

(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time

NIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU. General Atomics. Nimrod Project Review Meeting July 21 22, 1997

The Effect of Energetic Particles on Resistive Wall Mode Stability in MAST

GA A26242 COMPREHENSIVE CONTROL OF RESISTIVE WALL MODES IN DIII-D ADVANCED TOKAMAK PLASMAS

Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas

Effect of non-axisymmetric magnetic perturbations on divertor heat and particle flux profiles

Three Dimensional Effects in Tokamaks How Tokamaks Can Benefit From Stellarator Research

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Resistive wall mode stabilization by slow plasma rotation in DIII-D tokamak discharges with balanced neutral beam injection a

Disruption dynamics in NSTX. long-pulse discharges. Presented by J.E. Menard, PPPL. for the NSTX Research Team

Characterization of neo-classical tearing modes in high-performance I- mode plasmas with ICRF mode conversion flow drive on Alcator C-Mod

A simple model of the resistive wall mode in tokamaks

A New Resistive Response to 3-D Fields in Low Rotation H-modes

Comparison of Divertor Heat Flux Splitting by 3D Fields with Field Line Tracing Simulation in KSTAR

MHD. Jeff Freidberg MIT

DIII D. by F. Turco 1. New York, January 23 rd, 2015

GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIII-D

Performance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK

ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations

Study of chirping Toroidicity-induced Alfvén Eigenmodes in the National Spherical Torus Experiment

Kinetic calculation of the resistive wall mode and fishbone-like

RWM Control in FIRE and ITER

CSA$-9 k067s-~ STUDY OF THE RESISTIVE WALL MODE IN DIII-D GAMA22921 JULY 1998

D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706

HIGH PERFORMANCE EXPERIMENTS IN JT-60U REVERSED SHEAR DISCHARGES

Progressing Performance Tokamak Core Physics. Marco Wischmeier Max-Planck-Institut für Plasmaphysik Garching marco.wischmeier at ipp.mpg.

Supported by. The drift kinetic and rotational effects on determining and predicting the macroscopic MHD instability

Stabilization of sawteeth in tokamaks with toroidal flows

Non-inductive plasma startup and current profile modification in Pegasus spherical torus discharges

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science

Magnetohydrodynamic stability of negative central magnetic shear, high pressure ( pol 1) toroidal equilibria

Highlights from (3D) Modeling of Tokamak Disruptions

Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor

Modeling of resistive wall mode and its control in experiments and ITER a

W.A. HOULBERG Oak Ridge National Lab., Oak Ridge, TN USA. M.C. ZARNSTORFF Princeton Plasma Plasma Physics Lab., Princeton, NJ USA

Princeton Plasma Physics Laboratory. Multi-mode analysis of RWM feedback with the NMA Code

(Motivation) Reactor tokamaks have to run without disruptions

Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing

Effect of ideal kink instabilities on particle redistribution

arxiv: v1 [physics.plasm-ph] 24 Nov 2017

Observation of Neo-Classical Ion Pinch in the Electric Tokamak*

MHD limits and plasma response in high beta hybrid operations in ASDEX Upgrade

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER

THE DIII D PROGRAM THREE-YEAR PLAN

Divertor Heat Flux Reduction and Detachment in NSTX

Current Drive Experiments in the HIT-II Spherical Tokamak

Interaction of scrape-off layer currents with magnetohydrodynamical instabilities in tokamak plasmas

DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH

Macroscopic Stability of High β N MAST Plasmas

Characterization of Edge Stability and Ohmic H-mode in the PEGASUS Toroidal Experiment

The Role of Dynamo Fluctuations in Anomalous Ion Heating, Mode Locking, and Flow Generation

W.M. Solomon 1. Presented at the 54th Annual Meeting of the APS Division of Plasma Physics Providence, RI October 29-November 2, 2012

Rotation and Neoclassical Ripple Transport in ITER

Non-Solenoidal Plasma Startup in

Influence of Beta, Shape and Rotation on the H-mode Pedestal Height

Effects of stellarator transform on sawtooth oscillations in CTH. Jeffrey Herfindal

Advanced Tokamak Research in JT-60U and JT-60SA

Control of linear modes in cylindrical resistive MHD with a resistive wall, plasma rotation, and complex gain

Initial Investigations of H-mode Edge Dynamics in the PEGASUS Toroidal Experiment

GA A23114 DEPENDENCE OF HEAT AND PARTICLE TRANSPORT ON THE RATIO OF THE ION AND ELECTRON TEMPERATURES

Stability of a plasma confined in a dipole field

EXD/P3-13. Dependences of the divertor and midplane heat flux widths in NSTX

ELMs and Constraints on the H-Mode Pedestal:

D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706

INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL

Plasma Response Control Using Advanced Feedback Techniques

Is the Troyon limit a beta limit?

MHD Analysis of the Tokamak Edge Pedestal in the Low Collisionality Regime Thoughts on the Physics of ELM-free QH and RMP Discharges

STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK

Control of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform. Jeffrey Herfindal

GA A THERMAL ION ORBIT LOSS AND RADIAL ELECTRIC FIELD IN DIII-D by J.S. degrassie, J.A. BOEDO, B.A. GRIERSON, and R.J.

- Effect of Stochastic Field and Resonant Magnetic Perturbation on Global MHD Fluctuation -

MAGNETIC FIELD ERRORS: RECONCILING MEASUREMENT, MODELING AND EMPIRICAL CORRECTIONS ON DIII D

High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows

On the physics of shear flows in 3D geometry

Flow and dynamo measurements in the HIST double pulsing CHI experiment

DIII D Research in Support of ITER

Transcription:

1 EX/5-1 Advances in Global MHD Mode Stabilization Research on NSTX S.A. Sabbagh 1), J.W. Berkery 1), R.E. Bell 2), J.M. Bialek 1), S.P. Gerhardt 2), J.E. Menard 2), R. Betti 3), D.A. Gates 2), B. Hu 3), O.N. Katsuro-Hopkins 1), B.P. LeBlanc 2), F.M. Levinton ), J. Manickam 2), K. Tritz 5), H. Yuh ) 1) Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY, USA 2) Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ, USA 3) Laboratory for Laser Energetics, University of Rochester, Rochester, NY, USA ) Nova Photonics, Princeton, NJ, USA 5) Johns Hopkins University, Baltimore, MD, USA e-mail contact of main author: sabbagh@pppl.gov Abstract. Stabilizing modes that limit plasma beta and reducing their deleterious effects on plasma rotation are key goals for efficient operation of a fusion reactor. Passive stabilization and active control of global kink/ballooning modes and resistive wall modes (RWM) have been demonstrated on NSTX and research now advances to understanding the stabilization physics and reliably maintaining the high beta plasma for confident extrapolation to ITER and CTF. Active n = 1 control experiments with an expanded sensor set, combined with low levels of n = 3 field phased to reduce error fields, reduced resonant field amplification and maintained plasma rotation, exceeded normalized beta = 6, and produced record discharge durations limited by magnet system constraints. Details of RWM active control show the mode being converted to a rotating kink that decays, or saturates leading to tearing modes. Discharges with rotation reduced by n = 3 magnetic braking suffer beta collapse at normalized beta =.2 approaching the no-wall limit, while normalized beta greater than 5.5 has been reached in these plasmas with n = 1 active control, in agreement with single-mode RWM theory. Advanced state-space control algorithms proposed for RWM control in ITER theoretically yield significant stabilization improvements. Values of relative phase between the measured n = 1 mode and the applied correction field that experimentally produce stability/instability agree with theory. Experimental mode destabilization occurs over a large range of plasma rotation, challenging the notion of a simple scalar critical rotation speed defining marginal stability. Stability calculations including kinetic modifications to ideal theory are applied to marginally stable experimental equilibria. Plasma rotation and collisionality variations are examined in the calculations. Intermediate rotation levels are less stable, consistent with experimental observations. Trapped ion resonances play a key role in this result. Recent experiments have demonstrated magnetic braking by non-resonant n = 2 fields. The observed rotation damping profile is broader than found for n = 3 fields. Increased ion temperature in the region of maximum braking torque increases the observed rate of rotation damping, consistent with theory. 1. Introduction Stabilizing large scale magnetohydrodynamic (MHD) modes that limit plasma beta and reducing their deleterious effect on plasma rotation, ω φ, are key goals for efficient operation of a fusion reactor. Demonstrating steady-state high beta conditions in ITER advanced operational scenarios by passive or active means is a key step toward the high beta operation of DEMO. High beta operation is also important to efficiently achieve neutron fluence goals in a nuclear component test facility (CTF) and is required for a demonstration reactor based on the spherical torus (ST) concept. Research on the National Spherical Torus Experiment, NSTX, has demonstrated both passive stabilization [1] and active control [2] of global kink/ballooning modes and resistive wall modes (RWM), accessing high β t 2μ <p> / B 2 = 39% and β N 1 8 <β t >ab /I p = 7.2. Present experiments and analysis now advance toward understanding and optimizing the performance and reliability of active control to maintain the high β N plasma, understanding global mode stabilization physics, and testing present theories of plasma rotation damping due to non-axisymmetric fields, which are all important for confident extrapolation to future toroidal magnetic fusion systems.

2 EX/5-1 2. Resistive Wall Mode Active Control at Various Plasma Rotation An expanded set of up to 8 resistive wall mode (RWM) sensors was included in recent active n = 1 control experiments to determine optimal sensor configurations and spatial phase, and to increase the reliability of maintaining high β N greater than the ideal MHD no-wall beta limit, β no-wall N. When this control was combined with low levels of n = 3 field phased to reduce error fields, resonant field amplification was reduced and ω φ was maintained, producing record discharge durations in NSTX (up to 1.8s at I p =.8 MA) limited by magnet constraints. Without feedback control, high β N plasmas are more susceptible to RWM instability leading to disruption, even at high levels of ω φ, and to mode locking at reduced rotation. Fig. 1 compares a high β N discharge using n = 1 active mode control with RWM B p sensor arrays above and below the plasma midplane and n = 3 DC error field correction (EFC) to a discharge without active control that is terminated by RWM instability at β N > β no-wall N. In the former case, ω φ is maintained, β N reaches and exceeds β no-wall N = computed Without feedback With feedback by the DCON stability code and 6 remains above this value with β N β reaching 6 (β N /β no-wall N N = 1.5). In β 2 N > β no-wall N β no-wall N = (DCON) the latter case, plasma disruption 2 occurs with ω φ above 8 khz near ω 129283 φ /2π (khz) the q = 2 surface, significantly 1 larger than the critical value for 12967 ω φ maintained RWM stabilization of 3.8 khz 12 ΔB n=1 pu,l (G) observed in experiments using n 8 = 3 braking for ω φ control. Discharges both with and without 5 reduced rotation can suffer beta I A (A) no-wall collapse at β N =.., -5 n = 3 correction n = 1 feedback on..2..6.8 1. 1.2 Fig. 1: High β N discharge using n = 1 active mode control with B p sensor arrays above/below the plasma midplane and n = 3 DC error field correction (solid line) vs. a discharge without active control (dashed line) that is terminated by an RWM instability at β N >β N no-wall. typical for H-mode profiles, while β N > 5.5 has been reached at reduced ω φ with n = 1 control. Feedback control of n = 1 RFA and unstable RWMs was used as a routine tool for the first time in 28 with more than 2 shots taken in over ten experiments. Two quantifiable measures of success are the probability of reaching long pulse duration and maintenance of high β N. Restricting the database to standard high performance H-mode operation, high β N discharges without feedback have a greater susceptibility to RWM-induced disruptions, significantly decreasing the probability of pulse length durations greater than.6s (Fig. 2a). With n = 3 EFC and n = 1 feedback control, I p flat-top durations of up to 1.6s were achieved. Fig. 2b shows values <β N > and plasma internal inductance <l i > averaged over the I p flat-top interval for shots with.2s duration (> 6 RWM growth times) and longer. Highest <β N > is maintained with n = 3 EFC and n = 1 feedback. Details of n = 1 RWM feedback show that control occurs by converting the non-rotating RWM into a global kink which either spins up and stabilizes, or further leads to energy dissipation by tearing modes. This process provides a link between RWM destabilization and tearing mode onset. The conversion from the RWM to a rotating kink occurs on the eddy current decay time of the wall, τ w ~ 3-5ms for n = 1 modes. In rare cases, which are the most

3 EX/5-1 Fig. 2: (a) Frequency distribution of discharge duration for plasmas with standard high performance H-mode plasma shaping with no n=3 EFC and n = 1 feedback control (blue), and with correction and feedback control (red); (b) operational stability space showing level of <β N > for I p duration >.2s without feedback control (blue), and with control (red). ΔB n=1 p I A (ka) ΔB n=odd φ n=1 Bp.5 -.5-1. -1.5 2 1 3 2 1 1-1 RFA RWM Mode rotation Co-NBI direction 12896 RFA reduced.66.61.61 t (s).618 Fig. 3: Detail of RWM evolution and feedback current response under active control. core (arb) edge (a) 12896 RWM spin-up.66.68.61.612.61.616 t (s) kink (b) edge core Fig. : Ultra-soft X-ray measurement of mode activity during n = 1 feedback control: (a) chord integrated signals spanning from the plasma core to edge show RWM onset, unlock, and spin-up transition to global kink, (b) onset of tearing mode 3 ms after RWM onset, indicated by radial phase inversion of signals (circled). illustrative, the rotating kink can saturate (Fig. 3). In this case, resonant field amplification observed in the n = 1 amplitude of the RWM poloidal field sensors, ΔB n=1 p, precedes the onset of the RWM. As the mode grows to sufficiently high amplitude, increased drag by the toroidal plasma rotation on the mode causes the RWM to unlock and rotate (shown by the n = 1 RWM phase, Φ n=1 Bp ). Without feedback control, this signature leads to further RWM growth and disruption [1]. With feedback, the n = 1 control currents (representative current I A ) respond to the mode onset, the RWM unlocks from the wall and accelerates to frequency ~ ω φ, precluding the existence of the RWM which damps in Δt ~ τ w. The kink displacement is identified by ultrasoft X-ray emission (USXR) and is global (Fig. a). USXR emission shows a tearing mode displacement form about 3ms after the saturated kink forms (Fig. b). RWM onset time + 3 ms Δt = 1 ms Specific n = 1 control experiments were run to examine the effect of control system delay on RWM control in ITERrelevant low rotation plasmas. Control system

EX/5-1 response was slowed by applying a 75ms smoothing filter to the requested control field currents to preclude response to RWM typical growth times of 3-5ms. With the smoothing filter applied, mode growth led to disruption at β N > β N no-wall as ω φ was lowered, showing that slow EFC alone is insufficient for stabilization at low ω φ. With no filter, low global rotation, with ω φ near or below DIII-D balanced NBI levels [3], was produced at high β N = 5.3. Experimental active control performance at reduced plasma rotation was reproduced by the VALEN-3D code with an upgraded system model using the actual off-midplane sensor positions and applied control field compensation analogous to that used in the experimental feedback system. 1 Experimental β N reached Experiment (b) Discharges with ω φ γ (1/s) 1 3 1 2 1 1 1 1-1 (control off) (control on) (a) DCON no-wall limit active control passive active feedback (reached) withwall limit.5 5 5.5 6 6.5 7 7.5 β N adv. controller γ (1/s) reduced by n = 3 magnetic braking suffer beta collapse at β N =.2 while approaching the n = 1 ideal MHD beta limit, β no-wall N, as computed by DCON using kinetic/mse constrained plasma reconstructions (Fig. 5a). In contrast, β N > 5.5 has been reached in these plasmas with n = 1 active control on. This agrees with singlemode RWM theory using proportional gain. Advanced state-space control algorithms proposed for RWM control in ITER [] theoretically show potentially large improvements in stabilized β N < 6.7. The relative phase, Δφ f, between the measured n = 1 mode and the applied correction field theoretically determines if the mode is stabilized or is driven unstable. Values of Δφ f experimentally producing stability/instability agree with theory (Fig. 5b). 3. RWM Passive Stabilization Physics 1 8 6 2 unstable 5 mode 225 29 315 rotating mode locked passive unstable 6 12 18 2 3 36 Δφ f (deg) Passive stabilization of the resistive wall mode by plasma rotation was postulated in seminal theoretical RWM papers [5], however, a comprehensive physics model that unifies experimental RWM stabilization by rotation remains elusive. RWM passive stabilization has been investigated on NSTX, with results summarized in Ref. [6]. While simple RWM models with viscous energy dissipation [7] can qualitatively describe RWM dynamics, simple threshold models of the critical plasma rotation frequency, Ω crit, or disruption due to loss of torque balance by resonant fields [8] do not describe RWM marginal stability in NSTX which is more profound and shown to be at least related to the plasma rotation profile [6]. Recent experiments continue to show this greater complexity. The unstable plasma shown in Fig. 1 becomes unstable at relatively high ω φ > 8 khz near q = 2, while plasmas with ω φ lowered by n = 3 non-resonant magnetic braking reach marginal stability at ω φ = 3.8 khz. More recently, n = 3 magnetic braking has produced marginally stable plasmas that have ω φ = at the q = 2 surface. A more comprehensive stability model by Hu and Betti including kinetic stable Fig. 5: (a) Experimental β N reached with n = 1 active control with ω φ reduced by n = 3 magnetic braking compared to theoretical feedback performance, and computed feedback improvement with an advanced state-space controller; (b) Experimental feedback stability versus feedback phase compared to theory.

5 EX/5-1 modifications to ideal MHD theory [9] has been used evaluate the stability of NSTX plasmas with initial success. Modifications to Ideal Stability by Kinetic effects are computed with the MISK code, including the effect of trapped and circulating ions, trapped electrons, and Alfven damping [1]. The kinetic components of the perturbed pressure lead to a potential energy functional δw K. The modified RWM growth rate normalized to τ w is γτ w = ( δw + δwk )/( δwb + δwk ), where δw is δw computed with no stabilizing conducting structure and δw b is computed with a model of the experimental stabilizing conducting structure. The calculation involves an integration over energy of a frequency resonance term (Eq. 5 of Ref. 1). For the trapped ion component, ( ˆ ε 3 ) ω ω ω ω γ 5 * N + δ 2 * T + E i ˆ ε ˆ 2 W ε e d ˆ K ε ω D lωb iν eff ωe ω iγ, (1) + + 1.x1 where the mode frequency is ω, and εˆ 5 is the particle energy normalized to the ion temperature. The other six.5 frequencies are compared, as profiles of poloidal flux Ψ, for an RWM marginally. stable plasma in Fig. 6 where ω *N and ω *T are the density and temperature -.5 ω E ω b gradient components of the ion ω *T <ω D > diamagnetic frequency, ω E is the E B ω *N ν eff [s -1 ] frequency, ω D is the bounce-averaged precession drift frequency, ω b is the..2..6.8 1. bounce frequency, and ν eff is the Ψ/Ψ a effective collision frequency. Plasma Fig. 6: Frequency profiles for kinetic stabilization of toroidal rotation enters through ion force RWM marginally stable plasma 12183 at.75s. balance ω E = ω φ ω *N ω *T. The The collision frequency profile shown is for thermal resonances in Eq. 1 allow a more ions, bounce and diamagnetic frequency profiles are complex dependence of stability on the for thermal ions and zero pitch angle. rotation profile. In the plasma shown in.3 Fig. 6, n = 3 braking produces a peaked -.6 ω / φ /ω exp 2. φ.2 ω φ profile, with sufficiently small ω φ <. -..6 (ω *N + ω *T ) allowing increased δw K by.8 stabilizing resonance with the trapped.2 1. -.2 1.2 ion precession drift in the outer portion 1. γτ 1.6 w of the plasma (roughly q > 2). Variation =. 1.8 2. of ω φ from an experimental RWM.1 1. marginally stable equilibrium.2 reconstruction shows reduced stability at intermediate and low plasma rotation, unstable experiment rather than a simple critical rotation. threshold. This is shown by varying the ω (rad/s) ω [rad/s] -1.x1 5 Im(δW K ) Im( W K )..1.2.3. Re( W K ) Re(δW K ) Fig. 7: Effect of plasma rotation on n = 1 RWM stability for plasma 12183 at.75s. The markers indicate self-similarly varied ω φ profiles. experimental rotation profile ω exp φ and computing γτ w for the variations. A stability diagram is produced in which contours of constant γτ w are shown on a plot of Im(δW K ) vs. Re(δW K ) (Fig. 7). The region of γτ w > which is theoretically unstable to the RWM is highlighted. The rotation profile is varied self-similarly from.2 to 2 times the experimental profile. Stability decreases

6 EX/5-1 as plasma rotation decreases from ω φ /ω exp φ = 2. to the value 1., which represents the measured experimental profile. This point is experimentally observed to be on the verge of instability and is computed to be close to marginal (γτ w > -.1). As ω φ is reduced further, the plasma becomes more stable between 1. < ω φ /ω exp φ <.6. For low rotation (<.ω exp φ ), the plasma is predicted to be unstable. This is consistent with the general NSTX experimental no-wall observation that very low global rotation at β N > β N is typically unattainable without active n = 1 control. Further physical insight is attained by examining the components of δw K as rotation is varied. In Fig. 8 the real and imaginary parts of δw K are broken into the trapped ion, trapped electron, circulating ion and Alfven layer contributions. For ω φ /ω exp φ from.2 to.6, stability increases as the real and imaginary trapped ion components increase. From.6 to.8, the real part increases while the imaginary part decreases, leading to the turn in Fig. 7 back towards instability. Finally for ω φ > ω exp φ, the trapped ion contributions are nearly constant, but the increase in the circulating ion component leads to increased stability. Im(δW K ) Im( W K ) 5x1-2 K ) Re(δW K ) 3 2 1-1.3.2.1.2. γτ w.6 -.6 -. -.2 w =..2 unstable.8 1. 1.2 ω φ /ω exp φ.2 Trapped Ion Trapped Electron Circulating Ion Alfven Layer Re(dW_K) (Total) 1. 1.6 (TI) (CI) (TE) (AL) 1.8 2. ω φ /ω exp φ /.2..6.8 1. 1.2 1. 1.6 1.8 2. 5x1-2 K ) Im(δW K )...1.2.3. Re( W K ) Re(δW K ) 2. 3 2 1-1.2..6.8 Trapped Ion Trapped Electron Circulating Ion Alfven Layer Im(dW_K) (Total) 1. 1.2 ω φ /ω exp φ Fig. 8: Components of δw K for plasma 12183 at.75s vs. scaled rotation profile. Fig. 9: Effect of plasma collisionality on n = 1 RWM stability for plasma 12183 at.75s. Solid circles/dashed line indicates ω φ variation at increased v, while open squares/solid line show the ω φ variation at decreased v. 1. (TI) (TE) (CI) (AL) 1.6 1.8 (CI) (TI) (TE) (AL) Variation of plasma collisionality, ν, from the experimental equilibrium alters the dependence of stability on rotation. Fig. 9 shows analogous trajectories to Fig. 7, but with electron and ion temperatures halved/doubled while the densities are doubled/halved, producing collisionalities increased/ decreased by a factor of 5.7 while maintaining constant pressure. Increased ν simplifies the dependence of stability on ω φ, making the trajectory appear more like it has a single critical rotation profile for stability. This might be expected, as higher ν decreases the relative importance of the kinetic resonances. Lower ν produces the opposite effect. As ω φ is decreased, there is a broader reversal of γτ w from decreasing to increasing stability, marginal stability occurs at higher plasma rotation than in the actual experiment, and the trajectory shows stability at low ω φ. 2.

7 EX/5-1. Plasma Rotation Alteration by Even Parity Non-resonant Fields and T i Dependence Physics understanding of non-axisymmetric field-induced plasma viscosity is important to ITER, especially if magnetic ELM mitigation techniques will be used, and in CTF if ω φ profile alteration is desired for global mode control. Braking torque due to odd parity applied fields has been observed and quantitatively compared to neoclassical toroidal viscosity (NTV) theory in NSTX [11]. Recent experiments also demonstrate non-resonant braking by an n = 2 field configuration which has a significant n = component (66% of the n = 2 amplitude) as well as n = 8 and 1 components (each 15% of the n = 2 amplitude). The observed rotation damping profile is broader than for n = 3 fields (Fig. 1), which is theoretically expected due to the broader field spectrum and reduced radial falloff of the n = 2 field. ω φ (khz) A significant aspect of the NTV theory for ITER and a CTF is the strong increase in NTV torque that occurs as the ion collisionality drops below the trapped particle bounce frequency, ν *i < 1. In this ion collisionality regime, the nonambipolar flux that leads to NTV is decreasingly perturbed by collisions, which increases the NTV torque, causing a strong dependence on ion temperature, τ NTV ~ δb 2 ε 1.5 p i /ν i ~ T i 5/2 [11]. Recent experiments show that increased T i increases the plasma rotation damping rate during non-resonant magnetic braking. The experiments utilized lithium evaporation to pre-condition portions of the divertor and plasma facing component carbon tiles [12] to generate a significant increase in T i ΔB n=odd (G) T i (kev) ω φ (khz) I A (ka) 3 2 1. -. -.8 broader braking region than n = 3 case 12788.9 1. 1.1 1.2 1.3 1. 1.5 1.6 1372 13722 no lithium.5s.55s.65s.75s lithium t (s) a) b) broad braking region, peak change (1.3 < < 1.35 121 1. 1.1 1.2 1.3 1. 1.5 1.6 n = 2 braking 2 R = 1.37m..3.2 no lithium lithium.1 5-5 -1..5.6.7.525s.535s.55s.555s Fig. 1: Comparison of non-resonant magnetic braking by (a) n = 2 field configuration to (b) n = 3 configuration. Fig. 11: Comparison of n = 2 magnetic braking in plasmas with varied ion temperature. in the region of peak n = 2 magnetic braking. Fig. 11 compares two plasmas with equal magnitude of applied n = 2 braking currents, and with/without lithium wall preparation. The plasma with Li preparation shows higher initial T i (shown in the region of maximum change of ω φ, R = 1.37m) at the start of the magnetic braking pulse, with a ratio of T i between the two shots of T i (Li) / T i (no Li) = 1.32. This yields a theoretical increase in τ NTV of a factor of two in the plasma with higher T i and is consistent with the measured increase in dω φ /dt observed in the plasma. Fig. 12a and Fig. 12b show the plasma rotation profile evolution for these shots. A plasma using Li wall preparation was then run with the braking current reduced by 25% (~ 5% less braking torque) which allowed saturation of the ω φ profile and maintained RWM stability through the braking pulse (Fig. 12c).

8 EX/5-1 3 5. Summary and Discussion.615.615.615 With initial passive.625 2.625.625 (braking and active global 2 saturates) mode stabilization established on lithium 1 6A peak NSTX, research is no lithium lithium 1 coil current 8A peak 8A peak advancing toward (~ 5% less coil current coil current torque) understanding the 1372 13722 13723 physics of mode.9 1. 1.1 1.2 1.3 1. 1.5 1.6.9 1. 1.1 1.2 1.3 1. 1.5 1.6 1. 1.1 1.2 1.3 1. 1.5 stabilization and to maintaining high Fig. 12: Evolution of ω φ in discharges slowed by n = 2 magnetic braking. beta plasmas with high reliability. RWM instability occurs over a wide range of ω φ, and the combination of n = 3 static EFC and n = 1 feedback control greatly increases the probability of sustaining high β N over long pulses, producing record pulse duration in the device, limited by magnet system constraints. RWM dynamics during successful n = 1 feedback shows the mode convert to a rotating kink which damps, or saturates typically leading to tearing modes. Fast n = 1 feedback control response ~ O(τ w ) is necessary to avoid high β N plasma disruption at low ω φ. Simple models of a critical rotation threshold for RWM stabilization are too restrictive to explain the observed plasma stability vs. ω φ. A physics model including kinetic modifications to ideal stability shows regions of reduced passive stability at intermediate ω φ as well as low ω φ, which agree with experimental observations. Plasma collisionality alters the range of ω φ displaying reduced stability, which may unify the physics understanding of RWM stability between present STs and tokamaks. Plasma rotation alteration has been demonstrated by magnetic braking using non-resonant n = 2 fields. The observed rotation damping profile is broader than that found for n = 3 fields. Increasing ion temperature in the region of maximum braking torque increases the observed rate of rotation damping, consistent with the dominant scaling of the non-resonant NTV torque τ NTV ~ δb 2 ε 1.5 p i /ν i ~ T 5/2 i in the NSTX range of collisionality ν *i < 1. Research at lower ν i is important to determine at what collisionality the 1/ν i scaling will saturate and if τ NTV will eventually decrease as ν i is further reduced. ω φ (khz) 3 a) b).575.575 c).585.595.65.585.595.65 This research was supported by the U.S. Department of Energy under contracts DE-FG2-99ER552 and DE-AC2-76CH373. [1] S.A. SABBAGH, A.C. SONTAG, J.M. BIALEK, et al., Nucl. Fusion 6 (26) 635. [2] S.A. SABBAGH, R.E. BELL, J.E. MENARD, et al., Phys. Rev. Lett. 97 (26) 5. [3] H. REIMERDES, A.M. GAROFALO, G.L. JACKSON, et al., Phys. Rev. Lett. 98 (27) 551. [] O. KATSURO-HOPKINS, J.M. BIALEK, D.A. MAURER, et. al., Nucl. Fusion 7 (27) 1157. [5] A. BONDESON, and D.J. WARD, Phys. Rev. Lett. 72 (199) 279. [6] A.C. SONTAG, S.A. SABBAGH, W. ZHU, et al., Nucl. Fusion 7 (27) 15. [7] R. FITZPATRICK. Phys. Plasmas 9 (22) 359. [8] A. GAROFALO, G.L. JACKSON, R.J. LAHAYE, et al., Nucl. Fusion 7 (27) 1121. [9] B. HU and R. BETTI, Phys. Rev. Lett., 93 (2) 152. [1] B. HU, R. BETTI, J. MANICKAM, Phys. Plasmas 12 (25) 5731. [11] W. ZHU, S.A. SABBAGH, R.E. BELL, et al., Phys. Rev. Lett. 96 (26) 2252. [12] D.K. MANSFIELD, H.W KUGEL, R. MAINGI, et al. Transition to ELM-free improved H- mode by lithium deposition on NSTX graphite divertor surfaces, in proceedings of the 18th PSI meeting, Toledo, Spain May 26-3, 28, paper O-28; accepted for publication in Jour. Nucl. Mater..575.585.595.65