HT-托卡马克中低杂波和离子波恩斯坦波协同作用的数值研究.docVIP

HT-7托卡马克中低杂波和离子波恩斯坦波协同作用的数值研究 第23卷第3期 2OO6年5月 计算物理 CHINESEJOURNALOFCOMPUTATIONALPHYSICSVo1.23.No.3 Mav.2006 [文章编号]1001.246X《2006103.0325-10 HT.7托卡马克中低杂波和离子波恩斯坦波 协同作用的数值研究 张公让,匡光力 (1.中国科学院等离子体物理研究所,安徽合肥230031;2.合肥工业大学,安徽合肥230009) [摘要】用隐含5点模式对二维Fokker-Planck方程进行离散化,调用NAG程序库中子函数1303UAF求解该方 程,并将求解二维Fokker-Planck方程的code与波的射线轨迹code相结合,根据HT.7托卡马克装置参数,分别对 有或没有IBW情况下的等离子体电流的全波驱动和在低杂波和离子伯恩斯坦波协同作用下低杂波的功率沉积分 布,驱动电流分布等进行数值模拟研究.结果显示:有IBW情况下的全波驱动电流分布比没有IBW时全波驱动 的电流分布区域要广一些,而且驱动的总电流也大一点;在双波协同作用下,全波驱动的等离子体电流空间分 布随IBW的频率和n谱的变化而变化;在两波协同作用下低杂波驱动非感应电流的效率得到了很大的改善. [关键词]低杂波;离子伯恩斯坦波;双波协同;数值模拟;电子朗道阻尼;全波驱动 [中图分类号]0532.23[文献标识码]A 0引言 多年来在HT.7托卡马克装置上进行了多轮低杂波(LHW)与离子波恩斯坦波(IBW)的协同实验,目的是 为了探究低杂波与离子波恩斯坦波的局域协同作用对提高由低杂波驱动的等离子体电流,通过改变等离子 体电流分布从而形成内部输运垒最终改善等离子体的约束以及等离子体的边界参数涨落对低杂波在等离子 体中传播和吸收的影响等. 自Fisch提出射频波驱动电流理论以后,低杂波驱动电流(LHCD)的数值模拟也逐渐得到了发展和完 善.1978年,Fisch用一维Fokker.Planck方程求解了电子速度分布函数.1979年,Karney等人通过数值模拟 研究发现与一维Fokker.Planck方程相比,考虑二维速度空间效应以后,驱动电流增加40%左右.1996年, Shkaro~ky和Shoucfi通过解三维Fokker.Planck方程,完成了考虑空间扩散情形下的LHCD数值计算.以上模 拟均假设准线性扩散系数在一定的速度空间内保持为恒定值.1986年,Bonoli等通过数值模拟,得到了低杂 波电流驱动过程中的射线轨迹,功率沉积分布以及驱动电流的分布.与此同时,我所也开发了一套模拟计算 LHCD的程序,并针对具体装置在欧姆加热(OH)的等离子体背景下进行了相关的数值模拟与物理分 析卜.但是,对于完全由低杂波驱动电流和IBW与LHW协同作用情况下电流驱动的数值模拟计算,我所乃 至国内至今仍无人论及,尚属一个全新的课题. 1理论模型 低杂波驱动电流的物理机制是:通过天线辐射不对称谱的微波能量,使得入射的低杂波在环向传播的动 量具有单向性,这样具有单向动量的波在沿环向传播时,通过朗道阻尼与等离子体中的共振电子发生有效的 相互作用,将动量传递给电子,从而在电子速度分布的共振区形成平台,因此形成一个由高能电子携带的非 感应等离子体电流.IBW是离子回旋频率上的一支准静电波,它只有扰动电场,方向平行于波的传播方向,而 扰动磁场几乎可以忽略.由于LHW驱动非感应电流效应,使得电子的速度空间分布函数产生了高能尾部,而 IBW的电场与这些高能电子发生朗道共振,从而使得这种高能尾部得到了进一步地维持和提高,结果使得非 感应电流增大. [收稿日期]2004—12—20:[修回Et期]2005—06—27 [基金项目]国家自然科学基金资助项目 [作者简介]张公让(1966一),男,副教授,博士,从事等离子体物理和微波工程方面的研究,中科院合肥等离子体物理研究所1126信箱 326计算物理第23卷 实验中,等离子体电流全波驱动的实现通常是:首先,利用欧姆放电产生等离子体,并在欧姆场的作用下 形成等离子体电流.然后投入LHW或者相继投入LHW和IBW,让Tokamak中的等离子体电流缓慢地爬升至 平顶区,采用低杂波的反馈控制适时地调节欧姆场电压以便使等离子体电流较长时间地维持在一个恒定的 水平.最后,关断欧姆场,让等离子体电流完全由LHW或由LHW+IBW来驱动,直至等离子体熄灭.此时若 能长久维持等离子体电流不变,便实现了Tokamak的稳态运行,这正是我们的最终目的. 理论上,用带有准线性项的Fokker.Planck方程能够描述等离子体与射频波之间的相互作用.在用射频波 驱动电流的情况下,电子分布函数的演化可以利用Fokker.Planc