Search

Article

x

留言板

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

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

Effects of radiation from tungsten impurities on the thermal energy loss during the fast thermal quench stage of major disruption in tokamak plasmas

Zhang Qi-Fan Le Wen-Cheng Zhang Yu-Hao Ge Zhong-Xin Kuang Zhi-Qiang Xiao Sheng-Yang Wang Lu

Citation:

Effects of radiation from tungsten impurities on the thermal energy loss during the fast thermal quench stage of major disruption in tokamak plasmas

Zhang Qi-Fan, Le Wen-Cheng, Zhang Yu-Hao, Ge Zhong-Xin, Kuang Zhi-Qiang, Xiao Sheng-Yang, Wang Lu
cstr: 32037.14.aps.73.20240730
PDF
HTML
Get Citation
  • Recent studies based on the PLT, EAST, WEST, ASDEX-upgrade, JET and other tokamaks have shown that the accumulation of heavy impurities in the core regime is unavoidable, which may lead to the degradation of the plasma confinement and even trigger the major disruptions. The plasma thermal energy loss during the major disruptions mainly occurs during the fast thermal quench (TQ) stage. However, there is no comprehensive physical explanation for the scaling of the timescale of this stage. Tungsten as high Z impurity, which will be used as the wall material in International Thermonuclear Experimental Reactor (ITER), has strong radiation power, and may affect the thermal energy loss during the fast TQ. This work considers both the thermal diffusion induced by the stochastic magnetic fields and the radiation from tungsten impurities as the dominant thermal loss mechanisms in this stage, and construct a one-dimensional model of electron temperature evolution in tokamak plasmas. We numerically calculate and analyze the evolution of the electron temperature in this stage with the typical ITER-like parameters, and here are our main conclusions: 1) The order of magnitude of the fast TQ timescale is mainly determined by the level of thermal diffusion. However, the radiation from tungsten impurities can quantitively influence on the timescale of fast TQ and the electron temperature in the late phase of fast TQ. The higher the tungsten concentration, the shorter the TQ timescale and the lower the electron temperature it will lead to in the late phase. Both the numerical and analytical results show that the timescale is approximately linear with the tungsten impurity concentration. 2) Based on the evolution of the global energy loss and the global power loss during the fast TQ, it can be found that the global thermal energy loss via the radiation from tungsten impurities is much smaller than that via the thermal diffusion induced by the stochastic magnetic fields during the early phase of fast TQ stage. However, in the late phase of fast TQ stage, the global radiation power can be comparable to or even greater than that of the global thermal diffusion power. This is also the reason why the electron temperature in the late phase of fast TQ decreases as the concentration of tungsten impurities increases. Therefore, the contribution of the radiation from tungsten impurities to the thermal loss cannot be ignored in the late phase of fast TQ.
      Corresponding author: Xiao Sheng-Yang, d202380863@hust.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 12275097, 51821005).
    [1]

    Greenwald M, Terry J L, Wolfe S M, Ejima S, Bell M G, Kaye S M, Neilson G H 1988 Nucl. Fusion 28 2199Google Scholar

    [2]

    Troyon F, Gruber R, Sauremann H, Semenzato S, Succi S 1984 Plasma Phys. Control. Fusion 26 209Google Scholar

    [3]

    Hender T C, Wesley J C, Bialek J, et al. 2007 Nucl. Fusion 47 S128Google Scholar

    [4]

    ITER Physics Expert Group On Disruptions, Plasma Control, MHD, ITER Physics Basis Editors 1999 Nucl. Fusion 39 02251Google Scholar

    [5]

    Strait E J, Lao L L, Luxon J L, Reis E E 1991 Nucl. Fusion 31 527Google Scholar

    [6]

    Gill R D 1993 Nucl. Fusion 33 1613Google Scholar

    [7]

    Xia W, Zeng L, Tang T, et al. 2023 Plasma Phys. Control. Fusion 65 085011Google Scholar

    [8]

    Riccardo V, Loarte A 2005 Nucl. Fusion 45 1427Google Scholar

    [9]

    Sweeney R, Choi W, Austin M, et al. 2018 Nucl. Fusion 58 056022Google Scholar

    [10]

    Bondeson A, Parker R D, Hugon M, Smeulders P 1991 Nucl. Fusion 31 1695Google Scholar

    [11]

    Sheikh U A, Shiraki D, Sweeney R, et al. 2021 Nucl. Fusion 61 126043Google Scholar

    [12]

    Whyte D G, Jernigan T C, Humphreys D A, et al. 2003 J. Nucl. Mater. 313 1239Google Scholar

    [13]

    Hollmann E M, Jernigan T C, Groth M, et al. 2005 Nucl. Fusion 45 1046Google Scholar

    [14]

    Isler R C 1984 Nucl. Fusion 24 1599Google Scholar

    [15]

    Hinnov E, Mattioli M 1978 Phys. Lett. A 66 109Google Scholar

    [16]

    Wang F Q, Zha X J, Duan Y M, et al. 2018 Plasma Phys. Control. Fusion 60 125005Google Scholar

    [17]

    Yang X, Manas P, Bourdelle C, et al. 2020 Nucl. Fusion 60 086012Google Scholar

    [18]

    Neu R, Dux R, Geier A, et al. 2003 J. Nucl. Mater. 313 116Google Scholar

    [19]

    Köchl F, Loarte A, de la Luna E, et al. 2018 Plasma Phys. Control. Fusion 60 074008Google Scholar

    [20]

    Neu R 2006 Phys. Scr. T123 33Google Scholar

    [21]

    Noda N, Philipps V, Neu R 1997 J. Nucl. Mater. 241 227Google Scholar

    [22]

    Pütterich T, Neu R, Dux R, et al. 2010 Nucl. Fusion 50 025012Google Scholar

    [23]

    Krommes J A, Oberman C, Kleva R G 1983 J. Plasma Phys. 30 11Google Scholar

    [24]

    Xiao S Y, Wang L 2024 Phys. Plasmas 31 042511Google Scholar

    [25]

    Pütterich T, Fable E, Dux R, O'Mullane M, Neu R, Siccinio M 2019 Nucl. Fusion 59 056013Google Scholar

    [26]

    Kallenbach A, Bernert M, Dux R, et al. 2013 Plasma Phys. Control. Fusion 55 124041Google Scholar

    [27]

    Cheng F Y, Shi B R 2007 Chinese Phys. 16 3458Google Scholar

    [28]

    Abdullaev S S, Finken K H, Wongrach K, et al. 2015 J. Plasma Phys. 81 475810501Google Scholar

    [29]

    Militello F, Naulin V, Nielsen A H 2013 Plasma Phys. Control. Fusion 55 074010Google Scholar

    [30]

    Zhu B, Xu X Q, Tang X Z 2023 Nucl. Fusion 63 086027Google Scholar

    [31]

    Shiraki D, Commaux N, Baylor L R, et al. 2015 Nucl. Fusion 55 073029Google Scholar

    [32]

    Lehnen M, Gerasimov S N, Jachmich S, et al. 2015 Nucl. Fusion 55 123027Google Scholar

    [33]

    Tong R H, Chen Z Y, Jiang Z H, et al. 2018 Rev. Sci. Instrum. 89 10E113Google Scholar

  • 图 1  钨杂质冷却系数与电子温度的关系, 红色曲线表示线性拟合后的关系, 黑色圆点表示根据ADAS所提供数据计算的结果[25]

    Figure 1.  Dependence of the cooling factor of tungsten on the electron temperature. The red curve represents the fitting results we adopt here, and the black dots denote the calculated results based on the data provided by ADAS[25].

    图 2  典型的类ITER参数下, 电子密度剖面(黑色方形)、初始电子温度剖面(蓝色圆圈)和安全因子剖面(红色三角)

    Figure 2.  Profile of electron density (black squares), initial profile of electron temperature (blue circles) and the profile of safety factor (red triangles) with typical ITER-like parameters.

    图 3  对于典型的类ITER参数, 在不同钨浓度下快速TQ阶段不同时刻的电子温度剖面 (a) $ {c_{\text{w}}} = 0 $; (b) $ {c_{\text{w}}} = 6.7 \times {10^{ - 4}} $; (c) $ {c_{\text{w}}} = 1.6 \times {10^{ - 3}} $; (d) $ {c_{\text{w}}} = 3.3 \times {10^{ - 3}} $

    Figure 3.  Profile of $ {T_{\text{e}}} $ at different times during the fast TQ stage for typical ITER-like parameters with different tungsten concentrations: (a) $ {c_{\text{w}}} = 0 $; (b) $ {c_{\text{w}}} = 6.7 \times {10^{ - 4}} $; (c) $ {c_{\text{w}}} = 1.6 \times {10^{ - 3}} $; (d) $ {c_{\text{w}}} = 3.3 \times {10^{ - 3}} $.

    图 4  对于典型的类ITER参数, 钨杂质浓度分别为$ {c_{\text{w}}} = $$ 3.3 \times {10^{ - 4}} $ (红色圆圈), $ {c_{\text{w}}} = 6.7 \times {10^{ - 4}} $ (蓝色三角), $ {c_{\text{w}}} = $$ 1.6 \times {10^{ - 3}} $ (绿色方形), $ {c_{\text{w}}} = 3.3 \times {10^{ - 3}} $ (紫色菱形)和无钨杂质对照组$ {c_{\text{w}}} = 0 $(黑线)的芯部温度$ {T_{{\text{e, core}}}} $演化

    Figure 4.  Evolution of $ {T_{{\text{e, core}}}} $ with typical ITER-like parameters under $ {c_{\text{w}}} = 3.3 \times {10^{ - 4}} $ (red circles), $ {c_{\text{w}}} = 6.7 \times $$ {10^{ - 4}} $ (blue triangles), $ {c_{\text{w}}} = 1.6 \times {10^{ - 3}} $ (green squares), $ {c_{\text{w}}} = 3.3 \times {10^{ - 3}} $ (purple diamonds) and $ {c_{\text{w}}} = 0 $ (black line), respectively.

    图 5  $ {\tau _{90 \text{-} 20}} $对钨杂质浓度的依赖关系

    Figure 5.  Dependence of $ {\tau _{90 \text{-} 20}} $ on the tungsten impurity concentration.

    图 6  对于典型类ITER参数, 在r/a = 0.5时, 热扩散功率(红色点虚线)、钨杂质辐射功率(蓝色实线)和总损失功率(黑色短划线)随电子温度的演化

    Figure 6.  Dependence of the thermal diffusive power (red dash-dotted line), the radiation power (blue solid line), and the total power of the thermal loss (black dashed line) on the electron temperature at r/a = 0.5 for typical ITER-like parameters.

    图 7  对于典型类ITER参数, 快速TQ 阶段(a)损失的全局总热能(黑色三角)、全局辐射损失能量(红色方形)和全局热扩散损失能量(蓝色圆圈)随时间的演化; (b)全局总损失功率(黑色三角)、全局辐射损失功率(红色方形)和全局扩散损失功率(蓝色圆圈)随时间的演化

    Figure 7.  For typical ITER-like parameters, (a) the time evolution of the global total thermal energy loss (black triangles), the global energy loss through radiation (red squares) and diffusion (blue circles); (b) the time evolution of the global total power loss (black triangles), the global power loss through radiation (red squares) and diffusion (blue circles) in the tokamak during the fast TQ.

    表 1  本工作中用于定量计算电子温度演化的类ITER典型参数

    Table 1.  Typical ITER-like parameters used in this work to quantify the evolution of electron temperature.

    参数
    名称
    大半径
    R/m
    小半径
    a/m
    环向磁场
    B/T
    极向
    模数m
    归一化磁
    扰动幅度b
    参数值 6.2 2 5.3 30 2×10–3
    DownLoad: CSV

    表 2  对于典型的类ITER参数情况, 不同钨杂质浓度下的快速TQ时间尺度$ {\tau _{90 \text{-} 20}} $

    Table 2.  For the typical ITER-like parameters, the timescale of fast TQ, $ {\tau _{90 \text{-} 20}} $ with different tungsten impurity concentrations.

    $ {c_{\text{w}}} $ 0 3.3×10–4 6.7×10–4 1.6×10–3 3.3×10–3
    ${\tau _{90 \text{-} 20}}/{\rm ms} $ 1.709 1.646 1.589 1.458 1.284
    DownLoad: CSV
  • [1]

    Greenwald M, Terry J L, Wolfe S M, Ejima S, Bell M G, Kaye S M, Neilson G H 1988 Nucl. Fusion 28 2199Google Scholar

    [2]

    Troyon F, Gruber R, Sauremann H, Semenzato S, Succi S 1984 Plasma Phys. Control. Fusion 26 209Google Scholar

    [3]

    Hender T C, Wesley J C, Bialek J, et al. 2007 Nucl. Fusion 47 S128Google Scholar

    [4]

    ITER Physics Expert Group On Disruptions, Plasma Control, MHD, ITER Physics Basis Editors 1999 Nucl. Fusion 39 02251Google Scholar

    [5]

    Strait E J, Lao L L, Luxon J L, Reis E E 1991 Nucl. Fusion 31 527Google Scholar

    [6]

    Gill R D 1993 Nucl. Fusion 33 1613Google Scholar

    [7]

    Xia W, Zeng L, Tang T, et al. 2023 Plasma Phys. Control. Fusion 65 085011Google Scholar

    [8]

    Riccardo V, Loarte A 2005 Nucl. Fusion 45 1427Google Scholar

    [9]

    Sweeney R, Choi W, Austin M, et al. 2018 Nucl. Fusion 58 056022Google Scholar

    [10]

    Bondeson A, Parker R D, Hugon M, Smeulders P 1991 Nucl. Fusion 31 1695Google Scholar

    [11]

    Sheikh U A, Shiraki D, Sweeney R, et al. 2021 Nucl. Fusion 61 126043Google Scholar

    [12]

    Whyte D G, Jernigan T C, Humphreys D A, et al. 2003 J. Nucl. Mater. 313 1239Google Scholar

    [13]

    Hollmann E M, Jernigan T C, Groth M, et al. 2005 Nucl. Fusion 45 1046Google Scholar

    [14]

    Isler R C 1984 Nucl. Fusion 24 1599Google Scholar

    [15]

    Hinnov E, Mattioli M 1978 Phys. Lett. A 66 109Google Scholar

    [16]

    Wang F Q, Zha X J, Duan Y M, et al. 2018 Plasma Phys. Control. Fusion 60 125005Google Scholar

    [17]

    Yang X, Manas P, Bourdelle C, et al. 2020 Nucl. Fusion 60 086012Google Scholar

    [18]

    Neu R, Dux R, Geier A, et al. 2003 J. Nucl. Mater. 313 116Google Scholar

    [19]

    Köchl F, Loarte A, de la Luna E, et al. 2018 Plasma Phys. Control. Fusion 60 074008Google Scholar

    [20]

    Neu R 2006 Phys. Scr. T123 33Google Scholar

    [21]

    Noda N, Philipps V, Neu R 1997 J. Nucl. Mater. 241 227Google Scholar

    [22]

    Pütterich T, Neu R, Dux R, et al. 2010 Nucl. Fusion 50 025012Google Scholar

    [23]

    Krommes J A, Oberman C, Kleva R G 1983 J. Plasma Phys. 30 11Google Scholar

    [24]

    Xiao S Y, Wang L 2024 Phys. Plasmas 31 042511Google Scholar

    [25]

    Pütterich T, Fable E, Dux R, O'Mullane M, Neu R, Siccinio M 2019 Nucl. Fusion 59 056013Google Scholar

    [26]

    Kallenbach A, Bernert M, Dux R, et al. 2013 Plasma Phys. Control. Fusion 55 124041Google Scholar

    [27]

    Cheng F Y, Shi B R 2007 Chinese Phys. 16 3458Google Scholar

    [28]

    Abdullaev S S, Finken K H, Wongrach K, et al. 2015 J. Plasma Phys. 81 475810501Google Scholar

    [29]

    Militello F, Naulin V, Nielsen A H 2013 Plasma Phys. Control. Fusion 55 074010Google Scholar

    [30]

    Zhu B, Xu X Q, Tang X Z 2023 Nucl. Fusion 63 086027Google Scholar

    [31]

    Shiraki D, Commaux N, Baylor L R, et al. 2015 Nucl. Fusion 55 073029Google Scholar

    [32]

    Lehnen M, Gerasimov S N, Jachmich S, et al. 2015 Nucl. Fusion 55 123027Google Scholar

    [33]

    Tong R H, Chen Z Y, Jiang Z H, et al. 2018 Rev. Sci. Instrum. 89 10E113Google Scholar

  • [1] 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
    [2] 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
    [3] 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
    [4] 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
    [5] 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. Numerical simulation of synergistic effect of neoclassical tearing mode and toroidal field ripple on alpha particle loss in China Fusion Engineering Testing Reactor. Acta Physica Sinica, 2021, 70(11): 115201. doi: 10.7498/aps.70.20201972
    [6] 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
    [7] Zhang Tai-Yang, Chen Ran. A collisional-radiative model for lithium impurity in plasma boundary region of Experimental Advanced Superconducting Tokamak. Acta Physica Sinica, 2017, 66(12): 125201. doi: 10.7498/aps.66.125201
    [8] Huang Yan, Sun Ji-Zhong, Sang Chao-Feng, Ding Fang, Wang De-Zhen. Numerical study of the erosion of the EAST tungsten divertor targets caused by edge localized modes. Acta Physica Sinica, 2014, 63(3): 035204. doi: 10.7498/aps.63.035204
    [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] 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
    [19] 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
    [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:  882
  • PDF Downloads:  29
  • Cited By: 0
Publishing process
  • Received Date:  23 May 2024
  • Accepted Date:  19 July 2024
  • Available Online:  15 August 2024
  • Published Online:  20 September 2024

/

返回文章
返回