电子回旋波在D型截面等离子体中的传播.docVIP

  电子回旋波在 D 型截面等离子体中的传播# 李景春，龚学余，路兴强，彭晓炜，黄千红** （南华大学核科学技术学院，湖南 衡阳 421001） 5 摘要：利用能量原理，采用直接变分方法，开发出一种新的等离子体平衡方程求解程序。将 此程序与求解电子回旋（EC）波的波迹程序进行耦合，得到了 EC 波在 D 型截面托卡马克 等离子体中的传播轨迹。 关键词：等离子体物理；电子回旋波；D 型截面；波迹 中图分类号：TL61+4 10 The Propagation of Electron Cyclotron Waves in D-shaped Plasmas LI Jingchun, Gong Xueyu, LU Xingqiang, Peng Xiaowei, Huang Qianhong (School of Nuclear Science and Technology, University of South China, 15 HuNan HengYang 421001) Abstract: A plasma equilibrium code was developed by using direct variational method, the ray trajectories of EC wave in D-shaped Tokamak plasma were obtained by coupling the plasma equilibrium code to ray – tracing code. Keywords: Plasma Physics; electron cyclotron wave; D-shaped profile; ray trajectories 20 0 引言 电子回旋共振加热（ECRH）和电子回旋波电流驱动（ECCD）是磁约束聚变装置中等 离子体加热和电流驱动的重要方式。EC 波可以以很窄的波束传播,因此发射器不必很靠近等 离子体,另外它还具有高精确度的定域性能；随着长脉冲、高功率的回旋管源的显著发展[1-2]， 25 30 35 ECRH 和 ECCD 成为等离子体加热和电流驱动方案中的有力竞争者。在 ECRH 和 ECCD 理 论研究中，最基本的问题是计算出 EC 波在托卡马克等离子体中的传播。国外，文献[3-4] 应用等离子体色散关系，采用程函近似方法，追踪了电子回旋波在冷、热等离子体中的传播， 得到了与实验结果比较好的波迹。国内，彭晓炜、宋绍栋等人采用同样原理对 EC 波在等离 子体中的传播进行了研究[5-6]，并模拟了各装置上的 ECCD，但他们的模拟都是在圆截面的 近似下进行。然而真实的托卡马克为 D 型截面位形，D 型截面下的磁场、温度和密度分布 与圆截面不同,它们影响着 EC 波在等离子体中的传播。因此，数值计算 EC 波在 D 型截面等 离子体中的波迹便显得必要。 非圆截面等离子体的平衡位形一般通过数值求解平衡方程获得，数值求解 Grad-Shaf- ranov 方程可以直接求解（R,Z）坐标下的两维差分方程，也可以利用变分方法，通过求矩方 程给出近似解，前者需要更多的计算时间和更大的计算内存[7-8]。本文采用一种直接变分求 平衡的新方法，得到托卡马克的磁场、温度和密度分布，再将其插值到求解波迹方程的程序 中，得到 D 型截面下的 EC 波传播轨迹。 1 物理模型 射频波在托卡马克中的传播由射线轨迹方程描述。在满足几何光学近似下，射线轨迹方 基金项目：高等学校博士学科点专项科研基金资助课题（20094324120003） 作者简介：李景春，(1989-)，男，湖南永州人，硕士研究生，主要从事波与等离子体相互作用的研究。 通信联系人：龚学余，（1962-），男，教授，主要研究方向：核聚变与等离子体物理。 E-mail: gongxueyu@126.com -1-   40 程在环坐标 (r,? ,? ) 下，可以写成如下的形式:  d r d t  ??? ?D0??k r ?D0??? d?? d t  ??? ?D0??m ?D0??? d?? d t  ??? ?D0??n ?D0???  d k r d t  ? ?D0??r ?D0???  d m d t  ? ?D0??? ?D0???  d n d t  ? ?D0??? ?D0???  ? 0  （1） v 环 向 波 数 表 示 为 ： k??? m r ， k??? n R ， m 和 n 分 别 表 示 极 向 和 环 向 模 数 ，  45 R?? R0?? r cos? ，R 0 是托卡马克环的大半径。；D0 (?, k )