Analytical form of Jacobian (need cleaning up)

Consider the Jacobian of the form

$\displaystyle \mathcal{J}= \mu \left( \frac{R}{R_0} \right)^m \psi^n,$ (463)

where $ m$ and $ n$ are arbitrary integers which can be appropriately chosen by users, $ \mu_0$ and $ R_0$ are constants included for normalization. In the iterative metric method of solving fixed boundary equilibrium problem, we first construct a coordinates transformation $ (\theta, \psi) \rightarrow (R
(\theta, \psi), Z (\theta, \psi))$ (this transformation is arbitrary except for that surface $ \psi = 1$ coincides with the last closed flux surface), then solve the GS equation in $ (\theta, \psi)$ coordinate system to get the value of $ \Psi $ at grid points, and finally adjust the value of $ (R (\theta, \psi),
Z (\theta, \psi))$ to make surface $ \psi = \ensuremath {\operatorname {const}}$ lies on a magnetic surface. It is obvious the Jacobian of the final transformation we obtained usually does not satisfy the constraint given by Eq. (463) since we do not use any information of Eq. (463) in the above steps. Now comes the question: how to make the transformation obtained above satisfy the constraint Eq. (463) through adjusting the values of $ \theta $? To make the constraint Eq. (463) satisfied, $ \theta $ and $ \psi $ should satisfy the relation

$\displaystyle \mathcal{J} (\theta, \psi) = \mu \left[ \frac{R (\theta, \psi)}{R_0} \right]^m \psi^n,$ (464)

which

Using Eq. (), we obtain

$\displaystyle \int_0^{\theta} d \theta = \int_0^{\theta} d \theta
\frac{\mathcal{J}}{\psi^n} \frac{1}{\mu} \left( \frac{R_0}{R} \right)^m $

$\displaystyle \theta = \int_0^{\theta} d \theta \frac{\mathcal{J}}{\psi^n} \frac{1}{\mu} \left( \frac{R_0}{R} \right)^m$ (465)

normalized to $ 2 \pi $, the normalized $ \theta $ is written as

$\displaystyle \theta = 2 \pi \frac{\int_0^{\theta} d \theta \frac{\mathcal{J}}{...
...R} \right)^m}{\int_0^{2 \pi} d \theta \mathcal{J} \left( \frac{1}{R} \right)^m}$ (466)

yj 2018-03-09