基于三维Fokker-Planck方程的电子回旋波加热与电流驱动模拟

计算物理 ›› 2017, Vol. 34 ›› Issue (4): 387-393.

基于三维Fokker-Planck方程的电子回旋波加热与电流驱动模拟

郑平卫, 何丽华, 黄千红, 邓盛, 龚学余

南华大学, 湖南衡阳 421001

收稿日期:2016-07-06 修回日期:2016-09-21 出版日期:2017-07-25 发布日期:2017-07-25
通讯作者: 龚学余,男,主要从事核聚变等离子体物理研究,E-mail:gongxueyu@126.com
作者简介:郑平卫(1982-),男,博士生,主要从事射频波加热与电流驱动研究,E-mail:ghzpw@sina.com
基金资助:
国家自然科学基金（11405082；11375085）和湖南省教育厅项目（14C1001）资助

Electron Cyclotron Wave Heating and Current Drive Simulation Based on Three-dimensional Fokker-Planck Equation

ZHENG Pingwei, HE Lihua, HUANG Qianhong, DENG Sheng, GONG Xueyu

University of South China, Hengyang 421001, China

Received:2016-07-06 Revised:2016-09-21 Online:2017-07-25 Published:2017-07-25

摘要/Abstract

摘要： 利用反弹平均的三维Fokker-Planck方程，对电子回旋波加热和电流驱动进行数值模拟.考虑超热电子径向扩散对电流驱动的影响，在方程中加入径向扩散输运项，采用九点格式的中心差分对方程进行数值离散得到系数矩阵，采用不完全LU分解对系数矩阵进行预处理，利用双共轭梯度稳定法求解得到分布函数.在不考虑电子径向扩散输运条件下，得到电子回旋波驱动电流密度与功率沉积密度的分布；考虑径向扩散输运的计算结果与BANDIT3D进行比较，驱动电流分布的趋势基本一致.

关键词: Fokker-Planck方程, 分布函数, 电子回旋波, 加热与电流驱动

Abstract: Driven current profile and absorbed power profile are obtained with three-dimensional bounce-averaged Fokker-Planck equation in the presence of electron cyclotron wave. An additional term relate to radial transport is added to the equation accounting for influence of super-thermal electrons radial diffusion for current drive. Coefficient matrix is got by discretizing Fokker-Planck equation using nine-point center differential scheme. Incomplete LU decomposition is executed for precondition purpose. Making use of BiCGSTAB method, distribution function is calculated numerically, from which the drive current profile and absorbed power of electron cyclotron wave is calculated without radial transport. Driven current profile is compared with that of BANDIT3D with radial transport considered. It is shown that tendency of our result matched well with BANDIT3D.

Key words: Fokker-Planck equation, distribution function, electron cyclotron wave, heating and current drive

中图分类号:

O532

郑平卫, 何丽华, 黄千红, 邓盛, 龚学余. 基于三维Fokker-Planck方程的电子回旋波加热与电流驱动模拟[J]. 计算物理, 2017, 34(4): 387-393.

ZHENG Pingwei, HE Lihua, HUANG Qianhong, DENG Sheng, GONG Xueyu. Electron Cyclotron Wave Heating and Current Drive Simulation Based on Three-dimensional Fokker-Planck Equation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2017, 34(4): 387-393.

参考文献

[1] PEYSSON Y, BASIUK V, FOURMENT C, et al. Advances in plasma physics research(Vol.4)[M]. New York:Nova Science, 2003:1-29.
[2] HARVEY R W, SAUTER O, PRATER R, et al. Radial transport and electron cyclotron current drive in the TCV and DⅢ-D tokamaks[J].Physicals Review Letters, 2002, 88(20):205001(1)-(4).
[3] PRATER R. Heating and current drive by electron cyclotron waves[J].Physics of Plasmas, 2004, 11(5):2349-2376.
[4] FISCH N J, BOOZER A H. Creating an asymmetric plasma resistivity with waves[J].Physicals Review Letters, 1980, 45(9):720-722.
[5] OHKAWA T. Steady state operation of tokamaks by rf heating[R/OL]. General Atomics Report GA-A13847, 1976. http://www.ntis.gov/search/product.aspx?ABBR=PB2000108008.
[6] ZHENG P W, GONG X Y, YU J,et al. Fully implicit iterative solving method for the Fokker-Planck equation in tokamak plasmas[J].Plasma Science & Technology, 2014, 16(11):1000-1006.
[7] ROSA P R, ZIEBELL L F. Current drive by EC waves in the presence of magnetic islands and transport[J].Plasma Physics and Control Fusion, 2008, 50(9):095002(1)-095002(20).
[8] KENNEL C F, ENGELMANN F. Velocity space diffusion from weak plasma turbulence in a magnetic field[J].Physics of Fluids, 1966, 9(12):2377-2388.
[9] LERCHE I. Quasilinear theory of resonant diffusion in a magneto-active relativistic plasma[J].Physics of Fluids, 1968, 11(8):1720-1726.
[10] KARNEY C F F. Fokker-Planck and quasilinear codes[J]. Computer Physics Reports, 1986, 4(3-4):183-244.
[11] SHKAROFSKY I P, SHOUCRI M. Modelling of lower hybrid current drive in the presence of spatial radial diffusion[J].Nuclear Fusion, 1997, 37(4):539-547.
[12] KILLEEN J, MCCOY G D, MIRIN A A. Computational methods for kinetic models of magnetically confined plasmas[M]. New York:Springer, 1986:16.
[13] MCKENZIE J S, O'BRIEN M R, COX M. Solution of three-dimensional Fokker-Planck equations for tokamak plasmas using an operator splitting technique[J].Computer Physicals Communications, 1991, 66(2-3):194-206.
[14] PEYSSON Y, SHOUCRI M. An approcimate factorization procedure for solving nine-point elliptic difference equations application for a fast 2-D relativistic Fokker-Planck solver[J].Computer Physicals Communications, 1998, 109(1):55-80.
[15] 李荣华.偏微分方程数值解法[M].北京:高等教育出版社, 2005:210-214.