Core Magnetic Fluctuation Measurements In a High-Temperature Plasma - - PowerPoint PPT Presentation

core magnetic fluctuation measurements in a high
SMART_READER_LITE
LIVE PREVIEW

Core Magnetic Fluctuation Measurements In a High-Temperature Plasma - - PowerPoint PPT Presentation

Core Magnetic Fluctuation Measurements In a High-Temperature Plasma by Faraday Rotation B.H. Deng, W.X. Ding, D.L. Brower University of California, Los Angeles J.K. Anderson, D. Craig, G. Fiksel, C.B. Forest, D. Holly, V. Mirnov, S.C. Prager,


slide-1
SLIDE 1

B.H. Deng, W.X. Ding, D.L. Brower University of California, Los Angeles J.K. Anderson, D. Craig, G. Fiksel, C.B. Forest, D. Holly, V. Mirnov, S.C. Prager, J.S. Sarff, V. Svidzinski, and the MST Group University of Wisconsin-Madison Annual APS Meeting, Division of Plasma Physics Albuquerque, NM October 26-31, 2003 This work is supported by the U.S. DoE.

Core Magnetic Fluctuation Measurements In a High-Temperature Plasma by Faraday Rotation

slide-2
SLIDE 2

1. High-speed, far-infrared laser polarimetry for Faraday rotation measurements on MST – ∆t ~ 1 µs, ∆φ ~ 1 mrad 2. Current achievements on MST – Equilibrium magnetic field and plasma current profile – Comprehensive study of magnetic field and current density fluctuations

associated with tearing mode Characteristics: Amplitude, Frequency and Wave-number Spectrum, Spatial Distribution Dynamics:

– Magnetic fluctuation driven charge transport in MST – Fluctuation reduction in confinement-improved PPCD plasmas – Broadband magnetic and density fluctuations 3. Upgrading diagnostic system to three-wave polarimeter- interferometer for simultaneous measurements of ne and B.

Highlights

˜ B ⇒ J(r) ˜ B ⇒ τ E ,τ p B

θ (r,t ) ⇒ Jφ (r,t )

and ˜ B ω,k

( ) ⇒ ˜

J

slide-3
SLIDE 3

R-wave ω1 L-wave ω2

ER = ER cos(ω1t − kRz) EL = EL cos(ω 2t − kLz)

kL − kR = ω c (NL − NR) = ω 2

peωce

cω 2 ~ neBz

j ~ (ER + E L)2 ~ ER EL cos[(ω1 −ω2)t − (kR − kL)z]+ ...

Plasmas

Diode mixer

Faraday rotation is obtained by phase measurement

Faraday Rotation Measurement Method

slide-4
SLIDE 4

3-Wave Polarimeter

Reference Mixer Plasma L.O. Beam

ω2

λ/4 Plate λ/2 Plate FIR Lasers

ω2 ω1 ω1 ω3

Polarizer Probe Beams Lens Signal Mixer Lens Beam Splitter

F e z n e

c n B dz c ndz Ψ= Φ=

∫ ∫

Dodel and Kunz, Infrared Physics 18,773-776 (1978). Rommers and Howard, Plasma Phys. Control. Fusion 38,1805-1816(1996).

slide-5
SLIDE 5

MST R0 = 1.50 m a = 0.52 m Ip = 400 kA ne ~ 1019 m-3 B0 ~ 4 kG

3-Wave Polarimeter-Interferometer System

slide-6
SLIDE 6

I. System resolution

spatial: 11 discrete (vertical) chords temporal: ~1 µs time response phase: interferometer 50 mrad (8 x 1011 cm-2) polarimeter 1 mrad (10 G)

II. Measurement capabilities

Interferometer: Polarimeter: Determine both equilibrium: and fluctuating quantities:

III. At present, polarimeter and interferometer operated independently. The system is currently being upgraded to include a third FIR laser, forming a triple laser system, to provide high-speed simultaneous measurements of plasma density and Faraday rotation.

dl B n dl n

e e ||

~ ~

∫ ∫

Ψ φ ne (r,t) and Bθ (r,t) ⇒ Jφ (r,t) ˜ n k,ω

( )

˜ B k,ω

( )

˜ J

φ

FIR Polarimeter-Interferometer System

slide-7
SLIDE 7

Magnetic Fluctuation and Current Profile Play a Central Role in Self-Organized Magnetic Confinement

Magnetic reconnection Particle transport Momentum transport Energy transport

