基于TDLAS的固體推進劑裝藥羽流流速測量方法

姚 路，姚德龍，儀建華，闞瑞峰，楊燕京，許振宇，阮 俊，劉建國

(1. 中國科學院環境光學與技術重點實驗室，安徽光學精密機械研究所，安徽 合肥 230031；2. 中國兵器工業試驗測試研究院，陜西 華陰 714200；3.西安近代化學研究所燃燒與爆炸技術重點實驗室，陜西 西安 710065)

?

基于TDLAS的固體推進劑裝藥羽流流速測量方法

姚 路1，姚德龍2，儀建華3，闞瑞峰1，楊燕京3，許振宇1，阮 俊1，劉建國1

(1. 中國科學院環境光學與技術重點實驗室，安徽光學精密機械研究所，安徽 合肥 230031；2. 中國兵器工業試驗測試研究院，陜西 華陰 714200；3.西安近代化學研究所燃燒與爆炸技術重點實驗室，陜西 西安 710065)

介紹了一種基于可調諧半導體激光吸收光譜(TDLAS)的固體推進劑裝藥羽流流速測量新方法。選用推進劑燃燒產物H2O分子在1392nm處的吸收譜線，通過頻率標定、交叉布置光路，利用多普勒頻移效應可反演流速。在靜態試驗中對比了羽流流速特性曲線和發動機內壓強曲線，驗證了該技術用于固體推進劑裝藥羽流流速測量的可行性。研制了緊湊型測量樣機，采用鋰電池供電，響應時間30ms，質量2kg，功耗小于5W，在線反演流速并存儲結果。進行隨火箭橇飛行的動態測試，獲得了有效的羽流流速結果。結果表明，TDLAS技術可用于測量隨發動機飛行的動態羽流流速，是獲取推進劑裝藥實時在線燃燒狀態參數的有效途徑。

可調諧半導體激光吸收光譜；TDLAS；固體推進劑裝藥；羽流流速

引 言

固體推進劑裝藥在發動機燃燒室內燃燒生成的產物通過噴管排出，形成排氣羽流。排氣羽流場的特征與發動機及推進劑的性能密切相關，作為流場的基本物理量，研究流速及其分布對于深入認識流場的特性至關重要，可為發動機的結構設計和推進劑配方的優化提供參考依據[1]。

現有的推進劑排氣羽流流速測量技術可分為接觸式和非接觸式。由于發動機燃燒羽流高溫、高速等特點，測速管、熱線風速儀等傳統接觸式測量雖然結構可靠，但響應時間慢、空間分辨率低，且易于干擾流場，難以獲取真實的流場信息[2]。非接觸式測量主要以光學手段為主，如粒子圖像測速法(PIV)[3]、相位多普勒粒子分析儀(PDPA)[4]等可實現快速、高空間分辨的流場測量，但通常需要加入示蹤粒子，由于燃燒瞬間的不可控性，難以使示蹤粒子均勻分布于流場中，且粒子信號易被較強的火焰輻射信號湮沒，測量系統的復雜性也限制了其應用。可調諧半導體激光吸收光譜技術(TDLAS)通過測量燃燒產物中某種分子的吸收光譜特征，根據多普勒效應建立的光譜頻移和分子速度之間的關系來獲得氣流流速，具有非入侵、快速響應、高光譜分辨、結構緊湊等特點，并可以實現溫度、組分濃度、速度等多參數同時在線測量[5]。

