纳秒脉冲激光斜入射辐照铝靶等离子体羽流及微冲量特性研究
于程浩, 叶继飞*, 周伟静, 常浩*, 郭威
航天工程大学宇航科学与技术系, 激光推进及其应用国家重点实验室, 北京 101416
*通讯作者 e-mail: yjf1981@163.com; changhao5976911@163.com

作者简介: 于程浩, 1995年生, 航天工程大学宇航科学与技术系激光推进及其应用国家重点实验室博士研究生 e-mail: yuchenghao0536@163.com

摘要

为研究脉冲激光斜入射烧蚀铝靶冲量耦合机理, 直接测量其宏观冲量耦合特性是其中一种手段, 但激光烧蚀包含多种物理过程, 仅仅研究其宏观力学性能难以深入分析冲量形成机理, 脉冲激光烧蚀形成的等离子体羽流喷射是诱发力学效应的重要过程, 因此, 在研究宏观力学性能的基础上, 通过开展脉冲激光斜入射烧蚀铝靶等离子体羽流及发射光谱特性测量研究, 深入分析脉冲激光烧蚀冲量耦合机理。 围绕单脉冲1064nm激光斜入射烧蚀铝靶开展研究, 首先通过构建高速摄影测量系统和发射光谱测量系统, 获得了典型激光能量密度斜入射烧蚀铝靶产生的等离子体羽流图像、 等离子体光谱图像和等离子体发射光谱, 基于等离子体发射光谱, 利用Boltzmann作图法和Stark展宽法, 分别研究了脉冲激光多种斜入射角度下等离子体温度、 电子数密度随能量密度的变化关系; 通过搭建扭摆微冲量测量系统, 研究了脉冲激光多种斜入射角度下, 沿着激光入射方向的冲量耦合系数随能量密度的变化。 研究中遵循从羽流微尺度演化过程到冲量宏观力学性能测量分析的研究思路。 实验结果表明, 随着能量密度的增加, 等离子体羽流发光强度增强, 羽流离化程度增加, 等离子体温度、 电子数密度均先迅速增加, 冲量耦合系数也迅速增加; 当能量密度大于15 J·cm-2时, 由于等离子体屏蔽效应, 等离子体温度、 电子数密度均逐渐趋于饱和, 最终导致冲量耦合系数随着能量密度的增加而减小; 此外, 随着入射角度的增加, 等离子体温度、 电子数密度均逐渐减小, 导致冲量耦合系数也随之减小。 研究结果表明, 利用高速摄影和发射光谱可较好地分析脉冲激光烧蚀冲量耦合机理, 研究结果可为激光空间碎片清除、 空间微推力器、 空间非合作目标消旋等空间应用的关键参数优化提供参考。

关键词: 激光烧蚀; 斜入射; 等离子体羽流; 高速摄影; 发射光谱; 等离子体屏蔽; 冲量耦合
中图分类号:TN249 文献标志码:A
Characteristics of the Plasma Plume and Micro-Impulse Generated by Irradiating the Aluminum Target With a Nanosecond Laser Pulse at Oblique Incidence
YU Cheng-hao, YE Ji-fei*, ZHOU Wei-jing, CHANG Hao*, GUO Wei
State Key Laboratory of Laser Propulsion & Application, Department of Astronautical Science and Technology, Space Engineering University, Beijing 101416, China
*Corresponding authors
Abstract

