林家泉，李波，邱岳恒

(中國民航大學 電子信息與自動化學院，天津 300300)

據(jù)統(tǒng)計，民航飛機在起飛階段，由于燃料的不完全燃燒，會產(chǎn)生大量的NOx、SO2、CO、顆粒物等大氣污染物，這些污染物會長期積聚在跑道周圍，嚴重影響空氣質(zhì)量[1-2]。由于飛機在起飛前需要在跑道外等候一段時間，在這期間，大量的污染物會通過發(fā)動機引氣系統(tǒng)進入客艙，污染客艙內(nèi)空氣環(huán)境。因此，在一定的通風量下，提高空調(diào)系統(tǒng)的排污效率，快速地將污染物排出客艙外，減少對乘客的危害，顯得尤為重要。

目前，對于污染物在飛機客艙中的傳播規(guī)律已經(jīng)有了大量的研究，包括不同送風方式和速度下的污染物傳播規(guī)律研究[3-5]、不同種類的污染物傳播規(guī)律研究[6-9]。除了在現(xiàn)有送風形式上的研究之外，有部分研究提出將個性送風口安裝在座椅下方或扶手處等位置，通過個性送風來減少污染物傳播并保持熱舒適[10-13]，但是這些改進的新型送風方式需要對現(xiàn)有飛機上的空調(diào)系統(tǒng)及送風口位置進行大規(guī)模更改，因此更適合在未來飛機的設(shè)計上作為參考，而不能成為提高現(xiàn)有飛機客艙通風效果的選擇。

通過引氣系統(tǒng)增加新風量可以改善客艙內(nèi)空氣質(zhì)量，但這樣會導致更高的發(fā)動機燃油消耗，因此，探尋一種無需增加新風供給量就可以改善客艙空氣品質(zhì)的送風方式，有一定現(xiàn)實意義。非穩(wěn)態(tài)送風對建筑室內(nèi)空氣混合有積極效果，不少學者對此進行了研究。Mesenholler 等[14]對目前國內(nèi)外非穩(wěn)態(tài)送風研究進行了總結(jié)，主要包括間歇性氣流、正弦氣流和模擬自然風氣流，發(fā)現(xiàn)非穩(wěn)態(tài)送風在提高建筑室內(nèi)排污效果和通風質(zhì)量方面有著重要的研究意義。Fallenius 等[15]發(fā)現(xiàn)脈動送風可以產(chǎn)生更多數(shù)量的渦旋，減少停滯區(qū)產(chǎn)生，起到增強混合的積極作用。van Hooff 等[16]研究得出采用時間周期性供應(yīng)的正弦送風速度可以有效提高污染物去除效率。這些都證明了在建筑室內(nèi)環(huán)境中，非穩(wěn)態(tài)送風可以有效提高室內(nèi)污染物去除的效率。

相關(guān)人員在飛機客艙中對方波送風也進行了相關(guān)研究。Wu 和Ahmed[17]發(fā)現(xiàn)在飛機客艙環(huán)境中應(yīng)用方波信號的送風方式可以有效加強新鮮空氣和艙內(nèi)剩余空氣的混合；林家泉和戴仕卿[18]得出在方波信號下混合送風為客艙最優(yōu)的送風方式。但方波信號風速變化較為突然，氣流波動幅度大，易導致乘客的不適，因此，本文采用風速變化較為舒緩的正弦信號，以引氣污染物清除效果為研究對象，在保證通風量一定的情況下，對比不同送風模式下恒值信號送風和正弦信號送風的排污效果，以空氣齡為評價指標對空氣質(zhì)量進行評價，并結(jié)合吹風感指數(shù) (draft rating index，DR)評價指標，評價不同送風工況下乘客的熱舒適性，得出客艙空調(diào)最佳的送風工況。

1 本文研究方法

1.1 飛機客艙模型

