搜索

x

留言板

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

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

中子辐照诱导钨再结晶的模拟研究

张国帅 尹超 王兆繁 陈泽 毛世峰 叶民友

张国帅, 尹超, 王兆繁, 陈泽, 毛世峰, 叶民友. 中子辐照诱导钨再结晶的模拟研究. 物理学报, 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
引用本文: 张国帅, 尹超, 王兆繁, 陈泽, 毛世峰, 叶民友. 中子辐照诱导钨再结晶的模拟研究. 物理学报, 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
Zhang Guo-Shuai, Yin Chao, Wang Zhao-Fan, Chen Ze, Mao Shi-Feng, Ye Min-You. Simulation of neutron irradiation-induced recrystallization of tungsten. Acta Phys. Sin., 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531
Citation: Zhang Guo-Shuai, Yin Chao, Wang Zhao-Fan, Chen Ze, Mao Shi-Feng, Ye Min-You. Simulation of neutron irradiation-induced recrystallization of tungsten. Acta Phys. Sin., 2023, 72(16): 162801. doi: 10.7498/aps.72.20230531

中子辐照诱导钨再结晶的模拟研究

张国帅, 尹超, 王兆繁, 陈泽, 毛世峰, 叶民友

Simulation of neutron irradiation-induced recrystallization of tungsten

