Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Plasma in magnetic field in equilibrium

  1. Jan 23, 2012 #1
    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
     
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted



Similar Discussions: Plasma in magnetic field in equilibrium
Loading...