The Debye length, the plasma parameter, loss cone angle, mirror machines, Larmor radius, gyrofrequency, guiding center and gyro-period are all elements of basic plasma physics.
Below I show how example problems involving these parameters are solved. Most are taken from my book, Explorations in Mathematical Physics, Chapter 8,
1) A helium plasma is at a temperature of 10 6 K, and exhibits a number density of
n s = 10 10 /m3. Find the Debye length l D,s and the plasma parameter, L
Solution:
For helium plasma we have: T s= T He = 10 6 K
Debye length is:
l D =[kT / 4p n e e2 ] 1/2
Where k is Boltzmann’s constant and e is the electronic charge.
For this ion-species based Debye length we may write:
l D,s =[kT s / 4p Z s2 n s e2 ] 1/2
Where the subscript ‘s’ refers to the ion species. Here, for helium:
Z s= 2, e = 1.6 x 10 -19 C
n s = 10 10 /m3
l D,s =
[(1.38 x 10-23 (10 6 K) (8.85 x 10 -12 F/m) / (4p) 22 (10 10/m3) e2 ] ½
l D,s = 0.30 m
The plasma
parameter is defined:
L = n s l3 D
From the
values given (T =10 6 K, and n s = 10 10 /m3)
L = n s l3 D = ( 10 10 /m3) (0.30 m ) 3
L = 2.9 x
10 8
(2) Consider a plasma mirror machine of length 2L with a mirror ratio of 10 so that B(L) = B(-L) = 10 B(0). E.g.
Solution:
The mirror ratio is (Bmin / B max) = 10,
meaning that the induction strength at those end points will be ten times the
induction at the center point or apex of the magnetic loop or mirror machine.
We define what is called the “loss cone angle”:
sin (q L ) = ± Ö (Bmin
/ B max )
In the problem, Bmin = B(0) or the "zero
level" for the magnetic induction, say at position L = 0. This doesn't
mean the induction is zero at that point literally, however.
To do the problem, one must understand he's really being asked for the fraction
of electrons lost. A special condition obtains which applies to the angle - for
which the electrons will be TRAPPED only provided:
Θ (O) > (Θ)L
Thus, Θ (O) = (Θ)L
is said to be the "loss cone" of the system or machine.
If an isotropic particle distribution (in this case, electrons) is introduced
at a position L = 0, the fraction of particles that will be lost to the mirror
system is:
f(L) = 1/ 2 π ò o Θ(L) 2 π sin(Θ) d Θ
= f(L) = ò o Θ(L) sin(Θ) d Θ
= 1 - cos(Θ)L
From the problem, if N denotes the total of electrons released, then you will
have to find the fraction lost from :f(L) = N - [1 - B(0)/ B(L)]1/2
Bearing in mind, B(z=L) = B(z= -L)
Hence: f(L) = N - [1 - B(0)/ B(L)] 1/2
And,
B(z=L) = B(z= -L) = 10
(B(0))
This then yields:
f(L) = N - [1 - 1/ 10] 1/2
= N – [9/10] 1/2 = N – (0.9) 1/2
= N – 0.948
N > 1 the total released electrons, so the total fraction lost is:
f(L) = 1 – 0.948 = 0.052 (e.g. 5.2 % of
electrons lost)
(3) Make use of the revised force equation:
F z
= m ¶B z/
¶z
To show that: dm/ dt » 0
Solution:
Consider component of motion along Bz only then the revised force equation is:
F
z = m dvz/ dt = - m ¶B z/ ¶z
Multiply
by velocity vz on
the left and the equivalent dz/dt on right:
vz (m dvz/ dt) = - m dz/ dt (¶B z/ ¶z)
Where dz » ¶z
® m vz dvz/ dt
= - m ¶B z/ ¶t
Where ¶B z/ ¶t is a variation of Bz as seen by the particle.
Since the energy of the particle must be
conserved:
d/dt(½ m vz2 + ½ m v2⊥ ) = d/dt (m vz2 + mB) = 0
Thence:
m vz dvz/ dt
= d/dt (m vz2 ) = - m dB z/ dt
But: - m dB / dt +
d/dt (mB) = 0
Therefore: dm/ dt » 0
(4) If the perpendicular velocity component ( v⊥) is 105 m/s for the electron, find its Larmor radius and its gyro-period. Thence or otherwise obtain the gyration energy
Solution:
Larmor radius: r = m/ q [v⊥ / B] = v⊥/ (qB/ m e) = v⊥/ Ω e
r = (10 5 m/s) / 1.7 x 10 7 /s = 0.0056 m or: 0.56 cm
Gyro-period: T
= 2 p
/ Ω e
= 2 p / 1.7 x
10 7 /s = 3.5 x 10 -7
s
Gyration
energy E = m e
(v⊥)2/ 2 =
(9.1 x 10 -31 kg) (10 5 m/s) 2
/ 2
E = 4.5 x 10 -21 J
Or: in Electron volts: E= (4.5 x 10 -21 J )/ (1.6 x 10 -19 J/eV) = 0.028 eV
(5) A proton moves in a
uniform electric and magnetic field, with fields given by: E = 10 V/m (x^) and B = 0.0001 T (z^)where '^' denotes
vector direction.
a) Find the
gyrofrequency and the gyro-radius
b) Find the
proton's E X B drift speed
c)Find the gyration
speed v and compare to the drift speed
d)Find the
gyro-period, gyration energy and magnetic moment of the proton.
(a) The proton's gyro-frequency is just: Ω i = qB/ m i
Ω i =
qB/ m i = [(1.6
x 10 -19 C) (0.0001 T) ]/ 1.7 x 10 -27 kg
Then:
Ω i = 9.4 x 10 3 /s
The gyro radius r, is defined: r = v⊥ / (qB/m) = v⊥ / Ω
So we need to get v⊥ first. This can be obtained from the energy of gyration:
E = m(v⊥)2/ 2 = m m B where : m m = 8.5 x 10 -22 J/T
v⊥ = Ö (2 m m B/ m i ) =
Ö [2
(8.5 x 10 -22 J/T) (0.0001 T)/ 1.7
x 10 -27 kg] = 10 m/s
Then, the radius of gyration is: r = v⊥ / Ω i = 10 m/s/ 9.4 x 10 3 /s
And: r = 0.001 m = 0.1 cm or 1 mm
b) The proton's (E X B) drift speed is just the magnitude:
[E/B]^y = [ (10 V/m (x^)/ (0.0001 T (z^)) ] = 10 5 m/s ^y
c)The
gyration velocity is just: v⊥ = Ö (2 m m B/ m i ) = 10
m/s
d) The gyro-period is: 2 p / Ω i = 2 p /(9.4 x 10 3 /s ) = 6.6 x 10 -4 s
The gyration energy E = m m B where : m m = 8.5 x 10 -22 J/T
E = (8.5 x 10 -22 J/T) (0.0001 T) = 8.5 x 10 -26 J
The magnetic moment (m m ) of the proton for the problem was given but can also be worked out independently from:
E = m m B = m/2 (E/B) 2,
So: [E/B] = 10 5 m/s
And the magnetic moment of the proton is:
m m = m/2 (E/B) 2 B = m/2 (10 5 m/s ) 2 0.0001 T
= (1.7
x 10 -27 kg) (10 5 m/s ) 2 0.0001
T/ 2 = 8.5 x 10 -22 J/T
(6) Find the guiding center positions for the electron referenced in Problem 4 if t = T/4.
Solution:
The electron (or ion) moves
along the B-field referenced to a guiding
center, such that the following hold:
(a ) x – xo = - i v⊥ exp (i Ω t)/ W = r sin (Ω t)
and
(b) y – yo = ± v⊥ exp (i Ω t)/ W = r cos (Ω t)
The preceding equation pair (a, b) describes a circular orbit around
the guiding center (xo, yo), with the direction of
gyration always such that the magnetic field generated by the charged particle
is opposite to the externally imposed field.
The gyro-period is: T = 2 p / Ω
Bear in mind the gyration energy:
E = m m B = m/2 (E/B) 2,
Guiding
center positions:
= (0.0056 m) sin [(1.7 x 10 7
/s) (3.5 x 10 -7 s/ 4)] = 0.0056 m
y – yo = r cos
(Ω e t)
=
(0.0056 m) sin [(1.7 x 10 7 /s) (3.5 x 10 -7 s/ 4)] = 0
To see how these values can be, note that: (Ω e T) = 6.28 rad = 2p rad
So: (Ω e T/ 4) =1.57 rad = p / 2
But: sin (p /2) = 1 and (cos p /2) = 0
So
the value for x – xo is simply dictated by the value
for r
No comments:
Post a Comment