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

被引:22
作者
XIAO, Jianyuan [1 ]
QIN, Hong [2 ]
机构
[1] Univ Sci & Technol China, Sch Nucl Sci & Technol, Hefei 230026, Peoples R China
[2] Princeton Univ, Plasma Phys Lab, Princeton, NJ 08543 USA
基金
中国国家自然科学基金;
关键词
curvilinear orthogonal mesh; charge-conservative; particle-in-cell; symplectic algorithm; whole-device plasma simulation; MAXWELL-VLASOV EQUATIONS; SYMPLECTIC INTEGRATION; COLLECTIVE PROCESSES; CHARGE CONSERVATION; GYROKINETIC THEORY; PLASMA SIMULATION; TIME; IMPLICIT; SCHEMES; TRANSPORT;
D O I
10.1088/2058-6272/abf125
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s(-1), agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E x B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.
引用
收藏
页数:24
相关论文
共 164 条
[41]   MULTISCALE PARTICLE-IN-CELL PLASMA SIMULATION [J].
FRIEDMAN, A ;
PARKER, SE ;
RAY, SL ;
BIRDSALL, CK .
JOURNAL OF COMPUTATIONAL PHYSICS, 1991, 96 (01) :54-70
[42]   The Sunway TaihuLight supercomputer: system and applications [J].
Fu, Haohuan ;
Liao, Junfeng ;
Yang, Jinzhe ;
Wang, Lanning ;
Song, Zhenya ;
Huang, Xiaomeng ;
Yang, Chao ;
Xue, Wei ;
Liu, Fangfang ;
Qiao, Fangli ;
Zhao, Wei ;
Yin, Xunqiang ;
Hou, Chaofeng ;
Zhang, Chenglong ;
Ge, Wei ;
Zhang, Jian ;
Wang, Yangang ;
Zhou, Chunbo ;
Yang, Guangwen .
SCIENCE CHINA-INFORMATION SCIENCES, 2016, 59 (07)
[43]   The geometric theory of charge conservation in particle-in-cell simulations [J].
Glasser, Alexander S. ;
Qin, Hong .
JOURNAL OF PLASMA PHYSICS, 2020, 86 (03)
[44]  
Glasser AS., ARXIV190204395V1
[45]  
Glasser AS., ARXIV190204396V1
[46]   H mode confinement in Alcator C-Mod [J].
Greenwald, M ;
Boivin, RL ;
Bombarda, F ;
Bonoli, PT ;
Fiore, CL ;
Garnier, D ;
Goetz, JA ;
Golovato, SN ;
Graf, MA ;
Granetz, RS ;
Horne, S ;
Hubbard, A ;
Hutchinson, IH ;
Irby, JH ;
LaBombard, B ;
Lipschultz, B ;
Marmar, ES ;
May, MJ ;
McCracken, GM ;
OShea, P ;
Rice, JE ;
Schachter, J ;
Snipes, JA ;
Stek, PC ;
Takase, Y ;
Terry, JL ;
Wang, Y ;
Watterson, R ;
Welch, B ;
Wolfe, SM .
NUCLEAR FUSION, 1997, 37 (06) :793-807
[47]   On plasma rotation induced by waves in tokamaks [J].
Guan, Xiaoyin ;
Dodin, I. Y. ;
Qin, Hong ;
Liu, Jian ;
Fisch, N. J. .
PHYSICS OF PLASMAS, 2013, 20 (10)
[48]   On the toroidal plasma rotations induced by lower hybrid waves [J].
Guan, Xiaoyin ;
Qin, Hong ;
Liu, Jian ;
Fisch, Nathaniel J. .
PHYSICS OF PLASMAS, 2013, 20 (02)
[49]   NONLINEAR GYROKINETIC EQUATIONS FOR TOKAMAK MICROTURBULENCE [J].
HAHM, TS .
PHYSICS OF FLUIDS, 1988, 31 (09) :2670-2673
[50]  
Hairer E., 2002, Springer Series in Computational Mathematics