06.01.2015 Views

Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas

Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas

Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Magnetic</strong> fields <strong>and</strong> tokamak plasmas<br />

Alan Wootton<br />

2 ⎛ r ⎞<br />

b rms<br />

= b r ,rms<br />

= b θ ,rms<br />

= 0.5∑ b rmn r= rmn<br />

⎜ ⎟<br />

17.30<br />

m ,n ⎝ ⎠<br />

r mn<br />

−2(m +1)<br />

To go further <strong>and</strong> obtain a relationship between b θrms at r > a in terms of b rmn at r = r mn we must<br />

specify<br />

a given current density or safety factor profile q(r),<br />

the dependence of b rmn on m <strong>and</strong> n,<br />

locate all possible m <strong>and</strong> n pairs where q = m/n,<br />

specify upper <strong>and</strong> lower limits to either m or n.<br />

In general this must be done computationally. However, some progress can be made analytically.<br />

We can approximately solve Equation 17.30 <strong>for</strong> b θrms at r > a in terms of b rmn at r = r mn if we<br />

take b rmn independent of m <strong>and</strong> n over some range to be specified, <strong>and</strong> a q profile relevant to the<br />

2<br />

⎛ r<br />

plasma edge, <strong>for</strong> example q = q<br />

⎞<br />

a<br />

. Then we can relate the resonance location r<br />

⎝ a⎠<br />

mn to a given<br />

q,<br />

r mn<br />

a =<br />

⎛<br />

⎜<br />

⎝<br />

1<br />

q ⎞ 2 ⎛<br />

⎟<br />

m ⎞<br />

= ⎜ ⎟<br />

q a ⎠ ⎝ nq a ⎠<br />

1<br />

2<br />

17.31<br />

With ρ = r/a defining the measurement position <strong>and</strong> q = m/n, Equation 17.30 is written as<br />

⎛<br />

⎜<br />

⎝<br />

b rms<br />

b rmn<br />

⎞<br />

⎟<br />

⎠<br />

2<br />

ρ ≥1<br />

2(m+1)<br />

(m +1)<br />

⎛ 1⎞<br />

⎛ m ⎞<br />

= 0.5∑ ⎜ ⎟ ⎜ ⎟<br />

17.32<br />

m ,n ⎝ ρ⎠<br />

⎝ nq a ⎠<br />

ρ ≥ 1 is a constant, namely the location of the magnetic pick-up coils used to measure b rms .<br />

We now replace the summation by integrals. The space we are considering is shown in Figure<br />

17.4. The first integral, ∫dn, is from m/q a to m, allowing all n values in between. The second<br />

integral, ∫dm, is from m 1 to m 2 , with m 2 >> m 1 . In Figure 17.4 m 1 = 6 <strong>and</strong> m 2 = 12.<br />

Un<strong>for</strong>tunately our chosen profile allows q ≤ 1 (when r mn /a ≤ 1/√q a ). To restrict n ≥ 1 we should<br />

choose m 1 ≥ q a ≈ 3. The result of the integration is<br />

b rms<br />

b rmn<br />

=<br />

⎡ 1 ⎛ 1<br />

4q a<br />

ln( ρ)<br />

ρ 2( m 1 +1)<br />

− 1<br />

⎢ ⎜<br />

⎢ ⎝<br />

⎢<br />

1 ⎛<br />

⎢ −<br />

⎜<br />

⎣ ⎢ 2 ln( q a )+ 2 ln( ρ)<br />

⎝<br />

( )<br />

ρ 2( m 2 +1)<br />

⎞<br />

⎟<br />

⎠<br />

1<br />

( m<br />

q 1 +1 )<br />

a<br />

ρ 2 m 1 +1<br />

( ) − 1<br />

( m 2+1<br />

q ) a<br />

ρ 2( m 2+1)<br />

⎤<br />

⎥<br />

⎥<br />

⎞ ⎥<br />

⎟ ⎥<br />

⎠ ⎦ ⎥<br />

1<br />

2<br />

129

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!