周 鑫， 金 星， 王廣宇

（1.裝備學(xué)院 研究生管理大隊(duì)，北京101416； 2.裝備學(xué)院 航天裝備系，北京101416；3.裝備學(xué)院 基礎(chǔ)部，北京101416）

在研究引擎燃燒物理、化學(xué)過(guò)程時(shí)，需要對(duì)引擎中溫度、質(zhì)量流量等物理量進(jìn)行在線測(cè)量。此外，這些引擎和研究設(shè)備通常都處于高溫、高壓、高流速和多相流的惡劣物理環(huán)境下，這對(duì)引擎的傳感測(cè)量系統(tǒng)提出了可靠性高、重復(fù)性好、無(wú)氣體分子的交叉敏感等特殊要求，而傳統(tǒng)的測(cè)量系統(tǒng)無(wú)法或不能很好地滿足這些測(cè)量要求［1］。可調(diào)諧半導(dǎo)體激光吸收譜（TDLAS）分析測(cè)量技術(shù)能滿足高溫、高壓燈惡劣物理環(huán)境的要求，其TDLAS傳感系統(tǒng)具有體積小、重量輕、可靠性高和測(cè)量的特定性好等優(yōu)點(diǎn)［2］。

TDLAS技術(shù)的測(cè)量對(duì)象包括氣體溫度、組分濃度、氣體速度等物理量，其中溫度測(cè)量具有重要的地位。采用直接吸收譜技術(shù)中的兩線測(cè)溫技術(shù)對(duì)氣體溫度測(cè)量進(jìn)行建模與仿真，利用Matlab中的動(dòng)態(tài)仿真工具Simulink，建立了光源模型、氣室模型和數(shù)據(jù)檢測(cè)模型。在設(shè)定環(huán)境條件下，通過(guò)模型仿真得到測(cè)量的氣體溫度并進(jìn)行分析。本文旨在通過(guò)對(duì)氣體溫度測(cè)量進(jìn)行探索，為基于TDLAS技術(shù)的傳感器設(shè)計(jì)打好基礎(chǔ)。

1 基于TDLAS的溫度測(cè)量技術(shù)

TDLAS測(cè)量設(shè)備體積小，適于飛行測(cè)量，對(duì)于其他基于激光技術(shù)的測(cè)量方法，其非介入式、響應(yīng)快的優(yōu)勢(shì)也非常明顯［3－5］。

1.1 Beer－Lambert關(guān)系式

當(dāng)某一束頻率為υ的激光通過(guò)被測(cè)氣體，其入射光強(qiáng)和穿透光強(qiáng)滿足Beer－Lambert關(guān)系式

式 中：I0為 入 射 光 強(qiáng)；It為 穿 透 光 強(qiáng)；L為 吸 收 長(zhǎng)度（cm）；Kυ為吸收系數(shù)（cm－1），且

1.2 多普勒展寬、碰撞展寬

多普勒展寬是由分子無(wú)規(guī)則熱運(yùn)動(dòng)所產(chǎn)生。運(yùn)用統(tǒng)計(jì)學(xué)原理對(duì)稀薄氣體中的分子速度分布進(jìn)行研究，發(fā)現(xiàn)其分布符合麥克斯韋方程。因此，多普勒展寬可運(yùn)用高斯曲線來(lái)描述：

式（3）為高斯線型公式。ΔυD為高斯半高寬（cm－1），即多普勒寬度；υ0是中心頻率（cm－1）。

碰撞展寬是另一種重要的展寬機(jī)制，是由粒子間相互發(fā)射和吸收導(dǎo)致的碰撞所產(chǎn)生的。該機(jī)制基于2個(gè)假設(shè)：①碰撞是成對(duì)產(chǎn)生的；②相對(duì)于整個(gè)碰撞事件，碰撞持續(xù)時(shí)間短到可忽略。碰撞展寬通過(guò)洛倫茲曲線來(lái)反映，其方程如下：

式（4）為洛倫茲線型公式。ΔυC為洛倫茲半高寬（cm－1），即碰撞寬度。

多普勒展寬在低壓情況下非常明顯，碰撞展寬在高壓情況下占主導(dǎo)地位。整個(gè)展寬是碰撞展寬和多普勒展寬的組合。在展寬機(jī)制獨(dú)立的情況下，總的展寬機(jī)制可通過(guò)高斯線型公式和洛倫茲線型公式的卷積來(lái)得到，即Voigt線型公式為

1.3 溫度測(cè)量技術(shù)

典型的直接吸收實(shí)驗(yàn)結(jié)構(gòu)圖如圖1所示。

圖1 典型的掃描波長(zhǎng)直接吸收測(cè)量的結(jié)構(gòu)圖

在溫度測(cè)量方面主要使用兩線測(cè)量技術(shù)［6］。該技術(shù)使用2個(gè)具有不同溫度依賴的躍遷進(jìn)行比較，以得到氣體的溫度。線型強(qiáng)度S（T）可以從已知的關(guān)于T0的線型強(qiáng)度來(lái)得到