To study the impulse coupling mechanism of a pulsed laser ablation aluminum target, the direct measurement of its macroscopic impulse coupling characteristics is one of the means. However, laser ablation involves many physical processes. Therefore, analysing the impulse formation mechanism only by studying its macroscopic mechanical properties is difficult. Plasma plume ejection formed by pulsed laser ablation is an important process to induce the mechanical effect. Hence, based on studying the macroscopic mechanical properties, this paper deeply analyzes the impulse coupling mechanism of pulsed laser ablation by measuring the plasma plume and emission spectrum characteristics. In this paper, a single pulse laser with a wavelength of 1 064 nm is used to ablate aluminum targets. By constructing a fast photogrammetry system and optical emission spectroscopy measurement system, the plasma plume image, the plasma spectral image, and the plasma emission spectrum generated by laser oblique incident ablation of the aluminum target were obtained. Based on optical emission spectroscopy of the plasma plume, the Boltzmann plotting method and Stark broadening method were used to study the variation of plasma temperature and electron number density with the laser fluence at different incidence angles of a pulsed laser, respectively. Moreover, a torsion pendulum system was built to study the trend of the impulse coupling coefficient with the laser fluence along the direction of laser incident at various incident angles. The study follows the research ideas from the plume microscale evolution process to impulse macro mechanical properties analysis. The experimental results show that the luminescence intensity of the plasma plume strengthens with the laser fluence, accompanied by the rise of plume ionization degree. Moreover, the plasma temperature and electron number density increase rapidly, resulting in the impulse coupling coefficient heightening rapidly. When the laser fluence is greater than 15 J·cm-2, the plasma temperature and the electron number density are gradually saturated due to the plasma shielding effect. The change of plasma temperature and electron number density results in the decrease of impulse coupling coefficient with increased laser fluence. In addition, the plasma temperature and the electron number density decrease with the increase of incident angle, which reduces the impulse coupling coefficient. The results show that the coupling mechanism of the ablation impulse can be well analyzed using fast photography and optical emission spectroscopy. The results can provide a reference for optimising key parameters for space applications such as laser space debris removal, space micro-thruster and despinning non-cooperative targets in space.

Keyword: Laser ablation; Oblique incidence; Plasma plume; Fast photography; Optical emission spectroscopy; Plasma shielding; Impulse coupling
引言

激光烧蚀推进(laser ablation propulsion, LAP)可以远距离辐照空间目标, 使目标表面产生高温高压等离子体羽流, 羽流高速膨胀离开目标表面, 根据动量守恒定律可知, 目标会获得反冲冲量[1], 最终实现空间操控的目的。 激光烧蚀推进为空间碎片清除[2]、 空间非合作目标的消旋[3]等热点空间领域提供了解决方案。 激光烧蚀推进过程可以分为两个阶段[4]: 第一阶段, 脉冲激光辐照靶面, 部分激光能量被反射, 其余激光能量在靶材表面沉积, 导致辐照区温度升高; 第二阶段, 当辐照区靶材温度达到临界值时, 靶材熔化并产生靶蒸气, 靶蒸气进一步电离形成高温高压的等离子体, 形成的等离子体羽流还会与后续激光相互作用, 减弱激光能量与靶的耦合, 导致等离子体屏蔽效应[5]。 由激光与物质相互作用过程可以发现, 等离子体羽流对冲量耦合有着重要的影响, 然而, 关于等离子体羽流特性对激光烧蚀力学性能影响的研究依然很少。 脉冲激光烧蚀等离子体羽流特性受激光波长[6]、 激光能量[7]、 激光光斑尺寸[5]、 激光入射角度[8]、 激光脉冲数量[9]等因素影响。 针对激光烧蚀推进的空间应用, 国内外学者对脉冲激光烧蚀产生的冲量耦合特性进行了大量的研究, 但是以往的冲量耦合机理研究多集中在激光参数(波长、 脉宽等)[10]、 材料[11]等方面。 近年来, 随着对激光烧蚀推进的进一步研究, 激光入射角度的影响开始受到关注[1]

针对激光空间碎片清除等激光烧蚀推进应用背景, 基于高时空分辨率的羽流观测系统、 发射光谱测量系统, 开展激光斜入射辐照下等离子体羽流喷射过程和发射光谱特性研究, 通过建立扭摆微冲量实验测量系统, 研究激光斜入射辐照下的冲量耦合特性, 通过结合微观羽流喷射, 分析了激光斜入射辐照下的冲量耦合机理。

1 实验部分
1.1 实验系统

