INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL S.C. GUO* and M.S. CHU GENERAL ATOMICS *Permanent Address: Consorzio RFX, Padova, Italy **The authors acknowledge discussion with Dr. R. Fitzpatrick 328-99-1
Abstract Submitted for the DPP99 Meeting of The American Physical Society Sorting Category: 5.10(Theoretical) Interaction of an External Rotating Magnetic Field with the Plasma Tearing Mode Surrounded by a Resistive Wall 1 S.C. GUO, Consorzio RFX, Padova, Italy, M.S. CHU, General Atomics The effect of an externally rotating magnetic field on the plasma tearing mode surrounded by a resistive wall is studied. A pair of tearing mode evolution equations describing the magnetic energy and angular momentum balance across the magnetic island are used. The model is valid for both the RFP and the tokamak. The pair of equations is solved numerically to determine the equilibrium amplitude and phase of the tearing mode with respect to that of the external magnetic field and the phase stability of the combined system. When the external magnetic field amplitude is large, the tearing mode frequency is locked to that of the external field above minimum amplitude. Dependence of the critical unlocking amplitude and the phase stability on parameters relevant to present day experiments are obtained. In the opposite limit, when the amplitude is small, the external field is not sufficient to lock the tearing mode below a critical amplitude. Dependence of this critical locking amplitude on plasma characteristics and external wall distance is also obtained. Possible utilization of the external rotating field to stabilize the tearing mode is also discussed. 2 1 Supported by EURATOM and U.S. DOE Grant DE-AC03-95ER54309. 2 The authors acknowledge discussion with R. Fitzpatrick. X Prefer Oral Session Prefer Poster Session M.S. Chu chum@fusion.gat.com General Atomics Date printed: September 21, 1999 Electronic form version 1.4
OUTLINE Introduction Model Equations Steady state solutions and phase stability analysis Mode locked to the external rotating applied field; unlocking threshold Mode unlocked; locking threshold Intrinsically unstable modes and stable modes Summary and discussion 328-99-2
INTRODUCTION Key (MHD) issues in the operation of fusion devices: * Avoidance of locked mode induced by resistive wall and/or error field. * Active stabilization of low m resistive MHD modes. One of the possible approaches is to externally apply an rotating helical magnetic field, as some experiments have tested (DIII-D 1, JET 2, RFX 3...). The problem involves non-linear interaction of plasma modes with external magnetic perturbation and/or resistive wall. Various theoretical studies have been done on this field. 4 In this work, we derive and solve the tearing mode evolution equations (equilibrium solution), which involve the effects of both resistive wall and rotating external magnetic field. The model equation can be applied to both tokamaks and RFPs. 328-99-3
MODEL EQUATIONS Cylindrically symmetric plasma (r = a) surrounded by a resistive wall r w = b, and a perfect conducting shell r sh = c Plasma Resistive Wall Shell Vacuum Vacuum r s a b c r The equation governing the magnetic perturbation due to a tearing mode (in outer region) d f(r) dψ dr dr g(r)ψ =0 where ψ=rb r 328-99-4
Boundary Conditions MODEL EQUATIONS Thin shell approximation for the resistive wall δ b b << ωτ b < b δ b Externally applied rotating magnetic field ψ ex at with frequency ω ex r ex = c; Assumptions The plasma is assumed to rotate only in the toroidal direction 1 ω p 2 dω p dt << 1 Validity of Rutherford equation: w/δ VR (Fitzpatrick, 1999; Riconda et al., 1999) 1, Constant-ψ 328-99-5
TEARING MODE EVOLUTION EQUATIONS Without loss of generality, ψ m,n can be separated into three parts: Ψ s Ψ b Ψ ex ψ m,n =ψ s +ψ b +ψ ex =Ψ s ˆψ s (r) +Ψ b ˆψ b (r) +Ψ ex ˆψ ex (r) r s T m,n EM = 2π2 R 0 n µ 0 m 2 + n 2 2 ε I m s = s r s = r dψ m,n dr * ψ m,n ψ * m,n + dψ * m,n dr r dψ m,n dr ψ m,n rs r s + + r s r s 2 Ψ s Ψ s * Evolution equations d Ψ 4Iτ R dt 4 Ψ d ω dt = plasma + wall + ex = T plasma + T wall + T ex ω = ω 0 ω p 328-99-6
TEARING MODE EVOLUTION EQUATIONS plasma = b (0) 1 Ψ 0 Ψ 0 wall = E sb E bs E bb ω 2 p τb 2 + Ebb 2 ex = E sb E bex ω 2 ex τb 2 + Ebb 2 Ψ ex Ψ s cos( ϕ θ) T plasma = Gν τν ω 0 ω p ( ) T wall = Q E sb E bs ω p τ b E 2 bb +ωp 2 τb 2 T ex = Q E sb E bex E 2 bb +ω 2 2 ex τ Ψ s Ψ ex sin( ϕ θ) b E ij = r d ˆψ j dr i, j = s, b, ex i t φ = ω 0 ( ex ω p )d t +φ 0ex φ 0p θ=tan 1 ω ex τ b E bb 328-99-7
TEARING MODE EVOLUTION EQUATIONS (a) Resistive wall stabilization Ψ s =Ψ 0 1 E sb E bs E bb ω 2 p τ 2 2 b + E bb 1 b (0) τ b Ψ s =Ψ 0 ; τ b 0 Ψ s0 =Ψ 0 2 1 E sb E 2 bs E bb b (0) 2 =Ψ c (0) 0 b (0) 0.012 0.01 ψ 0 =0.00025 ψ 0 =0.00075 0.008 ψ s 0.006 0.004 0.002 0-20 -10 0 1 0 2 0 3 0 4 0 ω p τ b 328-99-8
PHASE STABILITY ANALYSIS Increment of δω p, δφ Unstable Steady state ω p, φ Perturbation δω p, δφ Variation of ψ s due to δω p, δφ Variation of torque balance Stable Variation of T EM due to δω p (or) δφ Decrement of δω p, δφ 328-99-9
STEADY STATE SOLUTION I: FOR Ψ ex = 0 plasma + wall = 0 T plasma + T wall = 0 10/31/99 01:43:36 0.150E+01 0.165E+02 0.315E+02 0.465E+02 0.615E+02 0.615E+02 0.465E+02 For ω 0 < ω c, continuous spectrum 0.315E+02 ω ψ < ω c, bifurcated spectrum 0.165E+02 0.150E+01 328-99-11
STEADY STATE SOLUTION II: LOCKED UNSTABLE MODE (a) Bifurcated solution ω p = ω ex (b) Continuous solution 10/31/99 fexv (psi, cm) 03:09:32 10/31/99 fexv (psi, cm) 06:50:32 0.125E+03 0.906E+03 0.169E+02 0.247E+02 0.325E+02 0.403E+02 0.481E+02 0.559E+02 0.638E+02 0.716E+02 0.875E+02 0.125E+03 0.263E+02 0.513E+02 0.763E+02 0.101E+01 0.297E+02 0.759E+02 0.231E+02 0.644E+02 0.528E+02 0.165E+02 0.413E+02 0.990E+01 0.297E+02 0.330E+01 0.182E+02 0.330E+01 0.660E+01 0.495E+01 0.990E+01 0.165E+02 0.165E+02 [m = 1, n = 0, ψ 0 = 0.25 10 3, b (0) = 0.3] 328-99-12
STEADY STATE SOLUTION II: LOCKED UNSTABLE MODE Stability Boundary [for case (a) lower frequency solution] 4 10-3 Unlock Threshold 0.6 Phase (on Threshold) ψ 3 10-3 ex 0.4 2 10-3 0.2 ( φ θ)/π 0 1 10-3 -0.2-0.4 0 10 0-15 -10-5 0 5 1 0 1 5 ω τb -0.6-15 -10-5 0 5 10 15 ω τb ϕ 0 θ π 2, ω p ~ ω p L ~0 ϕ 0 θ < ~ π 4, ω p ~ ω p H ~ ω 0 328-99-13
STEADY STATE SOLUTION II: LOCKED UNSTABLE MODE Stability Boundary [for case (a) lower frequency solution] 0.008 0.5 0.007 ψ ex 0.006 0.005-9.9 9.9 0.004 0.003 3.3 ( φ θ)/π 0 0.002 0.0 0.001 ωτ = 1.41 0 1 2 3 4 5 6 7 8 ψs (10-3 ) -0.5 2 3 4 5 6 7 8 (10-3 ) ψ s 328-99-14
(a) STEADY STATE SOLUTION III: UNLOCKED UNSTABLE MODE Analysis ( ω p ω ex ) Ψ s =Ψ s0 +δψ s Ψ s0 =Ψ 0 1+ wall b (0) 2 δψ s = 1 P Ψ ex 1 2Iτ R Ψ s0 1/2 cos( φ θ β) α 2 2 + (ω ex ω p ) δψ s oscillates in time, so the torque T ex also oscillates in time Assumption: Plasma is sufficiently viscous that it responds only to the steady components T ex = QP2 4Iτ R Ψ ex 2 Ψ s0 1/2 (ω ex ω p ) α 2 + (ω ex ω p ) 2 328-99-15
STEADY STATE SOLUTION III: UNLOCKED UNSTABLE MODE (b) Results (for ω H p > ω ex > w L p ) 10/31/99 psi_ext 08:20:29 0.125E+03 0.294E+02 0.575E+02 0.856E+02 0.114E+01 0.734E+02 The unlocked ω p solution exists when 0.549E+02 ω p < w L c, or ω p < w H c 0.363E+02 0.177E+02 0.825E+00 328-99-16
STEADY STATE SOLUTION III: UNLOCKED UNSTABLE MODE (c) Locking threshold ψ exc and corresponding ψ sc (10-3 ) 35 30 25 ω c τ b =3.3 ω ex τ b =26.2 ψ exc 20 ψ 15 10 5 0 ψ sc 2 2.5 3 ψ 3.5 o 4 4.5(10 ) Locking threshold ψ exc is much larger than unlocking threshold For ψ 0 = 2.5 10 4, ψ exc (unlock) = 2.4 10 4 (ω p τ b = 3.3), ψ exc (lock) = 4.8 10 3 328-99-17
STEADY STATE SOLUTION IV: LOCKED STABLE MODE For intrinsically stable plasma, if the external field is rotating in the plasma frame, there will still be a torque acting on the resonant surface due to inertia and dissipative effects. This torque acts to bring the tearing mode to rotate together with ω ex. Once the locking happens, a full reconnection occurs. 10/31/99 fexv (psi, om) 08:37:17 0.156E+03 0.104E+02 0.191E+02 0.279E+02 0.367E+02 0.206E+02 0.113E+02 0.206E+01 0.722E+01 0.165E+02 328-99-18
STEADY STATE SOLUTION IV: LOCKED STABLE MODE For intrinsically stable plasma, if the external field is rotating in the plasma frame, there will still be a torque acting on the resonant surface due to inertia and dissipative effects. This torque acts to bring the tearing mode to rotate together with ω ex. Once the locking happens, a full reconnection occurs. 11/09/99 fexv (psi, om) 10:27:48 0.200E-03 0.120E-02 0.220E-02 0.320E-02 0.420E-02 0.520E-02 0.620E-02 0.372E+02 0.252E+02 0.132E+02 0.120E+01 0.108E+02 0.228E+02 0.348E+02 328-99-18A
SUMMARY AND DISCUSSION A pair of tearing mode evolution equations, which describe the magnetic energy and angular momentum balance across the magnetic island are derived Steady state solutions are obtained numerically: (a) The locking solution (ω p =ω ex ) provides the required ψ ex and φ 0 When the external field is reduced beyond a critical amplitude (unlocked threshold), the locking is lost (b) The unlocking solution (ω p ω ex ) describes how the mode frequency and amplitude are influenced by the external field When the external field is increased beyond a critical amplitude (locking threshold), the tearing mode will then become locked The resistive wall introduces an extra phase shift Discussion: Is it possible to force a transition from ω p L to ω p H by applying a rotating external magnetic field? 328-99-19