基于电子束荧光的稀薄流场转动温度测量与数据处理方法与流程

1.本发明属于高超声速稀薄流风洞试验技术领域，具体涉及一种基于电子束荧光的稀薄流场转动温度测量与数据处理方法。


背景技术：

2.在高超声速稀薄流风洞流场中，由于空气中的氮气为双原子分子，在稀薄流动中会出现热力学非平衡效应，氮气气体分子的振动温度和转动温度不一致，需要测量氮气气体分子的振动温度和转动温度，用于热力学非平衡效应研究。目前，采用电子束荧光技术测量流场中气体的振动温度和转动温度，其中振动温度测量通过强度比方法和标定实现，转动温度测量通过测量氮气分子的电子束荧光转动光谱各枝强度拟合斜率计算获得。但是，由于氮气分子有两个转动自由度，在室温静态真空环境、高超声速稀薄空流场和激波后流场等不同的流场情况下，氮气分子的两个转动自由度会可能出现不一致，从平衡态到非平衡态，即整体转动自由度不服从玻尔兹曼平衡态分布，会导致转动谱线出现较为明显的两个斜率，难以给出一个合理的转动温度，仅采用同一个拟合斜率进行计算难以满足不同的流场情况下转动谱线的处理要求。
3.当前，亟需发展一种基于电子束荧光的稀薄流场转动温度测量与数据处理方法。


技术实现要素：

