-
为实现高效率激光退火,中国科学院近代物理研究所表面处理团队搭建了一套以kW级纳秒脉冲光纤激光器为基础的激光退火系统。
系统使用深圳水滴激光科技有限公司的 HL1000-RMP2 红外纳秒脉冲光纤激光器,固定参数激光光斑边长1 mm,脉宽25 ns,频率2~50 kHz内可调,单脉冲能量在100 mJ范围内可调,中心波长1 064 nm。参照实验条件(表1)对样品进行表面退火,模型设置如下。
参数 大小 单脉冲能量 100 mJ内可调 光斑长度 1 mm 脉宽 25 ns 频率 2~50 kHz 激光退火光路设计用到F200 mm复合镜片准直镜、动态反射镜、高速旋转电机、动态聚焦电机和动态反射镜电机等主要部件。铌腔视为轴向对称结构,将其横向分为五部分,激光头沿着铌腔中轴线移动时,需要激光处理区域的内径持续变化,同时待激光处理内腔表面沿中轴线的角度相对中轴线持续变化。因此,为保证激光光斑表面积在待处理内腔表面始终保持一致,采用上述三套电机。高速电机控制脉冲和脉冲重叠率,动态聚焦电机控制光斑直径保持光斑直径不变,动态反射电机保证激光与内腔待清洗区域始终保持90度直角,其光路图如图2所示。采用F200 mm的原因是其可在满足动态聚焦需求下保证光斑可以聚焦在1 mm以下。激光头由机械手夹持,沿着铌腔中轴线匀速从入口处向内推进,在推进过程中,动态聚焦镜电机与动态反射镜电机实时改变位置,实现三轴联动,确保激光照射在腔内部任意位置均为垂直角度,且光斑尺寸保持一致,针对性地解决了激光在复杂环境下覆盖不均匀的问题。于此同时,整个机构高速旋转,覆盖腔内轴向周圈内所有局部需照射部位。由于轴对称结构,处理完半边铌腔后,将激光头从铌腔另一入口伸入,重复以上激光处理步骤,分两步完成整腔激光处理。
-
针对铌的激光退火,本文考虑两种退火原理:
第一,激光加热使表面温度至铌再结晶温度以上,熔点以下,铌晶粒在受热过程中发生再结晶,由于初始位错、晶界较多,初始晶核易于在晶界上(也有可能在其他缺陷上)形成,再结晶将从这些位置开动,推测温度越高,再结晶粒径越大。退火结束后内应力消除,材料组织缺陷减少。第二,激光能量输出持续加大,表面温度升至铌熔点以上,使铌膜熔化重结晶,即发生相变重结晶退火。
两者相比,再结晶退火效果一般比重结晶的更好,所以在退火中要控制能量输出,尽量避免发生相变。但考虑到激光器的能量可调范围大,能量输出较高时可能会导致铌膜熔化,模拟过程中加入了相变模块。
如图3所示,激光作用形式可分为热吸收、热辐射、热扩散三部分,将热辐射和热扩散都看作能量损失,铌膜获得的能量仅来自热吸收,材料通过热吸收而不断升温。其中热吸收与材料的吸收率有关,常温下铌对1 064 nm波长的光吸收率约为0.2,且随温度不断变化。热辐射与表面辐射率有关,本文参考值为0.03。热扩散包括热传递和热交换,与材料属性和环境温度有关。其关系式将在2.3节中给出。
-
目前市场上的光纤激光普遍有单模与多模两种模式,本文采用的激光器激光类型为多模。单模激光的能量输出满足高斯分布,高斯分布中每个脉冲都存在明显上升沿和下降沿,最大值为脉冲峰值。而本文激光器出射激光(多模)光强分布可近似看作平顶光,脉冲较为稳定,能量分布均匀。
平顶光指光斑的能量在作用范围内均匀分布,其物理形式可以表示为
高斯光斑能量分布呈正态分布,其物理形式为
其中:
$ x_0 $ 和$ y_0 $ 表示光斑圆心坐标;x, y代表空间任意位置;$ \sigma_x $ 和$ \sigma_y $ 为标准差,为激光脉宽的一侧上升沿长度;$ A_0 $ 表示激光的能量振幅。高斯光处理时效率较低,易损伤材料,而平顶光能量输出稳定,使材料表面受热更加均匀,损伤性小,效率更高。 -
纯铌是一种熔点高、导热性差的金属,其密度、热容和导热系数受温度影响,即随退火进行上述材料属性不断改变,变化规律如表2所列。
参数 状态 固态 液态 密度
(kg·m−3)8 617.45−0.14T−8.08$ \times 10^{-5}T^2$ 8 580.66−0.23T 恒压热容
$C_{\rm{p}}$(J·kg·K−1)$\begin{array}{c} 236.95+0.106{{T} }-6.08 \times \\ 10^{-5}{T }^2+1.90 \times 10^{-8}T^3+\\2.35 \times 10^5T^{-2} \end{array}$ $\begin{array}{c} 360.28+5.38 \times 10^{-8}T-\\ 5.24 \times 10^{-13}T^2-4.89 \times \\ 10^{-16}T^3+ 624.28T^{-2} \end{array}$ 导热系数
$ k $/(W·m·K−1)48.67+0.016T−6.37$ \times 10^{-7}T^2$ 47.23+0.008T 热吸收过程的热传递方程可以表示为
上述表达式中:
$ ρ $ 表示材料密度;$C_{\rm{p}}$ 为材料恒压热容;$ u $ 表示表面熔化流动的速度场;$ T $ 表示环境温度;Q表示热源热量;Q$ _{{\rm{ted}}} $ 表示表面传递热量;$ q $ 是热传递矢量;可以写成:若材料表面温度升高至相变点
$T_{\rm{m}}$ 时将发生固液相变,由固态转变为熔融态,铌达到熔点后金属键作用会逐渐减小,使其无法保持晶体状态,该过程不断吸收能量,熔化并非瞬间完成,具有一定持续性,理论上此时温度保持不变,直到激光作用区材料完全熔化为止。该物理过程可表示为其中
$ \theta_1 $ 和$ \theta_2 $ 分别代表固相和液相占现有体系的百分比,$ \theta_1+\theta_2=1 $ ;$ L_{12} $ 是熔化潜热;$ k_1 $ 和$ k_2 $ 分别是固相和液相的热传递系数。在理论条件下,材料发生相变的过程中温度不变,即∆T = 0,而模拟软件计算为了数值收敛,∆T不能为0,且为了减小计算量,从较大∆T(如300 K→200 K→75 K)逐步递减使其无限趋近与0。表面辐射和热对流都看作能量损失,其方程表示为在表面辐射方程中:
$ \epsilon $ 表示铌材的环境辐射率;$ \sigma $ 是Stefan-Boltzmann常数;$ T_{{\rm{amb}}} $ 表示环境温度;$ T_{\rm{t}} $ 为材料内部温度,是与时间有关的函数;在热对流方程中,$ T_{{\rm{ext}}} $ 表示外部温度;h为传热系数;$ n $ 是常数;$ q $ 是热传递矢量;$ q_0 $ 为沿热传递矢量方向$ q $ 传递的总热量。
Simulation on Temperature Field of Laser-annealed Copper-based Niobium Thin Film by Finite Element Analysis Method
doi: 10.11804/NuclPhysRev.40.2022086
- Received Date: 2022-09-23
- Rev Recd Date: 2022-12-22
- Available Online: 2024-02-04
- Publish Date: 2023-12-20
-
Key words:
- copper-based niobium thin film /
- laser annealing /
- RF superconductivity /
- finite element analysis /
- fiber laser
Abstract: Copper-based niobium thin film cavity is one kind of superconducting radio frequency(SRF) cavities developed from bulk niobium cavity. This new type of cavity has similar SRF properties to the pure niobium cavities but is more economical, mechanically stable, and thermally stable, making it promising for industrial SRF applications. However, the inevitable surface defects of niobium film on the surface of Nb/Cu cavities induce extra power dissipation, and require post-treatment before cavity operation. Near-surface annealing of SRF cavity by nanosecond high-intensity pulsed laser is an emerging surface treatment technology for SRF cavity. Using high power and short pulsed laser with peak power flux exceeding 100 MW/cm2 to locally recrystallize the surface can reduce the surface roughness, increase the grain size and eliminates some surface defects, so as to obtain better surface RF properties. At present, the type of laser selected in this direction is mainly the solid-state laser with high peak power, which has a very low average power, not practical for the actual superconducting cavity with an m2 inner wall surface area. To solve this issue, a whole-cavity laser processing system based on a kW-level power nanosecond pulse fiber laser has been established by the Surface Treatment Research Team of the Institute of Modern Physics. In this study, finite element simulation was conducted with the parameters of the system to investigate the temperature distribution on the surface and inside the material under near-actual working conditions. Combined with the material properties, the laser surface treatment effect of pure niobium material was simulated. This study preliminarily confirms the feasibility of using nanosecond pulse laser for complete annealing of the inner wall of superconducting cavities, and provides an effective method for achieving better RF superconducting performance.
Citation: | Changlin WANG, Teng TAN, Pingran XIONG, Didi LUO, Tongtong ZHU, Qingwei CHU. Simulation on Temperature Field of Laser-annealed Copper-based Niobium Thin Film by Finite Element Analysis Method[J]. Nuclear Physics Review, 2023, 40(4): 541-546. doi: 10.11804/NuclPhysRev.40.2022086 |