基于自发辐射光谱的超声速流场测速技术
齐新华, 陈力*, 闫博, 母金河, 陈爽, 周江宁
中国空气动力研究与发展中心设备设计及测试技术研究所, 四川 绵阳 621000
*通讯作者 e-mail: chenli_03@163.com

作者简介: 齐新华, 1978年生, 中国空气动力研究与发展中心设备设计及测试技术研究所高级工程师 e-mail: qxh_78@sina.com

摘要

等离子体状态参数测量是研究等离子体特性, 开展等离子体模拟再入环境、 等离子体隐身、 等离子体减阻以及边界层控制等研究的重要基础。 利用等离子体射流的自发辐射光谱, 提出了一种基于光学多普勒频移效应的等离子体超声速射流测速的方法。 首先, 测量了等离子体中Ar原子产生的自发辐射光谱, 选择696.54 nm的特征谱线, 作为等离子体发生器测速实验的运动光源; 其次, 使用光谱仪、 传能光纤、 EMCCD相机和高光谱分辨法布里-珀罗(F-P)干涉仪, 设计了高温等离子体速度测量光路; 最后, 在氩壁稳电弧等离子体发生器上, 开展了超声速射流速度测量实验。 实验中, 同一测点的Ar原子产生的自发辐射光谱, 分别被与等离子体射流运动方向成49°和90°夹角的收集透镜收集进入光谱仪, 经光谱仪分光后仅保留特征谱线696.54 nm附近自发辐射光进入传能光纤, 从而消除其他波长的自发辐射光的影响; 光谱仪输出的特征辐射光谱, 经光纤传输及透镜整形成平行光后, 精细度30、 自由光谱范围6.6 GHz的F-P干涉仪, 形成多光束干涉圆环, 并由EMCCD相机采集, 实现对特征谱线的超高精度分辨; 根据多普勒原理, 不同角度收集的同一测点处Ar 696.54 nm特征谱线的频移将有所不同, EMCCD采集的干涉圆环半径也将不同, 通过测量同一级次不同收集方向特征谱线形成的干涉圆环半径改变量, 可测得高温等离子体射流流动速度。 针对同一喷管开展了两车对比试验, 实验测得两车射流轴向速度分别为791和783 m·s-1, 具有较好的重复性。 结果表明基于多普勒效应, 利用高温气体自发辐射光谱, 结合高光谱分辨F-P干涉仪, 能够实现高温等离子体射流速度的精确测量, 该方法属于非接触测量, 不干扰流场, 尤其适用于传统传感器难以应用的高温流场测量。

关键词: 自发辐射光谱; 瑞利散射; 速度
中图分类号:O433.1 文献标志码:A
Velocity Measurement Technology of Supersonic Flow Field Based on Spontaneous Emission Spectrum
QI Xin-hua, CHEN Li*, YAN Bo, MU Jin-he, CHEN Shuang, ZHOU Jiang-ning
Facility Design and Instrumentation Institute, China Aerodynamic Research and Development Center, Mianyang 621000, China
*Corresponding author
Abstract

The measurement of plasma state parameters is an important basis for studying plasma characteristics, including plasma simulated reentry environment, plasma stealth, plasma drag reduction, and boundary layer control. Based on the spontaneous emission spectrum of the plasma jet, a new method of plasma supersonic jet velocity is proposed in this paper. Firstly, the spontaneous emission spectrum generated by argon atoms in the plasma was measured, and the characteristic spectral line of 696.54 nm was selected as the moving light source for the speed measurement experiment of the plasma generator; secondly, the optical path of speed measurement was designed by using a spectrometer, energy transmission fiber, Electron-Multiplying CCD (EMCCD) camera and high spectral resolution Fabry-Perot (F-P) interferometer for high temperature plasma; finally, the velocity measurement experiment was carried out of supersonic jet on an argon wall stabilized arc plasma generator. In this experiment, the spontaneous emission spectrum of argon atom at the same measuring point was collected into a spectrometer by the collecting lens, which wasan angle of 49°and 90°between with the plasma jet motion direction, respectively. After being split by the spectrometer, only the characteristic line of 696.54 nm was retained into the energy transmitting optical fiber, to eliminate the influence of the spontaneous emission spectrum of other wavelengths; the characteristic emission spectrum from the spectrometer, which was transmitted by optical fiber and shaped into parallel light by the lens, irradiated the F-P interferometer with a fineness of 30 and a free spectral range of 6.6 GHz, then a multi-beam interference ring was formed and collected by an EMCCD camera, so as to realize the ultra-high precision resolution of characteristic spectral lines. According to the Doppler principle, the frequency shift of Ar 696.54 nm at the same measurement point collected at different angles was different, and the radius of the interference ring collected by EMCCD was also different. By measuring the radius of the interference ring formed by the characteristic spectral lines for the same level and different collection directions, the flow velocity of the high temperature plasma jet can be computed. The comparative experiments of two vehicles were carried out for the same nozzle, and the axial velocities of the two vehicles were 791 and 783 m·s-1, respectively, which had good repeatability. Based on the Doppler principle, the results show that, using the spontaneous emission spectrum of high temperature gas, combined with the high spectral resolution F-P interferometer, the high temperature plasma jet velocity can be accurately measured. This method belongs to non-contact measurement and does not interfere with the flow field and is especially suitable for the measurement of high temperature flow field, which is difficult to be applied by traditional sensors.

