9.2 Expression of metric elements of magnetic coordinates (ψ,𝜃,ϕ)

Metric elements of the (ψ,𝜃,ϕ) coordinates, e.g., ψ ⋅∇𝜃, are often needed in practical calculations. Next, we express these metric elements in terms of the cylindrical coordinates (R,Z) and their partial derivatives with respect to ψ and 𝜃. Note that, in this case, the coordinate system is (ψ,𝜃,ϕ) while R and Z are functions of ψ and 𝜃, i.e.,

R = R (ψ,𝜃),
(201)

Z = R (ψ,𝜃).
(202)

Then R and Z are written as

      ˆ
∇R  = R = Rψ∇ ψ + R𝜃∇ 𝜃,
(203)

∇Z  = ˆZ = Zψ∇ ψ + Z𝜃∇𝜃,
(204)

wehre Rψ ∂R∕∂ψ, etc. Equations (203) and (204) can be solved to give

           1
∇ψ =  R-Z--−-Z-R--(Z 𝜃Rˆ − R 𝜃ˆZ),
       ψ 𝜃    ψ 𝜃
(205)

     -----1------    ˆ     ˆ
∇𝜃 = ZψR 𝜃 − R ψZ𝜃(ZψR − R ψZ).
(206)

Using the above expressions, the Jacobian of (ψ,𝜃,ϕ) coordinates, 𝒥 , is written as

  −1
𝒥    = ∇ψ × ∇𝜃 ⋅∇ϕ
     = −-------1-------(Z  R ˆϕ − R Z  ˆϕ)× ϕˆ
        (RψZ 𝜃 − ZψR 𝜃)2 𝜃 ψ     𝜃 ψ     R
               1
     = −(Z𝜃R-ψ −-R-𝜃Z-ψ)R,                                (207)
i.e.,
𝒥 = R (R 𝜃Zψ − R ψZ𝜃).
(208)

Using this, Expressions (205) and (206) are written as

        R
∇ ψ = − 𝒥 (Z 𝜃Rˆ− R 𝜃ˆZ)
(209)

and

∇𝜃 = R-(ZψRˆ− R ψˆZ).
     𝒥
(210)

Then the elements of the metric matrix are written as

        R2
|∇ψ|2 = 𝒥2(Z2𝜃 + R2𝜃),
(211)

|∇ 𝜃|2 = R2-(Z2 +R2 ),
       𝒥 2  ψ    ψ
(212)

and

            2
∇ψ ⋅∇ 𝜃 = − R-(Z 𝜃Z ψ + R 𝜃Rψ).
           𝒥2
(213)

Equations (211), (212), and (213) are the expressions of the metric elements in terms of R, Rψ, R𝜃, Zψ, and Z𝜃. [Combining the above results, we obtain

∇ψ ⋅∇ 𝜃    Z Z  + R R
-----2-= − -𝜃--ψ2---𝜃2-ψ.
 |∇ ψ|         Z𝜃 + R 𝜃
(214)

Equation (213) is used in GTAW code. Using the above results, hαβ = 𝒥-
R2α ⋅∇β are written as

       𝒥         1
hψ ψ =--2|∇ψ|2 = -(Z2𝜃 + R2𝜃)
      R          𝒥
(215)

      𝒥         1
h𝜃𝜃 = -2|∇ 𝜃|2 = -(Z2ψ + R2ψ ),
      R         𝒥
(216)

      𝒥             1
hψ𝜃 = R2-∇ψ ⋅∇ 𝜃 = − 𝒥-(Z 𝜃Z ψ + R 𝜃Rψ)
(217)

As a side product of the above results, we can calculate the arc length in the poloidal plane along a constant ψ surface, dℓp, which is expressed as

     ∘ ----2------2
dℓp = ∘ (dR)-+-(dZ)------------------
   =   (Rψdψ + R𝜃d𝜃)2 + (Zψdψ + Z𝜃d𝜃)2.
Note that = 0 since we are considering the arc length along a constant ψ surface in (R,Z) plane. Then the above expression is reduced to
dℓ  = ∘(R--d𝜃)2-+(Z-d𝜃)2
  p   ∘---𝜃---    𝜃
    =  R2𝜃 + Z2𝜃d𝜃

    = |𝒥-∇-ψ|d𝜃,                                (218)
        R
which agrees with Eq. (168).]