21.04.2014 Views

Solutions to Chen's Plasma Physics

Solutions to Chen's Plasma Physics

Solutions to Chen's Plasma Physics

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.

Since B ∼ r −3 , we obtain<br />

So<br />

∇B = ∂B<br />

∂r ˆr = − 3 r 4 = −3B r<br />

v ∇B = 3E<br />

eBr =<br />

⇒<br />

3Er3<br />

eB 0 Rer 3 = 3Er2<br />

eB 0 Re<br />

3<br />

=<br />

| ∇B<br />

B | = 3 r<br />

3E(eV )r2<br />

B 0 R 3 e<br />

where I have used the fact that B = B 0 R 3 e/r 3 . Now we can plug in the numbers:<br />

(47)<br />

(48)<br />

v ∇B,e = 3 × 30 × 103 × (5 × 6 × 10 6 ) 2<br />

3 × 10 −5 × (6 × 10 6 ) 3 = 1.3 × 10 4 m/s ✷ (49)<br />

v ∇B,i = 3 × 1 × (5 × 6 × 106 ) 2<br />

3 × 10 −5 × (6 × 10 6 = 0.42 m/s ✷ (50)<br />

) 3<br />

b) The magnetic field is azimuthal, from north <strong>to</strong> south, i.e. the −ˆθ direction. The gradient of the<br />

magnetic field is clearly in the radial direction, so we have B × ∇B ∼ ˆr × ˆθ = ˆφ, which is eastward.<br />

This is for the electrons. The ions, which come without the minus sign <strong>to</strong> cancel the minus sign in<br />

equation 47, will go in - ˆφ, which is westward. ✷<br />

c) Well, it has <strong>to</strong> travel a distance L = 2π(5R e ), with the velocity v ∇B,e . So,<br />

T =<br />

L 2π × 5 × 6.4 × 106<br />

=<br />

v ∇B,e 1.3 × 10 4 ≈ 4.5 hours ✷ (51)<br />

d) The current density is given by J = nev, so, using the grad-B velocity in this expression, we get<br />

J = nev ∇B,e = 10 7 × 1.6 × 10 −19 × 1.3 × 10 4 = 2 × 10 −8 A/m 2 ✷ (52)<br />

2-10. A 20 − keV deuteron in a large mirror fusion device has a pitch angle θ = 45 o at<br />

the mid plane, where B = 0.7 T . Compute it’s Larmor radius.<br />

The Larmor radius is given by<br />

r L = mv ⊥<br />

(53)<br />

qB<br />

In natural units, a deuteron has m = 2 and q = +1. Furthermore, v ⊥ = vsin(θ). To find v, we<br />

convert the energy <strong>to</strong> velocity:<br />

√<br />

√<br />

r L = m 2E<br />

eB m sin(θ) = 2 × 1.6 × 10−27 2 × 20 × 10 3 × 1.6 × 10 −19<br />

1.6 × 10 −19 × 0.7 2 × 1.6 × 10 −27 sin(45 o ) = 0.3 m ✷ (54)<br />

2-12. A cosmic ray pro<strong>to</strong>n is trapped between two moving magnetic mirrors with<br />

R m = 5 and initially has W = 1 keV and v ⊥ = v || at the mid plane. Each mirror moves<br />

<strong>to</strong>ward the mid plane with a velocity v m = 10 km/sec (Fig. 2-10).<br />

a) Using the loss cone formula and the invariance of µ, find the energy <strong>to</strong> which the<br />

pro<strong>to</strong>n will be accelerated before it escapes.<br />

b) How long will take <strong>to</strong> reach that energy?<br />

1. Treat the mirrors as flat pis<strong>to</strong>ns and show that the velocity gained at each bounce<br />

is 2v m .<br />

2. Compute the number of bounces necessary.<br />

3. Compute the time T it takes <strong>to</strong> traverse L = 10 10 km that many times. Fac<strong>to</strong>r of<br />

two accuracy will suffice.<br />

a) The loss cone formula is sin 2 (θ m ) = 1/R m , where θ m is the angle of the magnetic mirror and<br />

R m is the mirror ratio B max /B 0 . We can also write down a formula for sin(θ m ) in terms of the<br />

parallel and perpendicular velocities:<br />

sin(θ m ) = 1 √<br />

Rm<br />

=<br />

v ⊥,f<br />

√<br />

v<br />

2<br />

⊥,f<br />

+ v 2 ||,f<br />

=<br />

1<br />

√<br />

1 + (<br />

v ||,f<br />

v ⊥,f<br />

) 2 (55)

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

Saved successfully!

Ooh no, something went wrong!