∇J// (r) δ r B , δ r J

Current Profile Control

< E >//=ηJ// + dynamo

Dynamo

Magnetic Fluctuation Induced Electromotive Force

ET

slide-8
SLIDE 8

Ip = 400 kA, 11 channel data, shot 101033036 Faraday Rotation

  • 4

4

Faraday Rotation (deg.)

80 60 40 20

Time (ms)

  • 32 cm
  • 17
  • 2

13 28 43

s1010330036

  • 4

4

Faraday Rotation (deg.)

80 60 40 20

Time (ms)

  • 24 cm
  • 9

6 21 36

1.6 1.2 0.8 0.4 0.0

nedl

70 60 50 40 30 20 10

Time (ms)

ne7 ne8 ne9 ne10 ne11

1.6 1.2 0.8 0.4 0.0

nedl

70 60 50 40 30 20 10

Time (ms)

ne1 ne2 ne3 ne4 ne5 ne6

Density: shot 1010330115

Lanier et al., Phys. Rev. Lett. 85,2120(2000); Phys.Plasmas 8,3402(2001)

Polarimeter-Interferometer Time History

slide-9
SLIDE 9

2.0 1.5 1.0 0.5 0.0 Jφ (MA/m

2)

0.4 0.2 0.0 r (m) before after Functional Fit 2.0 1.5 1.0 0.5 0.0 0.4 0.2 0.0 r (m) before after MSTFIT 4 2

  • 2
  • 4

Ψ (degrees) 2.0 1.5 1.0 R (m) before 4 2

  • 2
  • 4

2.0 1.5 1.0 R (m) after

Functional Fit

Jφ= J(0)[1-(r/a)2]γ,

Ψ = ne ∫ Bzdz

Equilibrium reconstruction: MSTFIT code

Brower, Ding, et al., PRL 88,185005-1(2002)

J(0) decreases at sawtooth crash

Faraday Rotation and Current Density Profiles

slide-10
SLIDE 10

J// Flattens at Sawteeth Crash

Equilibrium reconstruction code - solves Grad-Shafranov equation, fitting all experimental data (external magnetics, MSE, Faraday Rotation, pressure) At crash, plasma relaxes toward the Taylor minimum energy state

slide-11
SLIDE 11

0.20 0.15 0.10 0.05 0.00

  • 0.05

q 0.4 0.3 0.2 0.1 R-Ro (m)

  • .25 ms

.25 ms

1/6 1/71/8 1/9 1/10

q-Profile

slide-12
SLIDE 12

1.0 0.8 0.6 0.4 0.2 0.0 Coherence 80 60 40 20 Frequency [kHz] coherence statistical noise 12 kHz, m=1

crash (m=1,n=6)

˜ Ψ = ˜ n ∫ B0dl + n0 ∫ ˜ B dl ˜ n ∫ B0dl < rms noise, ⇒ ˜ Ψ ≈ n0 ∫ ˜ B dl

˜ B ≈ 33 G ~ 1%

( )

Core Magnetic Fluctuation Measurements

slide-13
SLIDE 13

→ Fluctuations approximately 1% before crash → Amplitude increases at crash → J(0) decreases

2.5 2.0 1.5 1.0 0.5 0.0 J(0) (MA/m2)

  • 2
  • 1

1 2 Time (ms) 400 300 200 100 Br (G)

  • 2
  • 1

1 2

B

~ r

B0 ≈ 1%

magnetic fluctuations act to redistribute current density

Magnetic Fluctuations Increase at Sawtooth Crash and Current Profile Flattens

slide-14
SLIDE 14

100 80 60 40 20 |Br|2 [G2/Hz] 100 80 60 40 20 Frequency [KHz]

5 6 7 8

10

2 3 4 5 6 7 8

100 |Br|2 [G2]

60 40 20

  • 20

Toroidal mode number

B r

2

~ n-5/ 3

400 300 200 100

  • 100

Phase (deg.) 80 60 40 20 frequency (kHz) 80 60 40 20

  • 20

toroidal mode number: n

Dispersion Relation

kφ = ∆φ d = n R

Frequency and Mode Number Spectra

slide-15
SLIDE 15

Current Fluctuation

Ampere's Law : δ r B • d r l = µ0

L

δI Faraday Rotation Fluctuation : δΨ = cF n0

δ r B • d r l ≈ cFn δ r B • d r l

x z r

δ r B • d r l ≈

L

δBz

dz

