程家磊,錢付平,魏舒婷,唐蓮花,姜榮賀,肖鵬程

(1.安徽工業(yè)大學(xué)建筑工程學(xué)院，安徽馬鞍山243032；2.河北億利橡塑集團有限公司，河北清河054800)

貨車排放的氮氧化物(NOx)和細(xì)顆粒物(PM)明顯高于客車，其中重型貨車是主要貢獻(xiàn)者，因此如何降低重型載貨汽車的污染物排放量，改善燃油經(jīng)濟性是目前節(jié)能減排的重要控制壞節(jié)?！皣鍢?biāo)準(zhǔn)”對重型載貨汽車生產(chǎn)廠家提出了新的較高要求，如何提高重型載貨汽車的高動力性、低排放性成為重型卡車行業(yè)面臨的主要任務(wù)。實現(xiàn)減排的重點在于降低NOx和PM的排放，而空氣濾清器是保證進入發(fā)動機氣缸內(nèi)部空氣為高潔凈度的主要部件[1-2]，可以為發(fā)動機提供清潔的空氣，使燃油燃燒更加充分。合理地組織燃燒是降低NOx和PM排放的主要措施，它依靠提高充氣效率及渦流比來實現(xiàn)燃料的完全燃燒，以此降低NOx和PM的排放。同時，空氣濾清器的空氣凈化能力對車輛發(fā)動機的動力性、使用壽命及經(jīng)濟性等有著直接影響。因此，如何在保證過濾效率的前提下，提高空氣濾清器系統(tǒng)內(nèi)部流體流動均勻性，降低其壓力損失是提升空氣濾清器性能的關(guān)鍵[3-4]。而作為汽車進氣系統(tǒng)一部分的前進氣道對空氣濾清器性能起著獨特的重要作用，它能夠進行預(yù)過濾，將大顆?；壹八值阮A(yù)先濾除，在一定程度上提高空氣濾清器的過濾效率和使用壽命。隨著計算機技術(shù)的發(fā)展，利用CFD軟件對汽車空氣濾清器及進氣道進行研究越來越多[5-8]，如韓青等[9]基于數(shù)值模擬對空氣濾清器結(jié)構(gòu)進行優(yōu)化，優(yōu)化后的結(jié)構(gòu)具有出口處的負(fù)壓較小和壓力分布均勻等優(yōu)點。袁志群等[10]對空氣濾清器的進氣阻力進行分析，得出了進氣阻力隨著流量的增加呈拋物線趨勢增加的規(guī)律。劉鵬飛等[11]通過實驗與數(shù)值模擬相結(jié)合的方法對干式空氣濾清器內(nèi)部阻力特性進行研究，得到了流阻特性隨流量的變化關(guān)系。Liu等[12]對汽車空氣濾清器聲學(xué)性能進行優(yōu)化設(shè)計，并用數(shù)值模擬方法對新型空氣濾清器的聲學(xué)特性進行預(yù)測，結(jié)果表明新型空氣濾清器聲學(xué)性能良好。He等[13]對某空氣濾清器及其管道開展數(shù)值模擬，并研究采用插入管結(jié)構(gòu)、內(nèi)置擋板結(jié)構(gòu)后的流動特性，得出了插入管結(jié)構(gòu)降低流動損失，而擋板結(jié)構(gòu)提高流動損失。Jiang等[14]對進氣道設(shè)計參數(shù)進行數(shù)值模擬，并與實驗結(jié)果進行對比，結(jié)果表明數(shù)值模擬方法可為氣道設(shè)計提供參考。Li等[15]通過數(shù)值模擬對汽車進氣道的空氣動力學(xué)噪聲進行分析，并通過與實驗對比驗證了數(shù)值模擬方法對汽車進氣道的空氣動力學(xué)噪聲預(yù)測的可行性。王治林等[16]通過實驗數(shù)據(jù)對空氣濾清器的數(shù)值模型進行驗證，結(jié)果表明模擬結(jié)果和實驗結(jié)果有很好的一致性，但該研究未對進氣管道的數(shù)值模型進行驗證。目前，對于汽車進氣系統(tǒng)的研究主要集中在空氣濾清器的結(jié)構(gòu)優(yōu)化以及壓力損失分析上，而對汽車進氣道的研究則集中在進氣道的參數(shù)設(shè)計以及進氣道的噪聲分析上，未考慮重型車空氣濾清器前進氣道的格柵角度及擋雨片結(jié)構(gòu)等對進氣道壓力損失的影響。

