A.3 Poloidal plasma current

As is discussed in Sec. 3.2, to satisfy the force balance, g RBϕ must be a magnetic surface function, i.e., g = g(Ψ). Using this, expression (50) and (51) are written

μ J  = − 1-dg-∂Ψ-= dg-B ,
 0 R     R dΨ ∂Z   dΨ  R
(527)

and

        1dg ∂Ψ    dg
μ0JZ = R-dΨ-∂R-=  dΨBz,
(528)

respectively. The above two equations imply that

JR-=  BR-,
JZ    BZ
(529)

which implies that the projections of B lines and J lines in the poloidal plane are identical to each other. This indicates that the J surfaces coincide with the magnetic surfaces.

The poloidal plasma current density is usually small (compared with the toroidal plasma current density) but is important for some cases of interest and thus could not be safely neglected. Many model equilibria (e.g., Solovev equilibria, DIII-D cyclone base cases) frequently used in simulations assume that g is a spatial constant, i.e., neglecting the poloidal plasmas current. The conclusions drawn from these simulations could be misleading.

Using this and ∇⋅ J = 0, and following the same steps in Sec. 1.7, we obtain

I   = 1-2π[g(Ψ ) − g(Ψ )],
 pol  μ0      2     1
(530)

where Ipol is the poloidal current enclosed by the two magnetic surfaces, the positive direction of Ipol is chosen to be in the clockwise direction when observers look along ˆϕ. Equation (530) indicates that the difference of g between two magnetic surface is proportional to the poloidal current. For this reason, g is usually call the “poloidal current function”.

In the above, we see that the relation of g with the poloidal electric current is similar to that of Ψ with the poloidal magnetic flux. This similarity is due to the following differential relations:

{
  B = ∇ × A
  Ψ = RA ϕ

{
 μ0J = ∇ × B
 g = RB ϕ

The poloidal plasma current density Jp can be further written as

Jp ≡ JR ˆR( + JZˆZ          )
     1--  1-∂g-ˆ   1-∂g-ˆ
  =  μ0  −R ∂Z R + R ∂R Z  .
     1
  =  μ-∇g × ∇ϕ                                   (531)
      0
Using g = g(Ψ), Eq. (531) can also be written as
Jp =-1-dg-∇Ψ × ∇ϕ
    μ0 dΨ
    -1-dg-
  = μ0 dΨ Bp.                               (532)