Development of Momentum Conserving Monte Carlo Simulation Code for ECCD Study in Helical Plasmas
Development of Momentum Conserving Monte Carlo Simulation Code for ECCD Study in Helical Plasmas
复制标题
螺旋等离子体 ECCD 研究动量守恒蒙特卡罗模拟代码的开发
DOI：
10.1051/epjconf/20158701010
复制
发表时间：
2015
影响因子：
--
通讯作者：
and Y. Moriya
中科院分区：
文献类型：
--
作者：
S. Murakami;S. Hasegawa;and Y. Moriya
为研究螺旋等离子体中的电子回旋电流驱动（ECCD），建立了GNET程序的并行动量守恒碰撞模型，在五维相空间中求解线性化的漂移动力学方程。为了保持平行动量守恒，我们在测试粒子碰撞项的基础上引入了场粒子碰撞项。两种类型的场粒子碰撞项被认为是。一种是高速度极限模型，其中动量守恒项不依赖于背景等离子体的速度，可以用简单的形式表示。另一种是直接从Fokker-Planck碰撞项导出的速度相关模型。在速度依赖模型中，场粒子算子可以用勒让德多项式表示，并引入Rosenbluth势，推导出每个勒让德多项式的场粒子项。在GNET代码中，我们引入了一个迭代过程来实现动量守恒碰撞算子。将高速极限模型应用于氦J等离子体的ECCD模拟。模拟结果表明，迭代格式具有良好的动量守恒性。
Parallel momentum conserving collision model is developed for GNET code, in which a linearized drift kinetic equation is solved in the five dimensional phase-space to study the electron cyclotron current drive (ECCD) in helical plasmas. In order to conserve the parallel momentum, we introduce a field particle collision term in addition to the test particle collision term. Two types of the field particle collision term are considered. One is the high speed limit model, where the momentum conserving term does not depend on the velocity of the background plasma and can be expressed in a simple form. The other is the velocity dependent model, which is derived from the Fokker–Planck collision term directly. In the velocity dependent model the field particle operator can be expressed using Legendre polynominals and, introducing the Rosenbluth potential, we derive the field particle term for each Legendre polynominals. In the GNET code, we introduce an iterative process to implement the momentum conserving collision operator. The high speed limit model is applied to the ECCD simulation of the heliotron-J plasma. The simulation results show a good conservation of the momentum with the iterative scheme.