鑒于此，本文采用商用計算流體力學(xué)軟件建立重型車空氣濾清器前進氣道數(shù)值計算模型，并通過實驗驗證數(shù)值模型的可行性，在此基礎(chǔ)上，同時改變進氣道內(nèi)格柵夾角和擋雨片大小及數(shù)量，研究其對前進氣道流場的影響。

1 數(shù)值計算模型及其驗證

1.1 三維模型的建立及網(wǎng)格生成

采用專業(yè)建模軟件SolidWorks建立空氣濾清器前進氣道幾何模型，前進氣道的氣流先經(jīng)氣道入口進入，再經(jīng)過格柵及其后面的擋雨片，最后由出口流出，其結(jié)構(gòu)如圖1所示，格柵的側(cè)面與坐標(biāo)平面XOY的空間夾角為格柵角度。本文采用ICEM的自動網(wǎng)格生成軟件對圖1所示的流場區(qū)域進行網(wǎng)格劃分，所生成的網(wǎng)格如圖2所示。

為了保證流場分析的精度和效率，在某些自動生成網(wǎng)格不理想的區(qū)域進行局部修正，最終得到的計算網(wǎng)格總數(shù)約為200萬。

1.2 數(shù)學(xué)模型

為便于計算，提出以下假設(shè)：氣體是不可壓縮流體；工作過程中溫度不變；入口處氣流速度分布均勻[17]。氣相控制方程(連續(xù)性方程及動量方程)為

圖1 前進氣道幾何模型Fig.1 Geometric model of forward intake port

圖2 前進氣道網(wǎng)格模型Fig.2 Grid model of the forward intake port

式中：p為靜壓，Pa；t為時間，s；u為氣相速度，m/s；ρ為氣相密度，kg/m3；x為方向向量；gi和Fi分別為i方向上重力體積力和外部體積力，N；τij為應(yīng)力張量，Pa。

1.3 邊界條件

模擬物性參數(shù)及其邊界條件如表1所示，其中前進氣道的進出口及壁面邊界條件的設(shè)置如圖3所示。

表1 邊界條件Tab.1 Boundary conditions

1.4 湍流模型及求解器設(shè)置

模擬時采用不可壓縮、三維穩(wěn)態(tài)流動數(shù)學(xué)模型和k-ε雙方程湍流模型[18]，收斂條件為所有殘差均低于10-3。穩(wěn)態(tài)湍流流場的計算采用SIMPLE算法。計算均在CPU為Intel?Xeon?E5-26702.8 GHz、內(nèi)核為8、內(nèi)存為64 GB的小型工作站上完成。

1.5 方案驗證

根據(jù)汽車行業(yè)標(biāo)準(zhǔn)“汽車用空氣濾清器試驗方法”[19]，在實驗臺上對空氣濾清器前進氣道進出口的壓力損失進行實驗測試，以驗證本文數(shù)值計算方法的可行性與準(zhǔn)確性。試驗臺主要由抽氣設(shè)備、空氣流量測量系統(tǒng)、空氣流量控制系統(tǒng)、靜壓和壓力降測量系統(tǒng)、空氣濾清器等組成，如圖4所示。

根據(jù)壓力損失的定義，當(dāng)空氣流過總成之時，壓力損失就是其上、下游規(guī)定測壓點處管道內(nèi)氣流的全壓差。因此，在確定壓力損失時，先確定被試總成上、下游測壓點之間的全壓值。

圖5為不同流量下前進氣道進出口壓差值的實驗值與基于上述數(shù)值計算模型的模擬值，由圖5可以看出，模擬計算結(jié)果與實驗結(jié)果的總體趨勢是一致的，最大誤差在15%以內(nèi)。因此可知模擬結(jié)果可信的，本文模擬計算方法可以用于預(yù)測空氣濾清器前進氣道的壓力損失。

圖3 前進氣道邊界條件的設(shè)置Fig.3 Boundary conditions of the forward intake port

圖4 進氣總成綜合實驗裝置Fig.4 Device of air intake assembly comprehensive experimental

圖5 不同流量下前進氣道壓力損失Fig.5 Pressure drop of forward intake port at different flow rates

2 實驗結(jié)果與分析

2.1 格柵角度和擋雨片對壓力損失的影響