計算流體力學(computational fluid dynamics，CFD)數(shù)值模擬是得到客艙環(huán)境中流場分布的有效方法。圖1 為客艙實驗平臺，圖2 為采用CFD 數(shù)值模擬建立的客艙仿真模型，模型中包含有天花板送風口2 個，側(cè)壁送風口10 個，個性化送風口30 個，底部排風口10 個，內(nèi)部包含有5 排，6 列的座椅。為保證計算精度，同時減少網(wǎng)格數(shù)量，采用局部精細化的方法對網(wǎng)格進行處理，采用四面體網(wǎng)格劃分方式，平均網(wǎng)格質(zhì)量均在0.7 以上，網(wǎng)格數(shù)量為534 萬。

圖1 客艙實驗平臺Fig.1 Cabin experiment platform

圖2 Boeing737 客艙仿真模型Fig.2 Boeing737 cabin simulation model

1.2 數(shù)值模型

在數(shù)值模擬中，根據(jù)飛機客艙內(nèi)部流場特點，在湍流計算中使用RNGk-ε湍流模型，可以得到較高的計算準確性[19]，其控制方程如下：

式中：Γ?,eff為廣義擴散系數(shù)；?ˉ為通用變量均值(速度、溫度、質(zhì)量濃度等)；uˉi為xi方向的平均速度，m/s；ρ為密度，kg/m3；S?為廣義源項。

求解壓力-速度耦合方程采用SIMPLE 算法，壓力插值為二階，控制方程的對流項和黏性項均采用二階離散格式。當所有的比例殘差變平并達到最小值時，計算獲得收斂。

1.3 客艙CFD 模型驗證

采用CFD 數(shù)值模擬和實驗驗證結(jié)合的方法，通過粒子圖像測速(particle image velocimetry，PIV)技術(shù)在實驗艙對所建立的數(shù)值模型進行驗證，利用煙霧發(fā)生器產(chǎn)生平均粒徑約為1.5 μm 的顆粒，從天花板和側(cè)壁送風口進入客艙，從底部排風口排出。將所需測量的截面劃分為4 個區(qū)域，使用1 600 萬像素的CCD 相機對第2 排和第3 排座椅中間的示蹤粒子圖像進行采集，單個采集區(qū)域面積為900 mm×600 mm，通過數(shù)據(jù)處理將各個子區(qū)域的流場圖拼接，得到整個截面的粒子運動時均場。圖3 為PIV實驗和仿真結(jié)果對比，客艙在建模過程中進行了簡化和近似，實驗結(jié)果和仿真結(jié)果中得到了基本一致的氣流運動形式。

圖3 PIV 實驗和仿真結(jié)果Fig.3 PIV experiment and simulation results

1.4 網(wǎng)格獨立性驗證

對客艙網(wǎng)格數(shù)量分別為327 萬、388 萬、470 萬、534 萬和739 萬的5 種情況進行模擬，計算采用穩(wěn)態(tài)方式，當變量曲線不隨迭代次數(shù)變化時，認為計算收斂。

圖4 為網(wǎng)格獨立性分析。取第1 排座椅上方10 個采樣點的空氣流速進行對比，可以發(fā)現(xiàn)，網(wǎng)格數(shù)量為534 萬和739 萬2 種情況下的采樣點的速度值較為接近，說明數(shù)量為534 萬的網(wǎng)格可以獲得較為準確的流場信息，增加網(wǎng)格對結(jié)果無明顯影響。因此，經(jīng)過網(wǎng)格獨立性驗證，采用534 萬數(shù)量的網(wǎng)格進行CFD 模擬。

圖4 網(wǎng)格獨立性分析Fig.4 Grid independence analysis

1.5 仿真邊界條件設(shè)定