4.本发明所要解决的一个技术问题是提供一种基于电子束荧光的稀薄流场转动温度测量方法，本发明所要解决的另一个技术问题是提供一种基于电子束荧光的稀薄流场转动温度测量方法的数据处理方法。
5.本发明的电子束荧光的稀薄流场转动温度测量方法使用的测量装置包括安装在高超声速低密度风洞的试验段上驻室的电子枪，安装在高超声速低密度风洞的试验段下驻室的法拉第杯，法拉第杯用于接收电子枪发射的电子束，产生的电子束荧光位于喷管和扩压器之间并垂直于喷管的轴线从上至下穿过试验段的中心空腔；
6.通过导线连接的光谱仪和ccd相机布置在试验段的观察窗外，光谱仪与观察窗之间还布置有凸透镜；
7.ccd相机与计算机之间通过网线连接；
8.观察窗采用光学玻璃；光谱仪的波长分辨率大于等于0.02nm，测量波长范围为388nm～391nm；
9.其特点是，所述的测量方法包括室温静态真空环境的转动温度测量方法、高超声速稀薄空流场的转动温度测量方法和激波后流场的转动温度测量方法，分别如下：
10.a.室温静态真空环境的转动温度测量方法；
11.先将高超声速低密度风洞的试验段抽真空至20pa以下，开启电子枪发射电子束，在试验段的测量区域产生的电子束荧光通过观察窗的光学玻璃被凸透镜收集至光谱仪的狭缝入口，电子束荧光在光谱仪内产生的电子束荧光转动谱线，通过ccd相机记录电子束荧
光转动谱线的波长和强度数据并传输至计算机存储；
12.b.高超声速稀薄空流场的转动温度测量方法；
13.先将高超声速低密度风洞的试验段抽真空至20pa以下，开启高超声速低密度风洞吹风，开启电子枪发射电子束，在试验段的测量区域产生的电子束荧光通过观察窗的光学玻璃被凸透镜收集至光谱仪的狭缝入口，电子束荧光在光谱仪内产生的电子束荧光转动谱线，通过ccd相机记录电子束荧光转动谱线的波长和强度数据并传输至计算机存储；
14.c.激波后流场的转动温度测量方法；
15.先将试验模型安装在高超声速低密度风洞的试验段内，试验段抽真空至20pa以下，开启高超声速低密度风洞吹风，开启电子枪发射电子束，调整电子束位置到试验模型的头部激波后方，在试验段的测量区域产生的电子束荧光通过观察窗的光学玻璃被凸透镜收集至光谱仪的狭缝入口，电子束荧光在光谱仪内产生的电子束荧光转动谱线，通过ccd相机记录电子束荧光转动谱线的波长和强度数据并传输至计算机存储。
16.本发明的基于电子束荧光的稀薄流场转动温度测量方法的数据处理方法，包括以下步骤：
17.s10.绘制电子束荧光转动谱线；
18.将计算机内存储的电子束荧光转动数据绘制成横坐标为波长、纵坐标为转动强度的电子束荧光转动谱线；
19.s20.标记电子束荧光转动谱线的激发态转动量子数j
′
；
20.激发态转动量子数j
′
跃迁至基态时转动量子数j
″
，取j
′‑j″
＝1的转动谱线的波长，为从波长391nm至波长388nm；利用计算机的数据处理程序，顺序标记电子束荧光转动谱线的波峰序号1～n为转动量子数j
′
，即j
′
＝1∶n；
21.s30.获得转动量子数j
′
对应的强度值ij′
和波长λj′
；
22.利用计算机的数据处理程序，读出激发态转动量子数j
′
对应的强度值ij′
和波长λj′
；其中转动量子数j
′
为偶数的，强度值ij′
取原读数的2倍；
23.s40.计算转动项g(j
′
，tr)；
24.假设转动温度tr＝200k，令特征温度θ＝2.878，计算转动项g(j
′
，tr)；
[0025][0026]
s50.获取数据点集{y(j
′
)}，{x(j
′
)}；
[0027]
由于：
[0028][0029]
其中，j
″
为基态时转动量子数；λj′j″
为从激发态转动量子数j
′
跃迁至基态时转动量子数j
″
的激发光波长；ij′j″
为波长λj′j″
的激发光强度；c为常数。
[0030]
令：
[0031]
[0032][0033]
得到数据点集{y(j
′
)}，{x(j
′
)}；
[0034]
s60.拟合数据点集{y(j
′
)}，{x(j
′
)}的一次函数；
[0035]
利用步骤s50获得的数据点集{y(j
′
)}，{x(j
′
)}，采用最小二乘法拟合一次函数：
[0036]
y(j
′
)＝a
×
x(j
′
) c
[0037]
则：
[0038][0039]
s70.通过一次函数，获得转动温度tr；
[0040]
取重复步骤s40～步骤s60，迭代计算转动温度tr，直至前后两次转动温度tr的差量δtr≤0.01％；
[0041]
s80.判断转动温度tr的合理性；
[0042]
利用计算机的数据处理程序，计算各转动量子数j
′
的强度值与一次函数y(j
′
)＝a
×
x(j
′
) c的偏离程度，统计偏离程度大于5％的强度值的百分比，如果百分比小于10％，则转动温度tr为最终值；否则，进入步骤s90；
[0043]
利用拟合的温度tr和常数c以及计算得到的偏离程度
[0044][0045]
s90.通过分段的一次函数，计算前段加权转动温度t
rc
和后段加权转动温度t
rh
；
[0046]
取j
′
＝1至5，重复步骤s40～步骤s60，获得前段一次函数拟合直线，计算得到前段加权转动温度t
rc
；取j
′
＝6至n，重复步骤s40～步骤s70，获得后段一次函数拟合直线，计算得到后段加权转动温度t
rh
；
[0047]
s100.计算权重系数获得转动温度tr；
[0048]
取j
′
为2、6，令
[0049]
本发明的基于电子束荧光的稀薄流场转动温度测量与数据处理方法简便可行，适用于室温静态真空环境、高超声速稀薄空流场和激波后流场的流场情况，具有工程应用价值。
附图说明
[0050]
图1为本发明的基于电子束荧光的稀薄流场转动温度测量方法使用的测量装置的示意图；
[0051]
图2为实施例1的电子束荧光转动谱线；
[0052]
图3为实施例1的一次函数拟合直线；
[0053]
图4为实施例2的电子束荧光转动谱线；
[0054]
图5为实施例2的一次函数拟合直线；
[0055]
图6为实施例2的前段一次函数拟合直线；
[0056]
图7为实施例2的后段一次函数拟合直线；
[0057]
图8为实施例3的电子束荧光转动谱线；
[0058]
图9为实施例3的一次函数拟合直线；
[0059]
图10为实施例3的前段一次函数拟合直线；
[0060]
图11为实施例3的后段一次函数拟合直线。
[0061]
图中，1.扩压器；2.试验段；3.电子枪；4.电子束荧光；5.观察窗；6.光谱仪；7.ccd相机；8.网线；9.喷管；10.计算机；11.法拉第杯。
具体实施方式
[0062]
下面结合附图和实施例详细说明本发明。
[0063]
本发明的基于电子束荧光的稀薄流场转动温度测量方法使用的测量装置见图1，包括安装在高超声速低密度风洞的试验段2上驻室的电子枪3，安装在高超声速低密度风洞的试验段2下驻室的法拉第杯11，法拉第杯11用于接收电子枪3发射的电子束，产生的电子束荧光4位于喷管9和扩压器1之间并垂直于喷管9的轴线从上至下穿过试验段2的中心空腔；
[0064]
通过导线连接的光谱仪6和ccd相机7布置在试验段2的观察窗5外，光谱仪6与观察窗5之间还布置有凸透镜；
[0065]
ccd相机7与计算机10之间通过网线8连接；
[0066]
观察窗5采用光学玻璃；光谱仪6的波长分辨率大于等于0.02nm，测量波长范围为388nm～391nm；
[0067]
所述的测量方法包括室温静态真空环境的转动温度测量方法、高超声速稀薄空流场的转动温度测量方法和激波后流场的转动温度测量方法，分别如下：
[0068]
a.室温静态真空环境的转动温度测量方法；
[0069]
先将高超声速低密度风洞的试验段2抽真空至20pa以下，开启电子枪3发射电子束，在试验段2的测量区域产生的电子束荧光4通过观察窗5的光学玻璃被凸透镜收集至光谱仪6的狭缝入口，电子束荧光4在光谱仪6内产生的电子束荧光转动谱线，通过ccd相机7记录电子束荧光转动谱线的波长和强度数据并传输至计算机10存储；
[0070]
b.高超声速稀薄空流场的转动温度测量方法；
[0071]
先将高超声速低密度风洞的试验段2抽真空至20pa以下，开启高超声速低密度风洞吹风，开启电子枪3发射电子束，在试验段2的测量区域产生的电子束荧光4通过观察窗5的光学玻璃被凸透镜收集至光谱仪6的狭缝入口，电子束荧光4在光谱仪6内产生的电子束荧光转动谱线，通过ccd相机7记录电子束荧光转动谱线的波长和强度数据并传输至计算机10存储；
[0072]
c.激波后流场的转动温度测量方法；
[0073]
先将试验模型安装在高超声速低密度风洞的试验段2内，试验段2抽真空至20pa以下，开启高超声速低密度风洞吹风，开启电子枪3发射电子束，调整电子束位置到试验模型的头部激波后方，在试验段2的测量区域产生的电子束荧光4通过观察窗5的光学玻璃被凸
透镜收集至光谱仪6的狭缝入口，电子束荧光4在光谱仪6内产生的电子束荧光转动谱线，通过ccd相机7记录电子束荧光转动谱线的波长和强度数据并传输至计算机10存储。
[0074]
本发明的基于电子束荧光的稀薄流场转动温度测量方法的数据处理方法，包括以下步骤：
[0075]
s10.绘制电子束荧光转动谱线；
[0076]
将计算机10内存储的电子束荧光转动数据绘制成横坐标为波长、纵坐标为转动强度的电子束荧光转动谱线；
[0077]
s20.标记电子束荧光转动谱线的激发态转动量子数j
′
；
[0078]
激发态转动量子数j
′
跃迁至基态时转动量子数j
″
，取j
′‑j″
＝1的转动谱线的波长，为从波长391nm至波长388nm；利用计算机10的数据处理程序，顺序标记电子束荧光转动谱线的波峰序号1～n为转动量子数j
′
，即j
′
＝1∶n；
[0079]
s30.获得转动量子数j
′
对应的强度值ij′
和波长λj′
；
[0080]
利用计算机10的数据处理程序，读出激发态转动量子数j
′
对应的强度值ij′
和波长λj′
；其中转动量子数j
′
为偶数的，强度值ij′
取原读数的2倍；
[0081]
s40.计算转动项g(j
′
，tr)；
[0082]
假设转动温度tr＝200k，令特征温度θ＝2.878，计算转动项g(j
′
，tr)；
[0083][0084]
s50.获取数据点集{y(j
′
)}，{x(j
′
)}；
[0085]
由于：
[0086][0087]
其中，j
″
为基态时转动量子数；λj′j″
为从激发态转动量子数j
′
跃迁至基态时转动量子数j
″
的激发光波长；ij′j″
为波长λj′j″
的激发光强度；c为常数。
[0088]
令：
[0089][0090][0091]
得到数据点集{y(j
′
)}，{x(j
′
)}；
[0092]
s60.拟合数据点集{y(j
′
)}，{x(j
′
)}的一次函数；
[0093]
利用步骤s50获得的数据点集{y(j
′
)}，{x(j
′
)}，采用最小二乘法拟合一次函数：
[0094]
y(j
′
)＝a
×
x(j
′
) c
[0095]
则：
[0096][0097]
s70.通过一次函数，获得转动温度tr；
[0098]
取重复步骤s40～步骤s60，迭代计算转动温度tr，直至前后两次转动温度tr的差量δtr≤0.01％；
[0099]
s80.判断转动温度tr的合理性；
[0100]
利用计算机10的数据处理程序，计算各转动量子数j
′
的强度值与一次函数y(j
′
)＝a
×
x(j
′
) c的偏离程度，统计偏离程度大于5％的强度值的百分比，如果百分比小于10％，则转动温度tr为最终值；否则，进入步骤s90；
[0101]
利用拟合的温度tr和常数c以及计算得到的偏离程度
[0102][0103]
s90.通过分段的一次函数，计算前段加权转动温度t
rc
和后段加权转动温度t
rh
；
[0104]
取j
′
＝1至5，重复步骤s40～步骤s60，获得前段一次函数拟合直线，计算得到前段加权转动温度t
rc
；取j
′
＝6至n，重复步骤s40～步骤s70，获得后段一次函数拟合直线，计算得到后段加权转动温度t
rh
；
[0105]
s100.计算权重系数获得转动温度tr；
[0106]
取j
′
为2、6，令
[0107]
实施例1
[0108]
本实施例用于室温静态真空环境的转动温度测量。
[0109]
获得的电子束荧光转动谱线如图2所示，n＝17，转动量子数j
′
＝1至17枝的强度值可明显分辨，转动温度tr＝300.96k。获得的一次函数拟合直线如图3所示，各转动量子数的点大多在一次函数拟合直线附近，只有j
′
＝2的偏离程度大于5％，转动谱线的两个斜率是一致的，转动温度tr＝300.96k为最终值。
[0110]
这与室温静态真空环境下，氮气分子的两个转动自由度是平衡的事实是一致的，测量当时的室温为25摄氏度，即298.15k，tr的偏差为0.94％，表明本发明的室温静态真空环境的转动温度测量方法合理，测量精度较高。
[0111]
实施例2
[0112]
本实施例用于高超声速稀薄空流场的转动温度测量。
[0113]
获得的电子束荧光转动谱线如图4所示，n＝9，转动量子数j
′
＝1至9枝的强度值可明显分辨，采用一次函数拟合的转动温度tr＝37.34k，一次函数拟合直线如图5所示，有j
′
＝1、2、5、6枝的点偏离直线超过范围，可见采用一次函数拟合并不合适。
[0114]
取j
′
＝1至5，计算前段加权转动温度t
rc
＝24.09k，获得的前段一次函数拟合直线如图6所示，j
′
＝1至5的点均在前段一次函数拟合直线附近，且转动谱线的两个斜率是一致的；
[0115]
取j
′
＝6至9，计算后段加权转动温度t
rh
＝55.67k，获得的后段一次函数拟合直线如图7所示，j
′
＝6至9的点均在后段一次函数拟合直线附近，且转动谱线的两个斜率是一致的；
[0116]
计算权重系数
[0117]
最终得到转动温度
[0118]
实施例3
[0119]
本实施例用于激波后流场的转动温度测量。
[0120]
获得的电子束荧光转动谱线如图8所示，n＝17，转动量子数j
′
＝1至17枝的强度值可明显分辨，采用一次函数拟合的转动温度tr＝221.03k，一次函数拟合直线如图9所示，有j
′
＝1-3、5-10、17枝的点偏离直线超过范围；可见采用一次函数拟合并不合适。
[0121]
取j
′
＝1至5，计算前段加权转动温度t
rc
＝40.98k，获得的前段一次函数拟合直线如图10所示，j
′
＝1至5的点均在前段一次函数拟合直线附近，且转动谱线的两个斜率是一致的；
[0122]
取j
′
＝6至17，计算后段加权转动温度t
rh
＝316.79k，获得的后段一次函数拟合直线如图11所示，j
′
＝6至17的点均在后段一次函数拟合直线附近，且转动谱线的两个斜率是一致的；
[0123]
权重系数
[0124]
最终得到转动温度
[0125]
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的技术领域。对于熟悉本领域的人员而言，在不脱离本发明原理的前提下，可容易地实现另外的改进和润饰，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。