Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas
Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas
Magnetic Fields and Magnetic Diagnostics for Tokamak Plasmas
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