脉冲激光斜入射辐照铝靶的实验系统示意图如图1所示。 由于扭摆的微小摆动会给光谱测量和羽流拍摄产生影响, 因此, 利用分束镜和反射镜将激光折转为两束, 其中一束进行扭摆微冲量测量, 另一束进行羽流和光谱的测量, 最后进行耦合分析。 利用DG645实现对激光器、 ICCD相机以及光谱仪的时序控制, 通过能量计实现对入射激光能量的实时监测。 实验过程中保持真空仓内压强为10-3 Pa, 使用Nd∶ YAG激光器作为烧蚀激光源辐照铝靶(Al 1050), 脉冲激光波长为1 064 nm, 脉宽为8 ns。 如图1中所示的扭摆结构局部放大图, 将铝靶固定在扭摆横梁顶部, 入射激光方向保持固定, 通过旋转靶片, 实现相对入射角度控制, 因此扭摆将测量得到烧蚀冲量沿着激光入射方向的分量。

图1 实验系统示意图Fig.1 Schematic diagram of experimental system

1.2 发射光谱方法

等离子体温度是研究等离子体羽流特性的重要参数之一, 当激光烧蚀等离子体羽流满足局部热力学平衡条件时, 可利用Boltzmann作图法获得等离子体温度, 选取同一电离态的多条特征谱线进行计算[7]

lnIkiλkgkAki=-EkkBT+lnN(T)U(T)(1)

式(1)中, Iki为谱线强度, λ 为跃迁波长, gk为能级统计权重, Aki为跃迁概率, Ek为能级能量, kB为Boltzmann常数, N(T)为总发射粒子数密度, U(T)为配分函数, T为等离子体温度。 在NIST光谱数据库中可获得gk, AkiEk等参数的值。 利用图2中拟合直线的斜率(-1/kBT), 即可计算获得等离子体温度。

图2 Al Ⅲ 特征谱线的典型Boltzmann斜线图Fig.2 Typical Boltzmann oblique plot of the Al Ⅲ spectral line

对等离子体羽流电子数密度的研究, 有助于分析脉冲激光烧蚀冲量耦合机制。 对于电子数密度的计算, 通常采用Stark展宽法, 基于原子谱线或单电荷离子谱线的Stark展宽Δ λ 1/2与电子数密度之间的简化关系为[6]

Δλ1/2=2ωNe1016(2)

式(2)中, ω 为电子碰撞参量, Ne为电子数密度。 Δ λ 1/2为Stark半高全宽, 可通过Lorentz函数拟合光谱实验数据获得, 典型特征谱线的Lorentz函数拟合如图3所示。

图3 466.31 nm特征谱线的Lorentz函数拟合Fig.3 The Lorentz function fitting of 466.31 nm spectral line

利用Boltzmann作图法进行等离子体温度求解的前提条件是等离子体羽流满足局部热力学平衡条件, 实际情况中, 往往先假设等离子体羽流满足局部热力学平衡条件, 在计算获得等离子体温度和电子数密度后, 再利用McWhirter标准判据验证假设是否成立[12]

ne1.6×1012T12ΔE3(3)

式(3)中, T是实验中计算获得的最大等离子体温度值, Δ E是计算等离子体温度时所用谱线的最大上下能级差。

1.3 微冲量测量

当高能脉冲激光辐照铝靶后, 基于扭摆微冲量测量原理[13], 利用传感器获得的扭摆横梁最大扭转角, 即可计算获得脉冲激光烧蚀产生的微冲量。 通常采用冲量耦合系数Cm对冲量耦合性能进行评估, 即冲量I与入射激光能量E的比值(Cm=I/E)。

2 结果与讨论
2.1 等离子体羽流高速摄影

实验发现, 等离子体羽流的持续时间为百纳秒左右, 因此, 将ICCD相机门宽设置为500 ns, 在脉冲激光辐照靶面后, 开始记录整个等离子体羽流喷射过程。 典型能量密度激光辐照铝靶后的自发光等离子体羽流图像如图4所示。 Φ 0表示在激光入射角度为0° 时的能量密度, 即垂直入射时对应的能量密度。 而当激光斜入射角度为θ 时, 对应的能量密度为Φ 0cosθ 。 在入射角度为0° 时, 自发光羽流的结构大致为半球形, 并且随着能量密度的增加, 发光强度明显增加, 这可能是由于靶材的烧蚀率[14]和电离程度[15]的增加导致的。 当入射角度为15° , 30° , 45° , 60° , 75° 时, 等离子体羽流强度均随能量密度的增加而增强。 在脉冲激光入射角度为0° , 15° 和30° 且激光能量密度较高时, 出现了强自发光区域— — 等离子体“ 核” [5], 这个区域是由于逆韧致吸收和多光子电离过程导致的, 该区域形成后被认为可进一步吸收后续入射激光能量, 从而形成等离子体屏蔽效应。 由图4还可以看出, 随着激光入射角度的增加, 羽流发光强度逐渐减弱, 羽流沿着靶面法向的喷射减弱, 而沿着平行靶面的方向喷射逐渐增强。

