搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

中国首台准环对称仿星器中离子温度梯度模的模拟研究

黄捷 李沫杉 覃程 王先驱

引用本文:
Citation:

中国首台准环对称仿星器中离子温度梯度模的模拟研究

黄捷, 李沫杉, 覃程, 王先驱

Simulation of ion temperature gradient mode in Chinese First Quasi-axisymmetric Stellarator

Huang Jie, Li Mo-Shan, Qin Cheng, Wang Xian-Qu
PDF
HTML
导出引用
  • 中国首台准环对称仿星器(CFQS)是目前世界上唯一在建的准环对称仿星器. 本文利用回旋弗拉索夫代码GKV开展了CFQS中离子温度梯度模(ITG)的模拟研究. 在静电绝热条件下, 模拟的结果给出了CFQS中纯的ITG与密度梯度和温度梯度间的依赖关系. ITG的激发存在温度梯度阈值, 此温度梯度阈值受到密度梯度的影响. ITG的增长率不仅与密度梯度的绝对值相关, 还取决于密度梯度的正负, 负密度梯度对ITG具有强的抑制作用. 非绝热的模拟结果表明, 捕获电子对ITG具有去稳作用, 电子温度梯度也对ITG具有去稳作用. 当考虑电磁条件时, 有限的等离子体比压会抑制ITG, 导致ITG向阿尔芬离子温度梯度模/动理学气球模(AITG/ KBM)的转化. 当密度和温度梯度都较大时, KBM的最大增长率与密度梯度和温度梯度近似成线性关系.
    The Chinese First Quasi-axisymmetric Stellarator (CFQS) is now the only quasi-axisymmetric stellarator under construction in the world. In this work, ion temperature gradient (ITG) mode in CFQS is studied by using gyrokinetic Vlasov code GKV. The basic characteristics of the eletrtostatic ITG are separately given under the adiabatic condition and the non-adiabatic condition. There is a critical temperature gradient for ITG. The growth rate of ITG is proportional to the temperature gradient. Furthermore, the growth rate depends on not only the absolute value of density gradient, but also the plus or minus sign of the density gradient. The negative density gradient can strongly suppress the ITG. The kinetic electron can destabilize the ITG and the electron temperature gradient can also destabilize the ITG. For electromagnetic condition, the ITG modes can be suppressed by the finite plasma beta, and then a transition from ITG to Alfvenic ion temperature gradient mode/kinetic ballooning mode (AITG/KBM) comes into being. The maximum growth rate of KBM is linearly proportional to density gradient and temperature gradient when both gradients are large.
      通信作者: 黄捷, jiehuang@swjtu.edu.cn
    • 基金项目: 国家磁约束核聚变能发展研究专项(批准号: 2022YFE03070000, 2022YFE03070001)、国家自然科学基金(批准号: 11820101004)和四川省国际科技创新合作项目(批准号: 2021YFH0066)资助的课题.
      Corresponding author: Huang Jie, jiehuang@swjtu.edu.cn
    • Funds: Project supported by the Chinese National Fusion Project for ITER (Grant Nos. 2022YFE03070000, 2022YFE03070001), the National Natural Science Foundation of China (Grant No. 11820101004), and the Sichuan Provincial International Science and Technology Innovation Cooperation Project, China (Grant No. 2021YFH0066).
    [1]

    Xu Y 2016 Matter Radiat. Extremes 1 192Google Scholar

    [2]

    Ho D D M 1987 Phys. Fluids 30 442Google Scholar

    [3]

    Boozer A H 1995 Plasma Phys. Controlled Fusion 37 A103Google Scholar

    [4]

    Subbotin A A, Mikhailov M I, Shafranov V D, Isaev M Yu, Nührenberg C, Nührenberg J, Zille R, Nemov V V, Kasilov SV, Kalyuzhnyj V N 2006 Nucl. Fusion 46 921Google Scholar

    [5]

    Garabedian P 1996 Phys. Plasmas 3 2483Google Scholar

    [6]

    Shimizu A, Liu H F, Isobe M, Okamura S, Nishimura S, Suzuki C, Xu Y, Zhang X, Liu B, Huang J, Wang X Q, Liu H, Tang C J, CFQS team 2018 Plasma Fusion Res. 13 3403123Google Scholar

    [7]

    Zarnstorff M C, Berry L A, Brooks A, Fredrickson1 E, Fu G Y, Hirshman S, Hudson, Ku L P, Lazarus E, Mikkelsen D, Monticello D, Neilson G H, Pomphrey N, Reiman A, Spong D, Strickler D, Boozer A, Cooper W A, Goldston R, Hatcher R, Isaev M, Kessel C, Ewandowski J L, Lyon J F, Merkel P, Mynick H, Nelson B E, Nuehrenberg C, Redi M, Reiersen W, Rutherford P, Sanchez R, Schmidt J, White R B 2001 Plasma Phys. Controlled Fusion 43 A237Google Scholar

    [8]

    Okamura S, Matsuoka K, Nishimura S, Isobe M, Nomura I, Suzuki C, Shimizu A, Murakami S, Nakajima N, Yokoyama M 2001 Nucl. Fusion 41 1865Google Scholar

    [9]

    Liu H F, Shimizu A, Isobe M, Okamura S, Nishimura S, Suzuki C, Xu Y, Zhang X, Liu B, Huang J, Wang X Q, Liu H, Tang C J, Yin D P, Wan Y, CFQS team 2018 Plasma Fusion Res. 13 3405067Google Scholar

    [10]

    Isobe M, Shimizu A, Liu H F, Liu H, Xiong G Z, Yin D P, Ogawa K, Yoshimura Y, Nakata M, Kinoshita S, Okamura S, tang C J, Xu Y, CFQS Team 2019 Plasma Fusion Res. 14 3402074Google Scholar

    [11]

    Liu H F, Shimizu A, Xu Y, Okamura S, Kinoshita S, Isobe M, Li Y B, Xiong G Z, Wang X Q, Huang J, Cheng J, Liu H, Zhang X, Yin D P, Wang Y, Murase T, Nakagawa S, Tang C J 2021 Nucl. Fusion 61 016014Google Scholar

    [12]

    Wang X Q, Xu Y, Shimizu A, Isobe M, Okamura S, Todo Y, Wang H, Liu H F, Huang J, Zhang X, Liu H, Cheng J, Tang C J, CFQS team 2021 Nucl. Fusion 61 036021Google Scholar

    [13]

    Horton W 1999 Rev. Mod. Phys. 71 735Google Scholar

    [14]

    Watanabe T H, Sugama H 2006 Nucl. Fusion 46 24Google Scholar

    [15]

    Nakata M, Nunami M, Sugama H 2017 Phys. Rev. Lett. 118 165002Google Scholar

    [16]

    Antonsen T M, Lane B 1980 Phys. Fluids 23 1205Google Scholar

    [17]

    Nakata M, Honda M, Yoshida M, Urano H, Nunami M, Maeyama S, Watanabe T H, Sugama H 2016 Nucl. Fusion 56 086010Google Scholar

    [18]

    Beer M A, Cowley S C, Hammett G W 1995 Phys. Plasmas 2 2687Google Scholar

    [19]

    Romanelli M, Bourdelle C, Dorland W 2004 Phys. Plasmas 11 3845Google Scholar

    [20]

    Du H R, Jhang H, Hahm T S, Dong J Q, Wang Z X 2017 Phys. Plasmas 24 122501Google Scholar

    [21]

    沈勇, 董家齐, 徐红兵 2018 物理学报 67 195203Google Scholar

    Shen Y, Dong J Q, Xu H B 2018 Acta Phys. Sin. 67 195203Google Scholar

    [22]

    Baumgaertel J A, Hammett G W, Mikkelsen D R, Nunami M, Xanthopoulos P 2012 Phys. Plasmas 19 122306Google Scholar

    [23]

    Dominguez R R, Waltz R E 1988 Phys. Fluids 31 3147Google Scholar

    [24]

    Nunami M, Watanabe T H, Sugama H, Tanaka K 2011 Plasma Fusion Res. 6 1403001Google Scholar

    [25]

    Alcusón J A, Xanthopoulos P, Plunk G G, Helander P, Wilms F, Turkin Y, Stechow A von, Grulke O 2020 Plasma Phys. Controlled Fusion 62 035005Google Scholar

    [26]

    罗一鸣, 王占辉, 陈佳乐, 吴雪科, 付彩龙, 何小雪, 刘亮, 杨曾辰, 李永高, 高金明, 杜华荣, 昆仑集成模拟设计组 2022 物理学报 71 075201Google Scholar

    Luo Y M, Wang Z H, Chen J L, Wu X K, Fu C L, He X X, Liu L, Yang Z C, Li Y G, Gao J M, Du H R, Kulun Integrated Simulation and Design Group 2022 Acta Phys. Sin. 71 075201Google Scholar

    [27]

    Mahmood M A, Rafiq T, Persson M, Weiland J 2009 Phys. Plasmas 16 022503Google Scholar

    [28]

    Dong J Q, Mahajan S M, Horton W 1997 Phys. Plasmas 4 755Google Scholar

    [29]

    Peeters A G, Angioni C, Apostoliceanu M, Jenko F, Ryter F, the ASDEX Upgrade team 2005 Phys. Plasmas 12 022505Google Scholar

    [30]

    Sandberg I, Isliker H, Pavlenko V P 2007 Phys. Plasmas 14 092504Google Scholar

    [31]

    Qi L, Kwon J, Hahm T S, Jo G 2016 Phys. Plasmas 23 062513Google Scholar

    [32]

    Malinov P, Zonca F 2005 J. Plasma Phys. 71 301Google Scholar

    [33]

    Kim J Y, Han H S 2017 Phys. Plasmas 24 072501Google Scholar

    [34]

    Pueschel M J, Jenko F 2010 Phys. Plasmas 17 062307Google Scholar

    [35]

    Xie H S, Lu Z X, Li B 2018 Phys. Plasmas 25 072106Google Scholar

    [36]

    Aleynikova K, Zocco A 2017 Phys. Plasmas 24 092106Google Scholar

    [37]

    Turnbull A D, Strait E J, Heidbrink W W, Chu M S, Duong H H, Greene J M, Lao L L, Taylor T S, Thompson S J 1993 Phys. Fluids B 5 2546Google Scholar

    [38]

    Dong J, Chen L, Zonca F 1999 Nucl. Fusion 39 1041Google Scholar

    [39]

    谢华生 2015 博士学位论文 (杭州: 浙江大学)

    Xie H S 2015 Ph. D. Dissertation (Hangzhou: Zhejiang University) (in Chinese)

    [40]

    Aleynikova K, Zocco1 A, Xanthopoulos1 P, Helander1 P, Nührenberg C 2018 J. Plasma Phys. 84 745840602Google Scholar

  • 图 1  不同密度梯度和温度梯度下绝热ITG的增长率波数谱 (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $

    Fig. 1.  Growth rate spectra of ITG for different density gradients and temperature gradients: (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $.

    图 2  不同密度梯度和温度梯度下绝热ITG的频率波数谱 (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $

    Fig. 2.  Real frequency spectra of ITG for different density gradients and temperature gradients: (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $.

    图 3  绝热ITG的最大增长率与密度梯度和温度梯度的关系

    Fig. 3.  Relationship of the maximum growth rate of adiabatic ITG to $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{T} $.

    图 4  考虑捕获电子效应后ITG的增长率(a)和频率(b)波数谱, 其中$ {R}_{0}/{L}_{n}=2 $, 电子温度梯度标长和离子温度梯度标长相等, 即$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $

    Fig. 4.  Growth rate (a) and real frequency spectra (b) of kinetic ITG for $ {R}_{0}/{L}_{n}=2 $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $.

    图 5  考虑捕获电子效应后ITG/TE-ITG最大增长率与密度梯度和温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $

    Fig. 5.  Contour map of the maximum growth rate of kinetic ITG/TE-ITG mode vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{T} $. Here, ${R}_{0}/ {L}_{{T}_{\mathrm{e}}} $$ ={R}_{0}/{L}_{{T}_{\mathrm{i}}}$.

    图 6  考虑捕获电子效应后ITG的增长率(a)和频率波数谱(b), 其中$ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

    Fig. 6.  Growth rate (a) and real frequency spectra (b) of kinetic ITG for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $.

    图 7  考虑捕获电子效应后ITG最大增长率与密度梯度和离子温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

    Fig. 7.  Contour map of the maximum growth rate of kinetic ITG vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}} $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $.

    图 8  当波数$ {k}_{y}{\rho }_{\mathrm{i}}=1.0 $时增长率(a)与频率(b)随比压的变化(其中, $ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}}=8 $), 图(a)中的箭头是从ITG转变为KBM的转变点

    Fig. 8.  Growth rates (a) and real frequencies (b) vs. $ \beta $ for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}}={R}_{0}/{{L}}_{{{T}}_{\mathrm{e}}}=8 $ at $ {k}_{y}{\rho }_{\mathrm{i}}=1.0 $. The arrow is plotted in panel (a) to point the transition point from ITG to KBM.

    图 9  $\beta =1{\text{%}}$时, KBM的增长率(a)和频率(b)波数谱, 其中$ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

    Fig. 9.  Growth rate (a) and real frequency spectra (b) of KBM for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $. Here, $\beta =1{\text{%}}$.

    图 10  $\beta =1{\text{%}}$时, KBM的最大增长率与密度梯度和离子温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

    Fig. 10.  Contour map of the maximum growth rate of KBM vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}} $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $ and $\beta =1{\text{%}}$.

  • [1]

    Xu Y 2016 Matter Radiat. Extremes 1 192Google Scholar

    [2]

    Ho D D M 1987 Phys. Fluids 30 442Google Scholar

    [3]

    Boozer A H 1995 Plasma Phys. Controlled Fusion 37 A103Google Scholar

    [4]

    Subbotin A A, Mikhailov M I, Shafranov V D, Isaev M Yu, Nührenberg C, Nührenberg J, Zille R, Nemov V V, Kasilov SV, Kalyuzhnyj V N 2006 Nucl. Fusion 46 921Google Scholar

    [5]

    Garabedian P 1996 Phys. Plasmas 3 2483Google Scholar

    [6]

    Shimizu A, Liu H F, Isobe M, Okamura S, Nishimura S, Suzuki C, Xu Y, Zhang X, Liu B, Huang J, Wang X Q, Liu H, Tang C J, CFQS team 2018 Plasma Fusion Res. 13 3403123Google Scholar

    [7]

    Zarnstorff M C, Berry L A, Brooks A, Fredrickson1 E, Fu G Y, Hirshman S, Hudson, Ku L P, Lazarus E, Mikkelsen D, Monticello D, Neilson G H, Pomphrey N, Reiman A, Spong D, Strickler D, Boozer A, Cooper W A, Goldston R, Hatcher R, Isaev M, Kessel C, Ewandowski J L, Lyon J F, Merkel P, Mynick H, Nelson B E, Nuehrenberg C, Redi M, Reiersen W, Rutherford P, Sanchez R, Schmidt J, White R B 2001 Plasma Phys. Controlled Fusion 43 A237Google Scholar

    [8]

    Okamura S, Matsuoka K, Nishimura S, Isobe M, Nomura I, Suzuki C, Shimizu A, Murakami S, Nakajima N, Yokoyama M 2001 Nucl. Fusion 41 1865Google Scholar

    [9]

    Liu H F, Shimizu A, Isobe M, Okamura S, Nishimura S, Suzuki C, Xu Y, Zhang X, Liu B, Huang J, Wang X Q, Liu H, Tang C J, Yin D P, Wan Y, CFQS team 2018 Plasma Fusion Res. 13 3405067Google Scholar

    [10]

    Isobe M, Shimizu A, Liu H F, Liu H, Xiong G Z, Yin D P, Ogawa K, Yoshimura Y, Nakata M, Kinoshita S, Okamura S, tang C J, Xu Y, CFQS Team 2019 Plasma Fusion Res. 14 3402074Google Scholar

    [11]

    Liu H F, Shimizu A, Xu Y, Okamura S, Kinoshita S, Isobe M, Li Y B, Xiong G Z, Wang X Q, Huang J, Cheng J, Liu H, Zhang X, Yin D P, Wang Y, Murase T, Nakagawa S, Tang C J 2021 Nucl. Fusion 61 016014Google Scholar

    [12]

    Wang X Q, Xu Y, Shimizu A, Isobe M, Okamura S, Todo Y, Wang H, Liu H F, Huang J, Zhang X, Liu H, Cheng J, Tang C J, CFQS team 2021 Nucl. Fusion 61 036021Google Scholar

    [13]

    Horton W 1999 Rev. Mod. Phys. 71 735Google Scholar

    [14]

    Watanabe T H, Sugama H 2006 Nucl. Fusion 46 24Google Scholar

    [15]

    Nakata M, Nunami M, Sugama H 2017 Phys. Rev. Lett. 118 165002Google Scholar

    [16]

    Antonsen T M, Lane B 1980 Phys. Fluids 23 1205Google Scholar

    [17]

    Nakata M, Honda M, Yoshida M, Urano H, Nunami M, Maeyama S, Watanabe T H, Sugama H 2016 Nucl. Fusion 56 086010Google Scholar

    [18]

    Beer M A, Cowley S C, Hammett G W 1995 Phys. Plasmas 2 2687Google Scholar

    [19]

    Romanelli M, Bourdelle C, Dorland W 2004 Phys. Plasmas 11 3845Google Scholar

    [20]

    Du H R, Jhang H, Hahm T S, Dong J Q, Wang Z X 2017 Phys. Plasmas 24 122501Google Scholar

    [21]

    沈勇, 董家齐, 徐红兵 2018 物理学报 67 195203Google Scholar

    Shen Y, Dong J Q, Xu H B 2018 Acta Phys. Sin. 67 195203Google Scholar

    [22]

    Baumgaertel J A, Hammett G W, Mikkelsen D R, Nunami M, Xanthopoulos P 2012 Phys. Plasmas 19 122306Google Scholar

    [23]

    Dominguez R R, Waltz R E 1988 Phys. Fluids 31 3147Google Scholar

    [24]

    Nunami M, Watanabe T H, Sugama H, Tanaka K 2011 Plasma Fusion Res. 6 1403001Google Scholar

    [25]

    Alcusón J A, Xanthopoulos P, Plunk G G, Helander P, Wilms F, Turkin Y, Stechow A von, Grulke O 2020 Plasma Phys. Controlled Fusion 62 035005Google Scholar

    [26]

    罗一鸣, 王占辉, 陈佳乐, 吴雪科, 付彩龙, 何小雪, 刘亮, 杨曾辰, 李永高, 高金明, 杜华荣, 昆仑集成模拟设计组 2022 物理学报 71 075201Google Scholar

    Luo Y M, Wang Z H, Chen J L, Wu X K, Fu C L, He X X, Liu L, Yang Z C, Li Y G, Gao J M, Du H R, Kulun Integrated Simulation and Design Group 2022 Acta Phys. Sin. 71 075201Google Scholar

    [27]

    Mahmood M A, Rafiq T, Persson M, Weiland J 2009 Phys. Plasmas 16 022503Google Scholar

    [28]

    Dong J Q, Mahajan S M, Horton W 1997 Phys. Plasmas 4 755Google Scholar

    [29]

    Peeters A G, Angioni C, Apostoliceanu M, Jenko F, Ryter F, the ASDEX Upgrade team 2005 Phys. Plasmas 12 022505Google Scholar

    [30]

    Sandberg I, Isliker H, Pavlenko V P 2007 Phys. Plasmas 14 092504Google Scholar

    [31]

    Qi L, Kwon J, Hahm T S, Jo G 2016 Phys. Plasmas 23 062513Google Scholar

    [32]

    Malinov P, Zonca F 2005 J. Plasma Phys. 71 301Google Scholar

    [33]

    Kim J Y, Han H S 2017 Phys. Plasmas 24 072501Google Scholar

    [34]

    Pueschel M J, Jenko F 2010 Phys. Plasmas 17 062307Google Scholar

    [35]

    Xie H S, Lu Z X, Li B 2018 Phys. Plasmas 25 072106Google Scholar

    [36]

    Aleynikova K, Zocco A 2017 Phys. Plasmas 24 092106Google Scholar

    [37]

    Turnbull A D, Strait E J, Heidbrink W W, Chu M S, Duong H H, Greene J M, Lao L L, Taylor T S, Thompson S J 1993 Phys. Fluids B 5 2546Google Scholar

    [38]

    Dong J, Chen L, Zonca F 1999 Nucl. Fusion 39 1041Google Scholar

    [39]

    谢华生 2015 博士学位论文 (杭州: 浙江大学)

    Xie H S 2015 Ph. D. Dissertation (Hangzhou: Zhejiang University) (in Chinese)

    [40]

    Aleynikova K, Zocco1 A, Xanthopoulos1 P, Helander1 P, Nührenberg C 2018 J. Plasma Phys. 84 745840602Google Scholar

  • [1] 周利娜, 胡汉卿, 刘钺强, 段萍, 陈龙, 张瀚予. 等离子体对共振磁扰动的流体和动理学响应的模拟研究. 物理学报, 2023, 72(7): 075202. doi: 10.7498/aps.72.20222196
    [2] 肖士妍, 贾大功, 聂安然, 余辉, 吉喆, 张红霞, 刘铁根. 开放式多通道多芯少模光纤表面等离子体共振生物传感器. 物理学报, 2020, 69(13): 137802. doi: 10.7498/aps.69.20200353
    [3] 沈勇, 董家齐, 徐红兵. 托卡马克离子温度梯度湍流输运同位素定标修正中杂质的影响. 物理学报, 2018, 67(19): 195203. doi: 10.7498/aps.67.20180703
    [4] 谢会乔, 谭熠, 刘阳青, 王文浩, 高喆. 中国联合球形托卡马克氦放电等离子体的碰撞辐射模型及其在谱线比法诊断的应用. 物理学报, 2014, 63(12): 125203. doi: 10.7498/aps.63.125203
    [5] 邹长林, 叶文华, 卢新培. 一维动理学数值模拟激光与等离子体的相互作用. 物理学报, 2014, 63(8): 085207. doi: 10.7498/aps.63.085207
    [6] 刘小龙, 黄建国, 雷开卓. 水下等离子体声源的冲击波负压特性. 物理学报, 2013, 62(20): 204301. doi: 10.7498/aps.62.204301
    [7] 陆赫林, 陈忠勇, 李跃勋, 杨恺. 磁场剪切对离子温度梯度模带状流产生的影响. 物理学报, 2011, 60(8): 085202. doi: 10.7498/aps.60.085202
    [8] 陆赫林, 王顺金. 离子温度梯度模湍流的带状流最小自由度模型. 物理学报, 2009, 58(1): 354-362. doi: 10.7498/aps.58.354
    [9] 简广德, 董家齐. 环形等离子体中电子温度梯度不稳定性的粒子模拟. 物理学报, 2003, 52(7): 1656-1662. doi: 10.7498/aps.52.1656
    [10] 黄朝松, 李钧. 等离子体交换不稳定性的模耦合理论. 物理学报, 1992, 41(5): 783-791. doi: 10.7498/aps.41.783
    [11] 张澄, 邓晓华, 霍裕平. 输运与撕裂模相互作用下等离子体的时空结构. 物理学报, 1990, 39(10): 1573-1582. doi: 10.7498/aps.39.1573
    [12] 汪茂泉. 流动等离子体对托卡马克中撕裂模的影响. 物理学报, 1986, 35(9): 1227-1232. doi: 10.7498/aps.35.1227
    [13] 潘传红, 丁厚昌, 吴灵桥. 耗散气球模动力理论. 物理学报, 1986, 35(11): 1411-1425. doi: 10.7498/aps.35.1411
    [14] 石秉仁, 隋国芳. 含第二稳定区的环流器微观气球模分析. 物理学报, 1984, 33(11): 1546-1555. doi: 10.7498/aps.33.1546
    [15] 董家齐. 有磁辫等离子体的双撕裂模研究. 物理学报, 1984, 33(10): 1341-1349. doi: 10.7498/aps.33.1341
    [16] 顾永年. 小环径比锐边界等离子体的扭曲模不稳定性. 物理学报, 1984, 33(4): 554-560. doi: 10.7498/aps.33.554
    [17] 石秉仁. 环流器等离子体高n气球模的第二稳定区. 物理学报, 1983, 32(11): 1398-1406. doi: 10.7498/aps.32.1398
    [18] 郭世宠, 沈解伍, 陈骝, 蔡诗东. 离子温度梯度不稳定性的解析理论. 物理学报, 1982, 31(1): 17-29. doi: 10.7498/aps.31.17
    [19] 石秉仁. 高比压环流器等离子体位形的稳定性. 物理学报, 1982, 31(10): 1308-1316. doi: 10.7498/aps.31.1308
    [20] 周玉美, 蔡诗东. 磁化等离子体里高频模和低频模耦合的参量不稳定性. 物理学报, 1980, 29(7): 916-926. doi: 10.7498/aps.29.916
