搜索

x

留言板

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

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

钨杂质辐射对托卡马克等离子体大破裂快速热猝灭阶段热能损失过程的影响

张启凡 乐文成 张羽昊 葛忠昕 邝志强 萧声扬 王璐

引用本文:
Citation:

钨杂质辐射对托卡马克等离子体大破裂快速热猝灭阶段热能损失过程的影响

张启凡, 乐文成, 张羽昊, 葛忠昕, 邝志强, 萧声扬, 王璐
cstr: 32037.14.aps.73.20240730

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
导出引用
  • 基于 PLT, EAST, WEST, ASDEX-Upgrade, JET等托卡马克装置开展的研究表明重杂质易产生聚芯现象, 这会导致等离子体约束性能降低甚至引发等离子体大破裂事件. 大破裂期间等离子体热能损失主要发生在快速热猝灭(thermal quench, TQ)阶段, 但目前对于这一阶段的时间尺度定标关系并没有较为全面的物理解释. 国际热核聚变实验堆(International Thermonuclear Experimental Reactor, ITER)将采用全钨壁材料, 而钨作为高Z杂质, 其较强的辐射能力将会对TQ过程的热能损失产生影响. 为研究钨杂质对快速TQ时间尺度的影响, 本工作通过同时考虑随机磁场导致的热扩散以及钨杂质辐射引起的热损失机制, 建立了托卡马克等离子体电子温度演化的一维模型, 并在典型类ITER参数下对该阶段的电子温度演化进行数值计算和分析. 主要结论为: 1)快速TQ时间尺度的量级由热扩散水平决定, 但钨杂质辐射可以定量上影响TQ时间尺度和TQ后期电子温度, 钨浓度越高TQ时间尺度越短、后期电子温度越低, 且数值与解析结果分析都表明该时间尺度与钨杂质浓度近似呈线性关系; 2)快速TQ阶段前期, 通过钨杂质辐射损失的能量远小于通过随机磁场引起热扩散损失的能量, 但在TQ后期, 钨杂质辐射功率量级可以接近甚至超过热扩散功率, 这也是导致TQ后期电子温度随钨浓度增大而降低的原因. 因此, 钨杂质辐射在TQ后期对热能损失的贡献不可忽略.
    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.
      通信作者: 萧声扬, d202380863@hust.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 12275097, 51821005)资助的课题.
      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]

    Fig. 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参数下, 电子密度剖面(黑色方形)、初始电子温度剖面(蓝色圆圈)和安全因子剖面(红色三角)

    Fig. 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}} $

    Fig. 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}}}} $演化

    Fig. 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}} $对钨杂质浓度的依赖关系

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

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

    Fig. 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)全局总损失功率(黑色三角)、全局辐射损失功率(红色方形)和全局扩散损失功率(蓝色圆圈)随时间的演化

    Fig. 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
    下载: 导出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
    下载: 导出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] 朱霄龙, 陈伟, 王丰, 王正汹. 托卡马克中低频磁流体不稳定性协同作用引起快粒子输运的混合模拟研究. 物理学报, 2023, 72(21): 215210. doi: 10.7498/aps.72.20230620
    [2] 刘冠男, 李新霞, 刘洪波, 孙爱萍. HL-2M托卡马克装置中螺旋波与低杂波的协同电流驱动. 物理学报, 2023, 72(24): 245202. doi: 10.7498/aps.72.20231077
    [3] 沈勇, 董家齐, 何宏达, 潘卫, 郝广周. 托卡马克理想导体壁与磁流体不稳定性. 物理学报, 2023, 72(3): 035203. doi: 10.7498/aps.72.20222043
    [4] 王福琼, 徐颖峰, 查学军, 钟方川. 托卡马克边界等离子体中钨杂质输运的多流体及动力学模拟. 物理学报, 2023, 72(21): 215213. doi: 10.7498/aps.72.20230991
    [5] 郝保龙, 陈伟, 李国强, 王晓静, 王兆亮, 吴斌, 臧庆, 揭银先, 林晓东, 高翔, CFETRTEAM. 中国聚变工程试验堆上新经典撕裂模和纵场波纹扰动叠加效应对alpha粒子损失影响的数值模拟. 物理学报, 2021, 70(11): 115201. doi: 10.7498/aps.70.20201972
    [6] 刘朝阳, 章扬忠, 谢涛, 刘阿娣, 周楚. 托卡马克无碰撞捕获电子模在时空表象中的群速度. 物理学报, 2021, 70(11): 115203. doi: 10.7498/aps.70.20202003
    [7] 章太阳, 陈冉. 东方超环(EAST)装置中等离子体边界锂杂质的碰撞-辐射模型. 物理学报, 2017, 66(12): 125201. doi: 10.7498/aps.66.125201
    [8] 黄艳, 孙继忠, 桑超峰, 丁芳, 王德真. 边界局域模对EAST钨偏滤器靶板腐蚀程度的数值模拟研究. 物理学报, 2014, 63(3): 035204. doi: 10.7498/aps.63.035204
    [9] 张重阳, 刘阿娣, 李弘, 陈志鹏, 李斌, 杨州军, 周楚, 谢锦林, 兰涛, 刘万东, 庄革, 俞昌旋. 双极化频率调制微波反射计在J-TEXT托卡马克上的应用. 物理学报, 2014, 63(12): 125204. doi: 10.7498/aps.63.125204
    [10] 杜海龙, 桑超峰, 王亮, 孙继忠, 刘少承, 汪惠乾, 张凌, 郭后扬, 王德真. 东方超环托卡马克高约束模式边界等离子体输运数值模拟研究. 物理学报, 2013, 62(24): 245206. doi: 10.7498/aps.62.245206
    [11] 卢洪伟, 查学军, 胡立群, 林士耀, 周瑞杰, 罗家融, 钟方川. HT-7托卡马克slide-away放电充气对等离子体行为的影响. 物理学报, 2012, 61(7): 075202. doi: 10.7498/aps.61.075202
    [12] 洪斌斌, 陈少永, 唐昌建, 张新军, 胡有俊. 托卡马克中电子回旋波与低杂波协同驱动的物理研究. 物理学报, 2012, 61(11): 115207. doi: 10.7498/aps.61.115207
    [13] 卢洪伟, 胡立群, 林士耀, 钟国强, 周瑞杰, 张继宗. HT-7托卡马克等离子体slide-away放电研究. 物理学报, 2010, 59(8): 5596-5601. doi: 10.7498/aps.59.5596
    [14] 徐强, 高翔, 单家方, 胡立群, 赵君煜. HT-7托卡马克大功率低混杂波电流驱动的实验研究. 物理学报, 2009, 58(12): 8448-8453. doi: 10.7498/aps.58.8448
    [15] 龚学余, 彭晓炜, 谢安平, 刘文艳. 托卡马克等离子体不同运行模式下的电子回旋波电流驱动. 物理学报, 2006, 55(3): 1307-1314. doi: 10.7498/aps.55.1307
    [16] 徐 伟, 万宝年, 谢纪康. HT-6M托卡马克装置杂质输运. 物理学报, 2003, 52(8): 1970-1978. doi: 10.7498/aps.52.1970
    [17] 王文浩, 俞昌旋, 许宇鸿, 闻一之, 凌必利, 宋梅, 万宝年. HT-7超导托卡马克边界等离子体参量及其涨落的实验研究. 物理学报, 2001, 50(8): 1521-1527. doi: 10.7498/aps.50.1521
    [18] 王文浩, 许宇鸿, 俞昌旋, 闻一之, 凌必利, 宋梅, 万宝年. HT-7超导托卡马克边缘涨落谱特征及湍流输运研究. 物理学报, 2001, 50(10): 1956-1963. doi: 10.7498/aps.50.1956
    [19] 张先梅, 万宝年, 阮怀林, 吴振伟. HT-7托卡马克等离子体欧姆放电时电子热扩散系数的研究. 物理学报, 2001, 50(4): 715-720. doi: 10.7498/aps.50.715
    [20] 石秉仁. 托卡马克低混杂波电流驱动实验中低混杂波传播的解析分析. 物理学报, 2000, 49(12): 2394-2398. doi: 10.7498/aps.49.2394
计量
  • 文章访问数:  592
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-23
  • 修回日期:  2024-07-19
  • 上网日期:  2024-08-15
  • 刊出日期:  2024-09-20

/

返回文章
返回