图4 典型入射角度下, 等离子体羽流图像随能量密度的变化Fig.4 The variation of the plasma plume image with the laser fluence at typical incident angles

2.2 等离子体发射光谱

为了进一步分析入射角度对羽流特性的影响, 开展了等离子体羽流参数特性研究。 当一束高能激光聚焦在目标靶片时, 电子碰撞激发促使一小部分原子进入激发态, 激发态具有不稳定性, 随后向低能级或基态跃迁并向外发射出特征波长的光子。 图5给出了350~650 nm范围内原子(Al Ⅰ )和离子(Al Ⅱ , Al Ⅲ )的发射光谱, 主要的波长分布与Shaim等[12]在相似实验工况下烧蚀铝靶获得的发射光谱结果基本一致。

图5 350~650 nm波长范围内的典型发射光谱Fig.5 Typical emission spectra in the wavelength range of 350~650 nm

基于图5获得的发射光谱实验结果, 我们选取451.26, 452.92, 569.66以及572.27 nm波长处的Al Ⅲ 分立谱线进行等离子体温度的计算, 计算过程中涉及到的相关参数列于表1中。 对于光学薄等离子体来说, 通常选用中性原子或者一价离子的Stark展宽进行电子数密度的计算。 谱线展宽主要由Stark展宽、 共振展宽、 Doppler展宽三种主要的展宽机制导致, 对于本实验中脉冲激光烧蚀铝靶产生的等离子体, Doppler展宽和共振展宽的影响与Stark展宽相比较小, 可以忽略不计。 在实验中, 选用Al Ⅱ (466.31 nm, 3s4p(1 P10)-3p2(1D2))特征谱线的Stark展宽进行电子数密度的计算。 综合考虑等离子体温度和电子数密度的计算需求, 我们选取450~470 nm波段和555~575 nm波段开展光谱测量。

表1 谱线相关参数 Table 1 Parameters related with spectral lines

实验中, 通过调整光纤端部的聚焦透镜到等离子体羽流之间的距离, 将整个等离子体羽流的辐射耦合到光纤中, 并传输至光谱仪。 使用1 200 lines· mm-1光栅进行光谱采集, 狭缝设置为10 μ m, 从脉冲激光出光时刻开始启动光谱仪工作, 门宽设置为1 μ s。 基于以上实验设置, 获得了不同激光入射角度辐照铝靶时, 发射光谱强度随能量密度的变化。 如图6所示, 典型入射角度下, 不同能量密度对应的发射光谱强度。 可以发现, 光谱强度随着能量密度的增加而增加, 但是能量密度大于15 J· cm-2后, 发射光谱强度增幅减小。 这是由于等离子体屏蔽效应导致的, 即激光脉冲前沿产生的等离子体减弱了靶材对后续激光脉冲的吸收[16], 使得光谱强度增幅减小。

图6 (a) 450~470 nm波段以及(b) 555~575 nm波段, 不同能量密度下的典型光谱强度Fig.6 Typical spectral intensities at different laser fluences in the wavelength range of 450~470 (a) and 555~575 nm (b)

脉冲激光斜入射辐照铝靶时, 电子数密度随能量密度Φ 0的变化趋势如图7所示。 可以看出, 电子数密度随着能量密度的增加逐渐增加。 这是因为随着能量密度的增加, 脉冲激光烧蚀羽流从气化机制向离化机制转变。 在低能量密度, 烧蚀羽流以靶蒸气为主, 随着能量密度的增加, 靶蒸气进一步吸收激光能量, 离化产生大量高温高密度等离子体, 因此随着能量密度的增加, 电子数密度也随之增加。 当能量密度大于15 J· cm-2时, 电子数密度的增加逐渐减缓, 原因是形成的等离子体逐渐增强, 等离子体与激光发生强相互作用, 通过逆韧致和多光子电离过程吸收入射激光能量[4], 激光能量不能完全与靶耦合, 导致与靶材的耦合效率降低, 因而形成的电子数密度增加变缓[7, 16]。 由图7还可以看出, 随着激光入射角度的增加, 电子数密度随之减小。