[ ]

x1 −

δBz

dz

[ ]

x2

≈ µ0δIφ = δΨ

1 −δΨ 2

cFn

Loop between polarimeter chords is equivalent to a Rogowski coil measurement Plasma

X1 X2

Ding,et al. PRL (2003)

slide-16
SLIDE 16

profile can be experimentally determined

δjϕ(r r , t)

δI5 δI4 δI3 δI2 δI

1

δj5 δj4 δj3 δj2 δj

1

δjϕ J0 ~ 6%

  • 15
  • 10
  • 5

5 10 15

δjϕ [ A/cm2 ]

  • 40
  • 20

20 40 r [cm]

Current Density Fluctuation Measurement

slide-17
SLIDE 17

∇×δB=µ0δJ , ∇•δJ = 0

In principle

ψ x = cF n0 ˜ B

z

dl ψ x

m =ψ x exp − cF ˜

n

1

r B

0 • d

r l

δjφ = ˜ j

0rexp(−(r − r s

w )

2)cosθ

χ

2 =

(ψi −ψi

m)2

σi

2 i=1,11

0.20 0.15 0.10 0.05 0.00 δΨ [deg.] 40 20

  • 20

X [cm]

Faraday Fluct. Fitting

Determining Magnetic Fluctuation Profile

slide-18
SLIDE 18

˜ j

1

J0 ≈ 4.5% J0 = 2 MA / m

2

˜ B

r

B ≈ 1%

  • 80
  • 60
  • 40
  • 20

20 40 60 80 Magnetic Field [Gauss] 1.0 0.8 0.6 0.4 0.2 0.0 r /a br bθ b

z

(1,6 ) mode

100 80 60 40 20 Current Density [kA/m

2]

1.0 0.8 0.6 0.4 0.2 0.0 r/a

w = 8 cm

rs = 17 cm

Global Local

Magnetic Fluctuation Spatial Structure

slide-19
SLIDE 19

2.5 2.0 1.5

  • 2
  • 1

1 2 Time [ms] Mean Current drops 20% 120 100 80 60 40 20

  • 2
  • 1

1 2 Time [ms] (1,6) mode current fluctutaion

Current Density Fluctuation over a Sawtooth

slide-20
SLIDE 20

− me e2ne ∂ r J ∂t + r E + r v × r B − 1 nee r J × r B + ∇P

e

nee = η r J

<E>// +<δr v ×δ r B >// −<δ r J ×δ r B >// /nee+...=η

// <J>//

The generalized Ohm’s law:

r J = r J

0 + δ

r J , r B = r B

0 +δ

r B , r v = r v

0 + δr

v

Mean field dynamics MHD dynamo Hall dynamo

Fluctuation Induced Electromotive Force

slide-21
SLIDE 21

< δJ × δB >// nee = 1 nee BP (rs) B(rs) 1+ BT BP      

2

        < (1 r ∂ ∂r rbθ )br > ≈ 1 nee BP B 1+ ( BT BP )2       < δjφbr >

δj is correlated with magnetic probe array data to determine phase difference of a specific mode (e.g.m/n=1/6).

∇•δB = 0, ∇ ×δB = µ0δJ

b

θ (r s) ~ 0

from measurement

Parallel Components of Hall Dynamo

slide-22
SLIDE 22

ηJ0 ≈ 0.5V /m

60 40 20

  • 2
  • 1

1 2 Time [ms] E //

<δJx δB> // /ne

1.7 V/m 0.50 V/m

30 20 10 0.8 0.6 0.4 0.2 r/a

< δJxδB>// /nee

0.3 0.2 0.1 0.0

  • 0.1
  • 2
  • 1

1 2 Time [ms]

Hall Dynamo and Inductive Electric Field

Hall dynamo effect is important near resonant surface!

slide-23
SLIDE 23

Current dynamics Magnetic Fluctuation Hall Electromotive Force (anti-current) Induced Electric Field Generation of Magnetic Field

100 50 δjϕ [kA/m2]

Current Fluctutaion

2.5 2.0 J [MA/m2]

Mean Current

60 40 20 V/m

Hall Dynamo

60 40 20 V/m

Electric Field

800 700 <BT> [Gs]

  • 2
  • 1

1 2 Time [ms]

Toroidal Flux

Magnetic Dynamics During Sawtooth Crash

slide-24
SLIDE 24
  • 6
  • 4
  • 2