Keyword: Spontaneous emission spectrum; Rayleigh scattering; Velocity
引言

随着科学技术的不断进步, 民用、 国防、 航空航天等领域的等离子体技术达到了一个新的高度, 其中等离子体再入、 等离子体减阻、 等离子体隐身、 等离子体助燃及边界层控制等是其典型应用的表现[1, 2]。 电弧等离子体是地面模拟超高速飞行器大气再入流动环境的主要工具之一, 对其流动特性的研究是高温空气动力学领域中一个非常重要的内容。 等离子体的性质是由其状态参数决定的, 如等离子体的温度、 电子密度及速度等, 对于等离子体温度和电子密度的测量试验研究相关报道较多[3, 4, 5], 而对于等离子体射流的速度测量, 一直缺乏有效的测量手段。 随着激光技术的不断发展, 将激光与流场速度测量联系起来, 先后发展了激光多普勒测速技术(laser Doppler velocimetry, LDV)、 干涉瑞利散射技术(interferometric Rayleigh scattering, IRS)等非接触测量技术, 实现了高速流场的速度测量, 在流体力学研究中获得了广泛应用[6, 7, 8, 9, 10, 11]。 但LDV技术需要在流场中加入示踪粒子, 而示踪粒子在高速流场中跟随性较差, 使其应用范围受到了一定限制; 而IRS技术虽不需外加示踪粒子, 但对超高速、 高超声速燃烧流场, 由于自身存在较强的自发辐射背景, 使得干涉瑞利散射信号受到很大的干扰, 信噪比很低。 因此, 提出了一种利用超高速或高超声速流场自身光谱辐射测量流场速度的方法, 基于氩壁稳电弧等离子体发生器的等离子体射流自发辐射光谱开展了试验研究。

1 基本原理

由多普勒效应[12]知, 当光源与探测器存在相对运动时, 探测器接收到的光频率(波长)将发生改变, 即产生多普勒频移, 其大小与相对运动的速度相关。 因此高速气流产生的自发辐射光相当于运动的光源, 被静止的探测器接收, 如图1所示, 假设光源的运动速度为v, 辐射光沿探测器收集方向的单位矢量为L, 根据多普勒效应, 多普勒频移Δ ν D可以表示为

ΔνD=1λvL=1λvcosθ(1)

图1 运动光源产生的多普勒效应Fig.1 Doppler effect produced by moving light source

式(1)中: λ 为发射光谱波长; θ 为流场速度与收集方向夹角; L为散射光沿探测器收集方向的单位矢量; v为流场速度。

由式(1)即可通过测量发射光谱的多普勒频移来获得流场的速度信息。

高速气流运动速度引起的多普勒频移非常小(MHz~GHz, 10-5 nm量级), 所以对光谱分析设备的分辨能力提出了非常高的要求, 常规光谱仪的分辨本领已远不能满足要求; 而基于多光束干涉原理的F-P(Fabry-Perot, F-P)干涉仪是利用多光束干涉产生细锐条纹的典型仪器, 所产生的干涉条纹具有十分清晰明锐的特点, 是研究光谱线超精细结构的强有力工具[13]

依据F-P干涉仪的工作原理知, 当波长为λ 的单色光经光学透镜形成平行光照射到F-P干涉仪后, 透过F-P干涉仪的多光束将在像平面形成非常明锐的干涉圆环, 形成亮环的条件为

δ=4πλnhcosi=2kπ(2)

式(2)中: δ 为相邻两束透射光之间的相位差; n为干涉仪介质折射率; h为两平面之间的距离; i为光束与光轴的夹角; k为整数。

图2给出了波长分别为λ λ λ 的单色光源在F-P干涉仪像平面形成的不同级次干涉圆环的示意图。

图2 F-P干涉仪形成的同心干涉圆环Fig.2 Concentric interference ring formed by F-P interferometer

在等离子体射流速度测量实验中, 等离子体气流中的某一特征辐射光波长λ 0, 相当于一运动的光源, 当成像透镜焦距为f, 流场速度与收集方向夹角为θ 1时, 由于多普勒效应, 静止的探测器接收到的辐射光波长将改变为λ 1, 其在F-P干涉仪像平面上形成的第k级亮环的半径可表示为