模型的格柵角度分為30°，45°，60°3種，模型分為加擋雨片和不加擋雨片2種方式，組合之后，總共有6種模型。根據(jù)空氣濾清器前進氣道的6種模型，分別對其進行數(shù)值計算，得到前進氣道的流場信息，其速度云圖如圖6。從圖6可以看出：不同格柵角度的前進氣道速度變化規(guī)律基本相同，且在格柵處速度較大；當(dāng)加擋雨片時，在擋雨片處速度也較大，主要是流通截面發(fā)生變化，導(dǎo)致流速的增加；在貼近前進氣道右側(cè)壁面處的速度較大，其主要是由于氣體在格柵處及擋雨片處速度突然增大，進入前進氣道后會形成氣流漩渦；此外，無論是否加擋雨片，在進氣道出口處右側(cè)的速度總是大于左側(cè)的速度，且在加了擋雨片之后更加明顯，主要是因為加了擋雨片之后，氣道的流通截面變小速度增大。速度的突然增大對前進氣道內(nèi)部的流場造成了擾動，因此未加擋雨片前進氣道的速度均勻性比加擋雨片的前進氣道速度的均勻性要好。

圖6 不同前進氣道的速度云圖Fig.6 Velocitydistributionofdifferentforwardintakeports

圖7為上述6種模型的壓力云圖。從圖7可以看出整個流場的壓力梯度變化較大，壓力損失主要集中在前進氣道的格柵和擋雨片兩個區(qū)域，這與速度場的分布相對應(yīng)，而這兩個區(qū)域的壓力損失主要是由于模型結(jié)構(gòu)的變化而導(dǎo)致氣流方向和流速的變化造成的。從圖中可知前進氣道的壓力梯度變化從上到下呈現(xiàn)由高到低的趨勢，主要原因是：在入口處，由于流場內(nèi)部的流通截面發(fā)生突變，導(dǎo)致流速激增、壓力降低，但隨著氣道的流通截面積增加，使得流體速度逐漸降低。從圖中還能看出未加擋雨片的前進氣道壓力分布比加擋雨片的前進氣道壓力梯度要小，且分布更為均勻，這與速度場的分析是一致的。加擋雨片之后，在前進氣道的出口管段存在負(fù)壓，主要是加了擋雨片之后，使氣流速度增大，所以在前進氣道出口局部區(qū)域出現(xiàn)了負(fù)壓。

圖7 不同前進氣道的壓力云圖Fig.7 Pressuredistributionofdifferentforwardintakeports

圖8是角度和擋雨片對前進氣道壓力損失的影響。由圖8可知，加擋雨片格柵角度為30°時壓力損失最大，為739 Pa；當(dāng)格柵角度從30°增加到60°時，前進氣道內(nèi)壓力損失顯著降低。由此可知，前進氣道格柵角度增大，降低了前進氣道內(nèi)的壓力損失，這主要是由于格柵角度增大，進氣口有效面積增大，氣流速度減小，因而壓力損失減小。對于未加擋雨片的前進氣道同樣遵循此規(guī)律，只是由于未加擋雨片，隨著格柵角度的增大前進氣道的壓力損失減小的幅度更大。同時可以看出無論是否加擋雨片，當(dāng)格柵角度從30°增加到45°時的壓力曲線的斜率要比格柵角度從45°增加到60°時壓力曲線的斜率大，這表明當(dāng)格柵角度從30°增加到45°時的壓力損失要比格柵角度從45°增加到60°時壓力損失減小的幅度更大。從圖中還可以看出對同一角度的進氣格柵，加擋雨片的前進氣道比未加擋雨片的前進氣道的壓力損失要大200 Pa左右，這主要是由于擋雨片改變了氣流的方向與流通截面，導(dǎo)致前進氣道的局部阻力增大。此外，擋雨片對前進氣道壓力損失的影響要比格柵角度對前進氣道壓力損失的影響更大，因為當(dāng)格柵角度減小時壓力損失增大200 Pa左右，而當(dāng)增設(shè)擋雨片時，壓力損失卻增大300 Pa左右，因此有必要對擋雨片的結(jié)構(gòu)進行優(yōu)化，來進一步降低前進氣道的壓力損失。

圖8 角度和擋雨片對前進氣道壓力損失的影響Fig.8 Influence of the angle and rain shield on pressure drop of forward intake ports

2.2 不同擋雨片結(jié)構(gòu)對壓力損失的影響

根據(jù)擋雨片的原始結(jié)構(gòu)A，對其大小及片數(shù)進行改變，得到3種不同結(jié)構(gòu)的擋雨片，如圖9所示。結(jié)構(gòu)B擋雨片的數(shù)量和A一致，而擋雨片的面積逐漸減小。結(jié)構(gòu)C是在結(jié)構(gòu)B的基礎(chǔ)上將最后一個擋雨片開一個矩形口。結(jié)構(gòu)D則是在結(jié)構(gòu)B的基礎(chǔ)上增加了擋雨片的個數(shù)，將擋雨片的片數(shù)由4片增加到了5片。

圖9 不同擋雨片結(jié)構(gòu)Fig.9 Different structures of the rain shield

