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

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

作　　者：郑平卫[1] 何丽华[1] 黄千红[1] 邓盛[1] 龚学余[1] 

机构地区：[1]南华大学,湖南衡阳421001

出　　处：《计算物理》2017年第4期387-393,共7页Chinese Journal of Computational Physics

基　　金：国家自然科学基金(11405082;11375085);湖南省教育厅项目(14C1001)资助

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

关 键 词：FOKKER-PLANCK方程 分布函数 电子回旋波 加热与电流驱动 

分 类 号：O532[理学—等离子体物理]