仿真過程中，選取NO2作為污染物，出口設(shè)置為壓力出口，送風口設(shè)置污染物NO2質(zhì)量濃度為5×10?5kg/m3，共釋放60 s，60 s 后開始送入新鮮空氣，送風溫度為295 K。在不同送風模式下，送風量均保持相同，經(jīng)計算，各送風口的送風速度如表1 所示。其中，正弦送風是以正弦函數(shù)為基準的時間周期型速度函數(shù)，周期為60 s，不同正弦函數(shù)的供氣速度通過UDF 導入FLUENT 中設(shè)置邊界條件進行模擬計算。

表1 送風速度Table 1 Air supply speed m/s

2 仿真結(jié)果分析

2.1 恒值信號送風客艙內(nèi)氣流特征

客艙內(nèi)空氣流動是影響氣體污染物傳播的重要原因。選取客艙內(nèi)第3 排區(qū)域截面的流場圖，觀察不同送風模式下在140，160，180 s 這3 個時刻的流場變化情況。

在恒值信號下，送風速度恒定，天花板送風模式下流場如圖5 所示。新鮮空氣沿艙壁向下運動，到達客艙底部，一部分通過排風口排出艙外，另一部分在客艙中部匯合并形成向上的運動氣流，在中間乘客位置處產(chǎn)生2 個主要的渦旋。側(cè)壁送風模式下流場如圖6 所示。側(cè)壁送風口送出的新鮮空氣，在康達效應(yīng)的影響下，沿行李架表面向上運動，在天花板處相遇后形成向下運動的氣流，在靠近過道乘客位置處形成一順一逆2 個渦旋?；旌纤惋L模式下流場如圖7 所示。從天花板和側(cè)壁2 個送風口送出的新鮮空氣，在行李架拐角處相遇后，產(chǎn)生向下運動的氣流，在中間乘客位置處形成較大的渦旋，向下的氣流撞擊地板后分別向客艙兩側(cè)運動，一部分氣流從排風口排出艙外，另一部分沿艙壁向上運動。通過觀察發(fā)現(xiàn)，隨著時間增加，天花板送風模式、側(cè)壁送風模式和混合送風模式下的渦旋位置均未發(fā)生明顯變化，由于渦旋的渦心處風速低，會形成大面積的空氣滯留區(qū)。

圖5 天花板送風模式下的流場（恒值信號）Fig.5 Flow field of ceiling air supply mode (constant value signal)

圖6 側(cè)壁送風模式下的流場（恒值信號）Fig.6 Flow field of side wall air supply mode (constant value signal)

圖7 混合送風模式下的流場（恒值信號）Fig.7 Flow field of mixed air supply mode (constant value signal)

2.2 正弦信號送風客艙內(nèi)氣流特征

正弦信號送風情況下，送風速度隨時間不斷變化，天花板送風模式下流場如圖8 所示。在t=140 s時刻，速度處于高速階段，渦旋處于中間乘客位置處；隨著送風速度減小，當t=160 s 時，渦旋逐漸向上偏移，且半徑逐漸減?。划攖=180 s 時，渦旋處于行李架附近。側(cè)壁送風模式下流場如圖9 所示。t=140 s 時在靠近過道乘客位置處存在2 個對稱的渦旋；隨著送風速度減小，當t=160 s 時，客艙中間向下的氣流動力減弱，渦旋半徑變小，并且向下發(fā)生偏移；t=180 s 時，在乘客呼吸區(qū)形成新的渦旋。混合送風模式下流場如圖10 所示，當t=140 s 時，送風速度大，在中間乘客頭頂位置形成2 個對稱的渦旋；隨著送風速度減小，當t=160 s 時，在客艙過道處向上的氣流動力減小，導致渦旋向下運動；當t=180 s 時，渦旋到達客艙底部，半徑變小，同時在客艙上部行李架附近形成新的渦旋，新的渦旋繼續(xù)向下移動，整個過程不斷循環(huán)。隨著時間增加，流場的變化導致渦旋位置發(fā)生變化，渦旋中心空氣滯留區(qū)被打破，客艙內(nèi)空氣流動性得到增強。

圖8 天花板送風模式下的流場（正弦信號）Fig.8 Flow field of ceiling air supply mode (sinusoidal signal)