目前，TDLAS技術已經在風洞、發動機進氣道、燃燒室出口流速測量方面得到了國內外研究機構的重視。斯坦福大學高溫氣體動力學實驗室 Hanson等[6-7]在TDLAS 技術的流場測量方面處于領先地位， 目前已經將 TDLAS 系統成功地應用到超燃沖壓發動機飛行實驗中,實現了進氣道質量流量捕獲測量和尾噴口流速、組分測量。Kent[8]將氣體流速和濃度測量相結合，實現了氧氣質量捕獲測量。Shawn[9]將TDLAS技術用于高超聲速氣流流速測量，利用H2O吸收線測得大于10馬赫的高速氣流。美國空軍WPAFB實驗室和Southwest Sciences，Zolo Technologies公司合作，將該技術用于高超聲速國際飛行研究實驗項目(HIFiRE)，并開發出多參數同時測量的TDLAS小型化便攜系統，在緊湊的封裝內可在線測量組分、溫度、流速等參數[10]。賈良權等[11]以760nm處的O2吸收譜線作為測量對象，實現高速流場的質量流量捕獲測量。李飛等[12]使用TDLAS技術研究了超燃沖壓發動機燃燒室入口和燃燒室內的流速、溫度、濃度等參數的測量。

取決其視線測量的原理，TDLAS技術需要光束傳輸路徑上無明顯干擾，如大顆粒遮擋光路、大折射率梯度引起光路偏折等因素，都會給測量結果帶來影響，故多用于純氣體測量環境。而固體推進劑裝藥燃燒產物中通常包含粉塵、顆粒等擋光因素。本研究首次探索了TDLAS技術在該場合的應用，在固體推進劑裝藥發動機出口處布置交叉光路，利用直接吸收技術，選取燃燒產物之一的H2O分子位于1392nm處的單根吸收譜線，進行了靜態和動態測量試驗，驗證了該技術用于固體推進劑裝藥羽流流速測量的可行性。

1 實 驗

1.1 測量原理

1.1.1 TDLAS直接吸收測量原理

當強度為I0、頻率為υ的光束通過一段氣體介質時，吸收會使光強發生衰減，根據Beer-Lambert定律，對于單成分待測氣體分子，定義吸光度A(υ)為:

(1)

=S(T)pXL

(2)

根據式(2)，積分后的吸光度不依賴于線型函數，根據積分面積即可反演出氣體的壓強p、濃度X或與溫度相關的S(T)，簡化了算法復雜度，提高了對測量環境的普適性。直接吸收原理如圖1所示[13]。

圖1 直接吸收技術示意圖 Fig.1 Schematic diagram of direct absorption technique

1.1.2 流速測量原理

根據多普勒頻移原理，當氣流流速在激光傳輸方向上存在速度分量時，運動的氣體分子接收到光子表觀頻率與實際激光發生頻率υ0之間有一定偏移，二者頻率之差即為多普勒頻移Δυ。通常采用雙光路測量方法，見圖2，當光路對稱布置時，兩條光路上探測器接收到的激光頻移大小通過式(3)計算

(3)

式中：υ0為吸收譜線中心頻率；u為氣體流速；c為光速，流速計算中取3×108m/s；θ為光束和流速垂直方向的夾角。

實際測量中，對原始光強信號做基線擬合扣除背景，得到兩條光路上的吸光度曲線，通過頻率定標和峰值擬合后，即可得到如圖3所示的譜線頻差。

圖2 雙光路流速測量示意圖Fig.2 Schematic diagram of double optical paths for flow velocity measuring

圖3 雙光束頻差示意圖Fig.3 Schematic diagram of double optical beam frequency difference

1.2 試驗裝置

試驗中使用的雙基和改性雙基推進劑的燃燒產物中均有水分子存在，故選取H2O分子在1392nm附近的吸收譜線用于流速測量。

TDLAS流速系統分為光學收發和測量控制兩部分，前者包括抗振動的光機收發結構，在燃燒區域實現雙光束的準直發射和聚焦接收，后者實現電子學控制功能，如激光器控制、探測器信號處理、數據采集和算法處理等，系統構成如圖4所示。光束夾角90°，光束交叉點距離發動機出口30cm，為減少燃燒區域兩側靜態空氣中的H2O吸收干擾，在不干擾流場的前提下光程需盡可能縮短，參考對燃燒羽流寬度的經驗值，光束沿流速垂直方向寬度為50cm。

