CFETR上alpha粒子波紋损失的初步计算草稿.docxVIP

CFETR上alpha粒子波纹损失的初步计算背景及程序介绍筹备中的中国聚变工程实验堆（CFETR）将要进行D-T燃烧反应，聚变功率会达到200MW，聚变反应产生的alpha粒子能量在3.5MeV左右，alpha粒子将把能量传递给热电子，其再通过碰撞加热等离子体中的离子，这种alpha粒子加热机制是实现点火和高聚变增益的关键。需要担心的是高能量粒子的约束问题，因为在装置中除alpha粒子外还有用于电流驱动和辅助加热的高能离子，这些高能量离子可通过波纹损失和反常输运损失掉，从而降低加热效率，而且长脉冲下对第一壁和偏滤器的热沉积是必须要考虑的问题。由于装置环向磁场线圈的分立性，CFETR的平衡磁场不是标准的轴对称位形，这对高能离子的约束是不利的。报告以此出发，据CFETR设计参数和现有装置上波纹损失的研究做一些初步计算，所用的模拟程序是导心轨道程序ORIBT。导心轨道程序ORBIT由美国普林斯顿大学White1984年开发，并完善至今。程序选取磁面坐标 ，将托卡马克平衡磁场定义为： ，磁面坐标分别表示极向磁通、极向角和环向角；)为磁场波纹； 为纵场线圈数。从平衡程序 EFIT 给出的托卡马克平衡磁场数值解中读取，由于该数值解中不包括EAST 纵场波纹，需在附属程序eqsub.f 中使用纵场波纹解析表达式叠加EAST 波纹场，从而得到ORBIT 所需的托卡马克平衡磁场。ORBIT 读取上述求解带电粒子导心运动方程所需的平衡量以及单个粒子的能量和速度 ，使用四阶的龙格-库塔算法求解粒子导心运动方程，并且追踪粒子的导心运动轨道，在二分之一能量慢化时间内，粒子轨道离开最后一个封闭磁面，定义为该粒子损失。ORBIT 中粒子的损失机理包括初始轨道损失、波纹损失(波纹磁阱捕获和随机波纹扩散)、螺距角散射效应以及能量慢化效应。ORBIT 所使用的描述带电粒子导心运动的Hamilton 方程组为： 公式（1） 式中， 为带电粒子在准稳态电场和磁场中的哈密顿量(动能+静电位能)； 为静电位能； 为安全因子；为扰动磁场；、分别为归一化的回旋半径和磁矩； 与磁面外的等离子体极向电流成正比； 与磁面内的环向总电流成正比， 。 工作内容概念设计中的CFETR主要参数如下：环向磁场BT（T）4.5-5等离子体电流IP（MA）8.5-10大半径R（m）5.7拉长比κ1.8-2小半径a（m）1.6三角形变δ0.4-0.8纵场线圈数N16辅助加热功率Paux（MW）65-70表格 1. CFETR主要物理参数 CFETR的磁体系统包括16个纵场磁体，中心螺线管和6个极向场磁体，为了满足实验的平衡需求，等离子体区域内（R=4100-7300mm）要求波纹度在0.5%以下，图1.是CFETR的线圈平面的环向磁场和波纹度分布。波纹度计算公式如下： 公式（2） 由于数值模拟的需要，还需得到波纹场 的解析表达式，对于非圆截面托克马克可以表示为： 公式（3） 其中。图表 1 . CFETR磁场(a)和波纹度(b)的分布根据CFETR波纹度分布使用最小二乘法拟合得出上述参数： 公式（4） 程序ORBIT所使用的是磁面坐标系下“平衡磁场+外加扰动场”的磁场位形，因此ORBIT中 由上面给出的 解析表达式经坐标转换得到。平衡磁场文件有平衡重建程序EFIT给出，报告分析了在等离子体电流IP=10MA下ITER-like、super-X和snowflake三种偏滤器平衡位形，如图表2.（a） (b) (c)图表 2 . 三种平衡磁场位形（a）ITER-like (b)super-X (c)snowflake程序中关于3.5MeV的Alpha粒子分布问题采用Monte-Carlo方法，设各向同性分布，剖面分布函数是：，a是小半径。计算中共模拟的10000个粒子，据alpha粒子慢化时间和磁轴上平行粒子绕行装置一周的时间，设置模拟时间是1/2能量慢化时间，即15000圈。调研CFETR文献，设置等离子体离子温度20KeV，另背景等离子体和杂质相关参数也根据典型的D-T燃烧实验设置。三个算例的参数文件除平衡位形外都一致，分别是ITER-like、super-X和snowflake，alpha粒子损失份额分别是14%、42%和60%，损失的alpha粒子的局域性很强都在外中平面以下，极向角分布范围分别是50o、65 o 和120 o，如图3。 (a) (b) (c)图表 3 . 损失alpha粒子极向角分布（a）ITER-like (b)super-X (c)snowflake 由上可知高能alpha粒子损失的局域性可能会对装置产生不利的影响，而长脉冲H-mode或hybrid-mode运行的CFETR的聚变功率将达到200MW，D-T燃料混合比例优化下以T作为工作气体注