Plasma in magnetic field in equilibrium

  • Thread starter vibe3
  • Start date
  • #1
46
1

Main Question or Discussion Point

Hello,

Starting from the equilibrum equation of a steady-state plasma in a magnetic field:
[tex]
(1) \nabla P = \vec{J} \times \vec{B}
[/tex]
where [itex]P = nkT[/itex] is the plasma pressure, [itex]\vec{J}[/itex] is the current and [itex]\vec{B}[/itex] is the magnetic field, I am trying to solve for the magnetic field [itex]\vec{B}[/itex] in terms of a known pressure [itex]P[/itex].

Using [itex]\vec{J} = \frac{1}{\mu_0} \nabla \times \vec{B}[/itex], and the identity [itex]\vec{B} \times (\nabla \times \vec{B}) + (\vec{B} \cdot \nabla) \vec{B} = \nabla( B^2/2)[/itex], the first equation becomes
[tex]
(2) \nabla \left( P + \frac{B^2}{2 \mu_0} \right) = \frac{1}{\mu_0} (\vec{B} \cdot \nabla) \vec{B}
[/tex]
Dotting both sides by a path length [itex]\vec{ds}[/itex] and integrating from the point of interest [itex]\vec{r}[/itex] to some other point [itex]p_2[/itex], and then using the gradient theorem gives
[tex]
(3) P(p_2) + \frac{B^2(p_2)}{2 \mu_0} - P(\vec{r}) - \frac{B^2(\vec{r})}{2 \mu_0} = \frac{1}{\mu_0} \int_{\vec{r}}^{p_2} \left[ (\vec{B} \cdot \nabla) \vec{B} \right] \cdot \vec{ds}
[/tex]
The integrand of (3) can be expressed as
[tex]
(4) \frac{1}{\mu_0} (\vec{B} \cdot \nabla) \vec{B} = \hat{b} \frac{\partial}{\partial l} \left( \frac{B^2}{2 \mu_0} \right) - \hat{n} \frac{B^2}{\mu_0 R_c}
[/tex]
where [itex]l[/itex] is a coordinate along a magnetic field line, [itex]\hat{b}[/itex] is a unit vector tangent to the field line, [itex]\hat{n}[/itex] is a unit vector normal to the field line and anti-radial, and [itex]R_c[/itex] is the radius of curvature of the field line (see for example Plasma physics and controlled nuclear fusion By Kenrō Miyamoto, eq 6.7).

By choosing the path of integration either in the [itex]\hat{b}[/itex] direction or [itex]\hat{n}[/itex] direction perhaps we can get a solution.

CHOICE 1: [itex]\vec{ds} = \hat{n} ds[/itex]

Here, we will choose [itex]p_2[/itex] to be infinity, where we have the advantage that all quantities will go to 0. So (3) becomes
[tex]
(5) P(\vec{r}) + \frac{B^2(\vec{r})}{2 \mu_0} = \frac{1}{\mu_0} \int_{\vec{r}}^{\infty} \frac{B^2}{R_c} ds
[/tex]
Unfortunately, it does not seem that there is a closed form solution to this integral.

CHOICE 2: [itex]\vec{ds} = \hat{b} dl[/itex]

Here, we are integrating along a field line, but there does not seem to be an obvious choice for [itex]p_2[/itex], but the integral will be trivial. Leaving [itex]p_2[/itex] in place for now, we find that (3) becomes
[tex]
(6) P(p_2) + \frac{B^2(p_2)}{2 \mu_0} - P(\vec{r}) - \frac{B^2(\vec{r})}{2 \mu_0} = \frac{B^2(p_2)}{2 \mu_0} - \frac{B^2(\vec{r})}{2 \mu_0}
[/tex]
which unfortunately leads to [itex]P(\vec{r}) = P(p_2)[/itex] which isn't very helpful. Interestingly, if we integrate the other way on the field line, so that [itex]\vec{ds} = -\hat{b} dl[/itex] we get a completely different solution:
[tex]
(7) B^2(\vec{r}) = B^2(p_2) - \mu_0 (P(\vec{r}) - P(p_2))
[/tex]
which has the unfortunate property that the solution depends on which point [itex]p_2[/itex] we pick. Can anyone shed some light on where I'm going wrong with these integrals? Many thanks
 

Answers and Replies

Related Threads on Plasma in magnetic field in equilibrium

  • Last Post
Replies
4
Views
3K
Replies
1
Views
2K
Replies
4
Views
1K
Replies
4
Views
565
  • Last Post
Replies
1
Views
2K
Replies
7
Views
522
Replies
3
Views
742
Replies
2
Views
600
  • Last Post
Replies
0
Views
2K
Top