試驗中選用NanoPlus DFB(Distributed Feed-Back)激光器，設置工作溫度為30℃，中心出光波長1391.7nm，調節激光器注入電流范圍，使其波長掃描范圍約為0.3nm，掃描重復頻率為1kHz。掃描波形為占空比99%的鋸齒，激光器平均出光功率為3mW。經由分束器分束后，分別經過自主設計的準直器出射。接收端利用感光面直徑3mm的InGaAs探測器接收，前級I-V轉換增益103V/A，帶寬2MHz，主控芯片選用MCU STM32F407，Cortex M4F內核，最高主頻168MHz，集成硬件浮點運算單元，利用片上ADC完成光譜數據的轉換采集，采樣率設置為800kHz，采樣精度12bits。每個掃描周期采集點數為800，取16個周期數據累加平均，理論可提高2位有效位數。在線進行光譜數據處理與計算，響應時間30ms，得到的流速結果存儲在板載Flash中，試驗結束后經由主控芯片將數據讀出。小型化系統使用四節18650鋰電池供電，功耗小于5W，采用割線方式觸發采集。

圖4 雙光束直接吸收法測速示意圖Fig.4 Schematic diagram of velocity measuring with double optical beam and direct absorption technique

2 結果與討論

2.1 標定結果及雙光路光譜信號

在反演流速時，需要得到兩條光路上吸收峰值在頻率域的偏移量，而由數據轉換采集到的波形數據為時域信號，需要對激光器在其出光范圍內進行波長標定，獲取每個采樣點上的相對頻率，從而進行流速反演。波長標定通常利用干涉計實現，獲取的原始干涉信號見圖5，其中相鄰干涉峰值之間的波長變化為干涉計的自由光譜范圍(FSR)，通過對干涉信號做背景扣除、尋峰、擬合后，得到采樣點數和相對頻率之間的關系如圖6所示。利用最小二乘法擬合，并將標定過程中采樣點數壓縮至與小型化樣機中采樣點數一致，選用高次擬合關系，得到相對波數υ(n)=( 8.36709879×10-5×n3-0.206195211×n2+ 28.0125375×n-477.432750)/29979.2，其中n為采樣點數，在擬合峰值位置后，其取值并不局限于整數，υ(n)單位cm-1。

圖5 波長標定中的原始干涉信號Fig.5 Raw interference signal of wavenumber calibrating

圖6 激光器出光波長的頻率標定曲線Fig.6 Frequency calibrating curve of the laser output

在實驗室環境中，對零流速進行了測量驗證，搭建光程近似相等的平行雙光路，在沒有流速時輸出各自光路上的原始信號和吸光度信號,如圖7所示。首先對吸光度信號做粗略峰值查找，定位至整數峰值點后，對該點左右小范圍內(20個點)做最小二乘擬合，擬合后的二次函數最值點作為最終吸收峰值位置，用于流速計算。

圖7 零流速時雙光束的原始信號和吸收峰相對位置Fig.7 Raw spectral signal and relative index of absorption peak under zero velocity

由圖7可知，在零流速狀態下，多普勒頻移量約為0.0015cm-1，換算得流速為45m/s，該頻移可作為系統本底偏移，在測量中進行扣除。

2.2 測量原理可行性驗證

首先，在火箭橇靜止狀態下，分別測量了雙基推進劑和改性雙基推進劑裝藥的羽流流速，并與同時測到的發動機工作壓強進行對比，結果見圖8。由圖8可知，測得的羽流流速曲線與壓強曲線變化趨勢一致，其中雙基推進劑流速在點火初始時刻產生一個峰值，隨后流速緩慢上升至位于0.9～1.39s的穩定區，整個點火過程中平均流速為831.8m/s；改性雙基推進劑裝藥流速在達到1100m/s后直接達到位于0.33～0.88s的穩定區，并且在該時段內緩慢下降，整個點火過程平均流速為1057.5m/s。由于流速測量的響應速度遠小于壓強采樣速率，部分快速變化的細節并沒有及時測量到。另外，由于流速測量位置位于發動機出口后30cm處，而壓強曲線變化直接反映發動機的內部情況，故從測得曲線上可以看出，流速響應滯后于壓強變化[1]。

