论文摘要
对托卡马克的理论和实验研究是当前受控核聚变研究的中心问题。对托卡马克中的电流驱动的数值计算研究也已成为当今数值模拟研究的一个焦点。特别是近几十年来,人们对离子回旋频率范围内的快波电流驱动的研究投入了大量的精力.本文主要对托卡马克中的快波电流驱动进行了较系统的理论研究,并确定了全波计算的数值计算方法。用全波方法研究环形对称托卡马克中的离子回旋频率范围内（ICRF）的快波电流驱动（FWCD）问题,研究中考虑了有限拉莫尔半径（FLR）效应和平行色散,建立起了全波计算的物理模型,得到了有限拉莫尔半径假设下的热等离子体介电张量和便于计算机编程计算的全波方程的具体形式,最终确定了研究托卡马克中离子回旋频率范围内的快波电流驱动的数值计算方法:在托卡马克位形下,采用Boozer坐标系即（ e? r ,e?η,e?ζ）, e?r沿径向, e?ζ沿磁场方向,（ e? r ,e?η,e?ζ）三者正交呈右手螺旋。选择此坐标系可以使得平行波数k //在磁面上有尽可能少的变量。对麦克斯韦全波方程应用迦辽金弱变分变换,使对变量两次求旋的二阶微分方程变为对变量一次求旋的一阶积分-微分弱变分全波方程。对积分-微分弱变分全波方程方程的求解,在径向r上采用三次哈密顿有限元方法,在极向θ上应用快速傅立叶变换和卷积定理,主要通过这两种数值计算方法把积分-微分弱变分全波方程化为一个关于电场强度的大型带状线性方程组。大型带状线性方程组的求解应用了多种先进的非对称求解理论,采用了动态内存加载以及符号预排序和分析等先进技术,运用这些方法编写的求解大型带状线性方程组的程序包可以较高效、准确地求解关于电场强度的线性方程组。综合运用上述各种数值计算方法,我们编写了研究托卡马克中快波电流驱动的相关程序,并得到等离子体与波的基本参数（安全因子、电子温度、密度等）的径向分布,及平行波数、色散函数的径向分布和平行方向电场强度的初步结果。以上所得计算结果,都与TORIC code所得到的结果进行了对比,发现输出结果的分布与参考资料基本一致,这说明我们目前研究快波电流驱动的思路和运用的数学处理方法是正确的。程序还在进一步升级和完善中,目前所做的工作为进一步研究快波电流驱动问题奠定了坚实的基础。