圖10為上述4種帶有不同擋雨片結(jié)構(gòu)的前進氣道內(nèi)的速度云圖。從圖中可以看出，氣體進入前進氣道后在入口處的速度場分布規(guī)律基本一致，在格柵和擋雨片處由于流通截面積的變化，導(dǎo)致局部速度的增大。從圖中還可以看出帶有結(jié)構(gòu)B擋雨片的前進氣道，其在出口處速度均勻性比其它三種要好，這主要是由于擋雨片面積依次減小，流通截面變大氣流速度減小，流速對進氣道內(nèi)流場的擾動減小，因此在出口處的流速也比較均勻。對于帶有結(jié)構(gòu)C擋雨片的前進氣道，由于在擋雨片上開口，改變了氣流的流向，對進氣道內(nèi)部流場的擾動增大，因此在前進氣道左側(cè)出口處流速局部增大。而對于帶有結(jié)構(gòu)D擋雨片的前進氣道，只有在格柵和擋雨片處速度較大，而在前進氣道的內(nèi)部速度較小，表明擋雨片的片數(shù)對前進氣道的內(nèi)部流速有著較大的影響。

圖10 不同前進氣道的速度云圖Fig.10 Velocity distribution of different forward intake ports

圖11為上述4種帶有不同擋雨片結(jié)構(gòu)的前進氣道內(nèi)的壓力云圖，從圖中可以看出，不同前進氣道的壓力變化規(guī)律基本相同，前進氣道的壓力損失主要集中在格柵和擋雨片處，且前進氣道上側(cè)和下側(cè)壓力梯度較大，在出口管段均出現(xiàn)了負(fù)壓，且在擋雨片處出現(xiàn)了局部負(fù)壓。帶有結(jié)構(gòu)C擋雨片的進氣道在擋雨片處局部負(fù)壓最明顯，主要是由于結(jié)構(gòu)C在擋雨片上開口，使局部流速增大，這與速度場的分析是一致的。

圖11 不同前進氣道的壓力云圖Fig.11 Pressure distribution of different forward intake ports

表2為擋雨片結(jié)構(gòu)對壓力損失的影響。由表2可知擋雨片結(jié)構(gòu)為B時前進氣道壓力損失最小，為727 Pa；而當(dāng)擋雨片結(jié)構(gòu)為D時，前進氣道壓力損失最大，為761 Pa。由此可知，擋雨片數(shù)量對壓力損失的影響比擋雨片大小對壓力損失的影響大。通過結(jié)構(gòu)B和結(jié)構(gòu)C的對比分析可知，在擋雨片上開口會增加進氣量，但相應(yīng)的壓力損失也會增大，主要是在擋雨片上開口之后改變了氣流的流向。通過結(jié)構(gòu)B和結(jié)構(gòu)D的對比分析可知，增加擋雨片的片數(shù)，會使壓力損失增大，主要增加了擋雨片的片數(shù)，使得氣道進氣口的有效面積減小，速度增大，導(dǎo)致局部壓力損失增大。通過結(jié)構(gòu)B和結(jié)構(gòu)A的對比分析可知，保持擋雨片的片數(shù)不變，依次減小擋雨片的大小可以減小壓力損失，主要原因是減小了擋雨片對氣流的阻擋作用。

表2 擋雨片結(jié)構(gòu)對壓力損失的影響Tab.2 Influence of the rain shield structure on pressure drop

3 結(jié) 論

基于CFD對空氣濾清器前進氣道內(nèi)氣相流動進行了數(shù)值模擬，并用實驗數(shù)據(jù)對模型進行了實驗驗證，同時分析了格柵角度和擋雨片對前進氣道壓力損失的影響，得出如下結(jié)論。

1)不同流量下空氣濾清器前進氣道的模擬值與實驗值吻合良好且最大誤差在15%以內(nèi)，表明所用模型和計算方法對空氣濾清器前進氣道的內(nèi)部流場具有良好的預(yù)測精度。

2)對于無論是否加擋雨片的前進氣道，隨著格柵角度的減小前進氣道壓力損失逐漸增大，其中加擋雨片、格柵角度為30°時壓力損失最大，為739 Pa。而對于格柵角度不變的前進氣道，加擋雨片的壓力損失要比未加擋雨片的壓力損失大200 Pa左右。

3)擋雨片結(jié)構(gòu)為B時前進氣道壓力損失最小，為727 Pa。而當(dāng)擋雨片結(jié)構(gòu)為D時，前進氣道壓力損失最大，為761 Pa，表明擋雨片數(shù)量對壓力損失的影響比擋雨片大小對壓力損失的影響大。在保持擋雨片的片數(shù)不變時，依次減小擋雨片的面積可以降低前進氣道壓力損失。