Abstract
In this research we measured the plasma boundary shift using array of magnetic pick up coils on the IR T1 tokamak. Also we approximated the magnetic surfaces by an equilibrium calculation. Firstly, four magnetic probes were designed, constructed, and installed on outer surface of the IR T1 tokamak chamber and then plasma boundary displacement measured from them. On the other hand, magnetic surfaces approximated by equilibrium calculation of the Grad Shafranov equation based on expansion of free functions as quadratic in flux function.
Tokamak; Plasma Boundary Shift; Magnetic Surfaces
Measurement of the plasma boundary shift and approximation of the magnetic surfaces on the IRT1 Tokamak
A. Salar Elahi^{I}; M. Ghoranneviss^{II}
^{I}Physics Department, Faculty of Engineering, Islamic Azad University, Saveh Branch, Saveh, Iran
^{II}Plasma Physics Research Center, Science and Research Branch, Islamic Azad University, P.O. Box 14665678, Tehran, Iran
ABSTRACT
In this research we measured the plasma boundary shift using array of magnetic pickup coils on the IRT1 tokamak. Also we approximated the magnetic surfaces by an equilibrium calculation. Firstly, four magnetic probes were designed, constructed, and installed on outer surface of the IRT1 tokamak chamber and then plasma boundary displacement measured from them. On the other hand, magnetic surfaces approximated by equilibrium calculation of the GradShafranov equation based on expansion of free functions as quadratic in flux function.
Keywords: Tokamak, Plasma Boundary Shift, Magnetic Surfaces.
1. INTRODUCTION
To a very good approximation in tokamaks, the problem of achieving toroidal equilibrium separates into two parts. First, the magnetic configuration must provide radial confinement, i.e., radial pressure balance in the poloidal plane so that the pressure contours form closed nested surfaces. Second, the configuration must compensate the radially outward expansion force inherent in all toroidal geometries, i.e., toroidal force balance. But, in the second one, plasma may intend to shift inward or outward, which is very dangerous for tokamak plasma. Therefore, plasma equilibrium study is one of the fundamental problems of the magnetically confined plasmas. There are many available global solutions of the steady state magnetohydrodynamics equations, in particular, the GradShafranov equation. Control of plasma position plays an important role in plasma confinement and to achieve optimized tokamak plasma operation. Accurate determination of the plasma position during confinement time is essential to transport it to a control system based on feedback. Over the years different methods have been developed to analyze the equilibrium problems and determination of both plasma boundary shift and magnetic surfaces [121].
In this paper we presented magnetic probes method for measurement of the plasma boundary shift and approximation of magnetic surface based on analytical solution of the GradShafranov equation on the IRT1 Tokamak, which is a small, air core, low beta and large aspect ratio tokamak with a circular cross section, (see Table 1). Details of the magnetic probes method for the determination of plasma boundary shift will be discussed in section 2. Approximation of magnetic surfaces based on analytical solution of the GradShafranov equation will be presented in section 3. Experimental result also will be presented in section 4. Summary and conclusion are also will be presented in section 5.
2. MAGNETIC PROBES TECHNIQUE FOR THE MEASUREMENT OF PLASMA BOUNDARY SHIFT
Because of dependence of the plasma position and plasma current distribution to the magnetic field distributions around the plasma, magnetic pickup coils can give us information about the plasma position or plasma boundary shift. Poloidal and normal components of the magnetic fields distribution in the quasicylindrical coordinates (r,θ,Φ) around the circular cross section plasma are [1] :
where Λ = β_{p} + li/2 1 is the Shafranov parameter and a, b, I_{p} , R_{0} , β_{p} , li, and ΔR are the plasma minor radius, chamber minor radius, plasma current, chamber major radius, poloidal beta, internal inductance, and plasma boundary horizontal shift, respectively. Therefore, by rearranging the above equations the following relation for the horizontal shift of the plasma boundary is obtained:
As observable, this relation depends only on the plasma current and magnetic fields distribution. Also the vertical displacement of the plasma boundary is [2] :
Equations (1) and (2) are accurate for the low beta, large aspect ratio, and circular cross section tokamaks as IRT1, also:
According to this approach, four magnetic probes were designed, constructed, and installed on the IRT1 tokamak, two magnetic probes were located on the circular contour Γ of the radius b = 16.5cm in angles of θ = 0 and θ = π to detect the tangential component of the magnetic field B_{θ} and two magnetic probes are also located above, θ = π/2, and below, θ = 3π/2, to detect the normal component of the magnetic field B_{r} (see the Fig. (1)).
By substituting the poloidal and normal components of the magnetic fields which are measured by the four magnetic pickup coils (after compensation and integration of their output) in Equations (3) and (4)(ΔB_{2}_{θ} measured using the poloidal Mirnov coils), horizontal and vertical displacements of the plasma boundary were determined. Experimental results will be presented in the section 4.
3. APPROXIMATION OF THE MAGNETIC SURFACES BASED ON ANALYTICAL SOLUTION OF THE GRADSHAFRANOV EQUATION
Maxwell's equations together with the force balance equation from MHD equations, in the cylindrical coordinates (R, Z) reduce to the twodimensional, nonlinear, elliptic partial differential equation, or GradShafranov equation [3] :
where
and p(ψ) and F(ψ) are two free functions, and where µ_{0} and J are the vacuum permeability and plasma current density respectively. Many authors solved the GradShafranov equation by expanding the free functions on different order in ψ. In this section, we presented quadratic order (which proposed by Guazzotto [4]), and approximated on the IRT1. If we choose the free functions to be quadratic in ψ as [4] :
where ψ_{0} , p_{0} , b_{0} are the values of ψ, p, and F on magnetic surfaces axis, R_{0} is the tokamak major radius, and B_{0} is the vacuum toroidal field.
The GradShafranov equation reduces to:
where β_{0} = 2µ_{0}p_{0} /. With normalizing variables as R^{2} /= x, and Z/a = y, Eq. (9) can then be written as:
where ε = a/R_{0} , is the inverse aspect ratio, and
The solution of Eq. (10) in cylindrical coordinates (R, Z) can be written as:
where x = i(ε/)ρ, and for updown symmetric case, Y_{m} obtained as:
Also X_{m} (ρ) can be written as:
where λ_{m} = i , W_{λ}_{m ,}_{µ} (ρ), and M_{λ}_{m ,}_{µ} (ρ) are the Whittaker functions, and in this model µ = 1/2.
Guazzotto proposed only three terms for m, and then ψ can be written as:
where a_{m} , b_{m} , and k_{m} are nine unknown coefficients which must be determined. The first six of them can be obtained from boundary conditions. The boundary conditions for the points of inner, outer, and top of the plasma cross section are (see Figure (2)):
where a is the plasma radius, and condition for right convexity on the inboard midplane is:
also two conditions (one defining the location of the magnetic axis (R_{axis}) and the other the normalization for ψ on magnetic axis) are:
These are seven boundary conditions for seven unknown coefficients (a_{m} , b_{m} , R_{axis}).
But for other three coefficients, by setting k_{1} = 0 (the simplest solution for GSE independent of Z), and introduce one approximate value for the γ (γ ≈ 2(q/ε)^{2} (δ B_{Φ} /B_{Φ} ), negative for diamagnetism plasma, and where q = aB_{Φ} /R_{0}B_{θ} and ε are the safety factor and inverse aspect ratio respectively), and assuming that k_{2} be imaginary and k_{3} be real, the values of α, k_{2} , k_{3} can be approximated by minimizing the error function between traditional plasma shape (R = R_{0} + a cosθ, Z = a sinθ), and analytical plasma shape. Appropriate error function between them defined as follow [4] :
where subscripts (A and T) indicates the analytical and traditional plasma shape parameters respectively, and sum in minimum include 3 angles (0, π /2, π) for our purpose (circular plasma).
In general by minimizing the error function (e.g. by Mathematica software) as possible to zero, and finding optimal values for α, k_{2} , k_{3} , and also solving seven equations for the boundary conditions (Eqs. (16), (17), and (18), six unknown coefficients (a_{m} , b_{m}), moreover R_{axis} can be find. Therefore the approximated magnetic surfaces can be plotted by substituting these nine coefficients and also input parameters as I_{p}, α, R_{0}, and e in Eq. (15). For example the magnetic surfaces at t = 15ms correspond to approximated γ ≈ 0.32 be shown in the Fig. (5).
4. EXPERIMENTAL RESULT
For the determination of plasma boundary shift using the magnetic probes method, we needed for determination of the magnetic fields distribution around the plasma. Therefore we designed, and constructed a four magnetic pickup coils, and installed them on outer surface of the IRT1 chamber, and then desired components of the magnetic fields were measured. From this measurement and using the Eq. (3), the plasma boundary shift measured. The Figure (4) shows the result.
5. SUMMARY AND CONCLUSION
In this research we measured the plasma boundary shift using array of magnetic pickup coils on the IRT1 tokamak. Also we approximated the magnetic surfaces by an equilibrium calculation. Firstly, four magnetic probes were designed, constructed, and installed on outer surface of the IRT1 tokamak chamber and then plasma boundary displacement determined from them. The possible minor errors in this research are because of (1) approximation in the measurement of magnetic fields distribution because of discrete probes measurement, (2) possible error in the compensation of error fields, (3) the approximate values chosen for γ, and (4) the errors do not become zero during minimizing the error function.
(Received on 15 February, 2010)
 [1] H. Niomiya and N. Suzuki, Jpn. J. Applied Phys., 21 (9), 1323 (1982).
 [2] M. Valovic, Czech. J. Phys., B38, 65 (1988).
 [3] J. P. Freidberg, Ideal MHD (Clarendon, Oxford, 1987).
 [4] L. Guazzotto and J. P. Freidberg, PSFC/JA0714, Phys. Plasmas, (2007).
 [5] C. V. Atanasiu et al., Phys. Plasmas, 11, 3510 (2004).
 [6] G. S. Lee et al., Nucl. Fusion 41, 1515 (2001).
 [7] S. H. Seo, Phys. Plasmas 16, 032501 (2009).
 [8] S. G. Lee et al., Rev. Sci. Instrum. 79, 10F117 (2008).
 [9] A. J. H. Donne, Nucl. Fusion 47, S337S384, (2007).
 [10] F. Najmabadi et al., Fusion Eng. Des. 65, 143 (2003).
 [11] L. Guazzotto et al., Phys. Plasmas 11, 604 (2004).
 [12] A. Salar Elahi et al., IEEE Trans. Plasma Science 38 (2), (2010).
 [13] A. Salar Elahi et al., J. Phys. Scripta 80, 045501 (2009).
 [14] A. Salar Elahi et al., J. Phys. Scripta 80, 055502 (2009).
 [15] M. Emami et al., J. Plasma Physics 76 (1), (2009).
 [16] A. Salar Elahi et al., J. Fusion Energy 28 (4), 346 (2009).
 [17] A. Salar Elahi et al., J. Fusion Energy 28 (4), 390 (2009).
 [18] A. Salar Elahi et al., J. Fusion Energy 29 (1), 83 (2010).
 [19] A. Salar Elahi et al., J. Fusion Energy 29 (1), 88 (2010).
 [20] A. Rahimi Rad et al., J. Fusion Energy 29 (1), 73 (2010).
 [21] A. Salar Elahi et al., J. Phys. Scripta 81 (5), 055501 (2010).
Publication Dates

Publication in this collection
27 Sept 2010 
Date of issue
Sept 2010
History

Received
15 Feb 2010