$r_{k}=\sqrt{2-\frac{k \lambda_{1}}{n h}} f$(3)

当流场速度与收集方向夹角改变为θ 2时, 由于多普勒效应, 静止的探测器接收到的辐射光波长将改变为λ 2, 其在F-P干涉仪像平面上形成的第k级亮环的半径同样可表示为

$r_{k}^{\prime}=\sqrt{2-\frac{k \lambda_{2}}{n h}} f$(4)

由式(3)和式(4), 通过分别测量夹角θ 1θ 2对应的干涉圆环半径rkr'k, 并结合式(1), 即可获气流的运动速度v

2 实验部分

基于光学多普勒效应, 利用氩壁稳电弧等离子体发生装置的自发辐射光谱, 开展高温流场速度实验, 首先需要确定合适的高温辐射特征光谱, 作为多普勒频移的运动光源; 其次通过高精度光谱分析设备, 实现对气流运动速度引起的多普勒频移量实现精确测量, 最后根据式(1)计算出气流速度。

2.1 装置

氩壁稳电弧等离子体发生器(图3)由控制电源、 发生器主体、 真空系统、 水电气保障系统及配套的试验测试系统构成, 工作气体为氩气。 其中, 在发生器与真空系统之间安装有拉瓦尔喷管, 可以产生亚音速到超音速气流(图4)。

图3 氩壁稳电弧等离子体发生器装置Fig.3 Argon wall stabilized arc plasma generator facility

图4 氩壁稳电弧等离子发生器形成的高温射流Fig.4 High temperature jet formed by argon wall stabilized arc plasma generator

2.2 辐射光谱筛选

为选择合适的氩原子特征辐射谱线作为多普勒频移测速的基准谱线, 首先开展了高温等离子体辐射光谱测量实验。 试验中, 光谱仪采用300 line· mm-1光栅, 氩气流量3.5 L· min-1, 获取的电流80~130 A状态下的氩原子光谱曲线(见图5); 随着电流的增大, 光谱强度明显增强; 在波长670~760 nm范围内, 有多条原子特征谱线, 其中696, 706, 738, 750和751 nm谱线强度相对较强。 综合后期测速实验中CCD相机的响应曲线, 以及F-P干涉仪的反射率曲线, 最终选择波长696.54 nm的特征谱线, 作为等离子体发生器测速实验的运动光源。

图5 不同电流状态下氩原子光谱曲线Fig.5 Atomic spectrum curves of Argon at different current states

2.3 测速系统搭建

根据测速原理, 设计了图6所示的实验光路。 实验中, 氩等离子体产生的原子谱线, 分别被与气流运动方向成不同夹角的收集透镜L1和L2收集后, 由光纤1传输至光谱仪, 实现光谱的筛选, 再经光纤2、 透镜L3形成平行光照射F-P干涉仪, 最后由透镜L4将干涉圆环成像至CCD相机。

图6 测速实验光路示意图Fig.6 Optical path diagram of velocity measurement experiment

实验中, 为保证较高的收光和传光效率, 光纤1和光纤2均采用的是芯径约0.4 mm的多模光纤; 光谱仪选用的分光光栅为300 line· mm-1; 透镜L1和L2的焦距f=300 mm, 口径为ϕ 100 mm; 透镜L3焦距f=150 mm, L4的焦距f=500 mm; 采用EMCCD(electron-multiplying CCD, EMCCD)相机记录辐射光谱形成的干涉圆环, 该相机具有较强的弱信号放大能力, 保证了采集的干涉图像的信噪比。 F-P干涉仪是整个测速光路中的核心分光元件, 其主要参数为反射率90%, 玻璃板通径50 mm, 厚度15 mm, 精细度30, 自由光谱范围6.6 GHz。

依据陈力等[7, 14]基于F-P干涉仪的流场速度测量准确度分析方法, 并考虑到本工作是利用与射流成49° 和90° 收集夹角的多普勒频移差来计算等离子体射流速度, 经分析该系统的速度测量准确度为7.6 m· s-1

3 结果与讨论

利用图6所示的实验光路, 在氩壁稳电弧等离子体开展了等离子体超射速射流速度测量实验。 针对2个不同的拉瓦尔喷管形成的超声速射流, 开展了速度测量。

图7为1#喷管形成的高温等离子体射流, 经过F-P干涉仪后形成的干涉圆环, 图7(a)信号收集方向与气流运动方向成90° 夹角, 图7(b)信号收集方向与气流运动方向成49° 夹角; 从图中可以看到, 不同的收集方向, 记录的干涉圆环的半径不同。

图7 1#喷管不同角度的干涉圆环Fig.7 Interference rings of nozzle No.1 at different angles

