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