圖8 雙基推進劑裝藥和改性雙基推進劑裝藥燃燒過程中羽流流速和發動機內壓強隨時間的變化曲線Fig.8 Change in plume velocity and pressure in the combustion chamber with time for double-base propellant charge and composite modified double-base propellant charge

通過對比可知，TDLAS技術在該測量場合能真實反映出流速的變化趨勢，從流速測量結果可以推測出壓強的變化趨勢，驗證了該技術在固體推進劑裝藥流速測量中的可行性。

2.3 在線測量試驗

利用前文提到的小型化測量樣機，搭載火箭橇，進行動態測量實驗。樣機集成鋰電池供電，可持續工作2h。試驗中，前2s先由助推火箭將火箭橇加速至240m/s，然后發動機開始工作0.5s，發動機燃燒結束后，借助軌道上的水袋使火箭橇減速至靜止。試驗中記錄的雙光路多普勒頻移和羽流流速測量結果見圖9，從圖9(b)可看出，在0～2s內，速度近似線性加速至240m/s，在第2s時刻，發動機開始工作，流速上升至穩定區，平均流速1249.8m/s；0.5s后，推進劑裝藥燃燒結束，流速降至助推末段的240m/s，并接近勻速保持到第4s到達水袋時刻，由于水袋破裂，大量液態水遮擋光路，第4s后數據無效。

圖9 試驗過程中的多普勒頻移和測得的羽流流速曲線Fig.9 The Doppler shift curve and measured plume velocity curve in experiment

3 結 論

(1)經過靜態和動態試驗驗證，TDLAS技術作為一種非接觸式光學測量手段，利用燃燒產物H2O分子位于1392nm處的吸收譜線，可有效測量固體推進劑裝藥羽流流速，且羽流流速變化和發動機壓強變化具有一致特征，為衡量發動機燃燒狀態提供了一項新的測量參數。

(2)TDLAS技術借助其結構緊湊、低功耗、易集成等特點，可用于開展隨發動機飛行的動態羽流流速測量，是獲取實時在線燃燒狀態參數的一個有效途徑，小型化羽流流速測量系統為后期在線測量方案提供了參考。

[1] 楊燕京, 趙鳳起, 孫美, 等. 基于TDLAS技術的雙基推進劑裝藥排氣羽流流速特征研究 [J]. 火炸藥學報, 2016, 39(3): 70-74.

YANG Yan-jing, ZHAO Feng-qi, SUN Mei, et al. Investigation on the exhaust plume velocity characteristics of double-base and CMDB propellants charge based on TDLAS technique[J]. Chinese Journal of Explosives & Propellants (Huozhayao Xuebao), 2016, 39(3): 70-74.

[2] 施衛平, 舒瑋. 慢周期性流動中畢托管和熱線測速的對比 [J]. 試驗力學, 1992, 7(4): 377-381.

SHI Wei-ping，SHU Wei. Comparison between pitot tube and hot-wire anemometer in slow periodic flows [J]. Journal of Experimental Mechanics, 1992, 7(4):377-381.

[3] 徐驚雷. PIV 技術在超及高超聲速流場測量中的研究進展 [J]. 力學進展, 2012, 42(1): 81-90.

XU Jing-lei. The development of the PIV experimental study of the super/hypersonic flowfield [J]. Advances in Mechanics， 2012，42(1):81-90.

[4] ZHU Hui-min，SU Ke. Investigation of spay characteristics by using phase-Doppler particle analyzer [J]. Journal of Propulsion Technology, 1998, 19(5):55-65.

