Search

Article

x

留言板

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

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

Numerical simulation of synergistic effect of neoclassical tearing mode and toroidal field ripple on alpha particle loss in China Fusion Engineering Testing Reactor

Hao Bao-Long Chen Wei Li Guo-Qiang Wang Xiao-Jing Wang Zhao-Liang Wu Bin Zang Qing Jie Yin-Xian Lin Xiao-Dong Gao Xiang CFETR TEAM

Citation:

Numerical simulation of synergistic effect of neoclassical tearing mode and toroidal field ripple on alpha particle loss in China Fusion Engineering Testing Reactor

Hao Bao-Long, Chen Wei, Li Guo-Qiang, Wang Xiao-Jing, Wang Zhao-Liang, Wu Bin, Zang Qing, Jie Yin-Xian, Lin Xiao-Dong, Gao Xiang, CFETR TEAM
PDF
HTML
Get Citation
  • 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_\zeta },\mu $), 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_\zeta },\mu $), 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.
      Corresponding author: Hao Bao-Long, blhao@ipp.ac.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11905142, 11875290)
    [1]

    Zhuang G, Li G Q, Li J, et al. 2019 Nucl. Fusion 59 112010Google Scholar

    [2]

    Wan Y X, Li J, Liu Y, et al. 2017 Nucl. Fusion 57 102009Google Scholar

    [3]

    Chen J, Jian X, Chan V, et al. 2017 Plasma Phys. Controlled Fusion 59 075005Google Scholar

    [4]

    Gorelenkov N N, Pinches S D, Toi K, et al. 2014 Nucl. Fusion 54 125001Google 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 021807Google Scholar

    [7]

    White R B, Rutherford P H, Colestock P, et al. 1988 Phys. Rev. Lett. 60 2038Google Scholar

    [8]

    Chapman I T 2011 Plasma Phys. Controlled Fusion 53 013001Google 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 032501Google Scholar

    [11]

    Gobbin M, Marrelli L, Fahrbach H, et al. 2009 Nucl. Fusion 49 095021Google Scholar

    [12]

    Strumberger E, Gunter S, Schwarz E, et al. 2008 New J. Phys. 10 023017Google Scholar

    [13]

    Garcia-Munoz M, Fahrbach H, Pinches S D, et al. 2009 Nucl. Fusion 49 085014Google Scholar

    [14]

    Garcia-Munoz M, Martin P, Fahrbach H, et al. 2007 Nucl. Fusion 47 L10Google Scholar

    [15]

    Mynick H E 1993 Phys. Fluids B 5 1471Google Scholar

    [16]

    郝保龙, 陈伟, 蔡辉山, 等 2020 中国科学: 物理学 力学 天文学 50 065201Google Scholar

    Hao B L, Chen W, Cai H S, et al. 2020 Sci. Sin. Phys. Mech. Astron. 50 065201Google 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 157Google Scholar

    [19]

    高翔, 万宝年, 宋云涛, 等 2019 中国科学: 物理学 力学 天文学 49 045202Google Scholar

    Gao X, Wan B N, Song Y T, et al. 2019 Sci. Sin.-Phys. Mech. Astron. 49 045202Google 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 016045Google Scholar

    [22]

    Hao B L, White R B, Gao X, et al. 2019 Nucl. Fusion 59 076040Google Scholar

    [23]

    Wu B, Hao B L, White R B, et al. 2017 Plasma Phys. Controlled Fusion 59 025004Google Scholar

    [24]

    Boozer A H, Kuo-Petravic G 1981 Phys. Fluids 24 851Google Scholar

    [25]

    Zhao R, Wang Z X, Wang F, et al. 2020 Plasma Phys. Controlled Fusion 62 115001Google Scholar

    [26]

    Carolipio E M, Heidbrink W W, Forest C B, et al. 2002 Nucl. Fusion 42 853Google Scholar

  • 图 1  CFETR混杂运行模式(v201806)中背景电子密度、温度和离子温度分布

    Figure 1.  Distributions of bulk electron density, electron temperature, and bulk ion temperature in CFETR hybrid scenario (v201806).

    图 2  CFETR初始alpha粒子密度分布

    Figure 2.  Density profile of alpha particle initial distribution in CFETR.

    图 3  CFETR纵场波纹度分布 (a) 解析函数实现值; (b) 工程设计值

    Figure 3.  Distributions of toroidal field ripple perturbation amplitude in CFETR: (a) Ripple data by analytical equation; (b) engineering data in design.

    图 4  CFETR混杂运行模式 (a) 2/1和3/2 NTM的平衡分量和安全因子分布; (b), (c)一阶分量和高斯拟合曲线

    Figure 4.  (a) Safety factor profile and equilibrium helical flux of 2/1 and 3/2 NTM in CFETR hybrid scenario; (b), (c) the first order harmonics and Gaussian fitting results.

    图 5  固定2/1和3/2模扰动幅度比例下扫描NTM扰动幅度下的庞加莱图 (a) ${\alpha _{2/1}} = 1.64 \times {10^{ - 4}}{R_0}$; (b) ${\alpha _{2/1}} = l $$ 6.00 \times {10^{ - 4}}{R_0}$

    Figure 5.  Poincare plot under different NTM perturbation amplitude with the ratio fixed ${{{\alpha _{2/1}}} / {{\alpha _{3/2}}}} = 4.38$: (a)${\alpha _{2/1}} = $$ 1.64 \times {10^{ - 4}}{R_0}$; (b)${\alpha _{2/1}} = 6.00 \times {10^{ - 4}}{R_0}$.

    图 6  在半个EP慢化时间内不同物理效应下计算得到的粒子损失份额

    Figure 6.  The particle loss fraction after following half a slowing down time under different perturbation and collision.

    图 7  纵场波纹扰动下损失alpha粒子信息 (a) 损失份额随时演化; (b) 损失粒子能量分布; (c) 损失粒子螺距角分布

    Figure 7.  Information of lost alpha particles under collision and toroidal field ripple: (a) Evolution of loss fraction; (b) energy distribution of lost particles; (c) pitch angle distribution of lost particles.

    图 8  一个慢化时间后初始分布alpha粒子波纹损失局域沉积在LCFS处得到的热负荷

    Figure 8.  The heat load at the last closed flux surface due to ripple loss of initial alpha particle distribution after a slowing down time.

    图 14  初始分布alpha粒子香蕉轨道漂移频率分布 (a) 极向回弹频率; (b) 环向进动频率

    Figure 14.  The drift frequency distributions of banana orbit of initial alpha particles: (a) Poloidal bounce frequency; (b) toroidal precession frequency.

    图 9  同向alpha粒子动理学庞加莱图 (a) 通行粒子; (b) 俘获粒子

    Figure 9.  The kinetic Poincare plot of co-current alpha particles: (a) Passing particle; (b) trapped particle.

    图 10  半个慢化时间内不同磁场波纹度和NTM扰动幅度下的初始分布alpha粒子损失份额

    Figure 10.  The loss fraction of initial alpha particles under different toroidal ripple and NTM perturbation amplitude.

    图 11  CFETR平衡位形中不同波纹度形成的波纹磁阱 (a) 工程设计值; (b) 5倍波纹度分布; (c) 9倍波纹度分布

    Figure 11.  The ripple well domain in CFETR equilibrium with different ripple amplitude: (a) Distribution with engineering design; (b) 5 times of ripple; (c) 9 times of ripple.

    图 12  (a) (${P_{\zeta} }, {{\mu {B_0}} / E}$)平面初始alpha粒子轨道类型分布; (b) 正值螺距角粒子分布; (c) 负值螺距角粒子分布

    Figure 12.  (a) Orbit classification in the plane of $({P_{\zeta} }, {{\mu {B_0}} / E})$ with initial alpha particles; (b) particles with positive pitch angle; (c) particles with negative pitch angle.

    图 13  不同倍数波纹度分布时的波纹磁阱和随机波纹扩散区域 (a), (b) 工程设计值; (c), (d) 5倍波纹度分布; (e), (f) 9倍波纹度分布

    Figure 13.  The ripple well trapping and stochastic ripple diffusion domain with different times of ripple distribution: (a), (b) Ripple in engineering design; (c), (d) 5 times of ripple; (e), (f) 9 times of ripple.

    表 1  CFETR与其他托卡马克装置主机参数对比

    Table 1.  Main parameters comparison of CFETR and other tokomak facilities.

    参数CFETRITEREAST
    磁轴场强${B_{{\rm{T}}0}}$/T6.55.32
    等体大半径${R_0}$/m7.26.21.9
    等体小半径a/m2.22.00.5
    等体电流${I_{\rm{p}}}$/MA14151
    纵场磁体柄数N161816
    DownLoad: CSV
    Baidu
  • [1]

    Zhuang G, Li G Q, Li J, et al. 2019 Nucl. Fusion 59 112010Google Scholar

    [2]

    Wan Y X, Li J, Liu Y, et al. 2017 Nucl. Fusion 57 102009Google Scholar

    [3]

    Chen J, Jian X, Chan V, et al. 2017 Plasma Phys. Controlled Fusion 59 075005Google Scholar

    [4]

    Gorelenkov N N, Pinches S D, Toi K, et al. 2014 Nucl. Fusion 54 125001Google 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 021807Google Scholar

    [7]

    White R B, Rutherford P H, Colestock P, et al. 1988 Phys. Rev. Lett. 60 2038Google Scholar

    [8]

    Chapman I T 2011 Plasma Phys. Controlled Fusion 53 013001Google 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 032501Google Scholar

    [11]

    Gobbin M, Marrelli L, Fahrbach H, et al. 2009 Nucl. Fusion 49 095021Google Scholar

    [12]

    Strumberger E, Gunter S, Schwarz E, et al. 2008 New J. Phys. 10 023017Google Scholar

    [13]

    Garcia-Munoz M, Fahrbach H, Pinches S D, et al. 2009 Nucl. Fusion 49 085014Google Scholar

    [14]

    Garcia-Munoz M, Martin P, Fahrbach H, et al. 2007 Nucl. Fusion 47 L10Google Scholar

    [15]

    Mynick H E 1993 Phys. Fluids B 5 1471Google Scholar

    [16]

    郝保龙, 陈伟, 蔡辉山, 等 2020 中国科学: 物理学 力学 天文学 50 065201Google Scholar

    Hao B L, Chen W, Cai H S, et al. 2020 Sci. Sin. Phys. Mech. Astron. 50 065201Google 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 157Google Scholar

    [19]

    高翔, 万宝年, 宋云涛, 等 2019 中国科学: 物理学 力学 天文学 49 045202Google Scholar

    Gao X, Wan B N, Song Y T, et al. 2019 Sci. Sin.-Phys. Mech. Astron. 49 045202Google 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 016045Google Scholar

    [22]

    Hao B L, White R B, Gao X, et al. 2019 Nucl. Fusion 59 076040Google Scholar

    [23]

    Wu B, Hao B L, White R B, et al. 2017 Plasma Phys. Controlled Fusion 59 025004Google Scholar

    [24]

    Boozer A H, Kuo-Petravic G 1981 Phys. Fluids 24 851Google Scholar

    [25]

    Zhao R, Wang Z X, Wang F, et al. 2020 Plasma Phys. Controlled Fusion 62 115001Google Scholar

    [26]

    Carolipio E M, Heidbrink W W, Forest C B, et al. 2002 Nucl. Fusion 42 853Google Scholar

  • [1] Zhang Qi-Fan, Le Wen-Cheng, Zhang Yu-Hao, Ge Zhong-Xin, Kuang Zhi-Qiang, Xiao Sheng-Yang, Wang Lu. Effects of radiation from tungsten impurities on the thermal energy loss during the fast thermal quench stage of major disruption in tokamak plasmas. Acta Physica Sinica, 2024, 73(18): 185201. doi: 10.7498/aps.73.20240730
    [2] Wang Fu-Qiong, Xu Ying-Feng, Zha Xue-Jun, Zhong Fang-Chuan. Multi-fluid and dynamic simulation of tungsten impurity in tokamak boundary plasma. Acta Physica Sinica, 2023, 72(21): 215213. doi: 10.7498/aps.72.20230991
    [3] Liu Guan-Nan, LI Xin-Xia, Liu Hong-Bo, Sun Ai-Ping. Synergistic current drive of helicon wave and lower hybrid wave in HL-2M. Acta Physica Sinica, 2023, 72(24): 245202. doi: 10.7498/aps.72.20231077
    [4] Wu Xiang-Feng, Wang Feng, Lin Zhan-Hong, Chen Luo-Yu, Yu Zhao-Ke, Wu Kai-Bang, Wang Zheng-Xiong. Numerical simulation of $\boldsymbol \alpha$ particle slowing-down process under CFETR scenario. Acta Physica Sinica, 2023, 72(21): 215209. doi: 10.7498/aps.72.20230700
    [5] Shen Yong, Dong Jia-Qi, He Hong-Da, Pan Wei, Hao Guang-Zhou. Ideal conductive wall and magnetohydrodynamic instability in Tokamak. Acta Physica Sinica, 2023, 72(3): 035203. doi: 10.7498/aps.72.20222043
    [6] Hao Bao-Long, Li Ying-Ying, Chen Wei, Hao Guang-Zhou, Gu Xiang, Sun Tian-Tian, Wang Yu-Min, Dong Jia-Qi, Yuan Bao-Shan, Peng Yuan-Kai, Shi Yue-Jiang, Xie Hua-Sheng, Liu Min-Sheng, ENN TEAM. Optimizing numerical simulation of beam ion loss due to toroidal field ripple on EXL-50U spherical torus. Acta Physica Sinica, 2023, 72(21): 215215. doi: 10.7498/aps.72.20230749
    [7] Zhu Xiao-Long, Chen Wei, Wang Feng, Wang Zheng-Xiong. Hybrid numerical simulation on fast particle transport induced by synergistic interaction of low- and medium-frequency magnetohydrodynamic instabilities in tokamak plasma. Acta Physica Sinica, 2023, 72(21): 215210. doi: 10.7498/aps.72.20230620
    [8] Liu Zhao-Yang, Zhang Yang-Zhong, Xie Tao, Liu A-Di, Zhou Chu. Group velocity in spatiotemporal representation of collisionless trapped electron mode in tokamak. Acta Physica Sinica, 2021, 70(11): 115203. doi: 10.7498/aps.70.20202003
    [9] Zhang Chong-Yang, Liu A-Di, Li Hong, Chen Zhi-Peng, Li Bin, Yang Zhou-Jun, Zhou Chu, Xie Jin-Lin, Lan Tao, Liu Wan-Dong, Zhuang Ge, Yu Chang-Xuan. Application of dual-polarization frequency-modulated microwave reflectometer to J-TEXT tokamak. Acta Physica Sinica, 2014, 63(12): 125204. doi: 10.7498/aps.63.125204
    [10] Du Hai-Long, Sang Chao-Feng, Wang Liang, Sun Ji-Zhong, Liu Shao-Cheng, Wang Hui-Qian, Zhang Ling, Guo Hou-Yang, Wang De-Zhen. Modelling of edge plasma transport during H-mode of EAST by SOLPS5.0. Acta Physica Sinica, 2013, 62(24): 245206. doi: 10.7498/aps.62.245206
    [11] Lu Hong-Wei, Zha Xue-Jun, Hu Li-Qun, Lin Shi-Yao, Zhou Rui-Jie, Luo Jia-Rong, Zhong Fang-Chuan. The effect of gas puffing on plasma during slide-away discharge in the HT-7 tokamak. Acta Physica Sinica, 2012, 61(7): 075202. doi: 10.7498/aps.61.075202
    [12] Hong Bin-Bin, Chen Shao-Yong, Tang Chang-Jian, Zhang Xin-Jun, Hu You-Jun. Study on synergy of electron-cyclotron and lower-hybrid current drive in Tokamak. Acta Physica Sinica, 2012, 61(11): 115207. doi: 10.7498/aps.61.115207
    [13] Lu Hong-Wei, Hu Li-Qun, Lin Shi-Yao, Zhong Guo-Qiang, Zhou Rui-Jie, Zhang Ji-Zong. Investigation of slide-away discharges in HT-7 tokamak. Acta Physica Sinica, 2010, 59(8): 5596-5601. doi: 10.7498/aps.59.5596
    [14] Xu Qiang, Gao Xiang, Shan Jia-Fang, Hu Li-Qun, Zhao Jun-Yu. Experimental study of large power lower hybrid current drive on HT-7 tokamak. Acta Physica Sinica, 2009, 58(12): 8448-8453. doi: 10.7498/aps.58.8448
    [15] Gong Xue-Yu, Peng Xiao-Wei, Xie An-Ping, Liu Wen-Yan. Electron cyclotron current drive under different operational regimes in tokamak plasma. Acta Physica Sinica, 2006, 55(3): 1307-1314. doi: 10.7498/aps.55.1307
    [16] Xu Wei, Wan Bao-Nian, Xie Ji-Kang. The impurity transport in HT-6M tokamak. Acta Physica Sinica, 2003, 52(8): 1970-1978. doi: 10.7498/aps.52.1970
    [17] WANG WEN-HAO, YU CHANG-XUAN, XU YU-HONG, WEN YI-ZHI, LING BI-LI, SONG MEI, WAN BAO-NIAN. MEASUREMENT OF EDGE PLASMA PARAMETERS AND THEIR ELECTROSTATIC FLUCTUATIONS ON THE HT-7 SUPERCONDUCTING TOKAMAK. Acta Physica Sinica, 2001, 50(8): 1521-1527. doi: 10.7498/aps.50.1521
    [18] ZHANG XIAN-MEI, WAN BAO-NIAN, RUAN HUAI-LIN, WU ZHEN-WEI. STUDY OF THE ELECTRON THERMAL CONDUCTIVITY OF THE OHMICALLY HEATED DISCHARGES IN THE HT-7 TOKAMAK. Acta Physica Sinica, 2001, 50(4): 715-720. doi: 10.7498/aps.50.715
    [19] WANG WEN-HAO, XU YU-HONG, YU CHANG-XUAN, WEN YI-ZHI, LING BI-LI, SONG MEI, WAN BAO-NIAN. ELECTROSTATIC FLUCTUATIONS AND TURBULENT TRANSPORT STUDIES IN THE HT-7 SUPERCONDUCTING TOKAMAK EDGE PLASMAS . Acta Physica Sinica, 2001, 50(10): 1956-1963. doi: 10.7498/aps.50.1956
    [20] SHI BING-REN. ANALYTIC STUDY OF LOWER HYBRID WAVE PROPAGATION IN TOKAMAK LHCD EXPERIMENTS. Acta Physica Sinica, 2000, 49(12): 2394-2398. doi: 10.7498/aps.49.2394
Metrics
  • Abstract views:  5122
  • PDF Downloads:  72
  • Cited By: 0
Publishing process
  • Received Date:  23 November 2020
  • Accepted Date:  16 January 2021
  • Available Online:  28 May 2021
  • Published Online:  05 June 2021

/

返回文章
返回
Baidu
map