PDRK:A General Kinetic Dispersion Relation Solver for Magnetized Plasma

作     者：谢华生 肖湧 

作者机构：Institute for Fusion Theory and Simulation and the Department of PhysicsZhejiang University 

出 版 物：《Plasma Science and Technology》 (等离子体科学和技术（英文版）)

年 卷 期：2016年第18卷第2期

页      面：97-107页

核心收录：

学科分类：070207[理学-光学] 07[理学] 070204[理学-等离子体物理] 0805[工学-材料科学与工程（可授工学、理学学位）] 0702[理学-物理学] 

基　　金：supported by the National Magnetic Confinement Fusion Science Program of China(Nos.2015GB110003,2011GB105001,2013GB111000) National Natural Science Foundation of China(No.91130031) the Recruitment Program of Global Youth Experts 

主　　题：plasma physics dispersion relation kinetic waves instabilities linear system matrix eigenvalue 

摘      要：A general,fast,and effective approach is developed for numerical calculation of kinetic plasma linear dispersion relations.The plasma dispersion function is approximated by J-pole expansion.Subsequently,the dispersion relation is transformed to a standard matrix eigenvalue problem of an equivalent linear system.Numerical solutions for the least damped or fastest growing modes using an 8-pole expansion are generally accurate;more strongly damped modes are less accurate,but are less likely to be of physical interest.In contrast to conventional approaches,such as Newton＇s iterative method,this approach can give either all the solutions in the system or a few solutions around the initial guess.It is also free from convergence problems.The approach is demonstrated for electrostatic dispersion equations with one-dimensional and twodimensional wavevectors,and for electromagnetic kinetic magnetized plasma dispersion relation for bi-Maxwellian distribution with relative parallel velocity flows between species.