式中：Q（T）是配分函數(shù)；h是普朗克常量（J·s）；c是光速（cm／s）；k是玻耳茲曼常數(shù)（J／K）；E″是低狀態(tài)能量（cm－1）。

配分函數(shù)Q（T）使用下面的多項(xiàng)式。多項(xiàng)式中常數(shù)可在HITRAN數(shù)據(jù)庫(kù)中查詢［7］。

因?yàn)楣饴窂介L(zhǎng)度和水的分壓都相同，所以這2個(gè)積分比R可化簡(jiǎn)為線型強(qiáng)度的比

式中：A為吸收系數(shù)Kυ的積分。氣體溫度可通過(guò)式（9）得到：

2 基于Simulink的TDLAS測(cè)量系統(tǒng)仿真模型

利用Matlab 7.4的Simulink環(huán)境中豐富的模塊庫(kù)，建立光源模型、氣室模型和數(shù)據(jù)檢測(cè)模型，測(cè)量并分析設(shè)定條件下的氣體溫度。

本仿真模型采用兩線測(cè)溫法，譜線1頻率為7 185.6cm－1，譜線2頻率為7 444.35cm－1；設(shè)定的水濃度為10%，光程為30cm；譜線1的低態(tài)能量為1 045.1cm－1，譜 線2的 低 態(tài) 能 量 為1 774.8cm－1；仿真模型的線型采用高斯線型，展寬機(jī)制為多普勒展寬，對(duì)應(yīng)該展寬的氣體壓強(qiáng)為幾個(gè)mPa，因而這里設(shè)定氣體壓強(qiáng)為0.002Pa。整個(gè)TDLAS仿真系統(tǒng)采用層次化結(jié)構(gòu)，由光源模塊、氣室模塊、數(shù)據(jù)檢測(cè)模塊3大部分構(gòu)成。

1）光源模塊。光源模塊如圖2所示。其中，調(diào)諧信號(hào)由鋸齒波發(fā)生器產(chǎn)生。調(diào)制后頻率的中心值設(shè)定為7 300cm－1。

2）氣室模塊。氣室模塊主要是根據(jù)待測(cè)氣體在制定中心波長(zhǎng)處譜線吸收線型函數(shù)模型來(lái)構(gòu)建，其功能是模擬激光在氣室中的吸收。氣室模塊建模關(guān)鍵在于譜線線型函數(shù)Φ（v）的選取。在本文中譜線線型函數(shù)選取高斯線型函數(shù)ΦD（v）。

高斯線型函數(shù)模型如圖3所示。氣室模型如圖4所示。

3）數(shù)據(jù)檢測(cè)模塊。數(shù)據(jù)檢測(cè)模塊主要是根據(jù)待測(cè)氣體在指定吸收譜線的吸收來(lái)構(gòu)建的，如圖5所示。

圖2 光源模塊

圖3 高斯線型函數(shù)模型

圖4 氣室模型示意圖

圖5 數(shù)據(jù)檢測(cè)模塊

數(shù)據(jù)檢測(cè)的關(guān)鍵在于測(cè)量方法的選取。兩線測(cè)溫法需要對(duì)穿透光強(qiáng)的吸收進(jìn)行積分，利用測(cè)溫法公式計(jì)算得到被測(cè)氣體的溫度。

3 仿真結(jié)果分析

通過(guò)已建立的測(cè)量系統(tǒng)仿真模型，結(jié)合給定的初始條件就能夠進(jìn)行測(cè)量系統(tǒng)的仿真。由于在氣壓、光程和初始光強(qiáng)不變的情況下，氣體濃度的變化會(huì)給氣體吸收系數(shù)帶來(lái)改變。在實(shí)際測(cè)量環(huán)境中，比如燃燒環(huán)境，水蒸氣的濃度會(huì)產(chǎn)生改變，從而帶來(lái)氣體吸收系數(shù)的改變。為了考慮這種改變對(duì)測(cè)量系統(tǒng)的影響，本項(xiàng)目首先仿真了氣體濃度改變對(duì)氣體吸收系數(shù)曲線的影響。然后給出了溫度測(cè)量的仿真過(guò)程圖，按照仿真過(guò)程步驟實(shí)現(xiàn)溫度測(cè)量的仿真。

3.1 氣體濃度改變對(duì)氣體吸收系數(shù)曲線的影響