Zhang Guo-Shuai, Yin Chao, Wang Zhao-Fan, Chen Ze, Mao Shi-Feng, Ye Min-You
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 钨作为未来聚变堆偏滤器靶板的候选材料, 需要长期服役在高温且受到高能中子辐照的严峻环境, 这将导致钨发生中子辐照诱导再结晶, 从而提高钨发生沿晶脆断的可能性, 威胁偏滤器的运行安全, 因此研究中子辐照诱导钨再结晶的物理机制具有重要意义. 然而, 与最近高通量同位素反应(HFIR)堆高温下中子辐照实验观察到的结果相比, 目前考虑辐照增强再结晶驱动力效应的模型低估了中子辐照对再结晶的影响, 结果表明仍有其他效应影响再结晶过程. 基于此, 本文在假设晶界迁移率与自体扩散系数成正比的前提下, 引入辐照增强晶界迁移因子(R), 建立了新的辐照诱导再结晶动力学模型. 模拟结果显示, 在综合考虑辐照增强再结晶驱动力和晶界迁移效应后, 模型计算出的850 ℃下达到一半再结晶分数所需要的时间(tX=0.5)和HFIR堆中子辐照实验结果相符, 这表明辐照增强晶界迁移效应是影响辐照诱导再结晶现象的重要因素之一. 另外, 模型研究了不同辐照温度下钨的tX=0.5. 结果表明辐照与未辐照的tX=0.5差别随温度升高而逐渐下降. 这是因为随着温度的升高, 辐照缺陷复合加剧, 辐照缺陷对再结晶驱动力的贡献下降, 且热激活扩散系数增大的幅度大于辐照下扩散系数的增大幅度, 所以热激活效应会逐渐主导再结晶过程.
    Tungsten is the candidate for divertor target material in future fusion reactors. The tungsten divertor target is expected to long serve in a harsh environment of high temperature and high-energy neutron irradiation. This can lead to neutron irradiation-induced recrystallization of tungsten, thereby increasing the possibility of intergranular brittle failure and compromising the safe operation of the divertor. Thus, clarifying the mechanism of neutron irradiation-induced tungsten recrystallization is important. However, the current model, which only considers the irradiation-enhanced effect on recrystallization driving force, underestimates the irradiation effect on recrystallization compared with the results observed in recent high-temperature neutron irradiation experiments in the HFIR reactor. It indicates that other irradiation effects can also influence the recrystallization process.In this study, we introduce the irradiation-enhanced grain boundary migration factor (R) into the established irradiation-induced recrystallization kinetic model, on the assumption that the grain boundary migration velocity is proportional to the self-diffusion coefficient. The simulation results show that after considering both irradiation-enhanced recrystallization driving force and grain boundary migration effect, the calculated half-recrystallization time (tX=0.5) at 850 ℃ from the model matches the one obtained in the neutron irradiation experiment in the HFIR reactor. This result indicates that the irradiation-enhanced grain boundary migration effect is one of the important factors affecting irradiation-induced recrystallization. In addition, the difference between irradiated and unirradiated tX=0.5 decreases with temperature increasing. This phenomenon is due to the fact that as the temperature increases, the contribution of irradiation defects to the driving force for recrystallization decreases owing to the irradiation defect recombination. Moreover, the increase of thermal activation diffusion coefficient is more significant than the increase of the irradiation-enhanced diffusion coefficient. These findings suggest that the thermal activation effect eventually dominates the recrystallization process over the irradiation effect as temperature increases.
      通信作者: 尹超, chaoyin@ustc.edu.cn ; 叶民友, yemy@ustc.edu.cn
    • 基金项目: 国家自然科学基金联合基金(批准号: U2267208)、中国博士后科学基金(批准号: 2021M703113)、中国科学院台湾青年人才计划(批准号: 2021TWGB0001)、安徽省高校协同创新项目(批准号: GXXT-2021-026)、中国科学院合肥大科学中心协同创新培育基金项目(批准号: 2022HSC-CIP010)和中央高校基本科研业务费专项资金(批准号: WK2140000015)资助的课题.
      Corresponding author: Yin Chao, chaoyin@ustc.edu.cn ; Ye Min-You, yemy@ustc.edu.cn
    • Funds: Project supported by the Joint Funds of the National Natural Science Foundation of China (Grant No. U2267208), the China Postdoctoral Science Foundation (Grant No. 2021M703113), the Chinese Academy of Sciences Taiwan Young Talent Program, China (Grant No. 2021TWGB0001), the University Synergy Innovation Program of Anhui Province, China (Grant No. GXXT-2021-026), the Collaborative Innovation Program of Hefei Science Center, Chinese Academy of Sciences, China (Grant No. 2022HSC-CIP010), and the Fundamental Research Funds for the Central Universities of China (Grant No. WK2140000015).

    钨由于具有高熔点、高热导、高溅射阈值、良好的抗热冲击性能及较低的氘氚滞留等优点, 被确定为国际热核聚变实验堆(ITER)的偏滤器靶板材料, 并且是聚变示范堆(DEMO)面向等离子体部件(plasma-facing component, PFC)中铠甲和靶板的候选材料[1-4]. 然而, 聚变堆严苛的服役环境, 会造成钨作为PFC服役性能的退化, 进而减短部件寿命[5]. 钨在聚变装置中的工作温度涵盖了500—1200 ℃的范围[6]. 研究表明塑性变形加工后的钨在高温下会发生退火现象, 包括回复、再结晶和晶粒生长过程[7-10]. 当钨发生再结晶时, 其转脆温度将从200 ℃上升至500 ℃, 从而在较高的工作温度下可能发生脆性断裂[11]. 因此, 分析并预测钨的再结晶过程对钨作为PFC的应用和评估具有重要意义. 根据未辐照等温退火实验外推预测, 钨在1000 ℃下服役两年后其再结晶分数约为50%[8,12]. 然而, 最近高通量同位素反应(high flux isotope reactor, HFIR)堆裂变中子辐照钨的实验显示, 钨在850 ℃下辐照24 d后的再结晶分数就已达到50%[13], 这表明中子辐照会诱导钨发生再结晶, 增加脆断风险.

    中子辐照会对钨造成辐照损伤, 由于缺陷之间的相互反应, 这些初级损伤产生的缺陷会进一步演化成尺寸更大的缺陷, 包括位错环和空洞[14-17]. 这些非热平衡的缺陷会增大钨内部储能[18], 从而提高了再结晶驱动力. Mannheim等[19,20]通过耦合团簇动力学(cluster dynamic, CD)模型和包含辐照增强再结晶驱动力效应的平均场再结晶模型, 从物理的角度模拟了再结晶过程中晶粒的成核和生长. 该模型预测在1000 ℃辐照温度下辐照近700 h后, 钨几乎不发生再结晶现象, 再结晶速率远低于HFIR堆实验结果, 这可能是该模型未考虑辐照增强晶界迁移效应所致.

    基于此, 本文在Johnson-Mehl-Avrami-Kolmogorov (JMAK)再结晶动力学模型的基础上, 同时考虑辐照增强再结晶驱动力效应和辐照增强晶界迁移效应, 建立了辐照诱导再结晶(irradiation-induced recrystallization, IIR)动力学模型.

    利用CD模型模拟钨在HFIR堆中子辐照下不同尺寸空位(vacancy, V)和自间隙原子(self-interstitial atom, I)团簇随时间的演化过程, 同时将模拟得到的缺陷尺寸与密度输入IIR模型以模拟高温中子辐照过程中钨再结晶分数的变化. 模拟的辐照温度区间为750—1300 ℃, 使用Matlab刚性求解器ODE15s进行数值求解, 相对容差设置为1×10–3, 并且设置每个解为非负值.

    如(1)式所示, 平均场CD模型基于扩散捕获速率理论, 通过一组常微分控制方程(ordinary differential equations, ODEs)来描述不同尺寸缺陷的密度随时间的演化[21]:

    dCθdt=Gθ+θJ+θ,θθJθ,θLθ,
    (1)

    其中Cθ为缺陷密度, 下标θ表示n个V组成的团簇(Vn)或n个I组成的团簇(In); θ为可移动缺陷, t为时间, Gθ为源项, J+θ,θ为吸收反应通量, 用于描述缺陷团簇吸收可移动缺陷的反应; Jθ,θ为发射反应通量, 用于描述缺陷团簇发射可移动缺陷的反应; Lθ为汇项, 用于描述缺陷阱对可移动缺陷的吸收. 除源项外, 吸收和发射反应项、汇项都是直接使用Mannheim等[19]建立和整理的反应式, 相应的参数在见补充材料中给出.

    为了获得源项, 首先通过SPECTRA-PKA程序计算特定中子能谱辐照钨所产生的初级反冲原子(primary knock-on atom, PKA)能谱[22,23], 再将PKA能谱叠加到PKA级联模拟结果[24]以确定级联存活损伤速率(Sdpa), 最后将Sdpa转换成单位体积内的级联存活损伤速率(Gtot):

    Gtot=SdpaρWNAMW,
    (2)

    其中ρW为钨的密度, NA为阿伏伽德罗常数, MW为钨的摩尔质量, 具体数值见补充材料. 另外, SPECTRA-PKA也给出了按照Norgett-Robinson-Torrens (NRT)模型得到的平均每离位原子速率(NRTdpa).

    表1为基于HFIR中子能谱[25]得到的NRTdpa, SdpaGtot. 这与文献[16]中根据NRT模型计算得出的平均每离位原子速率2.3×10–7 dpa/s基本吻合. 同时计算得到的级联存活效率(SdpaNRTdpa的比值)为0.30, 这与文献[26, 27]中MD模拟级联效率的结果一致. 由于级联内缺陷的相互反应, 使得级联内缺陷密度与尺寸呈幂律分布. 按照Vrielink等[28]的方法, 级联尺寸分布由下式给出:

    表 1  CD模型模拟HFIR堆中子辐照钨的源项相关参数
    Table 1.  Parameters related to source term of neutron irradiated tungsten in HFIR reactor simulated by CD model.
    参数NRTdpa/(10–7 dpa·s–1)Sdpa/(10–8 dpa·s–1)Gtot/(1021 m–3·s–1)
    数值2.166.414.06
    下载: 导出CSV 
    | 显示表格
    Gθn={(1fcθ)Gtot,n = 1,Aθ(1n1Nc)Sθ, 2nNc,0, otherwise,
    (3)

    其中fcθSθ分别为温度相关的缺陷团聚分数和幂律指数[19]. 钨低温自离子辐照实验结果表明级联产生的最大In尺寸约为1000[29], 并且MD级联模拟结果表明级联产生的最大Vn尺寸约为最大In尺寸的2倍[30]. 因此本文级联产生的最大Vn尺寸(Nc, V)取2000, 最大In尺寸(Nc, I)取1000.

    本文由于考虑的最大缺陷尺寸大于105, 为了提高计算效率, 在考虑大尺寸缺陷时, 采用Ghoniem和Sharafat[31]提出的Fokker-Plank近似方法. 此外, 测试的算例结果表明大尺寸缺陷(n > 100)的演化忽略I2反应项对缺陷的平均尺寸与密度影响不大, 这与Li等[32]的CD模拟结果一致. 因此本文中, 当缺陷尺寸大于100时, 通过二阶泰勒展开将(1)式转化成尺寸空间连续的形式[31]:

    C[x]t=G[x]x[FθC[x]x(DθC[x])]L[x],
    (4)

    其中, C[x]表示取x整数部分对应尺寸的缺陷密度, Fθ表示缺陷的净通量, Dθ表示缺陷的平均扩散通量, FθDθ分别可以写成[31]:

    Fθ=kθλθ,   Dθ=12(kθ+λθ),
    (5)

    其中kθλθ分别为相应缺陷团簇的生长速率和衰减速率[31].

    在尺寸空间x内按照Li等[32]的方法, 采取非均匀化的网格, 表达式如下:

    x1=1,    xi=xi1+Δxi,  2iNm,Δxi={1,2i100,Δxi1exp(ε),100iNm.
    (6)

    其中ε是远小于1的常数, Nm为考虑的最大缺陷尺寸, 且Nm需要满足边界条件CNm0或小于本研究考虑的缺陷的最大尺寸. 通过Fokker-Plank近似, 本文建立的CD模型可以考虑更大尺寸缺陷的演化, 其中最大的Vn尺寸约为3×105, 对应的空洞直径约20 nm, 最大的In尺寸约为6 × 103, 对应的位错环直径约20 nm.

    为了简化分析, 做出如下假设.

    1)认为整个空间缺陷密度分布是均匀的, 且缺陷密度足够低, 不影响基体钨的性质.

    2)认为只有V, I和I2是可移动的, 其他固定缺陷通过捕获可移动缺陷演化.

    3)所有Vn都是球体: 对于In, 小于7的是球体; 大于7的是位错环.

    4)无论是吸收反应还是反射反应, 只考虑吸收或发射一个可移动性缺陷, 即反应逐次进行.

    5)汇项只考虑位错对缺陷的吸收, 且位错不会因吸收缺陷而消失或改变密度.

    IIR模型通过未辐照下JMAK模型[9,33-35], 建立了中子辐照诱导再结晶模型, 假设中子辐照将会影响钨未发生再结晶的初始晶粒内再结晶驱动力P和晶界迁移率M. 温轧中等变形量钨的等温退火实验表明其再结晶过程符合位点饱和情况下的再结晶[7], 同时中子辐照产生的缺陷会提高钨自体扩散系数及内部储能, 应会加速新晶粒的成核过程, 因此本研究假设辐照下钨再结晶过程均为位点饱和情况. 对于位点饱和成核的JMAK模型, 再结晶分数(X)随时间的变化可以写成下列形式[33]:

    X=1exp[U(PMt)na],
    (7)

    其中U为成核数相关的材料常数; na为Avrami指数, 反映了晶粒的生长维度; 而中子辐照时的P可以写为

    P=EVold+EVolV+EVolI,
    (8)

    其中EVold表示位错造成的储能[33]:

    EVold=αρGb2,
    (9)

    其中α为常数, 取0.5; ρ表示位错线密度, G表示剪切模量, b是位错的伯格斯矢量; 具体数值见补充材料. 由于本研究假定初始晶粒的位错密度(7×1013 m–2[13])在再结晶过程中不变, 因此EVold为常数; EVolV表示空洞造成的储能; EVolI表示I位错环造成的储能. 钨的自离子辐照实验表明1/2111型I位错环占主导[36,37], 且1/2111型I位错环占比还会随温度的升高而增大[38]. 因此本研究中, 假设I位错环全部都是1/2111型. 根据Ma等[18]的模拟结果, 尺寸为n的空洞形成能(Ev(n))和尺寸为n的1/2111; 型I位错环形成能(El-1/2111(n))可以写成下列形式:

    Ev(n)=4πr2ˉγ=4π(3nΩ4π)2/3ˉγ,
    (10)
    El-1/2111(n)=2aE111nlnn+ε1n,
    (11)

    其中, a为晶格常数, r为空洞半径, n表示缺陷团簇包含的单体数, ˉγ表示空洞平均表面能密度, E111表示位错环平均线能量前对数因子, ε1取9.08 eV/[18](1 Å = 10–10 m), 具体数值见补充材料. 再将CD模型得到的缺陷密度信息结合(10)式和(11)式及小缺陷的形成能, 就可以求得任意时刻空位缺陷、自间隙缺陷的体积储能:

    EVolV=CVEfV+CV2EfV2+Nvmn=3CVnEv(n)TSv,
    (12)
    EVolI=6n=1CInEfIn+Nim7CInEl1/2111(n)TSi,
    (13)

    其中T为温度, SvSi分别是V缺陷和I缺陷的构型熵, 可利用Mannheim等[19]的方法得到.

    由于晶界迁移本质上是晶粒间原子的迁移, 因此假设晶界迁移率和自体扩散系数成正比. 在此定义辐照增强晶界迁移因子(R)为

    R=Drad/DradDtheDthe,
    (14)

    式中, R反映中子辐照增强晶界迁移效应, 其中DtheDrad分别是未辐照下热激活的自体扩散系数和中子辐照下的自体扩散系数. DtheDrad的表达式为[13,39]

    Dthe=DICeqI+DVCeqV,
    (15)
    Drad=DICI+DI2CI2+DVCV,
    (16)

    其中DI, DI2, DV分别是I, I2, V的扩散系数, 由(17)式给出. CeqICeqV分别是I和V的热平衡密度, 由(18)式给出. CI, CI2CV分别是辐照下I, I2和V的密度. D0θ, EmθEfθ分别为材料相关的扩散系数前置因子、可移动缺陷迁移能和可移动缺陷形成能, 具体数值见补充材料. Ω=a3/2为体心立方-钨的原子体积:

    Dθ=D0θexp(EmθkT),
    (17)
    Ceqθ=exp(EfθkT)/Ω.
    (18)

    未辐照钨等温退火实验数据拟合结果表明钨的晶界迁移率M满足阿伦尼乌斯的形式[7]. 通过引入辐照增强扩散因子, 中子辐照下晶界迁移率可以写成下列形式[40,41]:

    M=RM0Texp(QmkT),
    (19)

    其中M0为材料常数, Qm为再结晶激活能, T为温度, k为玻尔兹曼常数. 对(7)式进行微分得到辐照下再结晶分数随时间演化的微分表达式:

    dXdt=naURM0Texp[Qm/(kT)]naP(1X)×[ln(1X)](11/na).
    (20)

    为得到钨再结晶过程中的材料常数并与HFIR堆中子辐照钨(80%热轧变形量, 记为W80)实验结果进行比较, 通过对比文献[7, 8, 13]中W67(67%温轧变形量, 记为W67)和W80的显微组织和维氏硬度, 发现二者初始微观结构和力学性质相近, 因此这里可认为二者具有相似的再结晶动力学过程. 图1为使用JMAK模型拟合未辐照W67等温退火实验数据得到的结果[7]. 从图1可以看出, 除了再结晶的初始阶段, 拟合结果和实验结果基本符合. 拟合得到的材料常数如下:

    图 1 使用JMAK模型拟合钨等温退火实验再结晶分数的演变, 拟合用实验数据取自Lopez[7]\r\nFig. 1. Using the JMAK model to fit the evolution of the recrystallization fraction in the isothermal annealing experiment of pure tungsten, the experimental data used for the fitting were taken from Lopez[7].
    图 1  使用JMAK模型拟合钨等温退火实验再结晶分数的演变, 拟合用实验数据取自Lopez[7]
    Fig. 1.  Using the JMAK model to fit the evolution of the recrystallization fraction in the isothermal annealing experiment of pure tungsten, the experimental data used for the fitting were taken from Lopez[7].
    naUM0=5.02×1015Km3/J,
    EVold=6.04×105J/m3, Qm=6.08eV, na2.55.

    CD模型以HFIR堆中子辐照为例, 模拟了纯钨内部缺陷的演化. 图2为不同温度下缺陷演化的结果, 表明Vn和In的尺寸和密度对温度和时间存在不同的依赖关系. 从图2(a), (c), (e), (g), (i), (k)可以看出随着温度的上升, Vn的最大尺寸起初也随之增大, 在1100 ℃左右达到峰值, 之后最大尺寸急剧下降, 但其整体密度一直是呈降低的趋势. Vn随时间的演变与温度密切相关, 在低于1100 ℃时, Vn随时间一直生长; 高于1100 ℃, Vn很难生长起来. 这主要是由于V的迁移能(1.66 eV)较大的缘故. 根据扩散捕获速率理论, Vn与V发生吸收和发射反应的速率系数如下[32]:

    图 2 不同中子辐照温度下的缺陷团簇尺寸与密度随辐照时间的演变 (a), (c), (e), (g), (i), (k)分别为V团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变; (b), (d), (f), (h), (j), (l) 分别为I团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变\r\nFig. 2. Evolution of defect cluster size and density with irradiation time at different neutron irradiation temperatures: (a), (c), (e), (g), (i), (k) Evolution of V cluster at 750, 850, 950, 1100, 1200, 1300 ℃; (b), (d), (f), (h), (j), (l) I clusters evolution of cluster at 750, 850, 950, 1100, 1200, 1300 ℃.
    图 2  不同中子辐照温度下的缺陷团簇尺寸与密度随辐照时间的演变 (a), (c), (e), (g), (i), (k)分别为V团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变; (b), (d), (f), (h), (j), (l) 分别为I团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变
    Fig. 2.  Evolution of defect cluster size and density with irradiation time at different neutron irradiation temperatures: (a), (c), (e), (g), (i), (k) Evolution of V cluster at 750, 850, 950, 1100, 1200, 1300 ℃; (b), (d), (f), (h), (j), (l) I clusters evolution of cluster at 750, 850, 950, 1100, 1200, 1300 ℃.
    γ+n=4πrVnDV,
    (21)
    γn+1=γ+nΩexp(EbVn-VkT),
    (22)

    其中rVn为Vn的捕获半径; EbVn-V为Vn与V反应的结合能, 具体数值见补充材料.

    图3为由(21)式和(22)式得到的不同温度下Vn与V反应的速率系数, 而反应通量J±Vn, V正比于相应的速率系数. 从图3可知, 温度高于800 ℃, 吸收反应的速率系数γ+Vn就已经很大. 在750—1100 ℃之间, 随着温度的升高, Vn与V反应的吸收通量增大, 且吸收通量大于相应的发射反应通量, 导致Vn的生长速率与温度呈正相关. 由(22)式可知, 随着温度的升高, 发射反应的速率系数γVnγ+Vn增大得更快, 当温度超过1100 ℃后, Vn的发射反应通量占主导, 导致其平均尺寸和密度随温度的升高而急剧下降.

    图 3 不同温度下Vn与V发生反应的速率系数 (a)不同温度下Vn吸收V反应的速率系数; (b)不同温度下Vn+1发射V反应的速率系数\r\nFig. 3. Rate coefficients of the reaction between Vn and V at different temperatures: (a) Rate coefficients of Vn absorption V reactions at different temperatures; (b) rate coefficients of Vn+1 emission V reactions at different temperatures.
    图 3  不同温度下Vn与V发生反应的速率系数 (a)不同温度下Vn吸收V反应的速率系数; (b)不同温度下Vn+1发射V反应的速率系数
    Fig. 3.  Rate coefficients of the reaction between Vn and V at different temperatures: (a) Rate coefficients of Vn absorption V reactions at different temperatures; (b) rate coefficients of Vn+1 emission V reactions at different temperatures.

    有别于Vn, 图2(b), (d), (f), (h), (j)表明In的平均尺寸与整体密度都是随温度的上升而下降, 且在整个温度区间内其尺寸分布很快达到平衡且尺寸较小. 这主要是由于I的迁移能(0.013 eV)远低于V的迁移能. 在所考虑的温度区间内, I的扩散系数足够大, 导致In的分布很快达到平衡. 相较于In与I反应的吸收反应通量, 相应的发射反应通量随温度增大得更快, 导致其平均尺寸和密度随温度的升高而降低.

    上述结果与Klimenkov等[42]在裂变中子辐照钨实验上观察到的结果相比, 实验观察到空洞与析出物在1100 ℃时达到最大尺寸, 且空洞密度也基本随温度上升而下降, 这与本模拟结果一致; 但当辐照温度为1200 ℃时, 空洞密度不降反升, 另外1100 ℃时的位错环尺寸大于1000 ℃和1200 ℃, 这有别于本研究的模拟结果. 这可能是因为本研究没有考虑嬗变元素对缺陷演化的影响, 已有的实验和模拟研究表明不同密度的Re和Os会改变空洞与位错环的密度[17,43,44]. 此外, Vn并非只形成三维空洞, 部分Vn也会形成V位错环[18,45]. 因此, 模拟结果和实验有所出入.

    图4为基于3.1节CD模拟结果所得到的辐照增强P和辐照增强M, 其中考虑的最长时间为完全再结晶或辐照两年的时间. 图4(a)为不同温度下P随时间的演变, 可看出P与辐照温度呈负相关. 这是由于随着温度的上升, I和V扩散系数增大, 一方面使缺陷间复合增加, 另一方面使大尺寸Vn越容易生长, 进而吸收更多的点缺陷使得缺陷团簇密度下降, 而P与缺陷团簇密度成正相关, 因此当温度上升, 缺陷密度下降, P也跟着降低. 此外, 1100 ℃下P随着辐照时间成指数增大, 这与图2(g)中Vn的演变有直接关系. 为了说明这一点, 图4(b)展示了在所截取的时间点处, 不同缺陷类型对P的贡献占比. 可以明显看出低于1100 ℃下P主要由空洞贡献, 高于1100 ℃下P主要位错线贡献. 这是由于当辐照温度在1100 ℃以下, Vn一直生长, 所形成的空洞对P的贡献占比随时间增大; 而当辐照温度升高至1200 ℃及以上, Vn的平均尺寸和密度开始急剧下降, 且In的平均尺寸和密度也会降低, 因此P变成几乎全部都由位错线贡献. 值得一提的是, Vn在1100 ℃下生长得最快, 尤其是在演化后期出现更大尺寸的缺陷团簇, 虽然密度有所下降, 但尺寸相对于In高出2个数量级, 因此P在1100 ℃下的演化后期急剧上升.

    图 4 (a)不同辐照温度下P的演变; (b)在截取的时间点处不同辐照温度下P的占比; (c)不同辐照温度下R的演变; (d)不同辐照温度下P×M的演变\r\nFig. 4. (a) Evolution of driving force P at different irradiation temperatures; (b) proportion of defects contribution to driving force at different irradiation temperatures and time; (c) evolution of R at different irradiation temperatures; (d) evolution of product of driving force and grain boundary mobility (P×M) at different irradiation temperatures.
    图 4  (a)不同辐照温度下P的演变; (b)在截取的时间点处不同辐照温度下P的占比; (c)不同辐照温度下R的演变; (d)不同辐照温度下P×M的演变
    Fig. 4.  (a) Evolution of driving force P at different irradiation temperatures; (b) proportion of defects contribution to driving force at different irradiation temperatures and time; (c) evolution of R at different irradiation temperatures; (d) evolution of product of driving force and grain boundary mobility (P×M) at different irradiation temperatures.

    图4(c)所示, 与P随辐照温度变化的趋势相同, 中子辐照下的R与温度一直成负相关. 这是因为低温时点缺陷热平衡密度远低于辐照产生的点缺陷密度, 使得Dthe远低于Drad, 进而温度越低R越大. 因此可认为辐照温度越低, M越容易受到辐照增强自体扩散效应的影响[13]. 但随着温度逐渐升高, 点缺陷热平衡密度上升, 导致Dthe快速增大. 模拟结果表明辐照下Drad也会随温度的上升而上升, 但其上升的幅度比Dthe上升的幅度小, 造成R逐渐降低. 且在辐照温度升至1300 ℃时, R更是接近于1, 因此1300 ℃下辐照的M接近于未辐照的M.

    图4(d)为不同辐照温度下P×M的演变, 反映了二者对中子辐照诱导再结晶过程的综合影响. 钨的晶界迁移由于是热激活过程((19)式), 其随温度的上升而指数上升, 因此P×M基本上是随温度的上升而增大. 但在1100 ℃下, 由于P在演化后期急剧上升, 因此出现P×M > 1200 ℃的现象.

    图5为只考虑辐照增强P效应以及辐照增强P效应和M效应都考虑时再结晶分数随时间的演变曲线. 可以看出, 若只考虑辐照增强P, 当辐照温度低于850 ℃时, 两年内都不会发生明显的再结晶现象, 且1200 ℃下需要约171 h达到一半再结晶分数. 相比之下, 若中子辐照增强PM效应都考虑下, 850 ℃时纯钨的tX=0.5只有约503 h, 而1200 ℃下只需约27 h就能达到一半再结晶分数, 这比只考虑辐照增强P条件下快了144 h. 根据Lopez[7]的未辐照下等温退火实验结果, 1200 ℃纯钨的tX=0.5约为300 h, 因此, 可以认为中子辐照缺陷导致PM增大将使钨在1200 ℃下的tX=0.5缩短约273 h.

    图 5 只考虑辐照增强P及同时考虑辐照增强P和M (P+M)的再结晶分数演变曲线\r\nFig. 5. Recrystallization fraction (X) evolution curve considering only irradiation enhancement on driving force (P) and both irradiation enhancement driving force and grain boundary mobility (P+M).
    图 5  只考虑辐照增强P及同时考虑辐照增强PM (P+M)的再结晶分数演变曲线
    Fig. 5.  Recrystallization fraction (X) evolution curve considering only irradiation enhancement on driving force (P) and both irradiation enhancement driving force and grain boundary mobility (P+M).

    当将辐照增强PM都考虑下850 ℃的模拟结果与HFIR堆850 ℃下中子辐照钨24 d (576 h)达到一半再结晶分数的结果[13]进行对比可以发现, 模拟得到的tX=0.5 (503 h)与实验基本符合但略小些. 实验的tX=0.5略长的原因初步推测是HFIR堆高热中子通量辐照下产生的嬗变元素(Re和Os含量分别达到3.3%和1.9%原子分数[13]). 这些嬗变元素会形成纳米析出颗粒, 从会钉扎住晶界, 产生Zener阻力, 降低晶界迁移速率[33,46]. 因此, HFIR堆850 ℃下的实验结果相较于模拟结果而言有较长的tX=0.5.

    同时辐照增强PM时, 虽然M随着温度指数增长, 但1100 ℃和1200 ℃的完全再结晶时间却几乎相同. 这是由于在1100 ℃辐照温度下的演化后期, PM的乘积大幅增大(图4(d)), 导致后期再结晶速率加快. 另一个值得关注的现象是在1300 ℃的辐照温度下, 只考虑辐照增强P和辐照增强PM都考虑下的再结晶曲线几乎重合. 这是由于高温下热激活扩散效应占主导, 使R趋近为1, 因此二者的再结晶速率差别不大.

    图6为不同温度下基于热激活JMAK模型得到的tX=0.5和IIR模型得到的tX=0.5, 可以明显看出温度越高辐照与未辐照的tX=0.5量级相差越小. 相较于未辐照下的热激活再结晶过程, 750 ℃中子辐照下的tX=0.5缩短了8个量级, 1000 ℃下的缩短了5个量级, 而1300 ℃下的则相差无几. 这是由于随着温度的升高, 热激活主导了晶界迁移, 且温度高于1100 ℃后辐照缺陷密度和尺寸大幅下降, 使得辐照诱导再结晶现象基本消失.

    图 6 中子辐照下与未辐照下钨的半再结晶时间(${t_{X = 0.5}}$)随温度的演变\r\nFig. 6. Evolution of semi-recrystallization time (${t_{X = 0.5}}$) of tungsten under neutron irradiation and non-irradiation with temperature.
    图 6  中子辐照下与未辐照下钨的半再结晶时间(tX=0.5)随温度的演变
    Fig. 6.  Evolution of semi-recrystallization time (tX=0.5) of tungsten under neutron irradiation and non-irradiation with temperature.

    本研究通过建立CD模型和IIR模型模拟了不同温度下中子辐照钨缺陷团簇和再结晶分数的演化, 其中模拟的辐照缺陷尺寸和时间尺度分别达到十几纳米和年的量级, 结果包括以下3点.

    1) CD模拟结果表明不同温度下I团簇和V团簇随时间演变的趋势不同, 这主要是由于二者的迁移性质不同所导致的. 其中最突出的特征表现为V团簇尺寸在辐照温度为1100 ℃左右时达到极值, 这与中子辐照实验观察到的结果一致. 综合辐照增强PM的效应, 导致1100 ℃下完全再结晶时间和1200 ℃下的几乎相同.

    2) 本文在辐照诱导驱动力的基础上, 通过引入辐照增强晶界迁移因子R建立了新的IIR模型. 模拟结果与最近HFIR堆850 ℃中子辐照实验相符, 表明辐照增强晶界迁移是影响辐照诱导再结晶过程的重要因素. 模拟计算的850 ℃下中子辐照的半再结晶时间(tX=0.5)虽然与实验结果几乎相同, 但模拟预测的tX=0.5稍短于实验结果, 这可能是由于模型忽略了嬗变析出物对缺陷演化和晶界迁移影响的缘故. 由于是首次实验上观察到钨的辐照诱导再结晶现象, 模型需要进一步的实验数据进行完善.

    3) 通过IIR模型计算了不同辐照温度下钨再结晶分数曲线, 给出了偏滤器可能运行温度区间内辐照诱导再结晶动力学的完整图像. 同时, 模型比较了辐照和未辐照下半再结晶时间. 结果表明辐照增强效应和热激活效应对温度不同的依赖关系, 低温下辐照增强效应占主导, 高温下则是由热激活效应占主导.

    另外, 由于目前的IIR模型还并不完备, 在未来将会纳入以下因素以进一步的完善模型.

    1) 纳入钨微观结构的变化和嬗变元素对缺陷演化及晶界迁移的影响, 包括位错和晶界的演变以及晶界迁移过程中, 晶界与缺陷的相互作用.

    2) 将再结晶开始前的回复阶段纳入模型. 未辐照纯钨等温退火实验1200 ℃下都观察到了明显的再结晶成核孕育阶段.

    3) 本研究假设发生的再结晶过程是位点饱和成核, 即在短时间内再结晶晶核就已经饱和, 这在高温或辐照情况下可能是如此, 但还需要进一步探究成核过程以更确切的评估中子辐照下钨再结晶动力学过程.

    4) 改进CD模型, 考虑大尺寸位错环也具有一定的移动性以及位错环与位错线的相互作用.

    [1]

    Philipps V 2011 J. Nucl. Mater. 415 S2Google Scholar

    [2]

    Rieth M, Dudarev S L, Gonzalez de Vicente S M, et al. 2013 J. Nucl. Mater. 432 482Google Scholar

    [3]

    Norajitra P, Abdel-Khalik S I, Giancarli L M, Ihli T, Janeschitz G, Malang S, Mazul I V, Sardain P 2008 Fusion Eng. Des. 83 893Google Scholar

    [4]

    Abernethy R G 2017 J. Mater. Sci. Technol. 33 388Google Scholar

    [5]

    Coenen J W, Antusch S, Aumann M, et al. 2016 Phys. Scr. T 2016 014002Google Scholar

    [6]

    Hu X, Koyanagi T, Fukuda M, Katoh Y, Snead L L, Wirth B D 2016 J. Nucl. Mater. 480 235Google Scholar

    [7]

    Lopez A A 2015 Ph. D. Dissertation (Copenhagen: Technical University of Denmark)

    [8]

    Alfonso A, Jensen D J, Luo G N, Pantleon W 2014 J. Nucl. Mater. 455 591Google Scholar

    [9]

    Kang W A, Dr A, Xiang Z, Lla B, Xz B, Ywab C 2021 Mater. Sci. Eng. A 806 140828Google Scholar

    [10]

    Budaev V P, Martynenko Y V, Karpov A V, Belova N E, Zhitlukhin A M 2015 J. Nucl. Mater. 463 237Google Scholar

    [11]

    Bonnekoh C, Reiser J, Hartmaier A, Bonk S, Hoffmann A, Rieth M 2020 J. Mater. Sci. 55 12314Google Scholar

    [12]

    Ciucani U M, Thum A, Devos C, Pantleon W 2019 Nucl. Mater. Energy 20 100701Google Scholar

    [13]

    Gietl H, Koyanagi T, Hu X, Fukuda M, Hasegawa A, Katoh Y 2022 J. Alloys Compd. 901 163419Google Scholar

    [14]

    Duerrschnabel M, Klimenkov M, Jaentsch U, Rieth M, Schneider H C, Terentyev D 2021 Sci. Rep. 11 7572Google Scholar

    [15]

    Klimenkov M, Jaentsch U, Rieth M, Schneider H C, Armstrong D E J, Gibson J, Roberts S G 2016 Nucl. Mater. Energy 9 480Google Scholar

    [16]

    Fukuda M, Kumar N A P K, Koyanagi T, Garrison L M, Snead L L, Katoh Y, Hasegawa A 2016 J. Nucl. Mater. 479 249Google Scholar

    [17]

    Fukuda M, Tanno T, Nogami S, Hasegawa A 2012 Mater. Trans. 53 2145Google Scholar

    [18]

    Ma P W, Mason D R, Dudarev S L 2020 Phys. Rev. Mater. 4 103609Google Scholar

    [19]

    Mannheim A, van Dommelen J A W, Geers M G D 2018 Mech. Mater. 123 43Google Scholar

    [20]

    Mannheim A, van Dommelen J A W, Geers M G D 2019 Comput. Mater. Sci. 170 109146Google Scholar

    [21]

    Barbu A, Clouet E 2007 Solid State Phenom. 129 51Google Scholar

    [22]

    Gilbert M R, Marian J, Sublet J C 2015 J. Nucl. Mater. 467 121Google Scholar

    [23]

    Gilbert M R, Sublet J C 2018 J. Nucl. Mater. 504 101Google Scholar

    [24]

    Setyawan W, Nandipati G, Roche K J, Heinisch H L, Wirth B D, Kurtz R J 2015 J. Nucl. Mater. 462 329Google Scholar

    [25]

    Huang C H, Gilbert M R, Marian J 2018 J. Nucl. Mater. 499 204Google Scholar

    [26]

    Troev T, Nankov N, Yoshiie T 2011 Nucl. Instrum. Methods Phys. Res. B 269 566Google Scholar

    [27]

    Caturla M J, Rubia T, Victoria M, Corzine R K, Greene G A 2001 J. Nucl. Mater. 296 90Google Scholar

    [28]

    Vrielink M A O, Shah V, van Dommelen J A W, Geers M G D 2021 J. Nucl. Mater. 554 153068Google Scholar

    [29]

    Yi X, Sand A E, Mason D R, Kirk M A, Roberts S G, Nordlund K, Dudarev S L 2015 Epl 110 36001Google Scholar

    [30]

    Sand A E, Mason D R, De Backer A, Yi X, Dudarev S L, Nordlund K 2017 Mater. Res. Lett. 5 357Google Scholar

    [31]

    Ghoniem N M, Sharafat S 1980 J. Nucl. Mater. 92 121Google Scholar

    [32]

    Li Y G, Zhou W H, Ning R H, Huang L F, Zeng Z, Ju X 2012 Commun. Comput. Phys. 11 1547Google Scholar

    [33]

    Humphreys F J, Hatherly M 2004 Recrystallization and Related Annealing Phenomena (Oxford: Elsevier) pp232–242

    [34]

    Fanfoni M, Tomellini M 1998 Il Nuovo Cimento D 20 1171Google Scholar

    [35]

    Hallberg H 2011 Metals 1 16Google Scholar

    [36]

    Yi X, Jenkins M L, Hattar K, Edmondson P D, Roberts S G 2015 Acta Mater. 92 163Google Scholar

    [37]

    Yi X, Jenkins M L, Kirk M A, Zhou Z, Roberts S G 2016 Acta Mater. 112 105Google Scholar

    [38]

    Yi X 2014 Ph. D. Dissertation (Oxford: University of Oxford) pp207–234

    [39]

    Was G S 2017 Fundamentals of Radiation Materials Science (Berlin: Springer) pp191—203

    [40]

    Rollett A D, Gottstein G, Shvindlerman L S, Molodov D A 2004 Zeitschrift Fur. Metallkunde 95 226Google Scholar

    [41]

    Favre J, Fabregue D, Piot D, Tang N, Koizumi Y, Maire E, Chiba A 2013 Metall. Mater. Trans. A 44 5861Google Scholar

    [42]

    Klimenkov M, Duerrschnabel M, Jaentsch U, Lied P, Rieth M, Schneider H C, Terentyev D, Van Renterghem W 2022 J. Nucl. Mater. 572 154018Google Scholar

    [43]

    Li Y H, Zhou H B, Jin S, Zhang Y, Deng H, Lu G H 2017 Nucl. Fusion 57 046006Google Scholar

    [44]

    You Y W, Kong X S, Wu X, Liu C S, Fang Q F, Chen J L, Luo G N 2017 Nucl. Fusion 57 086006Google Scholar

    [45]

    Setyawan W, Selby A P, Juslin N, Stoller R E, Wirth B D, Kurtz R J 2015 J. Phys. Condens. Matter 27 225402Google Scholar

    [46]

    Nes E, Ryum N, Hunderi O 1985 Acta Metall. 33 11Google Scholar

    期刊类型引用(1)

    1. 李旭,廖晨伦,谢亚红,曹小岗,洪兵,梁立振,顾玉明,胡纯栋,孟献才. 基于全永磁微波离子源钛自成靶的D-D中子源设计研究. 真空与低温. 2024(03): 336-342 . 百度学术

    其他类型引用(0)

  • 图 1  使用JMAK模型拟合钨等温退火实验再结晶分数的演变, 拟合用实验数据取自Lopez[7]

    Fig. 1.  Using the JMAK model to fit the evolution of the recrystallization fraction in the isothermal annealing experiment of pure tungsten, the experimental data used for the fitting were taken from Lopez[7].

    图 2  不同中子辐照温度下的缺陷团簇尺寸与密度随辐照时间的演变 (a), (c), (e), (g), (i), (k)分别为V团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变; (b), (d), (f), (h), (j), (l) 分别为I团簇在750, 850, 950, 1100, 1200, 1300 ℃下的演变

    Fig. 2.  Evolution of defect cluster size and density with irradiation time at different neutron irradiation temperatures: (a), (c), (e), (g), (i), (k) Evolution of V cluster at 750, 850, 950, 1100, 1200, 1300 ℃; (b), (d), (f), (h), (j), (l) I clusters evolution of cluster at 750, 850, 950, 1100, 1200, 1300 ℃.

    图 3  不同温度下Vn与V发生反应的速率系数 (a)不同温度下Vn吸收V反应的速率系数; (b)不同温度下Vn+1发射V反应的速率系数

    Fig. 3.  Rate coefficients of the reaction between Vn and V at different temperatures: (a) Rate coefficients of Vn absorption V reactions at different temperatures; (b) rate coefficients of Vn+1 emission V reactions at different temperatures.

    图 4  (a)不同辐照温度下P的演变; (b)在截取的时间点处不同辐照温度下P的占比; (c)不同辐照温度下R的演变; (d)不同辐照温度下P×M的演变

    Fig. 4.  (a) Evolution of driving force P at different irradiation temperatures; (b) proportion of defects contribution to driving force at different irradiation temperatures and time; (c) evolution of R at different irradiation temperatures; (d) evolution of product of driving force and grain boundary mobility (P×M) at different irradiation temperatures.

    图 5  只考虑辐照增强P及同时考虑辐照增强PM (P+M)的再结晶分数演变曲线

    Fig. 5.  Recrystallization fraction (X) evolution curve considering only irradiation enhancement on driving force (P) and both irradiation enhancement driving force and grain boundary mobility (P+M).

    图 6  中子辐照下与未辐照下钨的半再结晶时间(tX=0.5)随温度的演变

    Fig. 6.  Evolution of semi-recrystallization time (tX=0.5) of tungsten under neutron irradiation and non-irradiation with temperature.

    表 1  CD模型模拟HFIR堆中子辐照钨的源项相关参数

    Table 1.  Parameters related to source term of neutron irradiated tungsten in HFIR reactor simulated by CD model.

    参数NRTdpa/(10–7 dpa·s–1)Sdpa/(10–8 dpa·s–1)Gtot/(1021 m–3·s–1)
    数值2.166.414.06
    下载: 导出CSV
  • [1]

    Philipps V 2011 J. Nucl. Mater. 415 S2Google Scholar

    [2]

    Rieth M, Dudarev S L, Gonzalez de Vicente S M, et al. 2013 J. Nucl. Mater. 432 482Google Scholar

    [3]

    Norajitra P, Abdel-Khalik S I, Giancarli L M, Ihli T, Janeschitz G, Malang S, Mazul I V, Sardain P 2008 Fusion Eng. Des. 83 893Google Scholar

    [4]

    Abernethy R G 2017 J. Mater. Sci. Technol. 33 388Google Scholar

    [5]

    Coenen J W, Antusch S, Aumann M, et al. 2016 Phys. Scr. T 2016 014002Google Scholar

    [6]

    Hu X, Koyanagi T, Fukuda M, Katoh Y, Snead L L, Wirth B D 2016 J. Nucl. Mater. 480 235Google Scholar

    [7]

    Lopez A A 2015 Ph. D. Dissertation (Copenhagen: Technical University of Denmark)

    [8]

    Alfonso A, Jensen D J, Luo G N, Pantleon W 2014 J. Nucl. Mater. 455 591Google Scholar

    [9]

    Kang W A, Dr A, Xiang Z, Lla B, Xz B, Ywab C 2021 Mater. Sci. Eng. A 806 140828Google Scholar

    [10]

    Budaev V P, Martynenko Y V, Karpov A V, Belova N E, Zhitlukhin A M 2015 J. Nucl. Mater. 463 237Google Scholar

    [11]

    Bonnekoh C, Reiser J, Hartmaier A, Bonk S, Hoffmann A, Rieth M 2020 J. Mater. Sci. 55 12314Google Scholar

    [12]

    Ciucani U M, Thum A, Devos C, Pantleon W 2019 Nucl. Mater. Energy 20 100701Google Scholar

    [13]

    Gietl H, Koyanagi T, Hu X, Fukuda M, Hasegawa A, Katoh Y 2022 J. Alloys Compd. 901 163419Google Scholar

    [14]

    Duerrschnabel M, Klimenkov M, Jaentsch U, Rieth M, Schneider H C, Terentyev D 2021 Sci. Rep. 11 7572Google Scholar

    [15]

    Klimenkov M, Jaentsch U, Rieth M, Schneider H C, Armstrong D E J, Gibson J, Roberts S G 2016 Nucl. Mater. Energy 9 480Google Scholar

    [16]

    Fukuda M, Kumar N A P K, Koyanagi T, Garrison L M, Snead L L, Katoh Y, Hasegawa A 2016 J. Nucl. Mater. 479 249Google Scholar

    [17]

    Fukuda M, Tanno T, Nogami S, Hasegawa A 2012 Mater. Trans. 53 2145Google Scholar

    [18]

    Ma P W, Mason D R, Dudarev S L 2020 Phys. Rev. Mater. 4 103609Google Scholar

    [19]

    Mannheim A, van Dommelen J A W, Geers M G D 2018 Mech. Mater. 123 43Google Scholar

    [20]

    Mannheim A, van Dommelen J A W, Geers M G D 2019 Comput. Mater. Sci. 170 109146Google Scholar

    [21]

    Barbu A, Clouet E 2007 Solid State Phenom. 129 51Google Scholar

    [22]

    Gilbert M R, Marian J, Sublet J C 2015 J. Nucl. Mater. 467 121Google Scholar

    [23]

    Gilbert M R, Sublet J C 2018 J. Nucl. Mater. 504 101Google Scholar

    [24]

    Setyawan W, Nandipati G, Roche K J, Heinisch H L, Wirth B D, Kurtz R J 2015 J. Nucl. Mater. 462 329Google Scholar

    [25]

    Huang C H, Gilbert M R, Marian J 2018 J. Nucl. Mater. 499 204Google Scholar

    [26]

    Troev T, Nankov N, Yoshiie T 2011 Nucl. Instrum. Methods Phys. Res. B 269 566Google Scholar

    [27]

    Caturla M J, Rubia T, Victoria M, Corzine R K, Greene G A 2001 J. Nucl. Mater. 296 90Google Scholar

    [28]

    Vrielink M A O, Shah V, van Dommelen J A W, Geers M G D 2021 J. Nucl. Mater. 554 153068Google Scholar

    [29]

    Yi X, Sand A E, Mason D R, Kirk M A, Roberts S G, Nordlund K, Dudarev S L 2015 Epl 110 36001Google Scholar

    [30]

    Sand A E, Mason D R, De Backer A, Yi X, Dudarev S L, Nordlund K 2017 Mater. Res. Lett. 5 357Google Scholar

    [31]

    Ghoniem N M, Sharafat S 1980 J. Nucl. Mater. 92 121Google Scholar

    [32]

    Li Y G, Zhou W H, Ning R H, Huang L F, Zeng Z, Ju X 2012 Commun. Comput. Phys. 11 1547Google Scholar

    [33]

    Humphreys F J, Hatherly M 2004 Recrystallization and Related Annealing Phenomena (Oxford: Elsevier) pp232–242

    [34]

    Fanfoni M, Tomellini M 1998 Il Nuovo Cimento D 20 1171Google Scholar

    [35]

    Hallberg H 2011 Metals 1 16Google Scholar

    [36]

    Yi X, Jenkins M L, Hattar K, Edmondson P D, Roberts S G 2015 Acta Mater. 92 163Google Scholar

    [37]

    Yi X, Jenkins M L, Kirk M A, Zhou Z, Roberts S G 2016 Acta Mater. 112 105Google Scholar

    [38]

    Yi X 2014 Ph. D. Dissertation (Oxford: University of Oxford) pp207–234

    [39]

    Was G S 2017 Fundamentals of Radiation Materials Science (Berlin: Springer) pp191—203

    [40]

    Rollett A D, Gottstein G, Shvindlerman L S, Molodov D A 2004 Zeitschrift Fur. Metallkunde 95 226Google Scholar

    [41]

    Favre J, Fabregue D, Piot D, Tang N, Koizumi Y, Maire E, Chiba A 2013 Metall. Mater. Trans. A 44 5861Google Scholar

    [42]

    Klimenkov M, Duerrschnabel M, Jaentsch U, Lied P, Rieth M, Schneider H C, Terentyev D, Van Renterghem W 2022 J. Nucl. Mater. 572 154018Google Scholar

    [43]

    Li Y H, Zhou H B, Jin S, Zhang Y, Deng H, Lu G H 2017 Nucl. Fusion 57 046006Google Scholar

    [44]

    You Y W, Kong X S, Wu X, Liu C S, Fang Q F, Chen J L, Luo G N 2017 Nucl. Fusion 57 086006Google Scholar

    [45]

    Setyawan W, Selby A P, Juslin N, Stoller R E, Wirth B D, Kurtz R J 2015 J. Phys. Condens. Matter 27 225402Google Scholar

    [46]

    Nes E, Ryum N, Hunderi O 1985 Acta Metall. 33 11Google Scholar

  • [1] 祁超, 马玉田, 齐艳飞, 肖善曲, 王波. 微观组织对叠片结构钨基面向等离子体材料的热疲劳效应的影响. 物理学报, 2024, 73(11): 112801. doi: 10.7498/aps.73.20240007
    [2] 彭超, 雷志锋, 张战刚, 何玉娟, 马腾, 蔡宗棋, 陈义强. 中子辐射导致的SiC功率器件漏电增加特性研究. 物理学报, 2023, 72(18): 186102. doi: 10.7498/aps.72.20230976
    [3] 徐驰, 万发荣. 聚变材料钨辐照后退火形成的位错环特性及inside-outside衬度分析. 物理学报, 2023, 72(5): 056801. doi: 10.7498/aps.72.20222124
    [4] 魏雯静, 高旭东, 吕亮亮, 许楠楠, 李公平. 中子对碲锌镉辐照损伤模拟研究. 物理学报, 2022, 71(22): 226102. doi: 10.7498/aps.71.20221195
    [5] 王凯, 孙靖雅, 潘昌基, 王飞飞, 张可, 陈治成. 飞秒激光辐照二硫化钨的超快动态响应及时域整形调制. 物理学报, 2021, 70(20): 205201. doi: 10.7498/aps.70.20210737
    [6] 黄文军, 乔珺威, 陈顺华, 王雪姣, 吴玉程. 含钨难熔高熵合金的制备、结构与性能. 物理学报, 2021, 70(10): 106201. doi: 10.7498/aps.70.20201986
    [7] 郝蕊静, 郭红霞, 潘霄宇, 吕玲, 雷志锋, 李波, 钟向丽, 欧阳晓平, 董世剑. AlGaN/GaN高电子迁移率晶体管器件中子位移损伤效应及机理. 物理学报, 2020, 69(20): 207301. doi: 10.7498/aps.69.20200714
    [8] 周良付, 张婧, 何文豪, 王栋, 苏雪, 杨冬燕, 李玉红. 氦泡在bcc钨中晶界处成核长大的分子动力学模拟. 物理学报, 2020, 69(4): 046103. doi: 10.7498/aps.69.20191069
    [9] 马玉田, 刘俊标, 韩立, 田利丰, 王学聪, 孟祥敏, 肖善曲, 王波. 氦离子显微镜对钨中氦行为的实验研究. 物理学报, 2019, 68(4): 040702. doi: 10.7498/aps.68.20181864
    [10] 袁伟, 彭海波, 杜鑫, 律鹏, 沈扬皓, 赵彦, 陈亮, 王铁山. 分子动力学模拟钠硼硅酸盐玻璃电子辐照诱导的结构演化效应. 物理学报, 2017, 66(10): 106102. doi: 10.7498/aps.66.106102
    [11] 孙凌涛, 郭朝中, 肖绪洋. Cu偏析诱导Co团簇结构及性质异常的动力学模拟. 物理学报, 2016, 65(12): 123601. doi: 10.7498/aps.65.123601
    [12] 郭洪燕, 夏敏, 燕青芝, 郭立平, 陈济红, 葛昌纯. 中能高浓度氦离子注入对钨微观结构的影响. 物理学报, 2016, 65(7): 077803. doi: 10.7498/aps.65.077803
    [13] 曾骏哲, 李豫东, 文林, 何承发, 郭旗, 汪波, 玛丽娅, 魏莹, 王海娇, 武大猷, 王帆, 周航. 质子与中子辐照对电荷耦合器件暗信号参数的影响及其效应分析. 物理学报, 2015, 64(19): 194208. doi: 10.7498/aps.64.194208
    [14] 王成龙, 王庆宇, 张跃, 李忠宇, 洪兵, 苏折, 董良. SiC/C界面辐照性能的分子动力学研究. 物理学报, 2014, 63(15): 153402. doi: 10.7498/aps.63.153402
    [15] 王欣欣, 张颖, 周洪波, 王金龙. 铌对钨中氦行为影响的第一性原理研究. 物理学报, 2014, 63(4): 046103. doi: 10.7498/aps.63.046103
    [16] 黄艳, 孙继忠, 桑超峰, 丁芳, 王德真. 边界局域模对EAST钨偏滤器靶板腐蚀程度的数值模拟研究. 物理学报, 2014, 63(3): 035204. doi: 10.7498/aps.63.035204
    [17] 谷文萍, 张林, 李清华, 邱彦章, 郝跃, 全思, 刘盼枝. 中子辐照对AlGaN/GaN高电子迁移率晶体管器件电特性的影响. 物理学报, 2014, 63(4): 047202. doi: 10.7498/aps.63.047202
    [18] 郭龙婷, 孙继忠, 黄艳, 刘升光, 王德真. 低能氢粒子沿不同角度轰击钨(001)表面的反射概率及入射深度分布的分子动力学研究. 物理学报, 2013, 62(22): 227901. doi: 10.7498/aps.62.227901
    [19] 冯锡淇, 林奇生, 满振勇, 廖晶莹, 胡关钦. 钨酸铅晶体的本征色心和辐照诱导色心. 物理学报, 2002, 51(2): 315-321. doi: 10.7498/aps.51.315
    [20] 李养贤, 刘何燕, 牛萍娟, 刘彩池, 徐岳生, 杨德仁, 阙端鳞. 中子辐照直拉硅中的本征吸除效应. 物理学报, 2002, 51(10): 2407-2410. doi: 10.7498/aps.51.2407
  • 期刊类型引用(1)

    1. 李旭,廖晨伦,谢亚红,曹小岗,洪兵,梁立振,顾玉明,胡纯栋,孟献才. 基于全永磁微波离子源钛自成靶的D-D中子源设计研究. 真空与低温. 2024(03): 336-342 . 百度学术

    其他类型引用(0)

计量
  • 文章访问数:  4769
  • PDF下载量:  132
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-04-04
  • 修回日期:  2023-05-31
  • 上网日期:  2023-06-20
  • 刊出日期:  2023-08-20

/

返回文章
返回