负三角形变托卡马克位形下高能量离子激发鱼骨模的模拟研究.pdf
《负三角形变托卡马克位形下高能量离子激发鱼骨模的模拟研究.pdf》由会员分享,可在线阅读,更多相关《负三角形变托卡马克位形下高能量离子激发鱼骨模的模拟研究.pdf(8页珍藏版)》请在咨信网上搜索。
1、专题:磁约束等离子体中的高能量粒子负三角形变托卡马克位形下高能量离子激发鱼骨模的模拟研究*任珍珍1)申伟2)1)(安徽大学物理与光电工程学院,合肥230601)2)(中国科学院等离子体物理研究所,合肥230031)(2023年 4月 23 日收到;2023年 7月 17 日收到修改稿)负三角形变位形下的托卡马克具有更低的湍流输运和更好的能量约束,被认为是未来聚变堆一个更好的选择.为了探索负三角形变位形下高能量粒子激发不稳定性的特征,使用动理学-磁流体混合模型程序M3D-K 开展了此位形下高能量离子激发鱼骨模的线性不稳定性和非线性演化的模拟研究.基于类 EAST 参数条件,模拟发现负三角形变解稳
2、理想内扭曲模不稳定性,但会致稳鱼骨模不稳定性.非线性模拟发现在没有磁流体非线性效应时,负三角形变位形下的鱼骨模更不容易饱和,可能的解释是相比于正三角形变位形,在负三角形变位形下的高能量离子轨道更接近与芯部,因而更容易驱动鱼骨模不稳定性.这些结果表明考虑高能量粒子激发的鱼骨模不稳定性后,负三角形变位形相比于正三角形变位形并没有明显优势.关键词:负三角形变,高能量离子,鱼骨模PACS:52.25.Xz,52.35.Bj,52.30.Gz,52.35.PyDOI:10.7498/aps.72.202306501引言N=2.7H98y2=1.2托卡马克是最有前景的磁约束聚变装置之一,其中等离子体截面的
3、形变对于其约束与稳定性具有显著的影响.因此,对于未来的聚变装置,形变参数是关键的设计参数之一14.近年在 TCV 托卡马克装置上的研究发现相比于传统的正三角形变,负三角形变可降低湍流进而改善等离子体约束5,6.随后,DIII-D 托卡马克装置上也采用了负三角形变位形在 L 模下放电,发现达到了较高的归一化比压值()以及相当于 H 模的约束性能(),并且没有边界局域模发生7.另外,负三角形变位形可以降低偏滤器的热负载,这些优点使其被认为是未来聚变堆的一个更好的选择8.高能量粒子物理对于未来聚变堆中的燃烧等离子体来说是一个很重要的研究方向9.这是由于在未来的磁约束核聚变堆中,实现自持燃烧的等离子体
4、中的氘氚聚变反应会产生大量的快阿尔法粒子.此外,辅助加热下如中性束注入、离子回旋加热等也会产生大量的快粒子.这些快粒子可以激发出各种不稳定性,这些不稳定性反过来会引起很强的快粒子输运并导致快粒子损失到装置的第一壁,大量的快阿尔法粒子输运会严重降低其对背景等离子体的加热效率.另外,快粒子可能对磁流体不稳定性如新经典撕裂模、内扭曲模、电阻壁模等有显著的影响,而这些磁流体不稳定性有可能严重降低等离子体约束性能甚至会导致托卡马克破裂.高能量粒子激发的不稳定性包括阿尔芬本征模不稳定性、快粒子模、鱼骨模等.其中,鱼骨模不*国家磁约束核聚变能发展研究专项(批准号:2019YFE03050002)、国家自然科
5、学基金(批准号:12005003,11975270)和等离子体所科学基金(批准号:DSJJ-2022-04)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.21(2023)215202215202-11稳定性最早在 PDX 装置上有中性束注入的条件下被观测到10,是环向模数和极向模数都为 的一种芯部不稳定性.鱼骨模不稳定性在实验中被观测到之后,有两种理论被用来解释鱼骨模不稳定性形成的物理机制.1984年 Chen 等11认为鱼骨模是捕获高能量离子通过进动频率共振激发的.另外,1
6、986 年,Coppi 和 Porcelli12在考虑热离子逆磁效应后认为鱼骨模是一种类磁流体模,其频率与热离子逆磁漂移频率相当.随后,鱼骨模在多个装置上被观测到1320,因此鱼骨模是一种常见的高能量粒子激发的不稳定性.然而,目前针对负三角形变位形的磁约束等离子体开展的高能量粒子物理相关的工作较少,DIII-D 装置上开展了负三角形变位形下的高能量粒子激发阿尔芬本征模不稳定性的实验,发现相比于正三角形变位形阿尔芬本征模引起的高能量粒子输运并没有明显减弱21,但是相关的工作没有理论模拟方面的分析,主要是实验观测结果.因此,本工作针对负三角形变位形等离子体,使用动理学-磁流体混合模型程序 M3D-
7、K开展了高能量粒子激发的鱼骨模不稳定性的模拟研究,线性模拟结果发现负三角形变位形解稳理想内扭曲模不稳定性,但会致稳鱼骨模不稳定性.然而,非线性模拟发现在没有磁流体非线性效应时,负三角形变位形下的鱼骨模更不容易饱和.这些结果表明考虑高能量粒子激发的鱼骨模不稳定性后,负三角形变位形相比于正三角形变位形并没有明显优势.本文分为以下几个部分,第 2 节简单介绍了M3D-K 程序所用的物理模型以及模拟所用的主要参数.第 3 节展示了负三角形变下鱼骨模不稳定性的线性模拟结果.第 4 节展示了负三角形变下鱼骨模不稳定性的非线性模拟结果.最后,第 5 节给出了本工作的结论.2M3D-K 程序模型以及参数设置2
8、.1 M3D-K 程序模型本文使用了动理学-磁流体混合模型程序M3D-K22,23.该程序使用电阻磁流体模型描述热等离子体,使用漂移动理学或者回旋动理学模型描述高能量粒子.磁流体方程组在极向截面采用有限元方法求解,在大环方向采用 4 阶有限差分或者拟谱法.由于 M3D-K 代码采用柱坐标系,在磁轴处不存在奇异点,因而适用于等离子体芯部磁流体不稳定性的数值模拟.对于高能离子部分,为了降低粒子噪声以及模拟所需粒子数量,程序采用了f 方法.M3D-K 程序通过在动量方程中包含高能离子压强张量,将高能离子与背景磁流体自洽耦合.对于背景等离子体部分,程序还可以包含双流体效应,以及包括热离子逆磁漂移、电子
9、压强等物理效应.M3D-K 程序已经成功被用来模拟各种磁流体不稳定性与高能离子的相互作用2435,是模拟高能离子不稳定性的有力工具.M3D-K 程序所用的方程组如下,高能量粒子效应通过压强耦合进入动量方程:dvdt=J B P Ph+2v.(1)连续性方程:t+(v)=0.(2)能量方程:dpdt=p v+p.(3)法拉第定律:Bt=E.(4)安培定律:0J=B.(5)欧姆定律:E+v B=J.(6),P,v,E,B,J,0,d/dt=/t+v 这里 分别表示质量密度、压强、流体速度、电场、磁感应强度、电流密度、绝热系数、真空磁导率、热传导系数张量、电阻率、黏滞系数.高能量粒子压强张量在忽略非
10、对角项的情况下使用 Chew-Goldberger-Low(CGL)形式:Ph=PI+(P/P)bb,(7)F(X,v/,)式中,I 表示单位矩阵;b 是磁场方向单位矢量;平行和垂直磁场方向的压强利用高能量粒子分布函数 在回旋中心坐标 中计算,P/(x)=mv2/(x X h)F(X,v/,)Bd3Xdv/dd,(8)物理学报ActaPhys.Sin.Vol.72,No.21(2023)215202215202-2P(x)=12mv2(x X h)F(X,v/,)Bd3Xdv/dd,(9)h=vb/eB/m(x,v)其中 是回旋半径矢量,是回旋频率,e 是电子电荷量,是高能量粒子相空间坐标,是
11、磁矩,是回旋角.F高能量粒子回旋中心分布函数 可以表示为F=F(X,v/,)=i(X Xi)(v/v/,i)(i),(10)其中 是狄拉克函数.高能量粒子的运动用回旋动理学或者漂移动理学来描述:dXdt=1Bv/B b0E1e(B0+B),(11)mdv/dt=eBBE 1e(B0+B),(12)EBb0B0BBB其中,是总的电场,是总的磁场,是沿平衡磁场方向的单位矢量,和 分别是平衡和扰动的磁场.表示回旋平均,如果是漂移动理学的话直接用粒子导心处的电磁场.变量 和 分别表示为B=B0+B+mv/e b0,(13)B=B b0.(14)2.2 初始参数与平衡剖面设置a=0.45 mR0=1.8
12、6 m=1.63=0.436B0=1.615 Tne0=5.52 1019m3vA=B0/(00)1/2=3.35 106m/sA=R0/vA=5.55 107stotal,0=4.61%0001q=1=0.171主要的平衡剖面和参数根据 EAST 托卡马克装置的#71320 炮设置36,主要参数如下:小半径,大半径 ,椭圆形变 ,三 角 形 变 ,磁 轴 处 磁 场 ,磁轴处电子密度 ,阿尔芬速度 ,阿尔芬时间 ,磁轴处包括热等离子体和快离子的总压强 .这里 是真空磁导率,是磁轴处的热等离子体质量密度.平衡的安全因子剖面和压强剖面如图1 中蓝线所示,其中 是归一化极向磁通,其在磁轴处值为 ,
13、在等离子体最外闭合磁面处为 .面的位置如图1 中红色虚线所示,在 处.在本工作中为了简化分析,热等离子体密度的径向分布取为常数.E0=58 keV考虑中性束注入产生的高能量离子,其入射能量为 .高能量离子在速度方向以及在实空间的径向分布表达式为f=cH(v0 v)v3+v3cexp(0)2/2 exp(/),(15)cHv0=2E0/mD=2.36 106m/smDB0/E0=0.8=0.2=0.27vc(0.75me/mi)1/3(2Te/me)1/2=1.93 106m/s其中 是常数,是阶跃函数,是高能量离子的入射速度,是氘离子质量,是高能量离子轨道平均的极向磁通,是临界速度,me是电子
14、质量,mi是离子质量,Te是电子温度.8642000.20.40.60.81.000.51.00-0.500.20.40.60.81.0BaselineFlat 图1安全因子与总压强平衡剖面Fig.1.Equilibriumprofilesofsafetyfactorandtotalpres-sure.0.03A0.020.01000.10.20.30.5hot,0total,00.4=0.436=0=-0.4360.06A0.040.02000.10.20.30.5hot,0total,00.4=0.436=0=-0.436Phot,0/Ptotal,0图2不同三角形变下模频率和线性增长率与
15、快离子压强比值 的关系Phot,0/Ptotal,0Fig.2.Modefrequency and linear growth rate as a func-tionofthefastionpressurefraction .物理学报ActaPhys.Sin.Vol.72,No.21(2023)215202215202-33线性模拟结果=0.436=0=0.436Phot,0/Ptotal,0Phot,0Ptotal,0Phot,0/Ptotal,0W首先分析在负三角形变条件下的鱼骨模的线性不稳定性结果.图 2 对比了正三角形变 ,无三角形变 ,负三角形变 位形下快离子激发内扭曲模以及鱼骨模的
16、情况,图中 是快离子压强与总压强比值,是磁轴处快离子压强值,是磁轴处总压强值.可以看出,当没有快离子效应时,三种位形下的理想内扭曲模都是不稳定的,并且负三角形变对理想内扭曲模不稳定性起解稳作用,正三角形变对理想内 扭 曲 模 起 致 稳 作 用.当 快 离 子 压 强 比 值 增大时,快离子的动理学效应会致稳内扭曲模3739.然而,当快离子压强增大到超过一定阈值时,鱼骨模不稳定性被快离子共振激发.另外,从图2 可以看出,无三角形变或者负三角形变下鱼骨模很难激发(激发阈值值比较高),但一旦被激发则增长更快;正三角形变下激发鱼骨模的阈值比较低,但被激发后却增长得慢一些.对于三角形变影响内扭曲模的分
17、析,之前 Eriksson 和 Wah-lberg40的工作发现正的椭圆形变和正的三角形变对理想内扭曲模起稳定作用.Martynov 等41发现三角形变对理想内扭曲模的势能 有重要的贡献,考虑到三角形变后解析结果与 KINX 程序结果一致.另外,他们使用 KINX 磁流体本征程序模拟发现一定参数条件下负三角形变对理想内扭曲模起解稳作用.这些理论分析结果与本文模拟的结果一致.Phot,0/Ptotal,0=00Phot,0/Ptotal,0=0.15Phot,0/Ptotal,0=0.4图 3 给出了不同高能量离子压强比值的模结构.在负三角形变位形下,当快离子压强比值 时,理想内扭曲模频率为 ,
18、且模结构是上下对称的(图3(a);当快离子压强比值 时,模结构变得扭曲并且有一定大小的模频率(图3(b);当快离子压强比值 时,模式变为鱼骨模不稳定性,并且有更高的模频率(图3(c).作为对比,图3(d)给出了正三角形变位形下的鱼骨模,可以看出模结构与负三角形变位形下的鱼骨模很相似.q=1另外,由于未来聚变堆运行时需要先稳定磁流体模式,选取了没有高能量粒子影响下理想内扭曲模稳定的算例开展研究.考虑到环效应对内扭曲模稳定性的影响时,当 面之内的压强梯度足够小时,内扭曲模会趋于稳定42.因此,选取了一个芯部压强剖面平坦并且磁轴处压强值保持不变的算例,重新计算了不同快粒子驱动时 3 个三角形变0.6
19、(a)0.40.2/m/m0-0.2-0.4-0.6-0.81.4 1.6 1.8 2.0 2.2 2.40.6(b)0.40.2/m0-0.2-0.4-0.6-0.81.4 1.6 1.8 2.0 2.2 2.40.6(c)0.40.2/m/m0-0.2-0.4-0.6-0.81.4 1.6 1.8 2.0 2.2 2.40.6(d)0.40.2/m0-0.2-0.4-0.6-0.81.4 1.6 1.8 2.0 2.2 2.4Phot,0/Ptotal,0U=0.436Phot,0/Ptotal,0=0=0.436Phot,0/Ptotal,0=0.15=0.436Phot,0/Ptota
20、l,0=0.4=0.436Phot,0/Ptotal,0=0.4图3不同高能量离子压强比值 下的流函数(a),;(b),;(c),;(d),UPhot,0/Ptotal,0=0.436Phot,0/Ptotal,0=0=0.436Phot,0/Ptotal,0=0.15=0.436Phot,0/Ptotal,0=0.4=0.436Phot,0/Ptotal,0=0.4Fig.3.Velocitystreamfunction atdifferentfastionpre-ssure fraction :(a),;(b),;(c),;(d),.0.015A0.0100.00500.40.50.6ho
21、t,0total,00.70.40.50.60.7=0.436=-0.436=0.436=-0.4360.054A0.0520.0500.0480.046hot,0total,0Phot,0/Ptotal,0图4芯部压强剖面平坦下模频率和线性增长率与快离子压强比值 的关系Phot,0/Ptotal,0Fig.4.Modefrequency and linear growth rate as a func-tionofthefastionpressurefraction withflatpressureprofile.物理学报ActaPhys.Sin.Vol.72,No.21(2023)2152
22、02215202-4参数下的内扭曲模的线性稳定性,压强剖面如图1的红线所示.新压强剖面下的鱼骨模增长率以及模频率如图4 所示.可以看出,在新压强剖面下鱼骨模的激发阈值更高,并且正三角形变下鱼骨模的激发阈值仍然低于负三角形变.4非线性模拟结果n=1Phot,0/Ptotal,0=0.4total,0=3.91%Phot,0/Ptotal,0=0.4A=0.00906本部分考虑没有磁流体非线性效应下的鱼骨模非线性演化,这里将热等离子体的磁流体响应通过只保留 环向扰动而限制为线性的.由于线性部分中 条件下负三角形变位形的鱼骨模比正三角形变位形的鱼骨模增长率高.为了更好地对比,选取了一个更低总比压值
23、的负三角形变位形下的鱼骨模算例,并且快离子压强比值仍保持为 ,如图5 所示.可以看出,图 5(a)中正三角形变鱼骨模与负三角形变鱼骨模的线性增长率基本一样,增长率均为 .但是正三角形变的鱼骨模可以非线性饱和,而负三角形变的鱼骨模不能饱和.图 5(b)给出了两种三角形变位形下鱼骨模向下扫频的特征,正三角形变位形下鱼骨模频率从=0.0413A=0.00707A=0.0336A0Phot,0/Ptotal,0=0.4,total,0=4.61%,=0.436n=1n=0 向下扫频到 ,负三角形变位形下鱼骨模频率从 向下扫频到 .图6 为 算例的动能和磁能的 分量随时间演化.可以看出,它们都随时间线性
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 三角 形变 马克 位形下 高能量 离子 激发 鱼骨 模拟 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。