隨著氣體濃度升高，氣體吸收系數(shù)的峰值不斷增大，而譜線產(chǎn)生吸收的頻率范圍（時(shí)間與頻率相對(duì)應(yīng)）沒(méi)有改變。在溫度測(cè)量中，利用了氣體吸收系數(shù)的比值，按照式（2），氣體濃度的改變不會(huì)影響溫度測(cè)量。氣體濃度為25%時(shí)氣體吸收系數(shù)峰值約為3.5E－3，而氣體濃度為5%時(shí)氣體吸收系數(shù)峰值約為6.5E－4，后者比前者小了一個(gè)數(shù)量級(jí)。在實(shí)際系統(tǒng)中，被測(cè)量物理量減小一個(gè)數(shù)量級(jí)，對(duì)應(yīng)的測(cè)量設(shè)備就要相應(yīng)的提高測(cè)量精度的數(shù)量級(jí)。

3.2 仿真過(guò)程圖與測(cè)量溫度

給定初始條件：氣體溫度1 601K；譜線1頻率7 185.6cm－1，譜線2頻率7 444.35cm－1；設(shè)定的氣體壓強(qiáng)0.002Pa，水濃度10%，光程30cm；譜線1的低態(tài)能量1 045.1cm－1，譜線2的低態(tài)能量1 774.8cm－1。仿真過(guò)程如圖6所示。

圖6 仿真過(guò)程圖

在仿真過(guò)程中，對(duì)吸收系數(shù)和穿透光強(qiáng)的數(shù)據(jù)進(jìn)行采集，得到穿透光強(qiáng)和氣體吸收曲線。譜線1穿透光強(qiáng)最大衰減發(fā)生在7 185.6cm－1，其對(duì)應(yīng)氣體吸收系數(shù)曲線峰值在7 185.6cm－1，二者對(duì)應(yīng)的波數(shù)一致，其內(nèi)在原理為譜線1在該處達(dá)到了氣體吸收曲線的峰值，氣體對(duì)穿過(guò)氣室的激光產(chǎn)生一個(gè)氣體吸收的峰值，光強(qiáng)隨之達(dá)到最大衰減，譜線2的情況同譜線1。同時(shí)，譜線1和譜線2的吸收峰值發(fā)生在不同的時(shí)間，這是由于激光頻率按照設(shè)定值在7 100～7 500cm－1之間掃描，2個(gè)譜線產(chǎn)生吸收峰值的時(shí)間所對(duì)應(yīng)的激光頻率正 好 是 譜 線 的 頻 率 中 心（7 185.6cm－1和7 444.35cm－1）。另外，譜線1的吸收峰值明顯大于譜線2的吸收峰值，譜線1光強(qiáng)衰減峰值約為5.8，譜線2光強(qiáng)衰減峰值約為4.7。

仿真計(jì)算中，根據(jù)氣體吸收系數(shù)數(shù)據(jù)，對(duì)氣體吸收系數(shù)曲線進(jìn)行積分，得到2個(gè)譜線的吸收系數(shù)積分比。使用該比值，結(jié)合式（9），得到測(cè)量溫度1 601.799K，誤差為0.05%。在實(shí)際的測(cè)量系統(tǒng)中，數(shù)據(jù)采集會(huì)有誤差，信號(hào)會(huì)有干擾，還有信噪比等影響，所以實(shí)際系統(tǒng)的誤差會(huì)較為明顯。

4 結(jié) 束 語(yǔ)

在條件給定的情況下，仿真模型能模擬出氣體的吸收系數(shù)曲線、穿透光強(qiáng)波形等仿真效果，實(shí)現(xiàn)對(duì)氣體溫度的測(cè)量，并能反映濃度等各參數(shù)在變化時(shí)對(duì)氣體的吸收系數(shù)曲線、穿透光強(qiáng)波形的影響。這對(duì)于TDLAS測(cè)量系統(tǒng)軟件、硬件的開(kāi)發(fā)都有一定的參考價(jià)值。

（

）

［1］LIU C，RIEKAER B.Near－infrared diode laser absorption diagnostic for temperature and water vapor in a scramjet combustor［J］.Optical Society of America，2005，31（6）：56－66.

［2］TUDOR I.Palaghita pattern factor sensing and control based on diode laser absorption［J］.Applied Physics Letters，2005，35（6）：26－28.

［3］DAVID S，ALAN C.Frequency modulation and wavelength modulation spectroscopies comparison of experimental methods using a lead－salt diode laser［J］.Applied Optics，1992，31（6）：124－132.

［4］JOEL S A.Frequency－modulation spectroscopy for trace species detection theory and comparison experimental methods［J］.Optical Society of America，1992，4（7）：26－28.

［5］LI H，F(xiàn)AROOQ A.Near－infrared diode laser absorption sensor for rapid measurements of temperature and water vapor in a shock tube［J］.Applied Physics，2007，89（2）：133－142.

［6］HANSON R K，JEFFRIES J B.Advances in Laser－based Sensors for Propulsion Systems［C］／／AIAA.24thAerodynamic Measurement Technology and Ground Testing Conference，Portland.Oregon：AIAA，2004：2476－2486.

［7］DHARAMSI N，BULLOCK M.Applications of wavelengthmodulation spectroscopy in resolution of pressure and modulation broadened spectra［J］.Applied Physics，1996，78（3）：283－292.