圖10 混合送風模式下的流場（正弦信號）Fig.10 Flow field of mixed air supply mode (sinusoidal signal)

2.3 客艙內(nèi)NO2 平均質(zhì)量濃度變化

為了分析NO2在客艙內(nèi)的變化情況，采集了不同時刻整個客艙內(nèi)的NO2平均質(zhì)量濃度，圖11 為2 個周期內(nèi)(60～180 s)6 種工況下的變化情況，正弦信號送風下的NO2質(zhì)量濃度下降更快。在t=180 s時，天花板送風模式中正弦信號送風比恒值信號送風下的NO2質(zhì)量濃度降低了7.95%，側(cè)壁送風模式降低了6.51%，混合送風模式降低了23.3%，可以看出，混合送風模式下NO2質(zhì)量濃度下降最多，說明混合送風模式中正弦信號送風排污效果最明顯。

圖11 客艙內(nèi)部NO2 質(zhì)量濃度Fig.11 Concentration of NO2 in cabin

2.4 乘客呼吸區(qū)截面NO2 分布

選取乘客呼吸區(qū)高度處的水平截面研究NO2分布情況，圖12 為t=180 s，風速為恒值信號時，天花板送風、側(cè)壁送風和混合送風3 種送風模式下的NO2分布。可以看出，不同送風模式下的NO2分布存在較大差異。天花板送風模式如圖12(a)所示，在第3 排右側(cè)區(qū)域存在小面積NO2聚集現(xiàn)象，兩側(cè)靠近壁面位置處NO2質(zhì)量濃度較低，且座椅與座椅之間的NO2質(zhì)量濃度高于過道區(qū)域；側(cè)壁送風模式如圖12(b)所示，整個截面內(nèi)NO2質(zhì)量濃度相差較小，說明在側(cè)壁送風模式下，NO2分布較為均勻；混合送風模式如圖12(c)所示，在第4 排和第5 排區(qū)域NO2質(zhì)量濃度較高，同時兩側(cè)靠近艙壁區(qū)域的NO2質(zhì)量濃度較低，這是因為新鮮氣流沿著客艙壁面向上流動，壁面周圍的NO2可以更快地得到稀釋。整體來看，在恒值信號下，天花板送風和混合送風模式都存在NO2聚集現(xiàn)象，這是由于客艙內(nèi)氣流循環(huán)基本處于穩(wěn)定狀態(tài)，存在大面積的空氣滯留區(qū)，滯留區(qū)內(nèi)NO2長期積聚，很難排出艙外。

圖12 t=180 s 時恒值信號下呼吸區(qū)截面NO2 分布Fig.12 Distribution of NO2 in respiratory area under constant signal at t=180 s

圖13 為t=180 s，風速為正弦信號時，3 種送風模式下的NO2分布。天花板送風模式如圖13(a)所示，在2 排座椅之間存在大量低質(zhì)量濃度區(qū)域，座椅附近區(qū)域NO2質(zhì)量濃度低于過道區(qū)域；側(cè)壁送風模式如圖13(b)所示，NO2分布較為均勻，過道處NO2質(zhì)量濃度低于座椅間質(zhì)量濃度；在混合送風模式下，如圖13(c)所示，過道側(cè)乘客附近NO2質(zhì)量濃度較低，大部分NO2主要集中于座艙尾部。

圖13 t=180 s 時正弦信號下呼吸區(qū)截面NO2 分布Fig.13 Distribution of NO2 in respiratory area under sinusoidal signal at t=180 s

對比圖12 和圖13，3 種送風模式呼吸區(qū)截面NO2分布情況，將恒值信號送風改為正弦信號送風后，NO2的質(zhì)量濃度均有不同程度的下降，其中混合送風模式下效果最為顯著，說明正弦信號送風更有利于為乘客創(chuàng)造健康的呼吸區(qū)環(huán)境。