2 4 6 Faraday Rotation (deg.) 30 25 20 15 10 5 Time (ms)

  • 24 cm
  • 9

6 21 36

1010326159

  • 6
  • 4
  • 2

2 4 6 Faraday Rotation (deg.) 30 25 20 15 10 5 Time (ms) PPCD ON PPCD OFF

  • 32 cm
  • 17
  • 2

13 28 43

Dynamo (sawteeth) Suppressed

Fast-Polarimetry During PPCD: Dynamo Suppressed

slide-25
SLIDE 25

Pulsed Poloidal Current Drive (PPCD) Changes Current Profile

7x10

6

6 5 4 3 2 1 λ= J// / B 0.8 0.6 0.4 0.2 r/a

No PPCD t=13 ms t=18 ms λ =constant

Brower,et al.PRL (2002)

Taylor minimum energy state

slide-26
SLIDE 26

Reduction of Magnetic Turbulence during PPCD

Broadband magnetic turbulence is reduced

80 60 40 20 P(f) [Gs2/kHz] 80 60 40 20 f [kHz] standard 400ka ppcd 400ka magnetic turbulence

Tearing Modes

slide-27
SLIDE 27

Broadband magnetic fluctuations above noise level are observed up to 250 kHz Noise can be further reduced using correlation techniques Better signal to noise ratio can be achieved with higher laser power Increased bandwidth allows us to study magnetic turbulence in high temperature plasma core

Measurement of Broadband Magnetic Fluctuations

2 3 4 5 6 7 8 9

10

  • 6

2 3

Faraday Rotation Power Spectrum (Deg.

2/Hz)

250 200 150 100 50 Frequency (kHz) Noise Level x = 6 cm

Standard Plasmas Ip = 340 kA ne = 1x10

13 cm

  • 13
slide-28
SLIDE 28

Broadband density fluctuations significantly above noise level are

  • bserved beyond 500 kHz at discharges with Ip=220 kA, ne~1x1019 m-3

Power falling off as ~f –α, α = 2.9 +/- 0.1

Measurement of Broadband Density Fluctuations

10

  • 13

10-12 10

  • 11

10

  • 10

10

  • 9

10

  • 8

Density Power Spectrum

5 6 7 8

10

2 3 4 5 6 7 8

100

2 3 4 5

Frequency (kHz)

Chan.1 Chan.2 Chan.3 Chan.4 Chan.5 Chan.6 Noise1 Noise2 Noise3 Noise4 Noise5 Noise6

~ f

  • α

10

  • 13

10

  • 12

10

  • 11

10

  • 10

10

  • 9

10

  • 8

Density Power Spectrum

5 6 7 8

10

2 3 4 5 6 7 8

100

2 3 4 5

Frequency (kHz)

Chan.6 Chan.7 Chan.8 Chan.9 Chan.10 Chan.11 Noise6 Noise7 Noise8 Noise9 Noise10 Noise11

~ f

  • α
slide-29
SLIDE 29

A third FIR laser will be added to achieve a three-wave polari- interferometer configuration for simultaneous Faraday rotation and density measurements A new high-power CO2 pump laser has been added (120 W CW) Precision bore Pyrex tubes will replace the metallic corrugated wave- guides of the FIR lasers for better performance at the desired wavelength (432.5 µm) The lasers will be housed in a laser safety room, which will also provide extra noise shielding as well as better IF frequency control

Diagnostic System Upgrade

slide-30
SLIDE 30

Diagnostic Technique directly applicable to Tokamaks! A high-speed far-infrared (FIR) laser polarimeter for Faraday rotation measurements in MST has achieved a time response of ~ 1 µsec, and a sensitivity of ~1 mrad. Equilibrium and fluctuating magnetic field and associated plasma current profile and fluctuations have been measured in the MST plasma core. Both low and high frequency fluctuations have been measured. Comprehensive study of the dynamics of plasma current density and magnetic field profiles and fluctuations in the core of MST plasmas (both standard and enhanced confinement discharges) has yield much new understanding on the behaviors of RFP plasmas: tearing mode instabilities, dynamo effects, and charge transport Plasma density profiles and fluctuations have also been independently studied in MST with the interferometer configuration. To enhance the power of the diagnostic, the system is being upgraded to a three-wave polari-interferometer for simultaneous high-speed magnetic field and density measurements.

Summary

slide-31
SLIDE 31

Laser Safety and Noise Reduction Room

Picture of the future site of the 3-wave interferometer/polarimeter