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 />

where the total current I flows inside the volume V, the line element dl is along the direction of<br />

the total current I, <strong>and</strong> r is the distance from the line element to the point of interest. A useful<br />

example <strong>for</strong> us is the vector potential of a circular filament. This is used to represent windings<br />

(vertical field, shaping, ohmic heating) on the tokamak, <strong>and</strong> elements of the plasma current itself.<br />

P<br />

r<br />

dl<br />

z<br />

Y<br />

dl<br />

R 0<br />

r<br />

R 0<br />

−φ<br />

φ<br />

R<br />

X<br />

Figure 1.6. The geometry used to evaluate the vector potential of a circular filament.<br />

Consider a circular filament of radius R 0 , with current I in the φ direction.<br />

A φ must be<br />

independent of φ, so choose the point of interest P in the (X,z) plane of Figure 1.6, where φ = 0.<br />

Pairing equidistant elements dl shown in thickened lines at ±φ we see the resultant is normal to<br />

(R,z). There<strong>for</strong>e only consider the component dl φ of dl in the direction normal to the plane (R,z);<br />

dl φ = R 0 cos(φ)dφ. The radial distance r from the point P to the element is given by<br />

r 2 = z 2 + R 2 0<br />

+ R 2 − 2R 0<br />

Rcos( φ). Then<br />

π<br />

A φ<br />

= µ Idl φ<br />

4π ∫ = µI R 0<br />

cos(φ )dφ<br />

r 2π ∫ 1.16<br />

0 R 2 0<br />

+ R 2 + z 2 − 2R 0<br />

Rcos(φ)<br />

( ) 1 2<br />

Far from the loop (i.e. a small loop) we have r 0 = (R 2 +z 2 ) 1/2 >> R 0 , <strong>and</strong> the integral becomes<br />

A φ<br />

= µI<br />

2π<br />

≈ RR 2<br />

0µI<br />

4r 3<br />

R 0<br />

cos(φ) ⎛<br />

⎜ 1 + RR cos(φ)<br />

0<br />

⎞<br />

⎟ dφ<br />

2<br />

r 0<br />

⎝ r 0<br />

⎠<br />

1.17<br />

= µ ( M × r)<br />

4πr 3<br />

π<br />

∫<br />

0<br />

Here we have written the magnetic moment of the loop M = πR 0 2 I , directed upwards.<br />

If the loop is not small, then let φ = π + 2θ, so dφ = 2 dθ <strong>and</strong> cos(φ) = 2 sin 2 (θ) - 1, <strong>and</strong> we obtain<br />

13

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

Saved successfully!

Ooh no, something went wrong!