綜合分析表明，在正弦信號送風的波動風速下，客艙內(nèi)空氣混合增強，波動的送風速度打破了恒值信號送風下存在的空氣滯留區(qū)，有利于NO2的排出，無論是對于整個客艙內(nèi)部，還是對于乘客的呼吸區(qū)，正弦送風都起到了很好的排污效果，其中正弦信號下混合送風模式的排污效果最好。

3 客艙空氣齡計算

為了更好地評價不同送風信號在客艙內(nèi)的通風換氣效果，使用空氣齡作為評價指標?？諝恺g指空氣進入房間的時間，某一點的空氣齡越小，說明該點的空氣越新鮮，空氣品質(zhì)就越好。客艙內(nèi)不同位置的平均空氣齡 τp計算表達式為[3]

式中：Cp(t)為所測點連續(xù)采集的污染物質(zhì)量濃度值；Cp(0)為初始時刻采集點處污染物的質(zhì)量濃度值。換氣效率 η可用于評價通風量的極限值，一般認為當空氣齡低于換氣效率 η時，可以認為此點的通風系統(tǒng)通風量充足。η表達式為[3]

式中：Vcabin為 客艙內(nèi)部體積；vventilation為客艙的送風量。

選取客艙第3 排的6 名乘客人嘴處為采樣點，計算空氣齡，圖14 為采樣點在6 種送風工況下的空氣齡。在相同送風模式下，采樣點處恒值信號下的平均空氣齡均高于正弦信號，說明具有波動風速的正弦信號，無論在哪種送風模式下，提供的新風都更容易到達乘客面前，為乘客提供更加新鮮的空氣。在正弦信號混合送風模式下的所有采樣點處的空氣齡均小于換氣效率，平均空氣齡約為118 s，說明這種送風工況下客艙內(nèi)通風充足。因此，在通風量相同的情況下，正弦信號混合送風模式的換氣效果最佳。

圖14 不同送風工況下的空氣齡Fig.14 Air age of different air supply modes

4 客艙熱舒適性分析

由于正弦送風會導致送風速度發(fā)生較大的變化，風速較大波動會造成人體不適，對不同送風工況下的吹風感指數(shù)進行計算，當吹風感指數(shù)DR<20%時，認為不會影響乘客舒適性體驗[18]，其表達式為

式中：ti為局部空氣溫度，℃；vi為局部空氣速度，m/s；Tu為局部湍流強度，%。

由于正弦送風速度在一個周期內(nèi)不斷變化，在正弦信號速度達到高速階段的情況下，計算乘客周圍的吹風感指數(shù)。截取5 排乘客周圍共5 個截面計算吹風感指數(shù)，圖15 為各個截面在不同工況下的吹風感指數(shù)。當采用正弦信號送風時，客艙內(nèi)乘客附近的吹風感指數(shù)均有明顯提高，其中天花板送風模式DR>20%；側(cè)壁送風模式和混合送風模式在正弦信號送風下的吹風感指數(shù)雖然高于恒值信號送風，但均小于20%，符合熱舒適性要求。

圖15 乘客周圍截面吹風感指數(shù)Fig.15 DR of cross section around passengers

5 結(jié) 論

1) 正弦送風有助于加強客艙內(nèi)空氣混合，波動的風速會導致渦旋的位置發(fā)生變化，打破渦旋內(nèi)的空氣滯留區(qū)，增強客艙內(nèi)空氣的流動性。

2) 在天花板送風、側(cè)壁送風和混合送風3 種送風模式下，正弦信號送風可以加快客艙內(nèi)NO2的排出，有效提高乘客附近空氣的新鮮程度。其中，正弦信號下的混合送風模式，具有最佳的排污效果。

3) 正弦信號下風速的波動變化會造成較大的吹風感，其中天花板送風模式會導致乘客吹風感明顯，但側(cè)壁送風和混合送風模式下吹風感指數(shù)均小于20%,符合熱舒適性要求。