[5] 張 亮，劉建國，闞瑞峰，等. 基于可調諧半導體激光吸收光譜技術的高速氣流流速測量方法研究[J]. 物理學報，2012，61(3)：034214-1-7.

ZHANG Liang，LIU Jian-guo，KAN Rui-feng，et al. On the methodology of measuring high-speed flows using tunable diode laser absorption spectroscopy [J]. Acta Physica Sinica, 2012, 61(3): 034214-1-7.

[6] Chang L S,Jeffries J B, Hanson R K. Mass flux sensing via tunable diode laser absorption of water vapor[J]. AIAA Journal, 2010, 48(11):2687-2693.

[7] Philippe L C, Hanson R K. Laser diode wavelength-modulation spectroscopy for simultaneous measurement of temperature, pressure, and velocity in shock-heated oxygen flows[J]. Applied Optics, 1993, 32(30): 6090-6103.

[8] Kent L. Development of a real-time diode laser mass flux sensor for simultaneous measurement of density and velocity of oxygen[D]. Stanford: Stanford University, 2005.

[9] Shawn D W. Development of a tunable diode laser probe for measurement in hypervelocity flows[D]. Stanford: Stanford University, 2000.[10] Todd B, Skip W. Development of an in-flight non-intrusive mass capture system[C]∥45nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Denver:American Institute of Aeronautics and Astronautics, 2009.

[11] 賈良權，劉文清，闞瑞峰，等. 波長調制-TDLAS技術測量風洞中氧氣流速方法研究[J]. 中國激光，2015，42(7)：304-311.

JIA Liang-quan, LIU Wen-qing, KAN Rui-feng, et al. Study on oxygen velocity measurement in wind tunnel by wavelength modulation-TDLAS technology[J]. Chinese Journal of Lasers, 2015, 42(7): 304-311.

[12] LI F, Yu X, Gu H, et al. Measurement of flow parameters in a scramjet combustor based on near-infrared absorption[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43:1-7.

[13] Chang L S. Development of a diode laser sensor for measurement of mass flux in supersonic flow[D]. Stanford: Stanford University, 2011.

Measurement Method of Plume Velocity for Solid Propellant Charge Based on TDLAS

YAO Lu1, YAO De-long2, YI Jian-hua3, KAN Rui-feng1, YANG Yan-jing3, XU Zhen-yu1, RUAN Jun1, LIU Jian-guo1

(1.Key Lab of Environment Optics and Technology, Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Science, Hefei 230031, China; 2. Test and Measurement Academy, China North Industries Group, Huayin Shaanxi 714200, China; 3.Science and Technology on Combustion and Explosion Laboratory, Xi′an Modern Chemistry Research Institute, Xi′an 710065, China)

A new method based on tunable diode laser absorption spectroscopy(TDLAS) was introduced to measure the plume velocity of solid propellant charge, which was achieved by selecting one absorption line near 1392nm of the major combustion product-H2O through frequency calibration, cross arrangement of optical path, and using the Doppler frequency shift effect. In the static state test, the characteristic curves of the plume velocity and the pressure curve in the combustion chamber were compared, the feasibility of the technique for the measurement of the plume velocity of solid propellant was verified. A compact measurement prototype was developed, which is powered by lithium battery, the response time is 30ms, the mass is 2kg, the power consumption is less than 5W, the velocity is measured online and the result is stored.The dynamic test flight on the rocket sled was carried out,and the available results of plume velocity were obtained.The results show that TDLAS can be used to measure the dynamic plume velocity and it is an effective process for obtaining the online combustion parameters.

tunable diode laser absorption spectroscopy; TDLAS; solid propellant charge; plume velocity

10.14077/j.issn.1007-7812.2016.05.005

2016-07-14；

2016-08-26

國家自然科學基金青年科學基金項目(No. 61205151)

姚路(1987-)，男，博士，從事激光光譜測量及診斷技術研究。E-mail: lyao@aiofm.ac.cn

TJ55;V512

A

1007-7812(2016)05-0035-05