Numerical Solution of Grad-Shafranov Equation for the Distribution of Magnetic Flux in Nuclear Fusion Devices

Authors : Selçuk Han AYDIN

Pages : 265-275

Doi:10.3906/sag-1301-67

View : 14 | Download : 4

Publication Date : 0000-00-00

Article Type : Research Paper

Abstract :The magnetohydrodynamic equilibrium in an axisymmetric plasma is described by the Grad-Shafranov insert ignore into journalissuearticles values(GS); equation in terms of the magnetic flux. Boundary element method insert ignore into journalissuearticles values(BEM); is suitable for plasma equilibrium since it requires the discretization of only the plasma boundary which changes shape during the operation of an actual fusion device. In this paper, numerical solutions of the GS equation are obtained by using the boundary element method, the finite element method insert ignore into journalissuearticles values(FEM); and the differential quadrature method insert ignore into journalissuearticles values(DQM); for a rectangular plasma when the source term insert ignore into journalissuearticles values(current density function); on the right hand side is assumed to be a monomial. Our aim is also to find the most applicable numerical procedure between those three methods for different plasma profiles. We transform the equation to the homogeneous one with a particular solution eliminating the domain integral in the BEM formulation. For the source term containing the magnetic flux insert ignore into journalissuearticles values(nonlinear right hand side); an iterative procedure is made use of in the BEM and FEM formulations. It is found that the FEM gives better accuracy for a D-shape tokamak plasma whereas the BEM is more suitable for a Solov`ev tokamak plasma. The solutions agree very well with the previously published numerical solutions for a rectangular plasma.
Keywords : Grad Shafranov equation, BEM, FEM, DQM