经过数据拟合, 精确获得了每个干涉圆环的干涉半径如图8所示。 当收集夹角为90° 时, 中心两级次干涉圆环半径分别为102.756和190.824 pixel; 当收集夹角为49° 时, 中心两级次干涉圆环半径分别为114.155和196.178 pixel。 通过干涉圆环的半径偏移量, 计算获得1#喷管形成的超射速射流轴向速度为670 m· s-1

图8 1#喷管不同角度干涉圆环拟合结果Fig.8 Fitting results of interference rings of nozzle No.1 at different angles

针对2#喷管, 共开展了两车速度测量实验, 图9和图10分别车次1和车次2高温射流辐射光谱经F-P干涉仪形成的干涉圆环及拟合结果, 进一步计算出两车实验射流轴向速度分别为791和783 m· s-1, 可见两车实验速度测量结果具有较好的重复性。

图9 2#喷管车次1形成的不同角度干涉圆环及拟合结果Fig.9 Interference rings of different angles formed by number 1 of 2# nozzle and fitting results

图10 2#喷管车次2形成的不同角度干涉圆环及拟合结果Fig.10 Interference rings of different angles formed by number 2 of 2# nozzle and fitting results

4 结论

针对高温等离子发生器高速射流速度测量需求, 发展了一种基于辐射光谱的多普勒测速技术。 利用高温等离子体自发辐射的原子特征谱线作为基准光源, 结合光谱仪、 传能光纤、 高光谱分辨F-P干涉仪和电子倍增CCD相机提升系统的光谱分辨能力, 设计了一套高温等离子体超声速射流速度测量光路系统。 针对两个不同的拉瓦尔喷管形成的超声速射流, 开展了高温流场速度测量的实验, 获得了射流轴向速度测量结果, 结果表明基于高温辐射光谱及多普勒效应, 开展高温等离子体射流速度测量是可行的。

参考文献
[1] SUN Zong-xiang(孙宗祥). Advances in Mechanics(力学进展), 2003, 33(1): 87. [本文引用:1]
[2] YANG Bo, BAI Xi-yao, BAI Min-dong, et al(杨波, 白希尧, 白敏冬, ). Modern Defence Technology(现代防御技术), 2005, 33(3): 25. [本文引用:1]
[3] LIU Yao-ze, YUAN Cheng-xun, GAO Rui-lin, et al(刘耀泽, 袁承勋, 高瑞林, ). Chinese Journal of Luminescence(发光学报), 2016, 37(10): 1299. [本文引用:1]
[4] YANG Wen-bin, ZHOU Jiang-ning, LI Bin-cheng, et al(杨文斌, 周江宁, 李斌成, ). Acta Physica Sinica(物理学报), 2017, 66(9): 262. [本文引用:1]
[5] ZHANG Wei, CHEN Lei, SONG Peng, et al(张维, 陈雷, 宋鹏, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2018, 38(12): 3678. [本文引用:1]
[6] YANG Fu-rong, CHEN Li, YAN Bo, et al(杨富荣, 陈力, 闫博, ). Journal of Experiment in Fluid Mechanics(实验流体力学), 2018, 32(3): 82. [本文引用:1]
[7] CHEN Li, YANG Fu-rong, SU Tie, et al(陈力, 杨富荣, 苏铁, ). Acta Photonica Sinica(光子学报), 2015, 44(1): 0112004-1. [本文引用:2]
[8] Bivolaru D, Danehy P M, Lee J W, et al. Single-Pulse Multi-Component Interferometric Rayliegh Scattering Velocimeter AIAA, 2006, AIAA-2006-0836, 44th Aerospace Sciences Meeting, Reno, NV, January 9-12, January. [本文引用:1]
[9] Mielke-Fagan A F, Clem M M, Elam K A. Rayleigh Scattering Measurements Using a Tunable Liquid Crystal Fabry-Perot Interferometer. AIAA, 2010AIAA-2010-4350, 27th Aerodynamic Measurement Technology and Ground Testing Conference, Chicago, Illinois, 2010. [本文引用:1]
[10] Doll U, Stockhausen G, Willert C. Optics Letters, 2017, 42(19): 3773. [本文引用:1]
[11] Estevadeordal J, Jiang N B, Cutler A D, et al. Applied Physics B, 2018, 124: 41. [本文引用:1]
[12] ZHAO Kai-hua, LUO Wei-yin(赵凯华, 罗蔚茵). New Concept Physics Tutorial: Mechanics(新概念物理教程: 力学) Beijing: Higher Education Press(北京: 高等教育出版社), 1995. 359. [本文引用:1]
[13] GUO Yong-kang(郭永康). Optics(光学). Beijing: Higher Education Press(北京: 高等教育出版社), 2005. 180. [本文引用:1]
[14] Chen Li, Yang Furong, Su Tie, et al. Chin. Phys. B, 2017, 26(2): 025205. [本文引用:1]