图7 不同入射角度对应的电子数密度随能量密度的变化关系Fig.7 The variation of the electron number density with the laser fluence at different incident angles

对任何光源的温度测量可以提供关于光源中物理和化学现象等重要信息, 研究激光烧蚀等离子体温度有助于更好地理解等离子体中复杂现象。 图8给出了等离子体温度随能量密度的变化。 从图8中可以看出, 在低能量密度区, 等离子体温度均比在高能量密度区增加得快。 在低能量密度区, 等离子体逆轫致辐射吸收系数较小, 随着能量密度的增加, 吸收系数迅速增大, 这也导致了该区域等离子体温度快速的增加; 而当能量密度较高时, 逆轫致辐射吸收系数达到饱和[16], 产生等离子体屏蔽效应, 使得等离子体温度趋于饱和[17]。 当能量密度Φ 0一定, 随着入射角度的增加, 等离子体温度逐渐降低。

图8 不同入射角度对应的电子温度随能量密度的变化关系Fig.8 The variation of the plasma temperature with the laser fluence at different incident angles

为了进一步分析等离子体羽流的离化程度, 开展了等离子体羽流光谱图像研究。 基于图5的特征谱线分布, 我们选择中心波长为450 nm、 半带宽为10 nm的带通窄带滤光片放置于ICCD相机前, 用于获得Al Ⅲ (447.99, 451.26和452.92 nm)的空间分布。 图9给出了典型能量密度辐照下的Al Ⅲ 组分空间分布。 可以发现, Al Ⅲ 羽流辐射强度随着能量密度的增加而增加, 即烧蚀羽流的电离程度随之增加[18]; 当能量密度大于15 J· cm-2时, Al Ⅲ 羽流辐射强度增加不明显, 这可能是因为羽流电离程度在能量密度较高时接近饱和[15]。 此外, 随着入射角度的增加, Al Ⅲ 羽流辐射强度逐渐降低。

图9 不同入射角度对应的Al Ⅲ 羽流随能量密度的变化关系Fig.9 The variation of the Al Ⅲ plume with the laser fluence at different incident angles

根据McWhirter标准验证发现[12], 所有工况下测量到的电子密度均高于McWhirter准则设定的临界值, 因此满足前期局部热力学平衡的假设。

2.3 烧蚀冲量

烧蚀冲量I与能量密度Φ 0的关系如图10所示。 从图10中可以明显看出, 当入射角度一定, 冲量随能量密度的变化分为两个线性阶段。 能量密度小于15 J· cm-2为第一阶段, 冲量随着能量密度的增加迅速增加。 能量密度大于15 J· cm-2为第二阶段, 等离子体屏蔽效应[19]导致冲量的增加相比于第一阶段有所减缓。

图10 不同入射角度对应的冲量随能量密度的变化关系Fig.10 The variation of the impulse with the laser fluence at different incident angles

冲量耦合系数Cm可以表征激光烧蚀推进的性能, 图11给出了冲量耦合系数Cm随能量密度Φ 0的变化。 可以看出, 随着能量密度的增加, Cm急剧增加, 并达到最大值。 这是由于随着能量密度的增加, 等离子体温度和电子数密度均迅速增加, 因此Cm急剧增加。 随着能量密度的进一步增加, 由于等离子体屏蔽效应[19], 等离子体温度和电子数密度的增加变缓甚至饱和, 进而导致Cm开始缓慢下降, 此时对应的冲量I增加减缓。 当能量密度Φ 0一定, 随着激光入射角度的增加, 等离子体温度和电子数密度均逐渐减小, 引起冲量耦合系数Cm的降低。

图11 不同入射角度对应的冲量耦合系数随能量密度的变化关系Fig.11 The variation of the impulse coupling coefficient with the laser fluence at different incident angles

