Hamada coordinates: Difference between revisions

m (Reverted edits by Otihizuv (Talk) to last revision by 130.206.40.141)
Line 2: Line 2:


== Form of the Jacobian for Hamada coordinates ==
== Form of the Jacobian for Hamada coordinates ==
In this section, following D'haseleer et al we will translate the condition of straight current density lines into one for the '''Hamada''' coordinates Jacobian. For that we will make use of the equilibrium equation <math>\mathbf{j}\times\mathbf{B} = p'\nabla\psi </math>, which written in a magnetic coordinate system reads
In this section, following D'haseleer et al we will translate the condition of straight current density lines into one for the '''Hamada''' coordinates Jacobian. For that we will make use of the equilibrium equation <math>\mathbf{j}\times\mathbf{B} = p'\nabla\psi </math>, which written in a general magnetic coordinate system reads
:<math>
:<math>
\frac{-I'_{tor}\Psi'_{pol} + I'_{pol}\Psi'_{tor}}{4\pi^2\sqrt{g_f}}
\frac{-I'_{tor}\Psi'_{pol} + I'_{pol}\Psi'_{tor}}{4\pi^2\sqrt{g_f}}
- \mathbf{B}\cdot\nabla\tilde{\eta} = p'~.
- \mathbf{B}\cdot\nabla\tilde{\eta} = p'~.
</math>
</math>
Taking the [[Flux coordinates#flux surface average|flux surface average]] <math>\langle\cdot\rangle</math> of this equation we find <math>(-{I}'_{tor}{\Psi}'_{pol} + {I}'_{pol}{\Psi}'_{tor})= 4\pi^2{p}'\langle(\sqrt{g_f})^{-1}\rangle^{-1}</math>, so that we have
Taking the [[Flux coordinates#flux surface average|flux surface average]] <math>\langle\cdot\rangle</math> of this equation we find the MHD equilibrium equation
:<math>(-{I}'_{tor}{\Psi}'_{pol} + {I}'_{pol}{\Psi}'_{tor})= 4\pi^2{p}'\langle(\sqrt{g_f})^{-1}\rangle^{-1} = p'V'~.
</math>
In the last identity we have used the general [[Flux coordinates#Useful properties of the FSA|property of the flux surface average]] <math>\langle\sqrt{g}^{-1}\rangle^{-1} = \frac{V'}{4\pi^2}</math>. Then, from the MHD equilibrium, we have
:<math> \mathbf{B}\cdot\nabla\tilde{\eta} = {p}'\left(\frac{\langle(\sqrt{g_f})^{-1}\rangle^{-1}}{\sqrt{g_f}}-1\right)
:<math> \mathbf{B}\cdot\nabla\tilde{\eta} = {p}'\left(\frac{{V'}/{4\pi^2}}{\sqrt{g_f}}-1\right)~,
</math>
</math>
where <math>\tilde{\eta}</math> and <math>\sqrt{g_f}</math> depend on our choice of coordinate system.


In a coordinate system where <math>\mathbf{j}</math> is straight <math>\tilde{\eta}</math> is a function of <math>\psi</math> only, and therefore LHS of this equation must be zero in such a system. It therefore follows that the Jacobian of the '''Hamada''' system must satisfy  
Now, in the '''Hamada''' magnetic coordinate system that concerns us here (that in which <math>\mathbf{j}</math> is straight) <math>\tilde{\eta}</math> is a function of <math>\psi</math> only, and therefore LHS of this equation must be zero in such a system. It follows that the Jacobian of the Hamada system must satisfy  
:<math>
:<math>
  \sqrt{g_H} = \langle\sqrt{g_H}^{-1}\rangle^{-1} = \frac{V'}{4\pi^2}~,
  \sqrt{g_H} = \frac{V'}{4\pi^2}~.
 
</math>
</math>
where the last identity follows from the [[Flux coordinates#Useful properties of the FSA|properties of the flux surface average]]. The Hamada angles are sometimes defined in 'turns' (i.e. <math>(\theta, \xi) \in [0,1)</math>) instead of radians (<math>(\theta, \xi) \in [0,2\pi)</math>)). This choice together with the choice of the volume <math>V</math> as radial coordinate makes the Jacobian equal to unity. Alternatively one can select <math>\psi = \frac{V}{4\pi^2}</math> as radial coordinate with the same effect.
 
The Hamada angles are sometimes defined in 'turns' (i.e. <math>(\theta, \xi) \in [0,1)</math>) instead of radians (<math>(\theta, \xi) \in [0,2\pi)</math>)). This choice together with the choice of the volume <math>V</math> as radial coordinate makes the Jacobian equal to unity. Alternatively one can select <math>\psi = \frac{V}{4\pi^2}</math> as radial coordinate with the same effect.


== Magnetic field and current density expressions in a Hamada vector basis ==
== Magnetic field and current density expressions in a Hamada vector basis ==
204

edits