Implementation of field-aligned coordinates in a semi-Lagrangian gyrokinetic code for tokamak turbulence simulation

Lei YE (叶磊); Xiaotao XIAO (肖小涛); Yingfeng XU (徐颖峰); Zongliang DAI (戴宗良); Shaojie WANG (王少杰)

doi:10.1088/2058-6272/aac013

Plasma Science and Technology > 2018 > 20(7): 74008-074008. > DOI: 10.1088/2058-6272/aac013

Lei YE (叶磊), Xiaotao XIAO (肖小涛), Yingfeng XU (徐颖峰), Zongliang DAI (戴宗良), Shaojie WANG (王少杰). Implementation of field-aligned coordinates in a semi-Lagrangian gyrokinetic code for tokamak turbulence simulation[J]. Plasma Science and Technology, 2018, 20(7): 74008-074008. DOI: 10.1088/2058-6272/aac013

Citation:

PDF (2104 KB)

Implementation of field-aligned coordinates in a semi-Lagrangian gyrokinetic code for tokamak turbulence simulation

¹ Institute of Plasma Physics, Chinese Academy of Science, Hefei 230031, People’s Republic of China
²Department of Modern Physics, University of Science and Technology of China, Hefei 230026, People’s Republic of China

Funds: This work is supported by National Natural Science Foun- dation of China under Grant Nos. 11505240, 11375196 and 11405174, and the National ITER program of China under Contract No. 2014GB113000.

More Information

Received Date: December 20, 2017

Graphical Abstract

Abstract

Abstract

Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al (2016 J. Comput. Phys. 316 180), to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field-aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.
- gyrokinetic simulation,
- field-aligned coordinates,
- quasi-neutrality equation,
- iterative algorithm,
- gyro-average operator

FullText(HTML)

References (39)

References

[1]	Garbet X, Idomura Y, Villard L and Watanabe T 2010 Nucl. Fusion 50 043002
[2]	Batchelor D et al 2007 Plasma Sci. Technol. 9 312
[3]	Lin Z, Tang W and Lee W 1995 Phys. Plasmas 2 2975
[4]	Watanabe T H, Sugama H and Ferrando-Margalet S 2008 Phys. Rev. Lett. 100 195002
[5]	Chen Y, Parker S E, Lang J and Fu G Y 2010 Phys. Plasmas 17 102504
[6]	Zhang W, Holod I, Lin Z and Xiao Y 2012 Phys. Plasmas 19 022507
[7]	Horton W 1999 Rev. Mod. Phys. 71 735
[8]	Beer M A, Cowley S and Hammett G 1995 Phys. Plasmas 2 2687
[9]	Chen Y and Parker S E 2003 J. Comput. Phys. 189 463
[10]	G?rler T, Lapillonne X, Brunner S, Dannert T, Jenko F, Merz F and Told D 2011 J. Comput. Phys. 230 7053
[11]	Hariri F and Ottaviani M 2013 Comput. Phys. Commun. 184 2419
[12]	Ye L, Xu Y, Xiao X, Dai Z and Wang S 2016 J. Comput. Phys. 316 180
[13]	Xiao X, Ye L, Xu Y and Wang S 2017 Commun. Comput. Phys. 22 789
[14]	Strang G 1968 SIAM J. Numer. Anal. 5 506
[15]	Ye L, Guo W, Xiao X and Wang S 2013 Phys. Plasmas 20 072501
[16]	White R B 2001 The Theory of Toroidally Confined Plasmas (Singapore: World Scientific)
[17]	Scott B 2001 Phys. Plasmas 8 447
[18]	Ottaviani M 2011 Phys. Lett. A 375 1677
[19]	Wang S 2012 Phys. Plasmas 19 062504
[20]	Wang S 2013 Phys. Rev. E 87 063103
[21]	Wang S 2014 Phys. Plasmas 21 072312
[22]	Wang S 2001 Phys. Rev. E 64 056404
[23]	Idomura Y, Tokuda S and Kishimoto Y 2003 Nucl. Fusion 43 234
[24]	Idomura Y, Ida M, Kano T, Aiba N and Tokuda S 2008 Comput. Phys. Commun. 179 391
[25]	Cheng C Z and Knorr G 1976 J. Comput. Phys. 22 330
[26]	Grandgirard V et al 2006 J. Comput. Phys. 217 395
[27]	Kwon J M, Yi D, Piao X and Kim P 2015 J. Comput. Phys. 283 518
[28]	Grandgirard V et al 2016 Comput. Phys. Commun. 207 35
[29]	Wang Z, Li J, Dong J and Kishimoto Y 2009 Phys. Rev. Lett. 103 015004
[30]	Crouseilles N, Mehrenberger M and Sellama H 2010 Commun. Comput. Phys. 8 484
[31]	Candy J and Waltz R 2003 J. Comput. Phys. 186 545
[32]	Lee W 1987 J. Comput. Phys. 72 243
[33]	Idomura Y, Wakatani M and Tokuda S 2000 Phys. Plasmas 7 3551
[34]	Kim C C and Parker S E 2000 Massively Parallel Three- Dimensional Toroidal Gyrokinetic Flux-Tube Turbulence Simulation (Academic)
[35]	Parker S E 2002 J. Comput. Phys. 178 520
[36]	Atkinson K E 2008 An Introduction to Numerical Analysis (New York: Wiley)
[37]	Dimits A et al 2000 Phys. Plasmas 7 969
[38]	Falchetto G et al 2008 Plasma Phys. Control. Fusion 50 124015
[39]	Xu Y, Ye L, Dai Z, Xiao X and Wang S 2017 Phys. Plasmas 24 082515