计量
  • 文章访问数:  903
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-17
  • 修回日期:  2022-06-28
  • 上网日期:  2022-09-08
  • 刊出日期:  2022-09-20

中国首台准环对称仿星器中离子温度梯度模的模拟研究

  • 西南交通大学物理科学与技术学院, 成都 610031
  • 通信作者: 黄捷, jiehuang@swjtu.edu.cn
    基金项目: 国家磁约束核聚变能发展研究专项(批准号: 2022YFE03070000, 2022YFE03070001)、国家自然科学基金(批准号: 11820101004)和四川省国际科技创新合作项目(批准号: 2021YFH0066)资助的课题.

摘要: 中国首台准环对称仿星器(CFQS)是目前世界上唯一在建的准环对称仿星器. 本文利用回旋弗拉索夫代码GKV开展了CFQS中离子温度梯度模(ITG)的模拟研究. 在静电绝热条件下, 模拟的结果给出了CFQS中纯的ITG与密度梯度和温度梯度间的依赖关系. ITG的激发存在温度梯度阈值, 此温度梯度阈值受到密度梯度的影响. ITG的增长率不仅与密度梯度的绝对值相关, 还取决于密度梯度的正负, 负密度梯度对ITG具有强的抑制作用. 非绝热的模拟结果表明, 捕获电子对ITG具有去稳作用, 电子温度梯度也对ITG具有去稳作用. 当考虑电磁条件时, 有限的等离子体比压会抑制ITG, 导致ITG向阿尔芬离子温度梯度模/动理学气球模(AITG/ KBM)的转化. 当密度和温度梯度都较大时, KBM的最大增长率与密度梯度和温度梯度近似成线性关系.

