-
Confinement of fusion born alpha particles in tokamak is the key issue to burning plasma. Apart from toroidal field ripple, instabilities can induce energetic particles to lose and be redistributed. Based on the parameters of China Fusion Engineering Testing Reactor (CFETT) hybrid scenario, alpha particle distribution and neoclassical tearing mode structure, the alpha particle loss induced under perturbation of ripple and neoclassical tearing mode (NTM) is calculated with the guiding center code ORBIT. The inputs have the initial distribution of alpha particles which is obtained with the TRANSP/NUBEAM code, the static NTM perturbation with different amplitudes which is obtained from TM1 code, and the ripple field from engineering design. The results show that the heat load on last closed flux surface is about 0.1 MW/m2, with ripple and collision included. The collisionless stochastic ripple diffusion is the main loss channel of initial alpha particle distribution in the CFETR, and the ripple perturbation has no influence on passing particles. The loss fraction does not increase with the NTM perturbation amplitude increasing, the synergistic effect is negligible. The scanning of ripple amplitude shows that the synergistic effect is slight. The monoenergetic initial distribution of alpha particles can give different types of orbits in the plane of (
Pζ,μ ), such as the domains of trapped particle and passing particle, lost particle and confined particle. The trapped fraction of initial alpha particles is about 27%, ripple loss region in phase space is narrow and away from the main trapped particle distribution. The increasing of ripple perturbation in simulation does enlarge the ripple loss domain in the phase space (Pζ,μ ), which is corresponding to a lager ripple loss fraction and has more trapped-passing boundaries. The NTM perturbation does enlarge the orbit excursions of trapped particles, and thus increasing the trapped passing transition near the boundary. The slight synergistic effect in calculation with larger ripple amplitude is explained by ripple loss region having more trapped-passing boundaries, not by the profile flattening of trapped particles. The NTM perturbation and finite collision can transit the passing particle to trapped particle near the boundary. With the help of kinetic Poincare plot, neither direct particle loss nor profile flattening of trapped particles is observed. The loss fraction enhancement can happen only when the profile flattening of trapped particles takes place within the ripple loss region, which is not the case in CFETR. The conclusion of this work contributes a lot to the design of CFETR and the study of alpha particle physics.-
Keywords:
- tokamak /
- alpha particle /
- magnetic ripple /
- neoclassical tearing mode
1. 引 言
磁约束核聚变研究在约70年的发展后到了建设聚变实验堆阶段, 国际热核聚变实验堆(ITER)的工程建设持续推进, 预期2025年开始放电. 我国在ITER物理设计和工程技术基础上发展的中国聚变工程试验堆(CFETR)作为ITER和聚变电站商用示范堆(DEMO)的过渡装置完成了初步的物理和工程设计[1-3]. 氘氚聚变产物alpha粒子的良好约束是获得高能量增益的前提, CFETR作为托卡马克类型的装置, 分立存在的纵场线圈引起的磁场波纹会引起一定份额的EP损失, 此外, 各种磁流体力学不稳定性不但引起直接的粒子损失, 还会改变EP分布函数, 聚变堆面向等离子体部件的热负荷限制和物理目标的实现只能允许小尺度的alpha粒子损失和再分布[4,5].
ITER和CFETR获得的燃烧等离子体让聚变研究有了新的参数范围和研究对象, 包括3.52 MeV高能量alpha粒子输运. 现有装置EP研究对象是中性束注入(NBI)和离子回旋加热引入的粒子分布, 由于背景参数和分布差异, 现有装置研究结论并不能直接外推到聚变堆级别参数[6].
理论研究发现alpha粒子会致稳锯齿模, 巨型锯齿可能激发新经典撕裂模(NTM), NTM限制等离子体比压并破坏EP约束, 需要额外的辅助手段控制NTM[7-9]. 在DIII-D, TFTR和ASDEX-U等装置上都观察到NTM引起了显著EP损失, 确认了粒子损失机制是相空间形成共振岛并被第一壁截断, 或共振岛相互重叠形成宽域的混沌区域[10-15]. 这些研究表明NTM对EP约束有显著影响, 包括粒子损失和分布剖面坍塌. 作者已发表工作也表明alpha粒子分布剖面外延、展平会引起波纹损失的显著增加[16], 原因是磁场波纹扰动幅度从等离子体芯部到边界呈自然指数增大, alpha粒子剖面坍塌会使更多的粒子进入波纹损失区. 显著的粒子损失不但会损伤第一壁, 且会降低聚变功率和辅助加热效率. 基于此, 本文利用CFETR混杂运行模式参数和NTM扰动数值结果, 研究NTM和磁场波纹扰动下的粒子损失叠加效应和损失机制. 本项研究对CFETR物理设计集成模拟流程中alpha粒子物理部分提供重要参考, 并有助于理解燃烧等离子体.
研究EP损失的主要方法是基于试探粒子模型的蒙特卡罗蒙特卡罗样本粒子跟踪, 本文的工作流程即在平衡位形叠加扰动后的总磁场中跟踪大量的样本粒子, 并在跟踪过程中考虑EP和背景的碰撞, 跟踪特定时间段后分析损失和未损失粒子信息. 跟踪电磁场中的带电粒子是经典研究, EP物理常需要研究长时间尺度下的粒子行为, 为了理论分析和数值计算简单高效, 对空间尺度远比拉莫尔回旋半径大且时间尺度远比回旋周期长的研究对象, 导心模型是有力的工具. 本文利用哈密顿导心轨道跟踪程序ORBIT[17], 在半个慢化时间内跟踪样本粒子分布, 初始alpha粒子分布来自输运程序TRANSP/NUBEAM[18], 平衡位形来自CFETR物理设计, 磁场波纹扰动来自工程设计, NTM扰动是基于初始值程序TM1物理分析结果的数值重建. ORBIT程序在直线磁面坐标系中求解哈密顿导心运动方程, 粒子的空间坐标为(
ψp,θ,ζ ), 分别是极向磁通, 极向角和广义环向角. 粒子的速度由能量E, 磁矩μ和螺距角正负号决定. 粒子导心的哈密顿量为H=ρ2//B2/2+μB+Φ ,ρ//=v///B 是归一化的平行速度, B是总磁场强度,μ=v⊥2/(2B) 是粒子磁矩, Φ是本文未引入的电势. 哈密顿运动方程为˙θ=∂H/∂Pθ ,˙Pθ=−∂H/∂θ ,˙ζ=∂H/∂Pζ ,˙Pζ=−∂H/∂ζ ,Pζ 和Pθ 是环向和极向正则角动量[17]. 由于托卡马克平衡位形的轴对称假设和粒子运动周期性质, 三个物理量(Pζ,E,μ) 即可决定EP轨道.论文第2节描述了CFETR混杂运行模式等离子体参数分布、alpha粒子初始分布、磁场波纹和NTM扰动的数值重建. 第3节是详细的粒子跟踪计算结果, 包括磁场波纹损失、NTM扰动叠加以及人为增大装置波纹度下的损失模拟. 第4节是结果讨论, 即粒子损失扰动叠加效应不显著的原因分析. 第5节是结论.
2. CFETR装置参数和全尺度粒子跟踪模拟输入
2.1 CFETR等离子体参数和初始alpha粒子分布
经过物理和工程设计迭代, 目前CFETR主机参数如表1所列. 作为实现稳态氚自持的托卡马克装置, CFETR物理实验分两个阶段实现200—1500 MW的聚变功率, 等离子体放电占空比约为0.3—0.5, 能量增益因子Q约为5—30[19]. 基于集成模拟流程提出了CFETR稳态运行模式和混杂运行模式两种, 本文基于混杂运行模式(v201806)参数, 利用输运程序TRANSP/NUBEAM计算得到了alpha粒子初始分布下的蒙特卡罗罗样本粒子. NUBEAM模块是分析托卡马克实验中中性束注入和聚变产物的含时集成工具. 混杂运行模式(v201806)的背景等离子体剖面见图1, 聚变功率约为
Pfus=1GW , 台基顶部的等离子体密度为格林极限的90%. 初始alpha粒子密度分布见图2, 初始alpha粒子产自氘氚聚变反应, 能量单值分布为3.52 MeV, 速度基本各向同性. 在没有不稳定性引起的显著输运和再分布时, alpha粒子经与背景等离子体库伦碰撞慢化并形成氦灰, 符合就地慢化模型. 平衡位形为下单零偏滤器位形, 纵场和等离子体电流方向都是俯视逆时针, 安全因子剖面见图4(a).表 1 CFETR与其他托卡马克装置主机参数对比Table 1. Main parameters comparison of CFETR and other tokomak facilities.参数 CFETR ITER EAST 磁轴场强BT0/T 6.5 5.3 2 等体大半径R0/m 7.2 6.2 1.9 等体小半径a/m 2.2 2.0 0.5 等体电流Ip/MA 14 15 1 纵场磁体柄数N 16 18 16 2.2 纵场波纹度分布和数值重建
CFETR的环向磁场由16柄D形线圈实现, 对比其他托卡马克装置, 柄数较少, 较大磁场波纹度对EP约束影响值得关注. 定义在极向截面的波纹度为
δ(r,θ)=(Bmax−Bmin)/(Bmax+Bmin) ,Bmax 和Bmin 分别是极向位置固定但在不同环向角处的磁场最大值和最小值. 工程设计中没考虑包层铁素体钢材料, 仅有纵场线圈时的波纹度分布见图3(b), 波纹度在等离子体外中平面附近最大, 约0.35%, 波纹扰动幅度从芯部到边界呈自然指数增大. 虽然CFETR满足托卡马克纵场磁体设计的一般要求, 即波纹度最大值小于0.5%, 但整体的波纹度较大, 芯部波纹度比ITER大一个量级. 在ORBIT中引入波纹扰动时, 磁面坐标系中叠加纵场波纹后的总磁场可表达为B(ψp,θ,ϕ)=B0(ψp,θ)[1+ δcos(Nϕ)] , 采用针对D形线圈波纹度几何分布的解析公式δ(R,Z)=δ0exp{[(R−Rrip)2+bripZ2]0.5/ wrip} 来拟合工程设计值. 拟合系数为δ0=1.57× 10−5 ,Rrip=6−0.0062Z2m ,brip=0.021 和wrip= 0.63m . 波纹扰动的解析函数实现结果和工程设计值对比见图3, 可知在等离子体约束区数值重建精确[16].2.3 NTM扰动函数分布
为了得到NTM扰动函数的模结构分布, 利用初始值程序TM1在柱坐标系中求解简化磁流体力学方程[20], 该坐标系中磁场表达为
B=B0t^et− (nr/mR)B0t^ep+∇ψh×^et , 即平衡磁场的环向和极向分量,再叠加螺旋磁通表达的扰动场. NTM的螺旋扰动场为δB=∇ψh×^et ,ψh 是螺旋磁通, 傅里叶展开为ψh= ∑iψi(r)exp[i(miθ−niϕ)] . 针对CFETR混杂运行模式, TM1程序计算了不同电子回旋电流驱动方案下NTM的饱和幅度, 得到主要的模数是2/1和3/2 NTM, 扰动函数剖面在不同饱和幅度下保持不变, 因此基于TM1程序数值结果通过高斯拟合在ORBIT中引入NTM模结构并扫描幅度. 图4给出的是2/1和3/2 NTM的平衡分量和一阶分量的高斯拟合结果, 平衡分量的极值点精确对应各自的有理面径向位置[21].ORBIT程序中磁面坐标系下扰动场表达为
δB=∇×αB0 ,α=αmn(ψp)sin(nζ−mθ−ωnt) 是时间和位置的标量函数. 柱坐标系中忽略高阶分量后的扰动磁场为δB=−∇×ψ1ˆϕ , 平衡场的共变基表达为B=∇ϕ+r2q(r)∇θ , 考虑到极向磁场远比环向场小, 磁面坐标系中的扰动函数可近似为−ψ1=α . 因此TM1程序中NTM计算结果的一阶分量用于表达ORBIT程序中扰动函数α[22]. 在给定平衡位形下, 2/1和3/2 NTM扰动幅度的比例固定, 本文算例为α2/1/α3/2=4.38 . 扰动函数α的量纲是长度, 在ORBIT程序中归一化单位为磁轴大半径R0 . 对于远比EP运动频率小的NTM, 假设为零频率磁扰动是合理且便于理论探讨的, 因此本文引入NTM为静止磁扰动. 为保持理论推导完整性, 下文保留扰动频率符号. 不同NTM扰动幅度引起的局域磁场拓扑变化见图5的庞加莱图, 磁岛的径向宽度正比于扰动幅度平方根, 较高的扰动幅度下可观察到磁岛相互重叠后的磁场无规化.3. 初始分布alpha粒子导心轨道跟踪模拟
3.1 alpha粒子磁场波纹损失
理论上任何破坏托卡马克环向轴对称性质的扰动都会降低EP约束, 如有限数目纵场线圈引起的磁场波纹. 纵场波纹扰动导致EP运动的正则环向角动量不守恒, 因此粒子轨道不闭合, 有额外的漂移, 通行粒子因磁力线回旋变换会平均掉该扰动, 俘获粒子受此扰动影响, 香蕉轨道转折点有一定扰动位移. 在磁场波纹扰动下, EP有两个损失通道, 局域波纹磁阱俘获和随机波纹扩散. 局域波纹磁阱是在波纹度较大时, 如外中平面附近, 形成二级磁阱, 平行速度较小的俘获粒子经过该区域时就俘获在该环向局域磁阱内, 只剩径向漂移快速损失掉或经碰撞和漂移逃出该损失区. 波纹磁阱俘获损失依赖于粒子速度、平衡位形、磁阱分布和碰撞率, 波纹磁阱俘获损失的时间尺度在几个极向回弹周期内. 当波纹度较小不足以形成局域磁阱时, 有限的波纹扰动也会使俘获粒子连续香蕉轨道之间解耦, 形成随机波纹扩散, 碰撞效应可降低该损失机制的随机阈值[23].
EP初始能量对应的碰撞频率分布由TRANSP/NBUEAM基于背景等离子体参数计算得到, 本文取典型代表值并在粒子跟踪过程中设定径向恒定, CFETR初始能量alpha粒子的典型能量慢化率为
νε0=1.8s−1 , 螺距角散射率为νpa0=0.16s−1 . 试探粒子程序时采用的碰撞修正模型涉及多个前提, 等离子体中EP产生以后, 分布函数在其损失或消失前经历复杂过程. EP分布函数改变的时间尺度远比其运动特征时间尺度长, EP密度远比背景等离子体密度低, 忽略频率较小的EP之间的碰撞, 仅考虑EP与背景的碰撞. EP与背景粒子之间小角度二体碰撞累计效应在托卡马克实验参数范围附近是主要的碰撞项, 考虑到等离子体德拜屏蔽, 忽略了频率远小于此的大角度碰撞, 虽然这种近距离碰撞可以引发核聚变反应. 基于此, ORBIT程序在跟踪粒子过程中对粒子即时能量Ei 和螺距角λi=v///v 在一个计算步长Δt 后分别修正为En=Ei−(2νεΔt)[Ei−(32+EνεdνεdE)T]±2[TEi(νεΔt)]0.5 和λn=λi(1−νpaΔt)±[(1−λi2)νpaΔt]0.5 ,T(eV) 是等离子体温度[24]. 高能alpha粒子的慢化主要来自与背景电子的碰撞, 电子慢化时间与EP能量无关, 因此对于短时间行为的EP能量慢化率可近似认为与能量无关νε≈νε0 , 螺距角散射率随能量慢化降低而增大νpa=νpa0(E0/E)1.5 .为有效和节省, 本文在半个慢化时间内跟踪初始alpha粒子样本, 只考虑纵场波纹扰动和碰撞效应的粒子损失份额约为1.4%, 如图6所示. 波纹损失粒子信息如图7, 经历和波纹扰动无关的初始轨道损失之后, 随着波纹扰动和碰撞效应的时间累计, 粒子损失份额在半个慢化时间后达到1.4%, 损失粒子能量有慢化成分, 损失粒子主要为俘获粒子. ORBIT程序在最后一个闭合磁面(LCFS)内跟踪粒子, 外出LCFS即定义为损失, 由于LCFS和装置第一壁接近, 因此可以根据该处的粒子沉积位置、数量和能量信息评估局域热负荷. 图8是CFETR混杂运行模式下初始分布alpha粒子在一个慢化时间内损失粒子在LCFS上形成的局域热斑, 峰值约0.1 WM/M2, 虽然在装置第一壁热负荷安全阈值范围内, 但是给其他损失通道预留安全空间较小, 考虑到各种损失通道的粒子并不沉积在同一位置, 相关设计还需要更进一步的优化和分析[25].
3.2 波纹叠加NTM扰动后的粒子损失
仅考虑NTM扰动下的初始分布alpha粒子的损失全部发生在几个极向回弹周期内, 即使粒子跟踪时间为半个慢化时间, 如图6所示. 经分析损失粒子能量、螺距角和损失时间分布后确认是初始轨道损失, 即由于产生在靠近边界的粒子有限轨道宽度引起, 扫描NTM扰动幅度不影响该损失份额. NTM扰动叠加碰撞后的粒子损失份额仅比初始轨道损失稍大, 是由于新经典效应引起. 考虑磁场波纹、NTM扰动以及碰撞效应下的粒子损失份额在1.4%附近, 扫描不同NTM扰动幅度, 损失份额基本不变. 俘获粒子轨道跨越宽度与扰动幅度成正比, 由图6计算结果可知, NTM扰动和EP的主要分布没有覆盖波纹损失区, 粒子进入损失区数目不受NTM扰动影响.
高频不稳定性引起EP损失的物理图像是变化磁场感生出电场, 在波粒共振时的相位耦合下, EP始终感受到恒向的径向电场漂移, 从而增大输运. 与此不同的是, 低频扰动或静态磁场扰动引起的EP损失起因是背景磁场拓扑结构的改变和EP轨道漂移. 在没有扰动时, EP相空间轨迹会形成完整的KAM(Kolmogorov Arnold Moser)曲面, 其阻止EP跨越输运, 低频或者零频率磁扰动引起通行粒子损失的物理图像是KAM面破缺形成给共振岛, 其与损失边界相交, 从而损失在第一壁上, 或者相邻共振岛相互重叠形成的混沌区域使粒子轨道随机[26]. NTM扰动下俘获粒子的共振损失是极向回弹频率和环向进动频率接近, 从而和较低的扰动频率达到共振条件, 形成磁场漂移. 扰动频率为零时形成的共振称为几何共振, 扰动和EP没有交换能量. CFETR初始alpha粒子分布的俘获轨道运动频率并不满足与NTM的共振条件, 见图14. 展示这些物理图像的有力工具是动理学庞加莱图.
动理学庞加莱图是记录EP轨迹点满足
nζ− ωt=2πk 时的平面坐标, k是整数, 平面坐标可为(ψp,θ ), (Pζ,θ )或(E,θ ). 跟踪的粒子分布需要满足磁矩μ和ωPζ−nE 为固定值, 且螺距角相同符号[17]. 通过不稳定性扰动下的EP动理学庞加莱图可以观察到KAM面的破缺与否、共振岛位置和混沌区域宽度等信息. 图9展示的是通行粒子和俘获粒子在扰动幅度较大时的共振岛分布, 通行粒子区域的3/2 和2/1 主共振岛和其间的高阶共振岛相互重叠形成宽域的混沌区域, EP因朗道相混而分布展平, 但是芯部和边界仍然存在完整的KAM面, 阻止了粒子直接损失. 俘获粒子区没观察到KAM面破缺, 即使增大动理学庞加莱图的记录分辨率仍如此. 值得指出的是香蕉轨道内外部分的漂移方向不同, 在记录轨迹点的庞加莱图(Pζ,θ )平面中相互重叠, 不容易分辨, 分析时需要记录香蕉轨道转折点, 以此观察进动岛结构. NTM扰动幅度较小时, 上述相空间的岛状结构相互分立, 仅有局域的KAM面破缺. 通过动理学庞加莱图分析可推知, 无论NTM扰动幅度大小都不会引起alpha粒子的显著损失, 仅有局域的通行粒子剖面展平.3.3 不同磁场波纹幅度下的粒子损失
为了验证图6中增大NTM扰动幅度而粒子损失份额不变的原因, 增大装置整体波纹度, 进行图6中磁场波纹叠加NTM扰动并考虑碰撞下的粒子损失计算, 结果如图10所示. 在CFETR装置实际波纹度分布下, alpha粒子的损失份额并不随NTM扰动幅度增大而增大, 但增大波纹度从工程设计值到其9倍时, 粒子损失份额随波纹度增大而显著增大. 在波纹度较大时, 可观察到一定程度的损失份额随NTM扰动幅度增大而增大. 波纹度分布整体增大引起波纹损失区扩大, 使波纹损失区朝芯部扩展并覆盖更多的粒子. 以波纹磁阱区域为例, 不同波纹度分布对应的磁阱区域如图11所示, 随着波纹损失区扩展, 逐渐靠近NTM扰动和EP主分布区, NTM扰动幅度增大开始影响进入损失区的粒子数.
4. 讨 论
为了系统地讨论CFETR上alpha粒子经NTM和波纹扰动的损失图像, 本节采用EP相空间损失区域和粒子分布密度来评估. 初始分布的alpha粒子能量单值分布, 即决定粒子轨道的三个物理量中固定E, 可以得到一个二维平面(
Pζ,μB0/E ),B0 是磁轴处磁场强度, 为ORBIT中磁场强度单位. 在轴对称位形和能量守恒关系下可以在(Pζ,μB0/E )平面内得到一系列曲线分割的不同轨道类型区域, 如图12(a)所示, 其中T-C指约束俘获粒子轨道, T-L指损失俘获粒子轨道,P±−LP±−L 指损失通行粒子轨道, p±-C指约束通行粒子轨道, 下标的正负号指螺距角符号[17]. 该平面中还有一些特殊的轨道, 如土豆轨道和极向静止轨道, 因在本文分析对象中占比过小, 及不可见没有标记出. 在不同平衡位形和EP能量时的该平面轨道类型区域占比不同, 但拓扑结构基本不变. 图12(b)和图12(c)是初始alpha粒子分布中不同螺距角正负号的样本粒子在(Pζ,μB0/E )平面的分布. CFETR混杂运行模式的初始分布alpha粒子中有28%份额的俘获粒子, 值得注意的是, 俘获粒子区域在不同螺距角正负号下的粒子轨道仅仅是香蕉轨道的内外分支, 是同一轨道.图13是(
Pζ,μB0/E )平面内不同波纹扰动幅度下的波纹损失区域, 其中波纹磁阱俘获区由磁阱形成条件ε|sinθ|<qNδ 确定, 无碰撞波纹随机扩散区的确定依据是简化判据GWB(Goldston, White and Boozer)δs=[ε/Nπq]1.5/(ρLq′) , ε是当地的纵横比倒数,q′=dq/dr 是安全因子径向导数. 图13中的波纹磁阱俘获损失区主要分布在粒子轨道的俘获-通行转换边界和外边界处,同图11所示一致. 图13展示了增大整体波纹扰动幅度时, 波纹损失区逐渐扩大, 占据大部分俘获粒子分布, 同图10纵轴方向的计算结果一致. 当波纹度增大倍数较大时, 在俘获-通行边界附近的粒子才可能进入损失区. NTM扰动和碰撞效应可以使粒子在俘获和通行状态间转换, NTM扰动幅度越大, 越多的粒子受其扰动进入俘获粒子区域, 从而经历波纹损失, 即图10观察到的波纹度较大时粒子损失份额随NTM扰动幅度增大而增大.EP与背景等离子体的碰撞是影响很多损失渠道的关键因素. 波纹磁阱俘获损失需要较高的波纹度以形成局域二级磁阱, 现在托卡马克装置已较难大范围满足, 相空间的损失区域占比也较小. 波纹度较小时只会引起俘获粒子香蕉轨道转折点位置的周期振荡, 并不会损失. 更大的波纹扰动或者碰撞会让香蕉轨道之间解耦, 前者即GWB判据的物理图像, 其对EP能量和波纹扰动幅度有一定要求. 对于能量固定的初始分布alpha粒子, 碰撞引起香蕉轨道之间的解耦对碰撞率有一定要求, 即
ν∗/(Nq)2−ν∗/(Nq)0.5 , 其中ν∗=εωb/(2π) . 螺距角散射引起香蕉粒子转折点相位随机, 对应径向位移随机震荡, 需要俘获粒子一次极向周期内经历足够多的碰撞次数, 经历一次碰撞期间经历了数个波纹扰动起伏, 连续轨道之间完全解耦[23]. 初始alpha粒子的香蕉轨道周期运动频率图如图14所示, 取极向回弹频率在ωb∼20 —40kHz , 判别频率ν∗/(Nq)2∼0.5 —4.0s−1 , 远高于初始分布alpha粒子的螺距角散射率, 因此无碰撞波纹随机扩散是CFETR初始分布alpha粒子的主要损失通道. 随着alpha粒子能量慢化, 粒子碰撞率增大, 碰撞效应在波纹损失过程中有更大贡献.5. 结 论
本文对CFETR上磁场波纹和NTM扰动叠加下的初始分布alpha粒子损失进行了详细的数值模拟, 结果显示粒子损失份额并不随NTM扰动幅度增大而增大, 与其他装置的结论不同. 原因是CFETR运行模式下的波纹损失区域较小, 没有覆盖EP和NTM扰动主要分布区. 人为增大装置波纹度才会观察到粒子损失份额随NTM扰动增大而增大, 原因是波纹损失区扩大靠近俘获-通行边界时, NTM扰动和有限碰撞引起边界附近粒子的轨道状态转换, 使粒子不断地从通行粒子区域进入俘获粒子区域. CFETR工程实际的波纹度较小, 波纹损失区离俘获-通行边界较远, 不足以发生粒子持续的跨过转换边界进入损失区.
通过动理学庞加莱图分析可知, NTM没有和俘获粒子形成有效共振, 因此俘获粒子剖面没有坍塌, 俘获粒子的分布剖面外延、展平会显著增大粒子损失, 即更多的粒子进入波纹损失区, 本文工作解除了这种担心. NTM和通行粒子形成了一系列的共振岛, 但远离边界并不会引起直接粒子损失, 只有局域的再分布. 值得指出的是, 本文结论只对低频率、局域且远离边界的不稳定性模式有效. 聚变堆上高频不稳定性下的alpha粒子约束研究需要更多工作, 扰动频率接近EP漂移运动频率时更容易形成共振.
[1] Zhuang G, Li G Q, Li J, et al. 2019 Nucl. Fusion 59 112010
Google Scholar
[2] Wan Y X, Li J, Liu Y, et al. 2017 Nucl. Fusion 57 102009
Google Scholar
[3] Chen J, Jian X, Chan V, et al. 2017 Plasma Phys. Controlled Fusion 59 075005
Google Scholar
[4] Gorelenkov N N, Pinches S D, Toi K, et al. 2014 Nucl. Fusion 54 125001
Google Scholar
[5] Fasoli A, Gormenzano C, Berk H L, et al. 2007 Nucl. Fusion 47 S267
[6] Pinches S D, Chapman I T, Lauber Ph W, et al. 2015 Phys. Plasmas. 22 021807
Google Scholar
[7] White R B, Rutherford P H, Colestock P, et al. 1988 Phys. Rev. Lett. 60 2038
Google Scholar
[8] Chapman I T 2011 Plasma Phys. Controlled Fusion 53 013001
Google Scholar
[9] Igochine V 2015 Active Control of Magneto-hydrodynamic Instabilities in Hot Plasmas (Berlin: Springer) p259
[10] Poli E, Garcia-Munoz M, Fahrbach H, et al. 2008 Phys. Plasmas 15 032501
Google Scholar
[11] Gobbin M, Marrelli L, Fahrbach H, et al. 2009 Nucl. Fusion 49 095021
Google Scholar
[12] Strumberger E, Gunter S, Schwarz E, et al. 2008 New J. Phys. 10 023017
Google Scholar
[13] Garcia-Munoz M, Fahrbach H, Pinches S D, et al. 2009 Nucl. Fusion 49 085014
Google Scholar
[14] Garcia-Munoz M, Martin P, Fahrbach H, et al. 2007 Nucl. Fusion 47 L10
Google Scholar
[15] Mynick H E 1993 Phys. Fluids B 5 1471
Google Scholar
[16] 郝保龙, 陈伟, 蔡辉山, 等 2020 中国科学: 物理学 力学 天文学 50 065201
Google Scholar
Hao B L, Chen W, Cai H S, et al. 2020 Sci. Sin. Phys. Mech. Astron. 50 065201
Google Scholar
[17] White R B 2014 The theory of toroidally confined plasmas (3th Ed.) (Singapore: World scientific publishing company) p73
[18] Pankin A, McCune D, Andre R, et al. 2004 Computer Physics Communications 159 157
Google Scholar
[19] 高翔, 万宝年, 宋云涛, 等 2019 中国科学: 物理学 力学 天文学 49 045202
Google Scholar
Gao X, Wan B N, Song Y T, et al. 2019 Sci. Sin.-Phys. Mech. Astron. 49 045202
Google Scholar
[20] Yu Q, Gunter S, Scott B D 2003 Phys. Plasmas. 10 798
[21] Wang X J, Yu Q, Zhang X D, et al. 2018 Nucl. Fusion 58 016045
Google Scholar
[22] Hao B L, White R B, Gao X, et al. 2019 Nucl. Fusion 59 076040
Google Scholar
[23] Wu B, Hao B L, White R B, et al. 2017 Plasma Phys. Controlled Fusion 59 025004
Google Scholar
[24] Boozer A H, Kuo-Petravic G 1981 Phys. Fluids 24 851
Google Scholar
[25] Zhao R, Wang Z X, Wang F, et al. 2020 Plasma Phys. Controlled Fusion 62 115001
Google Scholar
[26] Carolipio E M, Heidbrink W W, Forest C B, et al. 2002 Nucl. Fusion 42 853
Google Scholar
期刊类型引用(3)
1. 吴相凤,王丰,林展宏,陈罗玉,于召客,吴凯邦,王正汹. CFETR参数下α粒子慢化过程的数值模拟. 物理学报. 2023(21): 187-195 . 百度学术
2. 郝保龙,李颖颖,陈伟,郝广周,顾翔,孙恬恬,王嵎民,董家齐,袁保山,彭元凯,石跃江,谢华生,刘敏胜,ENN TEAM. EXL-50U球形环中快离子磁场波纹损失的优化模拟研究. 物理学报. 2023(21): 251-264 . 百度学术
3. Huishan Cai,Ding Li. Recent progress in the interaction between energetic particles and tearing modes. National Science Review. 2022(11): 240-261 . 必应学术
其他类型引用(1)
-
表 1 CFETR与其他托卡马克装置主机参数对比
Table 1. Main parameters comparison of CFETR and other tokomak facilities.
参数 CFETR ITER EAST 磁轴场强BT0/T 6.5 5.3 2 等体大半径R0/m 7.2 6.2 1.9 等体小半径a/m 2.2 2.0 0.5 等体电流Ip/MA 14 15 1 纵场磁体柄数N 16 18 16 -
[1] Zhuang G, Li G Q, Li J, et al. 2019 Nucl. Fusion 59 112010
Google Scholar
[2] Wan Y X, Li J, Liu Y, et al. 2017 Nucl. Fusion 57 102009
Google Scholar
[3] Chen J, Jian X, Chan V, et al. 2017 Plasma Phys. Controlled Fusion 59 075005
Google Scholar
[4] Gorelenkov N N, Pinches S D, Toi K, et al. 2014 Nucl. Fusion 54 125001
Google Scholar
[5] Fasoli A, Gormenzano C, Berk H L, et al. 2007 Nucl. Fusion 47 S267
[6] Pinches S D, Chapman I T, Lauber Ph W, et al. 2015 Phys. Plasmas. 22 021807
Google Scholar
[7] White R B, Rutherford P H, Colestock P, et al. 1988 Phys. Rev. Lett. 60 2038
Google Scholar
[8] Chapman I T 2011 Plasma Phys. Controlled Fusion 53 013001
Google Scholar
[9] Igochine V 2015 Active Control of Magneto-hydrodynamic Instabilities in Hot Plasmas (Berlin: Springer) p259
[10] Poli E, Garcia-Munoz M, Fahrbach H, et al. 2008 Phys. Plasmas 15 032501
Google Scholar
[11] Gobbin M, Marrelli L, Fahrbach H, et al. 2009 Nucl. Fusion 49 095021
Google Scholar
[12] Strumberger E, Gunter S, Schwarz E, et al. 2008 New J. Phys. 10 023017
Google Scholar
[13] Garcia-Munoz M, Fahrbach H, Pinches S D, et al. 2009 Nucl. Fusion 49 085014
Google Scholar
[14] Garcia-Munoz M, Martin P, Fahrbach H, et al. 2007 Nucl. Fusion 47 L10
Google Scholar
[15] Mynick H E 1993 Phys. Fluids B 5 1471
Google Scholar
[16] 郝保龙, 陈伟, 蔡辉山, 等 2020 中国科学: 物理学 力学 天文学 50 065201
Google Scholar
Hao B L, Chen W, Cai H S, et al. 2020 Sci. Sin. Phys. Mech. Astron. 50 065201
Google Scholar
[17] White R B 2014 The theory of toroidally confined plasmas (3th Ed.) (Singapore: World scientific publishing company) p73
[18] Pankin A, McCune D, Andre R, et al. 2004 Computer Physics Communications 159 157
Google Scholar
[19] 高翔, 万宝年, 宋云涛, 等 2019 中国科学: 物理学 力学 天文学 49 045202
Google Scholar
Gao X, Wan B N, Song Y T, et al. 2019 Sci. Sin.-Phys. Mech. Astron. 49 045202
Google Scholar
[20] Yu Q, Gunter S, Scott B D 2003 Phys. Plasmas. 10 798
[21] Wang X J, Yu Q, Zhang X D, et al. 2018 Nucl. Fusion 58 016045
Google Scholar
[22] Hao B L, White R B, Gao X, et al. 2019 Nucl. Fusion 59 076040
Google Scholar
[23] Wu B, Hao B L, White R B, et al. 2017 Plasma Phys. Controlled Fusion 59 025004
Google Scholar
[24] Boozer A H, Kuo-Petravic G 1981 Phys. Fluids 24 851
Google Scholar
[25] Zhao R, Wang Z X, Wang F, et al. 2020 Plasma Phys. Controlled Fusion 62 115001
Google Scholar
[26] Carolipio E M, Heidbrink W W, Forest C B, et al. 2002 Nucl. Fusion 42 853
Google Scholar
期刊类型引用(3)
1. 吴相凤,王丰,林展宏,陈罗玉,于召客,吴凯邦,王正汹. CFETR参数下α粒子慢化过程的数值模拟. 物理学报. 2023(21): 187-195 . 百度学术
2. 郝保龙,李颖颖,陈伟,郝广周,顾翔,孙恬恬,王嵎民,董家齐,袁保山,彭元凯,石跃江,谢华生,刘敏胜,ENN TEAM. EXL-50U球形环中快离子磁场波纹损失的优化模拟研究. 物理学报. 2023(21): 251-264 . 百度学术
3. Huishan Cai,Ding Li. Recent progress in the interaction between energetic particles and tearing modes. National Science Review. 2022(11): 240-261 . 必应学术
其他类型引用(1)
Catalog
Metrics
- Abstract views: 6258
- PDF Downloads: 75
- Cited By: 4