[1]	Debing ZHANG, Pengfei ZHAO, Yingfeng XU, Lei YE, Xianmei ZHANG. Gyrokinetic simulations of the kinetic electron effects on the electrostatic instabilities on the ITER baseline scenario[J]. Plasma Science and Technology, 2024, 26(9): 095101. DOI: 10.1088/2058-6272/ad4e78
[2]	Yifei LIU, Jiquan LI. Gyro-Landau-fluid simulations of impurity effects on ion temperature gradient driven turbulence transport[J]. Plasma Science and Technology, 2024, 26(1): 015101. DOI: 10.1088/2058-6272/ad0c9b
[3]	Xianzu GONG, Liqing XU, Jinping QIAN, Juan HUANG, Rui DING, A M GAROFALO, Annika EKEDAHL, Long ZENG, Erzhong LI, Bin ZHANG, Shiyao LIN, Biao SHEN, Mao WANG, Handong XU, Xinjun ZHANG, Miaohui LI, Guizhong ZUO, Qing ZANG, Haiqing LIU, Bo LYU, Liang WANG, Youwen SUN, Guosheng XU, Jiansheng HU, Damao YAO, Yu WU, Liqun HU, Bingjia XIAO, Nong XIANG, Kun LU, Yuntao SONG, Baonian WAN, Jiangang LI, the EAST Team. Realization of T_e0 > 10 keV long pulse operation over 100 s on EAST[J]. Plasma Science and Technology, 2023, 25(2): 022001. DOI: 10.1088/2058-6272/ac9cc6
[4]	Jianyuan XIAO (肖建元), Hong QIN (秦宏). Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics[J]. Plasma Science and Technology, 2021, 23(5): 55102-055102. DOI: 10.1088/2058-6272/abf125
[5]	JIANG Peng (江澎), LIN Zhihong (林志宏), Ihor HOLOD, XIAO Chijie (肖池阶). The Implementation of Magnetic Islands in Gyrokinetic Toroidal Code[J]. Plasma Science and Technology, 2016, 18(2): 126-130. DOI: 10.1088/1009-0630/18/2/05
[6]	ZHENG Pingwei(郑平卫), GONG Xueyu(龚学余), YU Jun(余俊), DU Dan(杜丹). Fully Implicit Iterative Solving Method for the Fokker-Planck Equation in Tokamak Plasmas[J]. Plasma Science and Technology, 2014, 16(11): 1000-1006. DOI: 10.1088/1009-0630/16/11/02
[7]	LI Dehui(李德徽), XIANG Nong(项农), LIN Yu(林郁), WANG Xueyi(汪学毅), YANG Cheng(杨程), MA Jun(马骏). Benchmark Simulations of Gyro-Kinetic Electron and Fully-Kinetic Ion Model for Lower Hybrid Waves in Linear Region[J]. Plasma Science and Technology, 2014, 16(9): 821-825. DOI: 10.1088/1009-0630/16/9/03
[8]	LIN Zhihong (林志宏), S. ETHIER, T. S. HAHM, W. M. TANG. Verification of Gyrokinetic Particle Simulation of Device Size Scaling of Turbulent Transport[J]. Plasma Science and Technology, 2012, 14(12): 1125-1126. DOI: 10.1088/1009-0630/14/12/17
[9]	T. TAKIZUKA. Development of the PARASOL Code and Full Particle Simulation of Tokamak Plasma with an Open-Field SOL-Divertor Region Using PARASOL[J]. Plasma Science and Technology, 2011, 13(3): 316-325.
[10]	GUO Wei, WANG Shaojie, LI Jiangang. Vacuum Poloidal Magnetic Field of Tokamak in Alternating-Current Operation[J]. Plasma Science and Technology, 2010, 12(6): 657-660.