English Abstract

    • 在磁约束聚变装置中, 湍流输运导致的等离子体粒子和能量损失严重影响等离子体的约束性能. 在托卡马克装置中, 湍流输运是最主要的等离子体输运损失方式. 但在仿星器中, 除湍流输运外, 新经典输运也是一种主要的等离子体输运损失方式[1]. 由于存在大的新经典输运, 长期以来仿星器的约束性能都不如托卡马克. 仿星器中大的新经典输运是由于仿星器中大的磁场波纹度导致的[2]. 因此, 可以通过优化仿星器的磁场位形, 降低其磁场波纹度, 进而提高对等离子的约束性能. 仿星器磁场位形的优化主要有三种方式: 准螺旋对称优化(如美国的HSX装置)、准力线优化(如德国的W7-X装置)和准环对称优化[3-5]. 利用前两种优化方式优化的仿星器都已建成并运行, 但在这两类仿星器上的研究结果表明, 以这两种优化方式优化的仿星器仍然具有大的新经典输运 [1]. 理论研究表明, 采用准环对称优化的仿星器具有与托卡马克相当的新经典输运水平, 这能大大地降低仿星器的新经典输运损失[6]. 准环对称优化的仿星器克服了传统仿星器新经典输运大的缺点, 同时保留了仿星器无等离子体电流的优点, 是一种比较理想的约束等离子体的磁场位形.

      到目前为止, 国际上尚无建成的基于准环对称优化的仿星器装置. 美国曾经计划建造一台准环对称仿星器装置NCSX[7], 日本也曾设计过一台准环对称仿星器装置CHS-qa[8]. 最近, 中国西南交通大学联合日本核融合科学研究所共同设计了一台准环对称仿星器—中国首台准环对称仿星器(CFQS), 并已开始装置的建造[6,9-12]. CFQS是目前世界上唯一在建的准环对称仿星器装置.

      虽然理论研究表明CFQS的新经典输运水平远低于现有的仿星器, 可以与托卡马克相比拟[6], 但CFQS对等离子体的整体约束性能不仅取决于其新经典输运水平, 还取决于其湍流输运水平. 因此, 要从理论上验证CFQS磁场位形的先进性, 还需要对其开展湍流输运的模拟研究. 等离子体中的湍流输运是由等离子体中的静电和电磁微观不稳定性驱动的, 如离子温度梯度模(ITG)、捕获电子模(TEM)等[13].

      本文利用回旋动理学弗拉索夫代码GKV开展CFQS中离子尺度的微观不稳定性的研究, 给出CFQS上ITG和阿尔芬离子温度梯度模/动理学气球模(AITG/KBM)的特征.

    • 本文利用回旋弗拉索夫代码GKV开展CFQS上静电和电磁ITG的模拟研究. GKV代码是基于$ \text{δ}f $的回旋动理学模型, 通过求解随时间演化的电子或离子扰动分布函数$ \text{δ}{f}_{s} $、静电势$\phi$和失势的平行分量$ {A}_{//} $来开展等离子体物理的研究, 其中, $ s=\mathrm{i} $$ s=\mathrm{e} $分别是离子和电子的扰动分布函数[14,15]. GKV代码的核心方程是在五维相空间里扰动分布函数的电磁回旋动理学方程[16,17], 方程的傅里叶展开形式如下:

      $ \begin{split} &\left(\frac{\partial }{\partial t}+{v}_{// }\boldsymbol{b}\cdot \nabla +\mathrm{i}{\omega }_{\mathrm{D}s}-\frac{\mu {b}\cdot \nabla B}{{m}_{s}}\frac{\partial }{\partial {v}_{// }}\right)\text{δ}{g}_{s{{\boldsymbol{k}}}_{\perp }}\\ &-\frac{c}{B}\sum _{\Delta }\boldsymbol{b}\cdot \left(\boldsymbol{k}'_{\perp }\times \boldsymbol{k}''_{\perp }\right)\text{δ}{\psi }_{s{\boldsymbol{k}}'_{\perp }}\text{δ}{g}_{s{\boldsymbol{k}}''_{\perp }} \\=\;&\frac{{e}_{s}{F}_{\mathrm{M}s}}{{T}_{s}}\left(\frac{\partial \text{δ}{\psi }_{{{\boldsymbol{k}}}_{\perp }}}{\partial t}+\mathrm{i}{\omega }_{*\mathrm{T}s}\text{δ}{\psi }_{{{\boldsymbol{k}}}_{\perp }}\right)+{C}_{s}\left(\text{δ}{g}_{s{{\boldsymbol{k}}}_{\perp }}\right)\text{, } \end{split}$

      其中$\text{δ}{g}_{s{{\boldsymbol{k}}}_{\perp }}$为扰动分布函数$ \text{δ}{f}_{s{{\boldsymbol{k}}}_{\perp }} $的非绝热(动理学)部分, 即$\text{δ}{g}_{s{{\boldsymbol{k}}}_{\perp }}=\text{δ}{f}_{s{{\boldsymbol{k}}}_{\perp }}+{e}_{s}{\mathrm{J}}_{0 s}\text{δ}{\phi }_{{{\boldsymbol{k}}}_{\perp }}{F}_{\mathrm{M}s}/{T}_{s}$, $ {\mathrm{J}}_{0 s} $为零阶贝塞尔函数; $ {F}_{\mathrm{M}s} $为分布函数的平衡部分; ${\boldsymbol{k}}_{\perp }$为垂直波数; $ {\omega }_{\mathrm{D}s} $$ {\omega }_{*\mathrm{T}s} $分别是磁漂移和抗磁漂移频率; $ \text{δ}{\psi }_{{{\boldsymbol{k}}}_{\perp }} $是回旋平均的势扰动, (1)式中其他符号的含义见参考文献[17].

      在GKV中采用通量管坐标系$ \left(x, y, z\right) $[18], 其中, $x = a\left(\rho -{\rho }_{0}\right)$, $y = a{\rho }_{0}{q\left({\rho }_{0}\right)}^{-1}\left[q\left(\rho \right)\theta - \zeta \right]$, $ z=\theta $. 这里$ a $为小半径, $ \rho =\sqrt{\psi /{\psi }_{\mathrm{e}\mathrm{d}\mathrm{g}\mathrm{e}}} $, $ \psi $为环向磁通, $ q\left({\rho }_{0}\right) $$ \rho ={\rho }_{0} $处的安全因子, $ \theta $为极向角. 在后面的模拟中, 在环向角$\zeta =0^{\circ}$、归一化小半径$ \rho =0.5 $处取一个通量管作为模拟区域. CFQS的大半径为1 m, 纵横比为4, 其准环对称磁场位形的环向周期数为2.

    • 离子温度梯度模是由离子温度梯度驱动的一种微观不稳定性, 除了离子温度梯度外, 等离子体的密度梯度、捕获电子、等离子体比压等都对ITG有影响. 本文主要考虑在绝热和非绝热情况下, 温度梯度和密度梯度对ITG的影响. 在绝热情况下, 等离子体的微观不稳定性是纯的ITG模, 这样可以在CFQS仿星器上研究纯的ITG模与温度梯度和密度梯度间的依赖关系.

    • 首先, 考虑静电绝热条件下CFQS中的ITG模. 在这种条件下, 可以研究纯的ITG与等离子体参数间的依赖关系, 为了解CFQS中ITG模的激发和稳定条件提供理论依据. 图1给出了静电绝热条件下ITG模的波数谱. 图1中的$ \gamma $为增长率; $ {k}_{y} $为波数; $ {\rho }_{\mathrm{i}} $为离子回旋半径; $ {R}_{0} $为装置大半径; $ {L}_{n} $$ {L}_{T} $分别为密度和温度梯度标长, 其定义为${L}_{n}= -{\left(\mathrm{d}\ln n/\mathrm{d}\, x\right)}^{-1},~ {L}_{T}=-{\left(\mathrm{d}\ln\,T /\mathrm{d}\,x\right)}^{-1}$, 其中$ n $$ T $分别为等离子体密度和温度. 由图1可以看出, 随着温度梯度增大, ITG的波数谱变得越来越宽, 增长率越来越大, 高波数的ITG变得更不稳定. 而随着密度梯度的增大, ITG的波数谱逐渐展宽, 最大增长率对应的波数往高波数区移动. 总体来看, 随密度梯度的增大, ITG增长率则先增大再减小. 由图1还可以看出, ITG的激发存在一个温度梯度阈值, 只有当温度梯度大于该阈值时, ITG才能被激发. 这个温度梯度阈值的大小强烈依赖于密度梯度. 鉴于实验上可能出现中空的密度剖面, 在CFQS位形下开展了负密度梯度对ITG影响的研究. 由图1(a)图1(b)可以看出, 虽然${R}_{0}/{L}_{n} = -2$${R}_{0}/{L}_{n} = 2$这两个密度梯度的绝对值相同, 但对应ITG的增长率却明显不同, 负密度梯度对ITG具有抑制作用. 这与托卡马克上的结果相同[19,20]. 因此, ITG的增长率不仅与密度梯度的绝对值相关, 还取决于密度梯度值的正负. 图2给出了图1条件中ITG的频率. ITG的传播方向为离子抗磁漂移方向, 在本文模型中, 以电子的抗磁漂移方向为正方向, 因此ITG的频率为负值. 由图2可以看出, 温度梯度对ITG频率的影响不大, 但密度梯度对ITG频率的影响很大. ITG的频率随密度梯度数值的增大而减小, 这跟增长率与密度梯度之间的依赖关系不同, ITG的频率只与密度梯度的数值大小有关, 密度梯度的数值越大, ITG的频率越小. 在托卡马克中也给出了相同的结果[20]. 结合图1, 对于同一个密度梯度, ITG越不稳定, 其频率越大. 对于固定的密度梯度和温度梯度, 频率与波数近似成线性关系, 即ITG的相速度近似为常数. 这与托卡马克和仿星器(如NCSX)中的结果一致[21,22]. ITG的相速度与温度梯度成正比, 与密度梯度成反比. 当密度梯度大时, 如图2中的$ {R}_{0}/{L}_{n}=8 $, 波数${k}_{y}{\rho }_{\mathrm{i}}\leqslant 1$时, ITG的频率趋近于0.

      图  1  不同密度梯度和温度梯度下绝热ITG的增长率波数谱 (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $

      Figure 1.  Growth rate spectra of ITG for different density gradients and temperature gradients: (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $.

      图  2  不同密度梯度和温度梯度下绝热ITG的频率波数谱 (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $

      Figure 2.  Real frequency spectra of ITG for different density gradients and temperature gradients: (a) $ {R}_{0}/{L}_{n}=-2 $; (b) $ {R}_{0}/{L}_{n}=2 $; (c) $ {R}_{0}/{L}_{n}=8 $.

      图3给出了CFQS中绝热ITG的最大增长率与密度梯度和温度梯度间的关系. 由图3可以看出, 当温度梯度最大、密度梯度在${R}_{0}/{L}_{n}\sim3$附近时, ITG最不稳定. 随着密度梯度的增大或者减小, ITG的增长率都减小. 这意味着有限的密度梯度既可以驱动ITG(当密度梯度较小时), 也可以抑制ITG(当密度梯度大时或者密度梯度为负值时). 由图3可以容易地看出, ITG的激发存在一个温度梯度阈值, 该阈值强烈受到密度梯度的影响. 对于CFQS, 激发ITG的温度梯度阈值$ {R}_{0}/{L}_{{T}_{\mathrm{c}}} > 3 $. 对于ITG, $ {\eta }_{\mathrm{i}}={L}_{n}/{L}_{T} $是一个重要的参量. 当密度梯度不变, 温度梯度增大($ {L}_{T} $减小)时, $ {\eta }_{\mathrm{i}} $随温度梯度增大而增大, 从图3可以看出, ITG的增长率也逐渐增大. 而若温度梯度不变, 密度梯度减小($ {L}_{n} $增大), 则$ {\eta }_{\mathrm{i}} $随密度梯度减小而增大. 从图3可以看出, 此时ITG的增长率随$ {\eta }_{\mathrm{i}} $增大不是单调变化的, 而是先增大再减小. 这与托卡马克和其他仿星器(如NCSX和LHD)中的结果相同[19,20,22-25]. 文献[20]给出了托卡马克等离子体中若固定密度梯度, 则ITG的增长率随着温度梯度增大(或$ {\eta }_{\mathrm{i}} $增大)而增大, 此外也给出了ITG增长率随密度梯度和温度梯度变化的等高线图, 这都与CFQS中的结果是一致的. 在仿星器中, 如NCSX上, 也给出了当固定温度梯度, ITG的增长率随密度梯度减小(或$ {\eta }_{\mathrm{i}} $增大)先增大再减小的结果[22].

      图  3  绝热ITG的最大增长率与密度梯度和温度梯度的关系

      Figure 3.  Relationship of the maximum growth rate of adiabatic ITG to $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{T} $.

      若在模型中考虑动理学电子效应, 则可以研究存在捕获电子时ITG的激发和稳定条件, 并研究捕获电子对TIG的影响. 但是, 考虑了捕获电子效应后, 可能会激发另一种微观不稳定性—捕获电子模(TEM). 通常, TEM的空间尺度比ITG小, 而本文研究低波数区的微观不稳定性, 因此微观不稳定性以ITG为主[26]. 但当ITG的驱动比较小时, 可能出现ITG和TEM共同存在的情况, 也可能出现耦合的TE-ITG模, 这个耦合的模在托卡马克和仿星器中都可以出现[27-31].

      存在捕获电子效应后, ITG的增长率和频率波数谱如图4所示. 由前面绝热的结果知, $ {R}_{0}/{L}_{n}= $2时ITG最不稳定, 故图4中等离子体的密度梯度取$ {R}_{0}/{L}_{n}=2 $. 与图1(b)图2(b)相比较, 发现在考虑捕获电子效应后, ITG的增长率增大. 这意味着捕获电子能使ITG变得更不稳定, 降低了激发ITG的温度梯度阈值[22,32]. 捕获电子对ITG的这种效应不依赖于装置的类型, 即与装置是托卡马克还是仿星器没有关系. 与绝热的结果一样, ITG的增长率越大, 其相速度也越大, 且ITG的相速度对同一密度梯度和温度梯度近似为常数. 但考虑了捕获电子后, ITG的频率和相速度都比绝热的结果小. 图5给出了当考虑捕获电子效应后ITG最大增长率与密度梯度和温度梯度间的关系. 由图5可看出: 在大的温度梯度区间, ITG最大增长率与密度和温度梯度的关系与绝热的结果(图3)相似, 只是增长率的数值更大一些; 但在低的温度梯度区, 其结果与绝热的结果存在一定的差异, 可能的原因是当温度梯度小的时候, 无法驱动出具有大增长率的ITG. 大的密度梯度一方面可以抑制ITG, 另一方面又可以驱动TEM, 产生所谓的密度梯度驱动的TEM[27]; 此外, 在该参数区间ITG和TEM的驱动都不强, 有可能形成耦合的TE-ITG模[28-31]. 故在该参数区间里ITG不是主导的不稳定性, 此时以耦合的TE-ITG模或混合的ITG和TEM为主. 这种耦合的TE-ITG模或混合的ITG和TEM在托卡马克和其他位形的仿星器中都能产生[27-31].

      图  4  考虑捕获电子效应后ITG的增长率(a)和频率(b)波数谱, 其中$ {R}_{0}/{L}_{n}=2 $, 电子温度梯度标长和离子温度梯度标长相等, 即$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $

      Figure 4.  Growth rate (a) and real frequency spectra (b) of kinetic ITG for $ {R}_{0}/{L}_{n}=2 $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $.

      图  5  考虑捕获电子效应后ITG/TE-ITG最大增长率与密度梯度和温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}} $

      Figure 5.  Contour map of the maximum growth rate of kinetic ITG/TE-ITG mode vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{T} $. Here, ${R}_{0}/ {L}_{{T}_{\mathrm{e}}} $$ ={R}_{0}/{L}_{{T}_{\mathrm{i}}}$.

      实际上, 在真实的磁约束等离子体中, 电子的温度容易受外界电子回旋共振加热(ECRH)的影响. CFQS未来的加热手段采用的是ECRH, 这导致CFQS等离子体中电子温度梯度和离子温度梯度的不同. 在前面的研究中, 电子温度梯度和离子温度梯度相等. 为研究不相同的电子温度梯度和离子温度梯度对CFQS中ITG的影响, 以及电子温度梯度如何影响ITG, 在下面的模拟中, 固定电子的温度梯度为$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $.

      图6给出了$ {R}_{0}/{L}_{n}=2 $时ITG增长率和频率的波数谱. 显然, 图6$ {R}_{0}/{L}_{{{T}}_{\mathrm{i}}}=8 $的曲线与图4$ {R}_{0}/{L}_{T} = 8 $曲线是完全相同的. 与图4相比, 当$ {R}_{0}/ {L}_{{T}_{\mathrm{i}}} > 8 $ 时, ITG的增长率略有降低, 这个降低主要体现在高波数的ITG上; 当$ {R}_{0}/{L}_{{T}_{\mathrm{i}}} < 8 $时, ITG的增长率略有增大; 而对于ITG的频率, 变化规律则相反. 这说明电子温度梯度可以影响ITG, 它对ITG具有一定的去稳作用, 特别是对高波数的ITG. 在托卡马克中也给出相同的结论[33]. 图7给出了当电子温度梯度固定时, ITG的最大增长率与密度梯度和离子温度梯度间的关系. 图7的结果与绝热的结果(图3)相似, 与图5的结果略有不同. 这主要是因为$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $能够驱动大的TEM, 故当离子温度梯度小时, 不稳定性以TEM为主; TEM的驱动大, 无法形成耦合的TE-ITG模[29].

      图  6  考虑捕获电子效应后ITG的增长率(a)和频率波数谱(b), 其中$ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

      Figure 6.  Growth rate (a) and real frequency spectra (b) of kinetic ITG for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $.

      图  7  考虑捕获电子效应后ITG最大增长率与密度梯度和离子温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

      Figure 7.  Contour map of the maximum growth rate of kinetic ITG vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}} $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $.

    • 在前面的模拟中, 设定等离子体比压 $\beta =0{\text{%}}$, 即模拟静电条件下ITG的特征. 如果在模型中考虑有限的比压, 则可以研究电磁条件下ITG的特征, 并获知在CFQS中有限比压对ITG的影响. 等离子体的比压定义为等离子体的热压与磁压的比值, 即$ \beta ={\mu }_{0}nT/{B}^{2} $, 其中$ {\mu }_{0} $是真空磁导率, $ B $是磁场. 图8给出了离子尺度($ {k}_{y}{\rho }_{\mathrm{i}}=1.0 $)的微观不稳定性的增长率和频率与等离子体比压间的关系. 由图8(a)可以看到, 增长率在$\beta =0.3{\text{%}}$处有个拐点, 在图中用蓝色箭头指出. 在拐点左侧, 不稳定性为ITG, 其增长率随比压的增大而减小, 而在拐点右侧, 不稳定性为动理学气球模(KBM), 其增长率随比压的增大而增大[34-36]. 离子尺度微观不稳定性的这种特征在托卡马克中同样也是存在的[34,35], 只是模拟的条件不同, 导致转换点处对应的比压值也不同. 其实, KBM也可以看成是电磁版本的ITG, 即AITG[37-39]. 考虑了电磁条件, $\beta \cong 0.3{\text{%}}$是ITG转变为AITG/KBM的转变点. 因此, 考虑有限比压效应后, 离子尺度的微观不稳定性以AITG/KBM为主.

      图  8  当波数$ {k}_{y}{\rho }_{\mathrm{i}}=1.0 $时增长率(a)与频率(b)随比压的变化(其中, $ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}={R}_{0}/{L}_{{T}_{\mathrm{i}}}=8 $), 图(a)中的箭头是从ITG转变为KBM的转变点

      Figure 8.  Growth rates (a) and real frequencies (b) vs. $ \beta $ for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}}={R}_{0}/{{L}}_{{{T}}_{\mathrm{e}}}=8 $ at $ {k}_{y}{\rho }_{\mathrm{i}}=1.0 $. The arrow is plotted in panel (a) to point the transition point from ITG to KBM.

      下面取$\beta =1{\text{%}}$, 由图8可知, 此时离子尺度的微观不稳定性以AITG/KBM为主. 图9给出了$ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $时KBM的增长率和频率波数谱. 与前面静电的模拟结果不同, 电磁条件下, 波数小的微观不稳定性具有更大的增长率, 增长率随波数近似成单调递减的关系. 但从图9也可以看出, 当波数$ {k}_{y}{\rho }_{\mathrm{i}}\lesssim 0.1 $时, 增长率随着波数减小而减小. 在托卡马克和其他位形的仿星器中也有类似的KBM增长率波数谱[35,40]. 而频率的大小与静电的结果类似, 随着波数的增大而增大, 但在电磁条件下, 频率与波数不再具有近似线性的关系, 即相速度不再近似为常数, 而是随着波数的增大而减小. 图10给出了考虑有限比压效应后KBM的最大增长率与密度梯度和温度梯度间的关系. 由于KBM的驱动源为压强梯度, 即大的密度梯度和大的温度梯度都能驱动大的KBM, 因此, 随着密度梯度和温度梯度的增大, KBM的增长率也逐渐增大. 由图10可以看出, 当密度梯度和温度梯度都较大时($ {R}_{0}/{L}_{n, {T}_{\mathrm{i}}}\gtrsim 4 $), KBM的最大增长率与密度梯度和温度梯度近似成线性关系. 与前面静电ITG的结果不同, 负密度梯度对KBM没有明显的抑制作用.

      图  9  取$\beta =1{\text{%}}$时, KBM的增长率(a)和频率(b)波数谱, 其中$ {R}_{0}/{L}_{n}=2 $, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

      Figure 9.  Growth rate (a) and real frequency spectra (b) of KBM for $ {R}_{0}/{L}_{n}=2 $ and $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $. Here, $\beta =1{\text{%}}$.

      图  10  取$\beta =1{\text{%}}$时, KBM的最大增长率与密度梯度和离子温度梯度的关系, 其中$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $

      Figure 10.  Contour map of the maximum growth rate of KBM vs. $ {R}_{0}/{L}_{n} $ and $ {R}_{0}/{L}_{{T}_{\mathrm{i}}} $. Here, $ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $ and $\beta =1{\text{%}}$.

    • 本文研究了CFQS中静电和电磁ITG, 给出了不同等离子体密度梯度和温度梯度条件下的ITG特征, 对CFQS上离子尺度的微观不稳定性有了一个初步的认识.

      在静电绝热条件下, 随着温度梯度的增大, ITG的增长率和频率都增大, ITG的波数谱越来越宽, 高波数的ITG变得不稳定. ITG的激发存在一个温度梯度阈值, 只有当温度梯度大于该阈值时, 才能出现ITG. 对于CFQS, 其温度梯度阈值$ {R}_{0}/{L}_{{T}_{\mathrm{c}}} > 3 $. 随着密度梯度的增大, ITG的增长率先增大再减小, 而ITG的频率单调减小. ITG的增长率不仅与密度梯度的绝对值相关, 还取决于密度梯度的正负, 负密度梯度对ITG具有强的抑制作用. 而ITG的频率只与密度梯度的数值大小有关, 密度梯度的数值越大, ITG的频率越小. 在CFQS中, 当密度梯度$ {R}_{0}/{L}_{n}\approx 2.5 $, 温度梯度取最大时, ITG具有最大的增长率. 当固定密度梯度和温度梯度, 不同波数的ITG的相速度近似为常数. ITG的相速度随温度梯度的增大而增大, 随密度梯度的增大而减小

      在考虑捕获电子效应后, 与绝热的结果相比, ITG的增长率增大, 意味着捕获电子对ITG具有去稳作用, 而ITG的频率略有降低, 相速度减小. 当温度梯度大时, ITG的最大增长率与密度梯度和温度梯度的依赖关系同绝热的结果相同. 但当温度梯度小时, 会出现耦合的TE-ITG模或者由大的密度梯度驱动的TEM, 导致在低温度梯度区以耦合的TE-ITG模或混合的ITG和TEM为主. 当固定电子温度梯度$ {R}_{0}/{L}_{{T}_{\mathrm{e}}}=8 $, 可以获知电子温度梯度和离子温度梯度对CFQS中ITG的影响. 电子温度梯度对ITG具有一定的去稳作用. 由于此时TEM的驱动较强, 无法形成耦合的TE-ITG模.

      当考虑电磁条件时, 有限的等离子体比压会抑制ITG, 随着等离子体比压的增大, 在$\beta \sim 0.3{\text{%}}$时出现了ITG向AITG/KBM的转化. 此后, 随着比压继续增大, KBM的增长率也增大. KBM增长率的波数谱与ITG的不同, KBM的增长率随着波数近似单调下降. 随着密度梯度和离子温度梯度的增大, KBM的增长率也逐渐增大. 当密度和温度梯度都较大时($ {R}_{0}/{L}_{n, {T}_{\mathrm{i}}}\gtrsim 4 $), KBM的最大增长率与密度梯度和温度梯度近似成线性关系.

      CFQS仿星器是三维平衡, 不同环向位置的磁场位形不同[11]. 取不同环向位置构建的通量管的磁场位形不同, 这对ITG的模拟结果会有影响. 在本文中, 在环向角为0°处取一个通量管作为模拟区域, 在此区域中ITG的增长率最大, 在半个环向周期内随着环向角增长, ITG的增长率逐渐减小.

参考文献 (40)

目录

    /

    返回文章
    返回