FEQ: a new flux coordinates based equilibrium solver including both magnetic axis and separatrix

doi:10.1088/2058-6272/ac363e

	FEQ: a new flux coordinates based equilibrium solver including both magnetic axis and separatrix
	Jiang, Xinhao 1,2; Sun, Youwen 1
刊名	PLASMA SCIENCE & TECHNOLOGY
	2022
卷号	24
关键词	equilibrium flux coordinates tokamak
ISSN号	1009-0630
DOI	10.1088/2058-6272/ac363e
通讯作者	Sun, Youwen(ywsun@ipp.ac.cn)
英文摘要	Accurate tokamak plasma equilibrium solution in flux coordinates is crucial for many stability and transport studies. Different approaches for dealing with singularities in solving the nonlinear Grad-Shafranov (GS) equation in flux coordinates or also known as straight field line coordinates are proposed in this paper. The GS equation is solved by iterating the position of grids directly in flux coordinates, and hence, no additional errors are introduced due to mapping process for a convergent solution. The singularity at magnetic axis in flux coordinates is removed by using a novel coordinate transform technique. Different from other techniques previously developed, no assumption in boundary condition at magnetic axis is used. This is consistent with the fact that there is no physical boundary at the magnetic axis. A flux coordinate system with poloidal coordinate chosen as the geometric poloidal angle is proposed. It conquers the difficulty in no definition of poloidal coordinate in flux coordinates at separatrix because of the singularity at x-point(s) in a divertor configuration. It also simplifies the process for computing poloidal flux coordinate during the iteration for solving the nonlinear GS equation. Non-uniform grids can be applied in both radial and poloidal coordinates, which allows it to increase the spacial resolution near x-point(s) in a divertor configuration. Based on the model proposed in this paper, a new Flux coordinates based EQuilibrium solver (FEQ) in tokamaks is developed. The numerical solutions from this code agree well with both the analytic Solov'ev solution and the numerical one from the EFIT code for a divertor configuration in the EAST tokamak. This code can be applied for simulating different equilibria with prescribed shape, pressure and current profiles, i.e. including both limiter and divertor configurations, positive triangularity and negative triangularity, different beta, arbitrary magnetic shear profile etc. It provides a powerful and convenient fixed-boundary inverse equilibrium solver including both magnetic axis and separatrix in the solution for tokamak researches.
资助项目	National Key R&D Program of China[2017YFE0301100] ; National Natural Science Foundation of China[11 875 292]
WOS关键词	MHD-EQUILIBRIUM ; STABILITY ; TRANSPORT
WOS研究方向	Physics
语种	英语
出版者	IOP Publishing Ltd
WOS记录号	WOS:000727697000001
资助机构	National Key R&D Program of China ; National Natural Science Foundation of China
内容类型	期刊论文
源URL	[http://ir.hfcas.ac.cn:8080/handle/334002/126810]
专题	中国科学院合肥物质科学研究院
通讯作者	Sun, Youwen
作者单位	1.Chinese Acad Sci, Inst Plasma Phys, HFIPS, POB 1126, Hefei 230031, Peoples R China 2.Univ Sci & Technol China, Hefei 230026, Peoples R China
推荐引用方式 GB/T 7714	Jiang, Xinhao,Sun, Youwen. FEQ: a new flux coordinates based equilibrium solver including both magnetic axis and separatrix[J]. PLASMA SCIENCE & TECHNOLOGY,2022,24.
APA	Jiang, Xinhao,&Sun, Youwen.(2022).FEQ: a new flux coordinates based equilibrium solver including both magnetic axis and separatrix.PLASMA SCIENCE & TECHNOLOGY,24.
MLA	Jiang, Xinhao,et al."FEQ: a new flux coordinates based equilibrium solver including both magnetic axis and separatrix".PLASMA SCIENCE & TECHNOLOGY 24(2022).