1. 引言
随着社会的发展和人类生活水平的提高,对能源的需求量逐渐增大,目前现存的化石能源已经不能满足人们对能源的需求量,因此寻找新能源已经成为能源工作的重点。而核聚变能因其产物清洁无污染且燃料储量丰富,被人们认为是最有可能“一劳永逸”解决能源问题的新能源。由于目前的大型核聚变反应堆,如国际热核巨变实验堆(International Thermonuclear Experimental Reactor, ITER),聚变商业示范堆(Demonstration, DEMO)等的燃烧等离子体都是由氢同位素组成,因此阐明氢同位素对离子温度梯度模(Ion temperature gradient mode, ITG)和捕获电子模(Trapped electron mode, TEM)等微观不稳定性,以及湍流输运和带状流的影响,对提高等离子体约束性能至关重要 [1] [2] [3] 。
目前在托卡马克(Tokamak)中已经有大量实验表明氢同位素离子质量对微观湍流和能量约束特性有很大的影响 [4] [5] [6] 。Bustos等人研究发现,同位素离子质量对托卡马克中TEM模的约束有明显的影响,而温度分布对同位素离子效应的影响取决于它们的梯度 [7] 。在托卡马克中,最近几项关于ITG-TEM驱动湍流模拟的研究表明,剪切旋转、快离子压力、电磁波动等现象对同位素离子质量有明显的依赖性 [6] [8] [9] 。
仿星器(Stellarator)是一个极具吸引力的聚变反应堆概念,它不需要内部的等离子体电流驱动,就可以有效地避免由电流引起的大破裂和电流驱动的不稳定性,从而实现稳态运行。近二十年来,随着仿星器装置设计水平和施工水平的提升,仿星器概念再一次受到了重视,并且从上世纪90年代开始,陆续设计建造了一些新的仿星器装置,如LHD (Large helical device),W7-X (Wendestein 7-X)等 [10] [11] 。特别是当2017年3月在LHD中首次进行了氘(D)实验之后 [12] [13] ,出现了一系列实验和数值模拟研究了氢同位素对微观不稳定性和湍流输运的影响。比如,Nakata等人使用GKV (Gyrokinetic Vlasov)代码研究了LHD中氢同位素等离子体对离子尺度和电子尺度的微观不稳定性以及TEM驱动的湍流和带状流的影响,发现ITG增长率表现出一种类似回旋玻姆标度的离子质量依赖性,TEM增长率对氢同位素离子质量有较强的依赖性,并且相比于氢(H)等离子体,氘等离子体中的湍流输运水平更低 [9] [14] 。同时,他们也发现相较于氢等离子体,氘等离子体的离子和电子热扩散系数更低,这表明同位素效应对等离子输运有显著的影响 [15] 。Ida等人发现同位素混合与湍流特性之间存在很强的相关性,并证明了等离子体中非混合剖面的形成以及向同位素混合状态的转变与离子沿抗磁方向传播的湍流水平的增加有关 [16] [17] 。Kobayashi等人首次清晰地观察到自组织内部输运垒(Internal transport barrier, ITB)结构中的同位素效应,发现氘等离子体中的ITB强度比氢等离子体中的ITB强度更强 [18] 。
从以上研究中可以发现在LHD中,氢同位素对微不稳定性和湍流输运的研究已经有了很大的进展,但在W7-X中很少有关于氢同位素影响的研究。之前对仿星器中微观不稳定性的非线性回旋动力学模拟都是在通量管和径向局部模区域中进行的 [19] [20] ,但是通量管模拟不能很好地处理仿星器内部三维磁场引起的多个环向模的耦合,因此在模拟中使用全局模拟是必要的。近几年,在托卡马克和仿星器中,GTC (Gyrokinetic toroidal code)代码已经被多次用于微观湍流不稳定性的全局模拟研究,其准确性已经得到了大量的验证 [21] [22] [23] 。因此,本文首次使用全局回旋动理学模拟研究了仿星器W7-X中氢同位素对微观湍流不稳定性的影响。在线性模拟中证明了ITG模归一化增长率符合回旋玻姆标度的离子质量依赖性,并且发现通量管模拟和全局模拟在模拟ITG增长率时几乎没有差异。非线性模拟中发现带状流对湍流具有抑制作用,且氢同位素离子质量越重,带状流对其湍流的抑制越强。
本文组织如下,第二部分介绍了W7-X仿星器三维磁场位型与模拟使用的回旋动理学模型,第三部分介绍了仿星器W7-X中氢同位素对微观湍流不稳定性影响的线性模拟结果,第四部分介绍了仿星器W7-X中氢同位素对微观湍流不稳定性影响的非线性模拟结果,第五部分为总结。
2. 模拟模型与仿星器位型
2.1. 模拟模型
在W7-X模拟中,等离子体的微观粒子在相空间中的分布函数满足的回旋动理学方程如下 [24] :
(1)
其中,
(2)
(3)
,
,
和t分别为离子回旋中心位置,磁矩,平行速度和时间。f是平衡分布函数,m是离子质量,Z为离子电荷。
和
分别为
漂移速度和磁漂移速度。B是平衡磁场,
,其中
,
是回旋频率。在本文中保留了ITG湍流产生的带状流,忽略了非轴对称系统中与新经典双极相关联的平衡径向电场。
模拟中采用
方法降低粒子噪声,提高计算效率 [25] 。离子回旋中心分布函数
可以分为两部分,第一部分为平衡部分
,第二部分为扰动部分
。式(1)可以简化为
,其中L为传播子,它也可以被分为平衡部分
和扰动部分
两部分。则式(1)可以表示为
(4)
其中
(5)
(6)
平衡扰动部分
满足
(7)
则扰动部分
满足
(8)
定义粒子权重为
,可以得到粒子权重方程
(9)
要求解式(9),需要
的解析式,但是新经典解
通常不可用,因此需要数值求解,我们将式(9)中
近似为局部麦克斯韦量。由于新经典项
没有出现在式(9)中,在本工作中不考虑新经典效应。
进一步对粒子扰动分布进行积分,即可得到粒子的扰动密度
,带入泊松方程式(10)求解,即可得到扰动静电势:
(10)
二次回旋平均电势
可以表示为
(11)
回旋平均电势
可以表示为
(12)
其中
是回旋半径频率。
在GTC模拟中,求解场方程与推动粒子都是在Boozer磁面坐标系
里完成的,平衡磁场满足
,在磁面坐标系下粒子的运动方程为 [21] [26]
(13)
(14)
(15)
(16)
以上各式中,
,极向流和环向流的径向导数分别为
,
。平行方向正则动量为
,
,
。值得注意的是,以上式子中含
的项为非线性项,只在非线性模拟中添加。
2.2. 仿星器位型和模拟设置
在GTC中一般采用数值方法构建复杂平衡位形,EFIT,VMEC,M3D-C1以及LR_eqMI是GTC中常用的几种构建平衡位型方法 [21] [22] [27] 。在本文中采用VMEC方法构建仿星器W7-X的三维平衡位形。VMEC方法所需要的平衡数据包括极向电流
,环向电流
,和磁场安全因子
,并且假设磁面是闭合的 [22] 。平衡磁场与平衡位形可以使用傅里叶级数在极向和环向的离散径向网格上展开,且在环向上通量等距。然后在二维极向平面上的离散网格点上,将这些左手坐标系中的三维数据转换为布尔坐标中的傅里叶级数,转换方式如下:
(17)
其中,
是Boozer坐标系中归一化的极向磁通量,
是极向,
环向环向角。
,
分别代表傅里叶余弦级数和正弦级数,n为环向模数。Boozer坐标
和柱坐标
也可以通过类似式(17)的方式进行转换,其中
是柱坐标中的环面角,R和Z是极向平面上的笛卡尔坐标。
仿星器W7-X的平衡磁场在环向方向的场周期为5,螺旋角
,这意味着所有的平衡量,包括磁场强度与度规张量均在沿环向旋转
后保持一致。本文中仿星器W7-X平衡磁场位型和等离子体分布参考Riemann的文章 [28] 。图1为仿星器W7-X中电子温度,离子温度,密度以及旋转变换
随径向位置的变化。其中红色实线表示旋转变换
,它通常用来表示磁场螺距,且
满足
,q是环形等离子体常用的安全因子,用于表示磁力线在磁面上单位极向角变化时所对应的环向角变化的比值。W7-X模拟中设置的径向区间为内边界
,外边界
。模拟中诊断磁面的径向位置为
,在诊断磁面处
,
,
,
,
,
,
。模拟中时间步长设置为
,其中
为声速,定义为
。在这里需要注意,模拟中使用的电子模型为绝热电子模型,离子使用回旋动理学离子,在本文中主要研究氢(H),氘(D),氚(T)三种氢同位素燃烧等离子体对ITG模的影响。
Figure 1. Radial profiles of equilibrium density, temperatures, and rotation transform
图1. 平衡密度,温度,以及旋转变换
的径向分布
在模拟中,电子温度,电子密度,与离子密度均假设为不随径向位置变化的常量。当离子碰撞到模拟设置的径向边界时,粒子仍旧被放回模拟区域内,但其比重被设置为0,这样处理可以保证运动离开模拟区域的粒子不会影响泊松方程的求解,与此同时,总的粒子数依旧保持守恒。在径向边界处,离子温度梯度依旧不为0,但在径向内边界处都有缓冲区,可以保证在边界处的扰动电势光滑地降为0。经数值收敛测试后,径向网格和极向网格数分别为110和2500,每个网格上的粒子数为150。
3. W7-X中ITG模的线性模拟
仿星器W7-X的环向场周期为5,这意味着在W7-X的模拟中存在五个线性本征模组,所对应每个本征模组的环向模数分别为
,k是正整数,
。每个本征模组由无数个本征模构成,例如,在
时,第一个本征模组由
等环向模构成。因为仿星器W7-X具有环向周期性,因此为了简便,我们在模拟中只模拟单个环周期,对应的环向角范围为
。
图2为H,D,T三种情况下,W7-X中ITG模的本征模结构在诊断磁面上的分布。在垂直于磁力线方向,本征模结构存在的区间比较窄,而在平行于磁力线方向,本征模则延展开来。H,D,T三种情况下在垂直于磁力线方向和平行于磁力线方向的本征模结构存在的区间基本一致。
Figure 2. The eigenmode structure of ITG mode in W7-X on the diagnostic magnetic surface for the H (a), D (b), and T (c) plasmas, respectively
图2. 当燃烧等离子体为H (a),D (b),T (c)时,W7-X中ITG模的本征模结构在诊断磁面上的分布
图3为H,D,T三种情况下,W7-X中ITG模的本征模结构在极向剖面上的分布。可以看出W7-X中ITG模的极向模结构只存在于较小范围的磁力线上,在极向剖面上的模结构也是由几块分散分布的结构组成,这与LHD的模结构有所不同,这样分布的模结构与W7-X环向类似磁镜的平衡磁场有关:平衡磁场在环向变化剧烈,这样会导致更多的环向模耦合在一起形成了本征模结构。同时还可以发现,ITG模结构对同位素离子质量有显著的依赖性,即氢同位素离子质量越大,ITG模结构的丝越宽,这主要是由于氢同位素离子质量越大,其极向波长越长造成的。
Figure 3. The eigenmode structure of ITG mode in W7-X on the
poloidalplane for the H (a), D (b), T (c) plasmas, respectively
图3. 当燃烧等离子体为H (a),D (b),T (c) 时,W7-X中ITG模的本征模结构在环向角
处极向剖面上的分布
为了更准确的分析W7-X中氢同位素对ITG模结构的影响,对图2中ITG模在诊断磁面上的本征模结构进行二维傅里叶变换,变换公式如下:
(18)
可以得到诊断磁面上与环向模数和极向模数相关的二维谱分布。
图4为H,D,T三种等离子体情况下,W7-X中ITG模的本征模结构在诊断磁面上的二维谱,环向谱以及极向谱分布。从图中可以看出,H,D,T等离子体的最不稳定环向谐波分别为
,
,和
,对应最不稳定的环向谐波,最不稳定的极向谐波分别为
,
,
,可以发现对于H,D,T等离子体,最不稳定的环向谐波n与其对应的最不稳定的极向谐波m的比值均近似等于
,即
。另外,我们可以发现,H等离子体的不稳定环向谐波范围为[100, 200],而D,T等离子体的不稳定的环向谐波范围分别为[60, 170],[50, 150],这说明ITG模环向耦合谐波的范围对氢同位素离子质量有明显的依赖关系,即氢同位素离子质量越重,其ITG模环向耦合谐波越多。
Figure 4. 2D spectrum
, toroidal spectrum
, and poloidal spectrum
on the diagnosed flux surface of ITG eigenmode, for H(a, d, g), D (b, e, h), and T (c, f, i) plasmas in the W7-X, respectively
图4. 当燃烧等离子体为H (a, d, g),D (b, e, h),T (c, f, i)时,W7-X中ITG模的模结构在诊断磁面上的二维谱分布
,环向谱分布
,极向谱分布
图5为H,D,T三种等离子体情况下,ITG模的归一化增长率
随环向谐波n的变化,需要注意的是在这里对增长率
和环向谐波n进行归一化时,使用的是每种同位素物种的离子质量,从图中可以看出,随环向谐波n的变化,三种等离子体的ITG模的归一化增长率
分布几乎完全相同,这符合回旋玻姆标度的离子质量依赖关系,即
。这与在通量管区域中的模拟结果基本一致,我们使用径向全局模拟再一次证明了ITG模归一化增长率类似回旋玻姆标度的离子质量依赖性,并且发现通量管模拟和径向全局模拟在模拟ITG增长率时几乎没有差异。
Figure 5. The normalized growth rate
of ITG mode as a function of toroidal harmonics n, for the H, D, and T plasmas in the LHD, respectively
图5. 当燃烧等离子体为H,D,T时,ITG模的归一化增长率
随环向谐波n的变化
4. W7-X中ITG模的非线性模拟
在非线性模拟中,使用与线性模拟相同的平衡位形,与等离子体剖面,但在非线性模拟中需在粒子导心运动方程中加入非线性项。其次为了降低粒子噪声,在非线性模拟中,每个格点上的粒子数增加到200,时间步长与其他格点数和线性模拟一致。非线性模拟中依旧只模拟单个环周期,对应的环向角范围为
。
图6给出了H,D,T三种等离子体情况下,仿星器W7-X中ITG模保留与抑制带状流时ITG模的本征模结构在极向剖面上的分布,右侧彩条代表扰动电势强度。将图3和图6比较可以发现,非线性阶
Figure 6. The eigenmode structure of ITG mode on the
poloidalplane, without and with the zonal flow, for the H (a, d), D (b, e), T (c, f) plasmasin W7-X, respectively
图6. 当燃烧等离子体为H (a, d),D (b, e),T (c, f)时,W7-X中抑制和保留带状流时ITG模的本征模结构在环向角
处极向剖面上的分布
段抑制带状流时,ITG模的本征模结构比线性阶段ITG模的本征模结构更宽,这主要是由于湍流先经过一段线性增长才达到非线性稳定的饱和阶段,而这段时间里湍流在不断地传播与扩散。将图6第一行与第二行比较可以发现,相比与抑制带状流,保留带状流时ITG模的本征模结构更宽,这主要是由于在非线性阶段并且保留带状流时,ITG模的本征模结构的径向分布被带状流切断,径向相关长度变短导致的。比较图6第一行,即抑制带状流时三种等离子体情况下的ITG模的扰动电势强度,发现他们的ITG模的扰动电势强度基本相同。但比较图6第二行,保留带状流时三种等离子体情况下的扰动电势强度,可以发现氢同位素离子质量较重时,其ITG模扰动电势强度较小,即W7-X中保留带状流时,ITG模扰动电势强度对氢同位素离子质量具有明显的依赖性。另外我们也可以发现,不管是保留还是抑制带状流时,ITG模的本征模结构在较重氢同位素离子质量情况下都会更宽。
图7为H,D,T三种等离子体情况下,仿星器W7-X中保留与抑制带状流时,ITG模模结构在诊断磁面上的环向谱分布,纵坐标为环向谱分布扰动电势的强度。比较图4和图7可以发现,相比于线性阶段,非线性阶段最不稳定的环向模耦合的谐波更多,即对应的环向不稳定范围更大。同时,在非线性阶段,W7-X中环向模的谱分布存在着逆级联现象(inverse cascade),即在线性阶段占主导地位的高环向模数与高极向模数的能量将在非线性阶段流入低环向模数。但是相比于仿星器LHD,W7-X中带状流的存在对逆级联的影响有所不同。在LHD中保留与抑制带状流时,环向谱分布既包含线性阶段占主导地位的高环向模,也包含着非线性阶段由逆级联产生的低环向模数。但在W7-X中,抑制带状流时与LHD相似,环向谱分布即包含线性阶段占主导的高环向模数,又包含非线性阶段产生的低环向模数,不同的是,在W7-X中保留带状流时,环向谱分布只包含着非线性阶段产生的低环向模数,高环向模扰动电势强度极低,这是仿星器W7-X的一个新的物理特性,产生这种现象的主要原因是W7-X中带状流在三维平衡磁场作用下通过线性耦合进一步提高了低环向模的幅度。比较图7第一行和第二行可以发现在H,D,T三种等离子体情况下,保留带状流时ITG模扰动电势强度都比抑制带状流时扰动电势强度更小,这说明带状流对湍流具有抑制作用。而分别比较图7第一行和第二行可以发现,抑制带状流时,H,D,T三种等离子体情况下,扰动电势强度基本相同,而保留带状流时,H,D,T三种等离子体情况下,扰动电势强
Figure 7. Toroidal spectrum
on the diagnosed flux surface of ITG eigenmode, without and with the zonal flow, for H (a, d), D (b, e), and T (c, f) plasmas in the W7-X, respectively
图7. 当燃烧等离子体为H (a, d),D (b, e),T (c, f)时,W7-X中抑制和保留带状流时ITG模的模结构在诊断磁面上的环向谱分布
度有很大差异,即氢同位素离子质量越重,其扰动电势强度越小,这主要是由于当氢同位素离子质量较重时,带状流对湍流的抑制越强而导致的,此结果与图6结果一致。
5. 总结
本文首次使用全局回旋动理学模拟研究了仿星器W7-X中氢同位素对ITG微观湍流不稳定性的影响。线性模拟表明,W7-X中ITG模的本征模结构在垂直于磁力线方向,本征模结构存在的区间比较窄,而在平行于磁力线方向,本征模则延展开来。H,D,T三种情况下在垂直于磁力线方向和平行于磁力线方向的本征模结构存在的区间基本一致,但ITG模结构的丝对氢同位素质量明显的依赖性。其次氢同位素离子质量越重,其ITG模耦合的环向谐波越少。同时我们使用全局回旋动理学模拟再一次证明了ITG模归一化增长率符合回旋玻姆标度的离子质量依赖关系,并且发现通量管模拟和径向全局模拟在模拟ITG增长率时几乎没有差异。非线性模拟表明保留和抑制带状流时,ITG模的本征模结构在较重氢同位素离子质量情况下都会更宽,且三种等离子体情况下,ITG模扰动电势强度在抑制带状流时基本相同,在保留带状时较重氢同位素离子,ITG模扰动电势强度会更小,表明带状流对湍流具有抑制作用,且氢同位素离子质量越重,对其湍流的抑制越强。
阐明氢同位素对微观湍流的影响对提高等离子体约束性能至关重要,本文主要阐明了W7-X仿星器中氢同位素对ITG微观湍流的影响,这是我们研究聚变装置中同位素效应这个关键物理问题迈出的第一步,之后我们将进一步研究氢同位素对捕获电子模和电子温度梯度模微观湍流的影响。
致谢
感谢美国加州大学林志宏教授以及GTC小组给予的帮助与支持。
参考文献
NOTES
*通讯作者。