CN112378616B - 基于波长调制吸收光谱的高速流场多参数测量系统及方法 - Google Patents

基于波长调制吸收光谱的高速流场多参数测量系统及方法 Download PDF

Info

Publication number
CN112378616B
CN112378616B CN202011278426.5A CN202011278426A CN112378616B CN 112378616 B CN112378616 B CN 112378616B CN 202011278426 A CN202011278426 A CN 202011278426A CN 112378616 B CN112378616 B CN 112378616B
Authority
CN
China
Prior art keywords
laser
collimator
flow field
water
measuring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011278426.5A
Other languages
English (en)
Other versions
CN112378616A (zh
Inventor
王兴平
金熠
彭冬
李佳胜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN202011278426.5A priority Critical patent/CN112378616B/zh
Publication of CN112378616A publication Critical patent/CN112378616A/zh
Application granted granted Critical
Publication of CN112378616B publication Critical patent/CN112378616B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

本发明涉及基于波长调制吸收光谱的高速流场多参数测量系统及方法,包括信号发生器、激光器、数据采集系统以及测量架,所述测量架及风洞试验舱外部均安装有一一对应的准直器形成三组准直器组;所述信号发生器的输出端分别连接有两激光器,两激光器的输出端均通过分束器连接形成四路激光,其中一路激光通过光纤传输到近红外标准具上用于激光频率标定、另外三路激光通过光纤用于分别测量水气温度、水气分压和水气速度;所述三组准直器组和近红外标准具的输出端分别与四个探测器一一连接,四个探测器的输出端与数据采集系统连接。本发明能够高精度、快速同步反演高速流场中的水气速度、温度和分压,适用于高超声速风洞中的流场质量评估。

Description

基于波长调制吸收光谱的高速流场多参数测量系统及方法
技术领域
本发明涉及可调谐激光气体吸收光谱探测技术领域,具体涉及一种基于波长调制吸收光谱技术的高速流场多参数快速测量系统及方法。
背景技术
随着先进航天飞行器的快速发展,与航天飞行器紧密相关的高超声速气动力学和气动加热问题受到高度重视,而进行这些研究首先要获取高速流场参数,包括气体的温度、分压及流速等。地面风洞模拟是一种效率高、成本低的研究方案,在与气流相关的流体力学领域中,已经得到了广泛的应用。
在风洞试验中,非接触测量技术是揭示流场特性的重要方法,这种技术不仅可以获得有关流动状态的直接现象及流动的发展过程,还能给出流动参数的定量结果。与其他试验方法的不同之处在于,非接触测量使得流场的某些特征可视化并且能获得整个流场中物理现象的信息。在风洞中常见的非接触测量可以分为以下6类:光学显示技术、光学测量技术、光谱诊断技术、微波测量技术、摄影测量及图像处理技术以及计算流动显示技术。
光谱诊断技术包括自发辐射光谱探测技术、激光诱导荧光技术、瑞利散射技术、拉曼光谱技术和可调谐二极管激光器吸收光谱法等。其中可调谐二极管激光器吸收光谱(TDLAS)诊断技术能实现待测流场的组分、浓度和温度信息,测量精度高、速度快,同时器件成熟,成本较低,是进行风洞非接触高速流场参数测量的有效方法之一。
当前利用TDLAS技术测量风洞流场关键参数多采用后期数据处理方法,处理时间长,多达半个月至一个月。也有部分研究人员采用了在线计算方法,但采用的是假设水气分压为常数,只求取流场水气温度一种关键参数,但该常数值通常由CFD模拟估计或经验判断,通常具有较大误差。鉴于此,有必要开发一种可以提高TDLAS技术使用效率、保证高测量精度,同时可实现多参数同步测量的流场关键参数快速测量TDLAS应用技术。
发明内容
本发明的目的在于提供一种基于波长调制吸收光谱技术的高速流场多参数快速测量系统及方法,该系统能够在线快速测量高速流场中的多种关键参数,包括流场中的水气温度、水气分压和水气速度。解决了目前TDLAS应用技术中数据处理时间长,难以实现快速在线计算的问题,同时保证待测流场中关键参数的计算精度。
为实现上述目的,本发明采用了以下技术方案:
基于波长调制吸收光谱的高速流场多参数测量系统,包括信号发生器、激光器、数据采集系统以及架设在风洞试验舱内部的测量架,所述测量架及风洞试验舱外部均安装有一一对应的准直器,并形成三组准直器组;所述信号发生器的输出端分别连接有第一激光器、第二激光器,所述第一激光器、第二激光器的输出端均通过分束器连接形成四路激光,其中一路激光通过光纤传输到近红外标准具上用于激光频率标定、另外三路激光通过光纤分别传输至相对应的一组准直器组上,用于分别测量水气温度、水气分压和水气速度;所述三组准直器组和近红外标准具的输出端分别与四个探测器一一连接,四个探测器的输出端与数据采集系统连接。
进一步方案,所述第一激光器与信号发生器之间连接有第一激光器控制器,所述第二激光器与信号发生器之间连接有第二激光器控制器;所述第一激光器控制器、第二激光器控制器分别控制第一激光器、第二激光器的电流和温度。
进一步方案,所述第一激光器、第二激光器的中心波长分别为1391.67nm、1395.00nm,激光器的调制信号为半锯齿波叠加高频调制信号;所述分束器为2×4光纤耦合分束器,使所述第一激光器、第二激光器中的激光等量输入分束器并分成四等份输出而形成四路激光。
进一步方案,所述测量架包括固定成十字形结构的固定柱、支撑柱,所述支撑柱上水平对称固定有两组肋板,所述肋板上对称安装有准直器。
优选的,所述肋板的中间垂直固定有L形支架,所述准直器的一端连接光纤头、另一端通过螺纹连接在L形支架的内侧壁;位于L形支架的外侧壁安装有保护劈。
进一步方案,所述固定柱、支撑柱、肋板、L型支架的迎风面均呈尖端,以减少激波的产生。
进一步方案,两组肋板中一组肋板为长肋板组、另一组为短肋板组,所述短肋板组位于长肋板组的两个肋板之间;所述长肋板组上对称安装有第一准直器A1、第一准直器A2,两者间形成第一束激光,所述第一束激光与流场运动方向呈90°夹角;
所述短肋板组上分别对称安装有第二准直器B1、第二准直器B2形成第二束激光以及第三准直器C1、第三准直器C2形成第三束激光;所述第二束激光、第三束激光中的一束激光与流场运动方向呈90°夹角、另一束激光与流场运动方向呈60°夹角。
更优选的,所述第一准直器A1、第二准直器B1、第三准直器C1为激光发射端,第一准直器A2、第二准直器B2、第三准直器C2为激光接收端。
本发明的另一个发明目的是提供一种基于波长调制吸收光谱的高速流场多参数测量方法,其包括以下步骤:
(1)根据待测流场的CFD模拟数据,初步确定待测流场中水气温度及水气分压的范围;
(2)根据免标定WMS-2f/1f理论公式
Figure BDA0002779913140000031
cos(2θ)dθ=f(T,Px),式中Px为水气分压,L为激光穿过流场的光程,S(T)为吸收线在温度T下的线强,i0为激光器光强调制系数,v0为激光器心频率,a为调制深度,φ为吸收线型函数,仿真得到中心波长分别为1391.67nm、1395.00nm的2f/1f信号值关于水气温度和水气分压的三维分布图;
(3)采用上述高速流场多参数快速测量系统分别采集通过风洞试验舱的三路激光的透射光光强信号,对其进行快速傅里叶变换,得到各自的2f/1f信号值;
(4)设步骤(3)中一路激光用于测量流场中的水气温度和水气分压,另外两路激光用于测量水气速度;
对于用于测量水气温度和水气分压的一路激光的两个不同波长激光的2f/1f信号值,利用线性插值算法,在步骤(2)中的三维分布图中查找同时满足前述两个2f/1f信号值的共解,该共解即为流场中待测的水气温度与水气分压;
对于用于测量水气速度的两路激光,利用吸收强度较大的、中心波长为1395.00nm的水气吸收线的两个2f/1f信号值的峰值频率差Δv,再按照式
Figure BDA0002779913140000032
进行流场水气速度计算,其中U为水气速度、Δv为两个2f/1f信号值的峰值频率差、为光速、v0为激光中心频率、β为用于测量水气速度的两路激光的角度差。
进一步方案,激光器光强调制系数i0、调制深度a均是通过标准具测量获得的,其中调制深度a的公式如下:
Figure BDA0002779913140000033
式中v0为激光中心频率,T为待测流场的估计温度。
本申请采用采用时分复用形式,两个激光器分别输出中心波长为1395.00nm及1391.67nm的波长调制扫描信号。两束激光通过2×4光纤耦合分束器,分为4路激光,每一路激光均含有25%的1391.67nm激光器输出的激光和25%的1395.00nm激光器输出的激光。其中第四路激光通过光纤传输到近红外标准具,进行激光频率标定,获取激光相对频率变化,包括激光频率扫描范围与调制深度a。第一路激光通过光纤传输至风洞试验舱内的测量架上的一组长肋板组上的准直器,接受激光,并通过光纤传输出风洞试验舱外,通过光电探测器接收激光信号,并将信号传输至数据采集系统,最后通过上位机计算流场参数,包括水气温度及水气分压。该路激光在风洞试验舱内,与流场运动方向呈90°。当激光方向与流场运动方向呈90°时,不存在流场运动引起的多普勒频移效应,即激光中心频率位置不会改变。
第二路和第三路激光通过光纤传输至测量架上的一组短肋板组上的准直器,接受激光并通过光纤传输出风洞试验舱外,通过光电探测器接收激光信号,并将信号传输至数据采集系统,最后通过上位机计算流场参数,即流场运动速度。该两路激光在风洞试验舱内,与流场运动方向分别呈90°与60°。当激光方向与流场运动方向呈60°时,因为激光在流场运动方向存在分量,因此存在流场运动引起的多普勒频移效应,即激光中心频率位置会改变。
其与90°激光的中心频率的位置差Δv,按下式计算可以得到流场速度信息。
Figure BDA0002779913140000041
其中c为光速,v0为激光中心频率,β为90°与60°的角度差30°。
其中激光方向与流场运动方向呈60°为一个较合适的角度,实验过程发现,角度大于60°时,观察到的激光的中心频率的位置差Δv过小,信噪比不高;而小于60°时,虽然观察到的激光的中心频率的位置差Δv更加明显,信噪比更高,但用于高速流场多参数测量的测量架机械结构会占用更大的实验空间。但如果实验现场条件允许,则可以考虑将角度由60°改为45°,甚至30°,以获得更准确的流场速度。
对采集得到的第一路激光透射光光强信号,利用快速傅里叶变换,得到2f/1f信号。对于1395及1392激光器,可以测得1395激光器及1392激光器的2f/1f信号。以测得的2f/1f信号为z坐标值,与数据库三维分布图中相交得到一条交线,该交线即为所选吸收线关于温度T及分压Px一组通解,两个激光器的通解中相同的解即为与待测流场参数一致的特解。具体在数据操作中,表现为2f/1f信号测得值与仿真三维分布图相交得到一条线,两个激光器的两条相交线在温度T及分压Px平面上投影的交点即为待测流场中的水气温度及水气分压。
对采集得到的第二路、第三路激光透射光光强信号,进行快速傅里叶变换,得到1395.00nm吸收线的2f/1f值。根据第二路与第三路2f/1f峰值信号的频率差Δv,按式
Figure BDA0002779913140000051
进行流场速度计算。
所以本发明能够实现高速流场中多参数快速测量,包括定量测量高速流场中的水气速度、水气温度和水气分压,可以为评估风洞流场状态提供关键参数依据,推动航空航天领域中风洞地面模拟技术发展。
附图说明
图1是本发明系统原理图,
图2是本发明中的测量架结构示意图,
图3是本发明测量流程示意图,
图4是本发明中的三维分布图,
图5是本发明中的第一路激光得到的激光透射信号,
图6是本发明中的第一路激光透射信号的2f/1f信号,
图7是本发明中的第一路激光2f/1f信号两解集在水气温度-水气分压上的投影,
图8是本发明中的第二路、第三路激光透射信号,
图9是本发明中的第二路、第三路激光透射信号的2f/1f信号,
图1、2中附图标记为:
1信号发生器,2-1第一激光器控制器,2-2第二激光器控制器,3-1第一激光器,3-2第二激光器,4分束器,5风洞试验舱,6标准具,7-1第一探测器,7-2第二探测器,7-3第三探测器,7-4第四探测器;8-1固定柱,8-2-支撑柱,8-3固定柱,8-4肋板,8-5L形支架,8-6保护劈;9-1第一束激光,9-2第二束激光,9-3第三束激光;10数据采集系统。
具体实施方式
实施例1:
如图1所示,基于波长调制吸收光谱的高速流场多参数测量系统,包括信号发生器1、激光器、数据采集系统10以及架设在风洞试验舱5内部的测量架,所述测量架及风洞试验舱5外部均安装有一一对应的准直器,并形成三组准直器组;所述信号发生器1的输出端分别连接有第一激光器3-1、第二激光器3-2,所述第一激光器3-1、第二激光器3-2的输出端均通过分束器4连接形成四路激光,其中一路激光通过光纤传输到近红外标准具6上用于激光频率标定、另外三路激光通过光纤分别传输至相对应的一组准直器组上,用于分别测量水气温度、水气分压和水气速度;所述三组准直器组和近红外标准具6的输出端分别与四个探测器一一连接,四个探测器的输出端与数据采集系统10连接。
进一步方案,所述第一激光器3-1与信号发生器1之间连接有第一激光器控制器2-1,所述第二激光器3-2与信号发生器1之间连接有第二激光器控制器2-2;所述第一激光器控制器2-1、第二激光器控制器2-2分别控制第一激光器3-1、第二激光器3-2的电流和温度。
进一步方案,所述第一激光器3-1、第二激光器3-2的中心波长分别为1391.67nm、1395.00nm,激光器的调制信号为半锯齿波叠加高频调制信号;所述分束器4为2×4光纤耦合分束器,使所述第一激光器3-1、第二激光器3-2中的激光等量输入分束器4并分成四等份输出而形成四路激光。
即本申请通过信号发生器产生低频锯齿波叠加高频正弦波的调制信号,输入到激光器控制器(美国ILX Lightwave,LDC-3724C)中,激光器控制器控制激光器(日本NEL,NLK1E5GAAA,中心波长1395nm、1392nm)的电流及温度(温度与电流的设置由激光器性能决定,在商用激光器的测试数据页中查询确定),使得激光器输出中心波长为所选谱线的激光。采用时分复用形式,两个激光器分别输出中心波长为1395.00nm及1391.67nm的波长调制扫描信号。两束激光通过2×4光纤耦合分束器,等分为4路激光,每一路激光均含有25%的1391.67nm激光器输出的激光和25%的1395.00nm激光器输出的激光。
如图2所示,测量架包括固定成十字形结构的固定柱8-1、支撑柱8-2,所述支撑柱8-2上水平对称固定有两组肋板8-4,所述肋板8-4上对称安装有准直器。
优选的,所述肋板8-4的中间垂直固定有L形支架8-5,所述L形支架8-5的水平段固接在肋板8-4的底端、水平段用于安装准直器(图中末画出,位于图中第一束激光9-1、第二束激光9-2、第三束激光9-3所指之处);所述准直器的一端连接光纤头、另一端通过螺纹连接在L形支架8-5的内侧壁;位于L形支架8-5的外侧壁安装有保护劈8-6。
进一步方案,所述固定柱8-1、支撑柱8-2、肋板8-4、L型支架8-5的迎风面均呈尖端,以减少激波的产生。
进一步方案,两组肋板8-4中一组肋板为长肋板组、另一组为短肋板组,所述短肋板组位于长肋板组的两个肋板之间;所述长肋板组上对称安装有第一准直器A1、第一准直器A2,两者间形成第一束激光9-1,所述第一束激光9-1与流场运动方向呈90°夹角;
所述短肋板组上分别对称安装有第二准直器B1、第二准直器B2形成第二束激光9-2以及第三准直器C1、第三准直器C2形成第三束激光9-3;所述第二束激光9-2、第三束激光9-3中的一束激光与流场运动方向呈90°夹角、另一束激光与流场运动方向呈60°夹角。
更优选的,所述第一准直器A1、第二准直器B1、第三准直器C1为激光发射端,第一准直器A2、第二准直器B2、第三准直器C2为激光接收端。
如图1所示,第一路激光经第一准直器A1、第一准直器A2、第一准直器A3,然后由第一探测器7-1接收激光信号;第二路激光经第二准直器B1、第二准直器B2、第二准直器B3,然后由第二探测器7-2接收激光信号;第三路激光经第三准直器C1、第三准直器C2、第三准直器C3,然后由第三探测器7-3接收激光信号;第四路激光经第四准直器D、远红外标准具后由第四探测器7-4接收激光信号。
实施例2:
如图3所示,一种基于波长调制吸收光谱的高速流场多参数测量方法,其包括以下步骤:
(1)根据待测流场的CFD模拟数据,初步确定待测流场中水气温度及水气分压的范围;
(2)根据免标定WMS-2f/1f理论公式
Figure BDA0002779913140000071
cos(2θ)dθ=f(T,Px),式中Px为水气分压,L为激光穿过流场的光程,S(T)为吸收线在温度T下的线强,i0为激光器光强调制系数,v0为激光器心频率,a为调制深度,φ为吸收线型函数,仿真得到中心波长分别为1391.67nm、1395.00nm的2f/1f信号值关于水气温度和水气分压的三维分布图;
其中激光器光强调制系数i0、调制深度a均是通过标准具测量获得的,其中调制深度a的公式如下:
Figure BDA0002779913140000072
式中v0为激光中心频率,T为待测流场的估计温度。
免标定WMS-2f/1f的理论公式为弱吸收假设下,根据Beer-Lambert定律,按二次谐波比一次谐波推到而得,属于已知内容。线强S(T)是关于温度T为变量的函数,可以通过公开的HITRAN高光谱数据库(http://www.hitran.org/)查得。线型函数φ是以温度T及分压Px为变量的函数,其理论公式已在许多文献上公开发表(如Goldenstein,Christopher(2014).Wavelength-modulation spectroscopy for determination of gas propertiesin hostile environments.)。其中激光器光强调制系数i0及调制深度a通过标准具测量获得。上式表明,2f/1f信号是关于温度T及分压Px为变量的函数,给定一定范围的温度T及分压Px,通过仿真,可以得到2f/1f信号关于温度T及分压Px的三维分布图。
(3)采用上述高速流场多参数快速测量系统分别采集通过风洞试验舱5的三路激光的透射光光强信号,对其进行快速傅里叶变换,得到各自的2f/1f信号值;
(4)设步骤(3)中一路激光用于测量流场中的水气温度和水气分压,另外两路激光用于测量水气速度;
对于用于测量水气温度和水气分压的第一路激光的两个不同波长激光的2f/1f信号值,利用线性插值算法,在步骤(2)中的三维分布图中查找同时满足前述两个2f/1f信号值的共解,该共解即为流场中待测的水气温度与水气分压;即以测得的2f/1f信号为z坐标值,与三维分布图中相交得到一条交线,该交线即为所选吸收线关于温度T及分压Px一组通解,两个激光器的通解中相同的解即为与待测流场参数一致的特解。具体在数据操作中,表现为2f/1f信号测得值与仿真三维分布图相交得到一条线,两个激光器的两条相交线在温度T及分压Px平面上投影的交点即为待测流场中的水气温度及水气分压;
对于用于测量水气速度的第二路激光、第三路激光,利用吸收强度较大的、中心频率为1395.00nm的水气吸收线的两个2f/1f信号值的峰值频率差Δv,再按照式
Figure BDA0002779913140000081
进行流场水气速度计算,其中U为水气速度、Δv为两个2f/1f信号值的峰值频率差、为光速、v0为激光中心频率、β为用于测量水气速度的两路激光的角度差。
实施例3:
针对某次高超声速风洞试验,如水气温度220K,水气分压5Pa,流场速度1500m/s,进行本发明的实施例详细说明。
在公开的网站(http://www.hitran.org/)HITRAN高光谱数据库中,查找1391.67nm及1395.00nm的水吸收线光谱参数,见下表。
Figure BDA0002779913140000091
调制信号:锯齿波频率、电压、相位、占空比、余弦波频率、幅值在信号发生器上设置,其中锯齿波频率的设置由风洞流场中待测参数的时间分辨率要求决定;锯齿波电压设置由激光器决定,其电压决定了扫描范围,通常扫描范围要求大于1cm-1;占空比与相位的设置可以实现两个激光器时分复用;余弦波频率参考锯齿波频率,通常大于100倍的锯齿波频率;幅值决定了调制深度,在标准具进行频率标定过程中可以观察到干涉条纹的个数来计算实际调制深度;调制深度的设置参考计算得到的理论最优的调制深度设置
Figure BDA0002779913140000092
其中系统的参数设置如下表:
Figure BDA0002779913140000093
利用下式确定最优的调制深度a,在实验中可以用过调节正弦信号的幅值改变调制深度到这一合适值。
Figure BDA0002779913140000094
由上式计算得到两条合适的吸收谱线调制深度分别为:a1395=0.02cm-1,a1392=0.02cm-1
搭建如图1所示的系统,其中信号发生器1与激光器控制器之间通过BNC线连接;激光器控制器与激光器通过RS232接口连接;激光器带尾纤输出,与2×4分束器通过光纤法兰连接;分束器其中三个输出端通过单模光纤接入风洞试验舱5中,分别与通过螺纹连接在测量架上的一个准直器(第一准直器A1、第二准直器B1、第三准直器C1)连接;相对侧一端分别通过第一准直器A2、第二准直器B2、第三准直器C2接收,并通过多模光纤接出风洞试验舱5,并通过第一准直器A3、第二准直器B3、第三准直器C3形成空间光,打入光电探测器(第一探测器7-1、第二探测器7-2、第三探测器7-3),形成投射信号。第一路激光通过光纤传输至测量架,用于测量水气温度和水气分压,第二路和第三路激光通过光纤传输至测量架,用于测量水气速度。
分束器剩下一路接第四准直器D,经过远红标准具6,再打入第四探测器7-4,形成标准具信号。四路信号通过16位数据采集系统采集得到数字信号。
风洞十字形测量架的布置如图2所示,第二路激光通过光纤传输至测量架,用于测量水气温度和水气分压,第三路和第四路激光通过光纤传输至测量架,用于测量水气速度。
即在风洞试验前,利用远红外标准具,通过最小二乘拟合得到激光光强调制系数i0、调制深度a和激光频率扫描范围。然后通过进入风洞试验舱内的三路激光分别测量水气温度、水气分压和测量水气速度。
针对待测风洞中流场的估计参数,设置数据库参数范围(如水气温度100-300K,水气分压0.1-10Pa),结合光谱参数与激光器调制参数(光强调制系数i0与调制深度a1395=0.02cm-1,a1392=0.02cm-1),利用2f/1f理论公式
Figure BDA0002779913140000101
进行数据库的建立。得到如图4所示的数据库三维分布图,其中R1395表示中心频率为1395.00nm的水气吸收线2f/1f值,R1392表示中心频率为1391.67nm的水气吸收线2f/1f值。
风洞试验时,第一路激光得到的激光透射信号如图5所示,对该透射信号进行快速傅里叶变换,获取2f/1f信号,如图6所示。获得两个吸收峰的峰值信号,分别在图4中的数据库三维分布图中获得与实验测得的峰值信号一致的水气温度与水气分压解集,即图4中显示的曲线对应的水气温度与水气分压。
即如图7所示的是图4中的曲线在水气温度-水气分压坐标面上的投影,两组解集在水气温度与水气分压做表面上的投影显示为两条曲线,其两线交点为两组解集的共解,即测得流场中水气温度为232K,水气分压为4.6Pa。
第二路与第三路激光得到的激光透射信号如图8所示,对第二路与第三路激光透射信号进行快速傅里叶变换得到图9所示的2f/1f信号,其中两路信号的峰值差为Δv=0.019cm-1。利用下式进行流场速度计算:
Figure BDA0002779913140000111
即测得待测流场的水气温度为232K,水气分压为4.6Pa、水气速度为1586.5m/s。
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。

Claims (8)

1.基于波长调制吸收光谱的高速流场多参数测量系统,包括信号发生器(1)、激光器、数据采集系统(10)以及架设在风洞试验舱(5)内部的测量架,其特征在于:所述测量架及风洞试验舱(5)外部均安装有一一对应的准直器,并形成三组准直器组;所述信号发生器(1)的输出端分别连接有第一激光器(3-1)、第二激光器(3-2),所述第一激光器(3-1)、第二激光器(3-2)的输出端均通过分束器(4)连接形成四路激光,其中一路激光通过光纤传输到近红外标准具(6)上用于激光频率标定、另外三路激光通过光纤分别传输至相对应的一组准直器组上,用于分别测量水气温度、水气分压和水气速度;所述三组准直器组和近红外标准具(6)的输出端分别与四个探测器一一连接,四个探测器的输出端与数据采集系统(10)连接;
所述测量架包括固定成十字形结构的固定柱(8-1)、支撑柱(8-2),所述支撑柱(8-2)上水平对称固定有两组肋板(8-4),所述肋板(8-4)上对称安装有准直器;两组肋板(8-4)中一组肋板为长肋板组、另一组为短肋板组,所述短肋板组位于长肋板组的两个肋板之间;
所述长肋板组上对称安装有第一准直器A1、第一准直器A2,两者间形成第一束激光(9-1),所述第一束激光(9-1)与流场运动方向呈90°夹角;
所述短肋板组上分别对称安装有第二准直器B1、第三准直器C1、第二准直器B2、第三准直器C2,其中第二准直器B1、第二准直器B2之间形成第二束激光(9-2),第三准直器C1、第三准直器C2之间形成第三束激光(9-3);所述第二束激光(9-2)、第三束激光(9-3)中的一束激光与流场运动方向呈90°夹角、另一束激光与流场运动方向呈60°夹角。
2.根据权利要求1所述的高速流场多参数测量系统,其特征在于:所述第一激光器(3-1)与信号发生器(1)之间连接有第一激光器控制器(2-1),所述第二激光器(3-2)与信号发生器(1)之间连接有第二激光器控制器(2-2);所述第一激光器控制器(2-1)、第二激光器控制器(2-2)分别控制第一激光器(3-1)、第二激光器(3-2)的电流和温度。
3.根据权利要求1所述的高速流场多参数测量系统,其特征在于:所述第一激光器(3-1)、第二激光器(3-2)的中心波长分别为1391.67nm、1395.00nm,激光器的调制信号为半锯齿波叠加高频调制信号;所述分束器(4)为2×4光纤耦合分束器,使所述第一激光器(3-1)、第二激光器(3-2)中的激光等量输入分束器(4)并分成四等份输出而形成四路激光。
4.根据权利要求1所述的高速流场多参数测量系统,其特征在于:所述肋板(8-4)的中间垂直固定有L形支架(8-5),所述准直器的一端连接光纤头、另一端通过螺纹连接在L形支架(8-5)的内侧壁;位于L形支架(8-5)的外侧壁安装有保护劈(8-6)。
5.根据权利要求1所述的高速流场多参数测量系统,其特征在于:所述固定柱(8-1)、支撑柱(8-2)、肋板(8-4)、L型支架(8-5)的迎风面均呈尖端,以减少激波的产生。
6.根据权利要求1所述的高速流场多参数测量系统,其特征在于:所述第一准直器A1、第二准直器B1、第三准直器C1为激光发射端,第一准直器A2、第二准直器B2、第三准直器C2为激光接收端。
7.基于波长调制吸收光谱的高速流场多参数测量方法,其特征在于:包括以下步骤:
(1)根据待测流场的CFD模拟数据,初步确定待测流场中水气温度及水气分压的范围;
(2)根据免标定WMS-2f/1f理论公式
Figure FDA0003259957380000021
式中Px为水气分压,L为激光穿过流场的光程,S(T)为吸收线在温度T下的线强,i0为激光器光强调制系数,v0为激光器心频率,a为调制深度,φ为吸收线型函数,仿真得到中心波长分别为1391.67nm、1395.00nm的2f/1f信号值关于水气温度和水气分压的三维分布图;
(3)采用如权利要求1-6任一项所述的高速流场多参数快速测量系统分别采集通过风洞试验舱的三路激光的透射光光强信号,对其进行快速傅里叶变换,得到各自的2f/1f信号值;
(4)设步骤(3)中一路激光用于测量流场中的水气温度和水气分压,另外两路激光用于测量水气速度;
对于用于测量水气温度和水气分压的一路激光的两个不同波长激光的2f/1f信号值,利用线性插值算法,在步骤(2)中的三维分布图中查找同时满足前述两个2f/1f信号值的共解,该共解即为流场中待测的水气温度与水气分压;
对于用于测量水气速度的两路激光,利用吸收强度较大的、中心波长为1395.00nm的水气吸收线的两个2f/1f信号值的峰值频率差Δv,再按照式
Figure FDA0003259957380000022
进行流场水气速度计算,其中U为水气速度、Δv为两个2f/1f信号值的峰值频率差、c为光速、v0为激光中心频率、β为用于测量水气速度的两路激光的角度差。
8.根据权利要求7所述的基于波长调制吸收光谱的高速流场多参数测量方法,其特征在于:激光器光强调制系数i0、调制深度a均是通过标准具测量获得的,其中调制深度a的公式如下:
Figure FDA0003259957380000023
(cm-1),式中v0为激光中心频率,T为待测流场的估计温度。
CN202011278426.5A 2020-11-16 2020-11-16 基于波长调制吸收光谱的高速流场多参数测量系统及方法 Active CN112378616B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011278426.5A CN112378616B (zh) 2020-11-16 2020-11-16 基于波长调制吸收光谱的高速流场多参数测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011278426.5A CN112378616B (zh) 2020-11-16 2020-11-16 基于波长调制吸收光谱的高速流场多参数测量系统及方法

Publications (2)

Publication Number Publication Date
CN112378616A CN112378616A (zh) 2021-02-19
CN112378616B true CN112378616B (zh) 2021-12-10

Family

ID=74585383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011278426.5A Active CN112378616B (zh) 2020-11-16 2020-11-16 基于波长调制吸收光谱的高速流场多参数测量系统及方法

Country Status (1)

Country Link
CN (1) CN112378616B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112729417B (zh) * 2021-03-29 2021-06-22 湖南久钰电子有限公司 基于感知数据融合的多参数智能监测系统及方法
CN113654757B (zh) * 2021-08-27 2022-10-28 中国科学技术大学 一种模拟喷管中高超声速凝结过程的装置及诊断方法
CN114184348B (zh) * 2021-10-26 2022-07-05 中国航天空气动力技术研究院 一种高焓流场光电特性辨识装置及方法
CN114518230B (zh) * 2022-03-01 2022-10-28 上海交通大学 一种发动机羽流场速度温度同步测量系统
CN115683228B (zh) * 2022-11-16 2023-08-11 西安交通大学 一种高温热管内部气相工质多参数测量系统及工作方法
CN116718344B (zh) * 2023-08-10 2023-10-27 中国空气动力研究与发展中心高速空气动力研究所 一种多参数的推力矢量喷流光学标定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108008146A (zh) * 2017-10-31 2018-05-08 天津工业大学 一种提高气体流速测量精度的方法和装置
CN111795829A (zh) * 2020-07-07 2020-10-20 西安热工研究院有限公司 一种基于tdlas技术燃气轮机燃烧状态参数监测系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102135547B (zh) * 2010-12-28 2012-07-11 中国科学院安徽光学精密机械研究所 基于可调谐半导体激光吸收光谱技术的气体流速测量仪
JP6057430B2 (ja) * 2013-08-21 2017-01-11 国立大学法人徳島大学 レーザ光を用いたガス分析装置及びガス分析方法
CN104749396A (zh) * 2015-03-25 2015-07-01 中国科学院安徽光学精密机械研究所 一种基于波长调制-tdlas的高速气流流速检测系统
CN106017725B (zh) * 2016-05-26 2019-07-09 中国人民解放军战略支援部队航天工程大学 一种适用于燃烧流场气体二维重建的测量装置
CN111474138B (zh) * 2020-04-20 2023-03-28 东南大学 一种基于高频参考光频分复用技术的气体浓度测量装置及测量方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108008146A (zh) * 2017-10-31 2018-05-08 天津工业大学 一种提高气体流速测量精度的方法和装置
CN111795829A (zh) * 2020-07-07 2020-10-20 西安热工研究院有限公司 一种基于tdlas技术燃气轮机燃烧状态参数监测系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
利用波长调制光谱技术测量高温环境中的气体压强、温度和H2O组分浓度;屈东胜;《光谱学与光谱分析》;20170531;第37卷(第5期);第1339-1344页 *

Also Published As

Publication number Publication date
CN112378616A (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
CN112378616B (zh) 基于波长调制吸收光谱的高速流场多参数测量系统及方法
CN108981953B (zh) 一种基于干涉调制原理的激光吸收光谱测温方法和系统
CN101813623B (zh) 散射物质多维光谱测量装置与方法
CN107872274B (zh) 一种光纤色散系数的测量方法
CN103557886B (zh) 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时测量的方法
CN113447458B (zh) 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法
CN105891434B (zh) 一种海水盐度的在线检测方法及其装置
CN111077110B (zh) 一种基于双光梳光谱的温度场和浓度场测量系统和方法
CN108760681A (zh) 一种基于波形分解的路径平均温度测量系统与方法
CN111077109A (zh) 一种基于双光梳光谱技术的温度和浓度测量系统和方法
CN103048050B (zh) 一种测量目标真实温度的辐射测温方法和仪器
CN107894327B (zh) 一种光纤色散系数的测量装置
CN104864979A (zh) 一种分布式拉曼光纤测温系统测量误差的修正方法
CN204556093U (zh) 一种低噪声微悬臂梁热振动信号测量装置
CN110207733B (zh) 基于扫频激光的光纤干涉仪臂长差测量装置及方法
CN108020527A (zh) 一种用于火焰中气体浓度及温度的检测方法和装置
CN103364371A (zh) 同轴式光热干涉的大气气溶胶吸收系数差分测量新方法
CN108801496B (zh) 一种基于交叠吸收谱的路径温度直方图测量系统与方法
CN112729591B (zh) 一种基于低带宽光电探测器的激光吸收光谱温度测量方法
CN104655029B (zh) 一种位相增强型薄膜厚度测量方法和系统
CN110361359B (zh) 一种基于偶次谐波的吸收率函数复现方法
CN104535199A (zh) 一种相干测量太赫兹波频率的方法
CN108507679A (zh) 一种宽谱段高分辨光谱干涉测量方法及装置
CN114548172A (zh) 一种波长调制激光气体吸收光谱的参数提取方法与系统
CN208012761U (zh) 一种宽谱段高分辨光谱干涉测量装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant