丁水汀，鄧長春，邱天，*，李江涵，單曉明，賀宜紅

（1.北京航空航天大學(xué) 能源與動力工程學(xué)院，北京100083； 2.航空發(fā)動機氣動熱力國家級重點實驗室，北京100083；3.中國航發(fā)湖南動力機械研究所，株洲412000）

航空發(fā)動機空氣系統(tǒng)是處于發(fā)動機主流道內(nèi)側(cè)、外側(cè)的轉(zhuǎn)子與靜子或轉(zhuǎn)、靜子之間的聯(lián)系各腔、間隙、流孔的空氣流路的總稱[1]。廣泛應(yīng)用流體網(wǎng)絡(luò)法[2-12]計算空氣系統(tǒng)壓力分布、流量分配等參數(shù)。但目前能夠計算瞬態(tài)換熱特性的流體網(wǎng)絡(luò)法在公開文獻中鮮見報道。為彌補這一不足，需開展空氣系統(tǒng)元件瞬態(tài)換熱試驗。容腔元件作為典型瞬態(tài)元件，是空氣系統(tǒng)瞬態(tài)換熱特性試驗的基礎(chǔ)研究元件。

國內(nèi)外相關(guān)學(xué)者對容腔充放氣過程進行過一系列研究。Dutton和Robert[13]開展了單孔容腔充、放氣試驗，對比了絕熱模型和等熵模型的區(qū)別，但未研究充、放氣過程中的換熱特性。Thorncroft等[14]針對單孔容腔充、放氣過程，利用模型與試驗對比的方法確定多變過程中多變指數(shù)，建立了可以預(yù)測容腔充、放氣過程中腔內(nèi)壓力與溫度隨時間變化的模型。李世武[15]、馬崇揚[16]等對方腔內(nèi)自然對流換熱展開過數(shù)值模擬的研究。楊麗紅[17]對容腔放氣過程進行過數(shù)值模擬研究，并建立了熱力學(xué)模型。郭鐘華等[18]針對氣動真空系統(tǒng)中容腔，利用“停止法”得出了真空抽吸過程中的換熱系數(shù)。但此方法對于試驗系統(tǒng)控制能力要求較為苛刻，并不適用于大多數(shù)類型的試驗。丁水汀、于航、邱天[19-20]等利用一種非絕熱單孔容腔瞬態(tài)響應(yīng)的零維建模方法，進一步表明非絕熱容腔模型較絕熱模型能夠更精確地反映容腔的真實響應(yīng)規(guī)律，且在10 Hz的數(shù)采頻率下對非絕熱單孔容腔元件充氣過程換熱特性開展過試驗研究，獲得了容腔元件的初步換熱規(guī)律。

由于真實航空發(fā)動機空氣系統(tǒng)中的瞬態(tài)響應(yīng)時間在毫秒量級，試驗中為捕捉更準確的容腔瞬態(tài)換熱特性，需要提高數(shù)據(jù)采集頻率。在高頻狀態(tài)下，某些試驗測試誤差會被數(shù)據(jù)處理的差分過程放大。因此了解試驗中換熱特性對各種測試參數(shù)隨機誤差的敏感性并掌握有效的誤差抑制方法是試驗成功的關(guān)鍵。為指導(dǎo)相關(guān)試驗的開展，本文利用瞬態(tài)數(shù)值計算結(jié)果分析50 Hz數(shù)據(jù)采集頻率下，測試隨機誤差對試驗結(jié)果的影響。將利用計算流體力學(xué)（Computational Fluid Dynamics，CFD）技術(shù)計算所得的非絕熱單孔容腔充氣過程中滿足各守恒定律的數(shù)據(jù)作為整個充氣過程中真實數(shù)據(jù)。結(jié)合試驗方案和現(xiàn)有技術(shù)水平，對相關(guān)數(shù)據(jù)附加噪聲以模擬試驗測試隨機誤差。將利用無噪聲數(shù)據(jù)與有噪聲數(shù)據(jù)算出的換熱特性進行對比，分析換熱特性誤差對不同參數(shù)測試隨機誤差的敏感性。并利用改進經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）[21-23]算法對實驗 數(shù)據(jù)進行光順處理，降低測試隨機誤差的影響。本文研究內(nèi)容對空氣系統(tǒng)中其他類型容腔元件的充放氣過程換熱特性試驗的開展均有借鑒意義。

1 相關(guān)理論

1.1 物理模型

航空發(fā)動機空氣系統(tǒng)中的一些腔室結(jié)構(gòu)可抽象為單孔容腔結(jié)構(gòu)[24]。本文將實際航空發(fā)動機空氣系統(tǒng)內(nèi)部復(fù)雜單孔容腔結(jié)構(gòu)簡化為單孔方形腔，其物理模型如圖1所示。方形腔邊長0.3 m，進氣口是一收縮狀的噴嘴，處于壁面中心位置。噴嘴進口直徑為30 mm，出口直徑為10 mm。

圖1 容腔物理模型Fig.1 Physical model of cavity

1.2 試驗原理及測試方法

單孔容腔充氣過程換熱特性試驗應(yīng)用的假設(shè)與文獻[19]中提出的假設(shè)相同。實際試驗中根據(jù)能量守恒方程獲取充氣過程中的換熱特性。單孔容腔在整個充氣過程中任意時間段內(nèi)滿足開口系統(tǒng)能量守恒方程：

腔內(nèi)氣體的內(nèi)能變化量d U等于流進的空氣帶來的焓值增量d H和腔內(nèi)氣體與容腔壁面的換熱量δQ之和。在任一Δt時間段內(nèi)，d U和d H可以用向后差分的方法近似表示為

式中：cV為空氣比定容熱容；T（t）為腔內(nèi)平均溫度；cPin為進口空氣比定壓熱容；Tin為進口氣流溫度；m（t）為腔內(nèi)t時刻氣體質(zhì)量；Δt=0.02 s。腔內(nèi)壁面與腔內(nèi)氣體交換的熱流密度q表示為

在計算熱流密度的過程中，忽略噴嘴進口面積，即容腔表面積A與邊長L的關(guān)系為：A=6L2。對流換熱系數(shù)h及平均努塞爾數(shù)Nu表示為

式中：Tw為容腔壁面溫度；λ為空氣導(dǎo)熱系數(shù)。試 驗 中 所 需 測 量 的 物 理 量 包 括L、Tin、T（t）、m（t）和Tw。其 中L、Tin與Tw在試 驗中保持 不變，屬于穩(wěn)態(tài)測量參數(shù)。T（t）與m（t）是瞬態(tài)測量參數(shù)。

針對T（t）的測量，可以采用響應(yīng)時間在毫秒量級的高頻響熱電偶[25-27]。利用流量計測量進口流量隨時間變化關(guān)系，再進行積分處理可測得m（t）。但在目前的氣體流量測量技術(shù)中，毫秒級的熱式氣體質(zhì)量流量計往往在某些特定的流量下具備快速響應(yīng)特性，流量測試范圍也存在一定的局限[28]。因此試驗中難以實現(xiàn)m（t）的直接測量。一種簡單的方法是根據(jù)理想氣體狀態(tài)方程反推m（t），如式（7）所示：

式中：P（t）為腔內(nèi)平均壓力；Rg為氣體常數(shù)；V為容腔體積。此方法要求測量P（t）。文獻[29-30]介紹的壓阻型壓力傳感器可以滿足本文試驗頻率要求。

1.3 數(shù)據(jù)樣本的獲取

樣本數(shù)據(jù)來源于瞬態(tài)CFD結(jié)果。工質(zhì)為理想氣體。腔內(nèi)溫度變化不大時，忽略物性參數(shù)隨溫度變化的影響。物性參數(shù)如表1所示。進口氣流溫度為295 K，進口絕對壓力為191 504 Pa，壁面溫度為400 K，進口噴嘴壁面絕熱。腔內(nèi)初始壓力為101 325 Pa，初始速度為0，初始溫度為295 K。

CFD結(jié)果如圖2所示。充氣初始階段，P（t）與T（t）快速增高。充氣階段結(jié)束后不再有外部氣體對腔內(nèi)進行能量補充，但腔內(nèi)氣體溫度低于腔壁面溫度，氣體與壁面存在換熱，T（t）緩慢升高，P（t）變化不明顯。根據(jù)CFD結(jié)果可以確定在此工況下，充氣過程大約持續(xù)0.68 s，故取0.68 s前的計算數(shù)據(jù)作為樣本數(shù)據(jù)。

表1 理想氣體物性參數(shù)Table 1 Physical property parameters of ideal gas

圖2 CFD結(jié)果Fig.2 CFD results

2 誤差分析

2.1 測試隨機誤差給定

隨機誤差是試驗中難以消除的誤差，大小取決于現(xiàn)有測試技術(shù)水平。本文試驗中隨機誤差來源于試驗件尺寸、氣流壓力與溫度的測量。根據(jù)文獻[31-39]介紹的高精度加工及測試技術(shù)，試驗中各測試數(shù)據(jù)的誤差給定如表2所示。穩(wěn)態(tài)測試誤差 δL、δTin與 δTw給定為常數(shù)，瞬態(tài)測試誤差δP（t）與 δT（t）按正態(tài)分布3σ原則[40]給定。

表2 試驗隨機誤差Table 2 Random error in experiment

2.2 分析結(jié)果

定義Nuer為含有誤差的努賽爾數(shù)，換熱特性相對誤差百分比δNu為

表2給定的常數(shù)誤差引起的結(jié)果誤差如圖3所示。幾何尺寸的測量誤差引起的換熱特性誤差為常數(shù)，相對誤差只有0.07%。在充氣初始階段的0.02 s內(nèi)，δTin引起的誤差存在誤差峰值。因為這一時間段內(nèi)腔內(nèi)參數(shù)不嚴格符合集總假設(shè)，按照差分的方法算得的d H相對誤差較大。在整個充氣過程中由于非集總性引起的峰值較小，因此本文不對此進行詳細研究。整個充氣過程中常數(shù)誤差引起的換熱特性偏差先增后降。這是因為隨著充氣過程的進行，腔內(nèi)溫度與壓力不斷升高，對換熱誤差產(chǎn)生放大效應(yīng)。到充氣階段后期，腔內(nèi)壓力與溫度變化緩慢，進氣流量不斷下降直至為0。進氣流量的減小對誤差產(chǎn)生縮小效應(yīng)。在包含全部單向偏差的狀態(tài)下，換熱特性相對誤差在1.80%以內(nèi)。

圖3 穩(wěn)態(tài)誤差對Nu的影響Fig.3 Influence of steady error on Nu

圖4給出腔內(nèi)瞬態(tài)相對壓力prel測試偏差對換熱特性的影響。在充氣過程的大部分時間點上，prel測試偏差對換熱特性的影響不大。prel偏差在38.78 Pa以內(nèi)的情況下，換熱特性最大相對偏差不到2.00%。

圖5給出腔內(nèi)瞬態(tài)溫度測試誤差對換熱特性的影響。腔內(nèi)溫度測試誤差對試驗結(jié)果影響顯著，所得換熱特性折線圖存在大幅波動。溫度誤差在0.464℃以內(nèi)的情況下?lián)Q熱特性相對誤差最高可達130.22%。

圖6展示了腔內(nèi)溫度與壓力誤差同時存在時的換熱特性結(jié)果。換熱特性波動狀態(tài)與圖5（b）相似，波動幅值較大，最高誤差達到129.07%，時均誤差為25.24%?？梢姡粌?nèi)瞬態(tài)溫度測試隨機誤差對試驗結(jié)果的影響最為顯著，使得試驗結(jié)果不能準確反映換熱規(guī)律。

在直接測量的瞬態(tài)誤差無明顯波動狀態(tài)下，換熱特性的大幅度波動表明在數(shù)據(jù)處理過程中存在誤差放大環(huán)節(jié)。誤差放大的環(huán)節(jié)存在于腔內(nèi)氣體內(nèi)能增量d U與質(zhì)量流量m（t）-m（t-Δt）的計算中。波動來源于 ΔP與 Δ（P/T），表達式如下：

圖4 壓力偏差對Nu的影響Fig.4 Influence of pressure error on Nu

圖5 溫度誤差對Nu的影響Fig.5 Influence of temperature error on Nu

圖6 全部瞬態(tài)誤差對Nu的影響Fig.6 Influence of all transient errors on Nu

瞬態(tài)測試誤差同時存在的情況下 ΔP 與Δ（P/T）的計算結(jié)果如圖7所示?？梢园l(fā)現(xiàn)，ΔP的計算對瞬態(tài)測試誤差的放大作用不明顯，Δ（P/T）影響顯著，結(jié)果折線圖存在較大波動現(xiàn)象。這一波動現(xiàn)象在后續(xù)數(shù)據(jù)處理中會被進一步放大，導(dǎo)致試驗結(jié)果難以反映實際的換熱特性。因此，除了提高測試精度外，抑制 Δ（P/T）的波動也是試驗的關(guān)鍵環(huán)節(jié)。

圖7 ΔP與 Δ（P/T）計算結(jié)果Fig.7 ΔP andΔ（P/T）calculation results

3 EMD光順處理

3.1 EMD算法應(yīng)用

EMD算法對于非線性、非平穩(wěn)信號的處理具有獨特的優(yōu)越性，其基本原理及信號處理過程在文獻[21-23]中有詳細介紹。同時此算法對于離散數(shù)字曲線具有良好的光順作用[22]，因此適用于試驗數(shù)據(jù)Δ（P/T）的光順處理中。

應(yīng)用EMD算法進行濾波及光順處理時，處理效果會受到端點效應(yīng)[16-18]的影響。本文基于極值平移法[18]的基本思想，根據(jù)信號端點附近的極值特征，在信號端點后增加兩對極值點。具體方法是取0.62～0.68 s的數(shù)據(jù)，每個數(shù)據(jù)減去一常數(shù)后作為延續(xù)信號直接放在原始數(shù)據(jù)后。極值點的數(shù)量、選擇方法及被減常數(shù)取值不同會對處理結(jié)果產(chǎn)生一定影響，本文不重點研究。

3.2 光順處理效果

利用EMD算法對 Δ（P/T）進行光順處理的效果如圖8所示。應(yīng)用傳統(tǒng)EMD算法對 Δ（P/T）進行處理時，端點效應(yīng)的影響較為明顯。主要體現(xiàn)在充氣過程接近尾聲時處理后數(shù)據(jù)存在反?，F(xiàn)象。應(yīng)用改進EMD算法進行處理，端點效應(yīng)得到有效抑制，Δ（P/T）的波動狀態(tài)得到了明顯的改善。

應(yīng)用改進EMD算法進行光順處理后，得到的換熱特性折線圖如圖9所示。改進EMD算法的應(yīng)用對最終換熱特性的結(jié)果產(chǎn)生顯著改善，換熱特性波動幅值明顯縮小。在充氣過程接近尾聲時的換熱特性誤差較大，最大誤差從129.07%降到63.62%，時均誤差從25.24%降到8.12%。精度大幅提高，換熱特性折線圖可反映容腔充氣過程近似真實的換熱特性變化規(guī)律。

圖8 EMD算法對 Δ（P/T）處理效果對比Fig.8 Comparison ofΔ（P/T）treatment results by EMD algorithm

圖9 改進EMD算法的應(yīng)用效果Fig.9 Application effect of improved EMD algorithm

4 結(jié) 論

針對容腔充氣過程換熱試驗的數(shù)據(jù)處理過程中差分處理方法放大隨機誤差的問題，分析了誤差放大效應(yīng)的主要環(huán)節(jié)并采用改進EMD算法進行誤差抑制，得出以下結(jié)論：

1）誤差放大環(huán)節(jié)主要存在于進氣流量的計算中，此環(huán)節(jié)將腔內(nèi)平均溫度測試誤差放大，影響試驗結(jié)果。

2）腔內(nèi)平均溫度以外的測試參數(shù)隨機誤差對試驗結(jié)果影響較小。

3）改進EMD算法的應(yīng)用可以顯著降低試驗測試隨機誤差的影響，在本文研究的工況中，此算法使試驗最大誤差從129.07%降到63.62%，時均誤差從25.24%降到8.12%。