反磁剪切托卡马克等离子体中低频剪切阿尔芬波的理论研究.pdf
《反磁剪切托卡马克等离子体中低频剪切阿尔芬波的理论研究.pdf》由会员分享,可在线阅读,更多相关《反磁剪切托卡马克等离子体中低频剪切阿尔芬波的理论研究.pdf(16页珍藏版)》请在文库网上搜索。
1、专题:磁约束等离子体中的高能量粒子反磁剪切托卡马克等离子体中低频剪切阿尔芬波的理论研究*马瑞瑞1)陈骝2)3)仇志勇2)1)(核工业西南物理研究院,成都610041)2)(浙江大学物理系,杭州310027)3)(加州大学欧文分校物理与天文学系,加州CA92697-4574)(2023年 2月 21 日收到;2023年 6月 14 日收到修改稿)位于动理学热离子带隙附近的低频阿尔芬扰动因可以与高能量粒子或背景热粒子发生相互作用而引起了广泛的关注.本文在一般鱼骨模色散关系的理论框架下,针对反磁剪切托卡马克等离子体中观测到的由高能量粒子或者背景热粒子激发的低频剪切阿尔芬波的线性特性进行了一系列的理论
2、研究.由于这些低频剪切阿尔芬波与 2019 年 DIII-D 开展的专门研究高能量离子驱动低频不稳定的实验密切相关,因此本文通过采用 DIII-D 具有代表性实验的平衡参数,证明了实验上观测的低频模和比压阿尔芬本征模分别是以阿尔芬极化为主的反应型和耗散型不稳定模,因此,将前者称为低频阿尔芬模更准确.由于受逆磁和捕获粒子动理学效应的影响,低频阿尔芬模既可以在低频区(频率小于热离子的渡越或反弹频率)与比压阿尔芬声模耦合,又可以在高频区(频率大于或近似等于热离子的渡越频率)与比压阿尔芬本征模耦合.此外,由于受到不同激发机制的影响,与局域在安全因子最小值有理面附近的低频阿尔芬模相比,驱动比压阿尔芬本征
3、模的高能量离子的压强梯度在偏离安全因子最小值的有理面时达到最大值,相应的比压阿尔芬本征模的本征函数在高能量离子驱动最大的径向位置处出现峰值.通过改变安全因子最小值理论上重现了实验上观测的比压阿尔芬本征模和低频阿尔芬模的上升频谱特征.研究还表明,比压阿尔芬声模由于受到强烈的朗道阻尼,因而很难被高能量粒子激发,这与基于第一性原理的理论预测和模拟结果一致.本文证明了一般鱼骨模色散关系在解释和预测实验和数值模拟结果方面的强大能力.关键词:低频剪切阿尔芬波,高能量粒子,不稳定性,回旋动理学PACS:52.35.Bj,52.55.Pi,52.35.Py,52.30.GzDOI:10.7498/aps.72
4、.202302551引言高能量粒子(energeticparticle,EP)在聚变等离子体中的约束特性对决定当前和未来接近或在反应堆参数区间运行的托卡马克的性能具有重要意义.例如,为了点燃氘-氚(D-T)等离子体,必须对在 D-T 聚变反应中产生的阿尔法粒子进行良好的约束.同样地,在射频或中性束注入实验中产生的高能量离子也必须受到良好的约束才能成功地实现等离子体加热和/或电流驱动.尽管目前基于库仑碰撞理论对 EP 约束性的评估是可行的,但这一理论仍然存在严重的问题,问题主要来自于由 EP 激发的集体振荡所引起的 EP 的反常输运损失.这是由于携带着足够强的来自压强梯度自由能的 EP,其特征频
5、率(例如渡越/反弹频率)与剪*国家磁约束核聚变发展研究(批准号:2022YFE03040002,2018YFE0304103)、国家自然科学基金(批准号:12261131622,12175053)和四川省自然科学基金(批准号:2022NSFSC1814)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.21(2023)215207215207-1切阿尔芬波(shearAlfvnwave,SAW)的频率在同一量级(均在 MHz 范围内),并且 SAW 的群速度与磁场平行,因此通过波
6、-粒共振相互作用,可以直接激发多种 SAW,并引起 EP 的反常输运.在聚变装置中,由于磁场和等离子体的不均匀性,使得各种阿尔芬本征模(Alfvneigenmode,AE)存在于 SAW 连续谱的带隙中.这些 AE 因不受连续谱阻尼的影响,很容易被 EP 激发,因此通常被认为是 EP 输运的主要因素1.pitibi在这些由不同效应引起的 AE 中,其中由有限比压引起的阿尔芬本征模(beta-inducedAlfvneigenmode,BAE)存在于由热等离子体可压缩效应产生的间隙中2,3.自首次观测到 BAE 以来2,3,动理学热离子带隙频率4范围内的低频阿尔芬波谱一直是人们研究的热点.这些模
7、的频率()与热离子逆磁频率()同一量级,且与热离子渡越()和/或反弹()频率相当,因而可以与背景热粒子和 EP 发生相互作用1,59.由这种相互作用引起的有限电磁扰动和带状场结构会对相应粒子的输运过程产生正面或负面的影响1,4,10.这类低频 SAW 扰动包括但不局限于:考虑了有限逆磁漂移效应的动理学气球模(kineticballooningmode,KBM)1113、由热离子可压缩性而引起的BAE2,3、由动理学热离子可压缩性及波-粒共振效应引起的阿尔芬的离子温度梯度(Alfvniciontemperaturegradient,AITG)模6,14、以及由两个边带声模相互耦合再与 SAW 耦
8、合而产生的比压阿尔芬声模(beta-induced Alfvn acoustic eigenmode,BAAE)15,16.对于这些低频 SAW,由于热等离子体的可压缩效应会对其色散关系产生重要的修正,因此,低频SAW 的线性性质往往与已被广泛研究的环向阿尔芬本征模(toroidalAlfvneigenmode,TAE)17和反磁剪切阿尔芬本征模(reversedshearAlfvneigenmode,RSAE)18,19等高频阿尔芬本征模的线性性质不同.一般地,包含任何热粒子动理学效应的由高能量粒子激发的 Alfvn 波的不稳定性,均可以在由 Zonca 和 Chen20,21发展的一般类鱼
9、骨模色散关系(generalfishbone-likedispersionrelation,GFLDR)的理论框架下进行统一的描述.GFLDR 是类比磁流体能量原理的推导方法,将描述波色散关系的式子表示为类能量形式的方程.它首先要求扰动量在空间上存在两个径向尺度,然后采用模分解22和 WKB(Wentzel-Kramers-n=3m=48|bi|A|A0|mLPE;rel m本文从不考虑47以及考虑 EP48的效应来研究低频 SAW 的线性特征.在这种情况下,不同种类粒子的动理学效应将通过其压强特征尺度线性地进入低频 SAW 的色散关系并影响其行为.对于DIII-D 第 178631 次放电,
10、图 3 所示为热粒子和高能量粒子的压强特征尺度(和 ),以及在弱和/或零磁剪切情况下()模的宽度()的径向依赖关系.这里,r 为沿着小半径方向的径向位置.具体而言,EP 的压强分布由“弛豫的(relaxed)”EP 剖面和“经典的(classical)”EP 剖面这两种极限情况给出.“弛豫的”EP 剖面由 EFIT49重建提供,其中快离子的压强是通过计算平衡压强与热粒子压强之差得到的.另一种由TRANSP/NUBEAM50给出,是指在考虑由不稳定性引起的快离子输运的情况下得到的“经典的”EP 分布.对于这两种极限情况,EP 的压强特征尺度分别由 和 表示.当模失稳时,真实的 EP分布可能介于这
11、两个极限之间,且最接近基于 EFIT获得的结果.在弱和/或零磁剪切区域,对于给定的环向和极向模数 ,其归一化平行波矢,模的径向度可以用 45,46来估算.这里,表示在 处的平行波矢,其中 q 剖面具有最小值并由 给出,表示 q 在径向上的二次导数,R 是托卡马克的大半径.由图 3 可知,在弱和/或零磁剪切的区域,.这一关系通常对应于将全局问题简化为求解模的局域色散关系的前提条件.因此,对于与高能量粒子激发无关的 LFAM 的分析,可以通过求解低频 SAW 的局域色散关系来处理47.然而,对于高能量离子驱动的 BAE,需要讨论两种不同的情况48:一种是在适中的 EP 压强梯度下,即仍有 的关系,
12、这种情况可近似通过求解局域的 GFLDR5,20,21,25,4346来研究 EP激发的 BAE 的性质;另一种是当系统中存在强的EP 压强梯度时,则有 ,对于这种情况,需要求解低频 SAW的全局色散关系.关于这部分的讨论将在理论模型中具体展开.80010001200140020050100150Frequency/kHzTime/ms#178631RSAEBAELFM(2,1)(5,3)(6,4)(3,2)(13,9)(14,10)(15,11)(14,11)(13,10)(12,9)(11,8)(10,7)(9,6)(7,5)(8,6)(9,7)(10,8)(11,9)(12,10)(13
13、,11)(7,5)(8,6)(9,7)(10,8)(8,7)(7,6)(6,5)(5,4)(4,3)(a)8001000120014001.61.22.0minTime/ms(b)qminm/n图2Heidbrink 等34对 DIII-D 参考炮数#178631 分析的实验结果(a)布局在大半径 R=192201cm 之间的ECE 通道的互功率谱图;(b)EFIT 重建得到的 与时间的关系.其中频谱图上所示的不同符号分别对应不同 值的不同模式:RSAE(),BAE()及 LFAM()qminm/nFig.2.The DIII-D experimental results from Ref.3
14、4 byHeidbrinketal.:(a)Cross-powerspectrogramintherefer-enceshotforECEchannelsbetween192201cm;(b)mea-sured fromEFITreconstructionsvs.time.TheRSAE(),BAE(),andLFAM()symbolsrepresenttheval-uesof shownonthespectrogram.物理学报ActaPhys.Sin.Vol.72,No.21(2023)215207215207-4LPm本文的结构如下:首先,第 2 节讨论弱磁剪切和/或零磁剪切附近低频 S
15、AW 的局域和全局色散关系,这些色散关系的选取取决于 和 的相对大小;紧接着,第 3 节详细讨论不包含及包含EP 效应的低频 SAW 的数值结果和相应的理论分析,并进一步与实验结果进行比较;最后第 4 节对本文进行总结和讨论.2理论模型LPthLPEmLPth mLPE;rel mLPE;clm本节将讨论具有弱反磁剪切 DIII-D 放电中低频 SAW 的色散关系.正如第 1 节所讨论的,我们将依据 及 与 的相对大小,讨论以下两种情况来研究低频 SAW 的稳定性.情况:不考虑 EP 的作用()以及考虑适中 EP 梯度()情形下,描述低频 SAW 的局域GFLDR 模型;情况:考虑强EP 压强
16、梯度()情形下,适用于描述低频SAW 的全局GFLDR模型.2.1 局域模型LPth mLPE;rel m对于情况:不考虑 EP 的作用()以及考虑适中 EP 梯度()情形下,背景热离子的压强特征尺度及 EP 的特征压强尺度均可以与模结构的径向宽度很好地分开,因此描述 SAW 的涡量方程1,5,20,21,43,44可以简化为通常局域极限下的 SAW 色散关系.关于这一点,已有大量文献进行了详细的推导和讨论1,20,21,25,43,44,此|bir0s=0S (r/q)q1/2处将不再赘述.我们只需注意,针对 DIII-D 的实验,在应用局域 GFLDR 时需要考虑两个关键点.首先,对于描述
17、 的低频 SAW,需要考虑热粒子的可压缩效应,即考虑包括波-粒子(包含通行和捕获的热粒子)的相互作用8;其次,还需处理 DIII-D 芯部具有反磁剪切的位形特征.因此,对于局域在 处,磁剪切 但 具有有限值的低频 SAW 的局域 GFLDR 可表示为20,21,25,46,51iS(2n k2/n0q2minR2)1/2(1/n)1/2k/n0qminR i(2n k2/n0q2minR2)1/21/2=Wnf+Wnk(),(1)n()n()WnfWnk|2n|1WnfWnf其中,为广义的惯性项,它既包括逆磁效应,也包括通行热粒子和捕获热粒子的动理学效应.关于 的具体推导请参考文献8.这里将主
18、要结果总结于附录 A 中以方便读者阅读.此外,方程(1)的右边包含了对理想区域势能的“流体”()和“动能”()贡献.对于低频模(),与频率无关,考虑圆截面近似并采用(s,)模型52,则 的表达方程可表示为Wnf4(S2k/0qminRn322S?k/0qminRn?1/2+9324),(2)=c+Ec=Rq2mindc/drE=12Rq2mind(E/+E)/drcE/EWnkF0E=/F0E=B0E(r)2522mEb(1 0B0)3/2(0)mE0E(r)8PE(r)/B20PEB0(x)=2/2 bb其中,这里 为背景热粒子的比压,和 分别为平行于磁场方向和垂直于磁场方向的高能量离子的比
19、压.方程(2)中包括高能量粒子的绝热和对流响应的贡献53.反映了共振和非共振的波-粒相互作用,因而是模频率的函数.简单起见(但仍与 DIII-D 情况相关),我们假设高能量离子的平衡分布函数 遵从具有单一俯仰角()的慢化分布,即 .其中,是高能量离子的质量,是初始俯仰角,是热压()和磁压的比值,是磁轴上平衡磁场,是狄拉克函数,是磁矩,是单位质量的 EP 初始能量.那么,高能量通行离子的非绝热贡献的表达方程可写为43,440.200.240.280.3200.51.01.5th(m)E;rel(m)E;cl(m)(m)LPthLPEs=rq/qm图3热粒子和高能量粒子压强特征尺度(和 ),以及在
20、弱和/或零磁剪切()情况下模宽度()的径向依赖关系LPthLPEmFig.3.Radial dependences of the typical scale lengths ofthermalandenergeticparticlepressure(and ),magnetic shear(s)as well as the estimated radial modewidth().物理学报ActaPhys.Sin.Vol.72,No.21(2023)215207215207-5WnkuE25/2(1 0B0/2)2 ln(+1 1),(3)=/tEmtEm2b/qR其中 ,是粒子能量最大时 E
21、P 的渡越频率.r=r0k/n0qminRSf(iE/k/)a.c./d.c.a.c.d.c.Sf|Sf|Sf|这里需要注意的是,方程(1)的中括号的第一项代表在 处有限的 ,它来自有限的场线弯曲效应,并在方程中起着重要的稳定作用20,21,25.此外,通过 来定义模的极化.这里的 和 分别对应于平行电场、波矢和标势扰动的正弦分量和几乎恒定的分量.的详细表达方程也在附录 A 中给出.根据文献 8,28,51 中的讨论,的值远远小于 1,表示模具有 SAW 的极化特征;而 的值与 1 可比拟或大于 1 则代表模具有离子声波的极化特征.当粒子压强特征尺度与模的径向宽度可以很好区分开时,我们得到了描
22、述低频 SAW 的局域GFLDR,即方程(1).进一步,可以借助数值方法研究方程(1)中涉及的各项,来描述实验观测和模拟中低频 SAW 的基本物理.然而,在强 EP压强梯度存在的情况下,如上述讨论的情况,此时描述低频 SAW 的局域 GFLDR 将不再适用,需要借助全局的 GFLDR 来处理.2.2 全局模型LPE;cl mLPE;clm针对情况,即 ,此时两个特征尺度 和 不再能很好地区分,即 EP 的特征平衡参数相对于模的宽度不再是缓变量,因此,模在径向局域的位置将受高能量离子分布的影响,因而需要全局的色散关系来描述低频 SAW 的性质.要解决这个问题,需要采用傅里叶逆变换将涡量方程变回至
23、实空间46,相应地,描述由单个环向和极向模数(m,n)主导的低频 SAW 的全局色散关系为(e er)2 2A,m(1+x2A,m+x442A,m)(e er)m(F+K)m=0,(4)k/k=(e er)erex2=nqmin(r r0)2(i/n1/2)S(/x)m其中 ,这里的 和 分别为径向和极向的单位矢量,是标量场扰动的 m 次极向谐波.值得注意的是,不同极向谐波之间的环向耦合对于局域在反磁剪切区域的模方程通常并不重要.此外,方程(4)中的 F 和 K 项分别表示类流体粒子和高能量离子的贡献,其表达方程为F D2S 42DS+2D2S(+1)+23,K 2q2Eq2R2mEc22dE
24、QF0E2tE 2=2Wnku,(5)DS=SA,m/nqEdE=(2E/2+2E/)/cER tE=E/qR QF0E=(+E)F0E EF0E=1cE(kb)F0EcE=qEB/mEc()=d3()/b其中 ,为高能量离子的电荷.,代表对速度空间 的积分,下标 和 分别表示相对于平衡磁场 的平行和垂直分量.E()=c11 (c2)2/c23c1=0.7099c2=0.3018c3=0.2944=c2=0.3018qminE求解常微分方程(4)一般需要借助数值方法.然而,对于 DIII-D 的情形,如图 4 中的黑色曲线所示的“经典”分布下高能量离子归一化压强梯度的径向关系,它可以用红色曲线
25、所示的解析方程 来很好地拟合.这里,.这一方法使得我们可以解析得到描述低频 SAW 的全局色散关系.需要注意的是,高能量离子的最大驱动位于 附近,它偏离了 的径向位置.此外,方程(3)中 为 r 的函数:E(r)=aE01(r r0+b)22cL2PE;cl,(6)a=c1/E0b=r0 c2ac=c3a/LPE;clE0LPE;clr=r0其中 ,a 为等离子小半径,和 表示在 的值.x=r r0=z b进一步,引入变量代换 ,则方程(4)可进一步简化为0.200.250.300.350.620.640.660.680.700.72E;cl=/E;cl fitPrediction bound
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 剪切 马克 等离子体 低频 阿尔芬波 理论研究