- #1

- 46

- 1

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