3 结论

采用1 064 nm脉冲激光斜入射辐照铝靶产生等离子体, 利用高速摄影研究了等离子体羽流特性; 利用发射光谱法研究了等离子体羽流的电子数密度和等离子体温度特性; 基于扭摆微冲量测量原理研究了沿着激光入射方向的冲量耦合特性。 研究发现, 随着能量密度的增加, 等离子体羽流始终沿着靶面法向喷射, 等离子体羽流发光强度增强, 等离子体温度和电子数密度迅速增加, 导致冲量耦合系数Cm迅速增加; 随后由于等离子体屏蔽效应的影响, 冲量耦合系数Cm缓慢下降, 此时, 等离子体羽流贴近靶面处出现了强辐射区域, 等离子体温度和电子数密度增加变缓甚至饱和。 此外, 当能量密度Φ 0一定时, 随着入射角度的增加, 等离子体羽流强度逐渐减弱, 并且等离子体羽流沿着靶面法向的喷射减弱, 沿着靶面平行方向的喷射增强, 等离子体温度和电子数密度均减少, 冲量耦合系数Cm也随之降低。

参考文献
[1] Wang B. Applied Physics Letters, 2017, 110(1): 014101. [本文引用:2]
[2] Ebisuzaki T, Quinn M N, Wada S, et al. Acta Astronautica, 2015, 112: 102. [本文引用:1]
[3] Kumar R, Sedwick R J. Journal of Spacecraft and Rockets, 2015, 52(4): 1129. [本文引用:1]
[4] Marla D, Bhand arkar U V, Joshi S S. Journal of Applied Physics, 2011, 109(2): 021101. [本文引用:2]
[5] Li X, Wei W, Wu J, et al. Journal of Applied Physics, 2013, 113(24): 243304. [本文引用:3]
[6] Hussein A E, Diwakar P K, Harilal S S, et al. Journal of Applied Physics, 2013, 113(14): 143305. [本文引用:2]
[7] ZHENG Pei-chao, LIU Hong-di, WANG Jin-mei, et al(郑培超, 刘红弟, 王金梅, ). Acta Photonica Sinica(光子学报), 2014, 43(9): 82. [本文引用:3]
[8] Chang H, Jin X, Zhou W. Optik, 2014, 125(12): 2923. [本文引用:1]
[9] Smijesh N, Rao K H, Chetty D, et al. Optical Letters, 2017, 43: 6081. [本文引用:1]
[10] WEN Ming, LI Nan-lei, WU Jie(文明, 李南雷, 吴洁). Infrared and Laser Engineering(红外与激光工程), 2017, 46(S1): 1. [本文引用:1]
[11] DU Bang-deng, XING Bao-yu, YE Ji-fei, et al(杜邦登, 邢宝玉, 叶继飞, ). Infrared and Laser Engineering(红外与激光工程), 2020, 49(S2): 147. [本文引用:1]
[12] Shaim Md H A, Elsayed-Ali H E. Journal of Applied Physics, 2017, 122(20): 203301. [本文引用:3]
[13] Zhou W, Chang H, Ye J, et al. Laser Physics, 2020, 30(6): 066002. [本文引用:1]
[14] Vladoiu I, Stafe M, Popescu I M. University Politehnica of Bucharest Scientific Bulletin-Series A: Applied Mathematics and Physics, 2007, 69(3): 81. [本文引用:1]
[15] CHANG Hao, JIN Xing, WEN Ming, et al(常浩, 金星, 文明, ). High Power Laser and Particle Beams(强激光与粒子束), 2013, 25(5): 1110. [本文引用:2]
[16] Haq S U, Ahmat L, Mumtaz M, et al. Physics of Plasmas, 2015, 22(8): 083504. [本文引用:3]
[17] Amin S, Bashir S, Anjum S, et al. Physics of Plasmas, 2017, 24(8): 083112. [本文引用:1]
[18] Sinko J E, Phipps C R. Applied Physics Letters, 2009, 95(13): 131105. [本文引用:1]
[19] Phipps C R, Luke J R, Lippert T, et al. Applied Physics A, 2004, 79(4-6): 1385. [本文引用:2]