韓云龍,胡永梅,錢付平

(安徽工業(yè)大學(xué)建筑工程學(xué)院,安徽馬鞍山243002)

室內(nèi)空氣品質(zhì)問題,如CO2濃度高或新風不足、揮發(fā)性有機氣體的排放、細(病)菌等更易獲得大家的關(guān)注,尤其對于那些大部分時間在封閉的室內(nèi)工作和生活的人來說,更是如此,可是大家往往忽略了空氣中可吸入顆粒物濃度的影響。對于中央空調(diào)系統(tǒng)而言,如果通風系統(tǒng)只是裝備了低效或中效過濾器,那么室內(nèi)可吸入顆粒物濃度將會受到影響。室外空氣中的灰塵會通過過濾器進入送風系統(tǒng),其中一部分會沉積在管道內(nèi),而送風管道內(nèi)的溫濕度非常適合病菌、微生物的生存繁殖,在空調(diào)系統(tǒng)啟動過程中或管道振動則可能會促使氣流將這些附有病菌的積塵重新?lián)P起并送至室內(nèi)。若人們長期暴露于此環(huán)境下,呼吸道疾病或病態(tài)建筑綜合癥將不可避免。因此,諸多研究者[1-5]對管道內(nèi)顆粒沉積進行了相關(guān)的研究工作,他們的研究結(jié)果表明風速、顆粒粒徑、管壁表面粗糙度等對顆粒物沉積具有重要的影響,一般而言,隨風速、粒徑的增加,顆粒沉積速率提高。對于水平的通風管道,氣流中的顆粒物在重力沉積作用下部分沉積于管道底部,同時慣性和湍流擴散也會加速顆粒向管道周壁的沉積。

對于空調(diào)房間,為了維持室內(nèi)一定的溫、濕度環(huán)境,要進行一個熱濕交換過程,經(jīng)過熱濕處理后的空氣經(jīng)由通風管道送往空調(diào)房間。因此,若管道保溫措施不當,則氣流與管壁間的溫差產(chǎn)生的熱泳力會加速顆粒的沉積,一些研究者[6-8]認為即使氣流與管壁之間存在一個小的溫差也會加速顆粒沉積速率。然而,前人[9-12]的研究多基于小斷面通風管道內(nèi)的顆粒物沉積,同時空氣濕度變化對通風管道內(nèi)顆粒物沉積的影響還未有報道。該文采用雷諾應(yīng)力模型(RSM)模擬管道內(nèi)湍流,在接近實際空調(diào)通風管道尺寸斷面的管道內(nèi)應(yīng)用拉格朗日隨機軌道模型對完全發(fā)展湍流中顆粒物沉積進行模擬研究,對于管道內(nèi)氣流與管壁溫差所產(chǎn)生的熱泳力對顆粒物沉積的影響進行探討,同時研究空氣相對濕度對顆粒物沉積的影響。

1 模型

1.1 流體與顆粒運動方程

通風管道內(nèi)空氣與顆粒之間的運動屬于典型的氣固兩相流動問題,目前,對此問題的解決方法主要就是歐拉法與拉格朗日法。拉格朗日法可追蹤單個顆粒的運動軌跡,但是需要花費較多的計算時間。一般而言,空調(diào)通風管道內(nèi)的灰塵濃度較低,由于過濾器的裝備,粒徑也較小,因此,空氣與灰塵可分別處理為連續(xù)相與分散相。通常來說,由于灰塵粒徑較小及濃度較低,連續(xù)相對分散相有重要的影響,而灰塵顆粒對連續(xù)相的影響可以被忽略,所以采用拉格朗日法的單相耦合計算方法。

空調(diào)通風管道內(nèi)湍流流動的模擬計算應(yīng)用商業(yè)軟件FLUENT 6.1,一些研究者[13-15]采用雷諾平均的Navier-Stokes(RANS)方程,即k-ε雙方程模擬不可壓縮湍流流動。Tian和Ahmadi[16]于湍流管道內(nèi)納米和微米級的顆粒物沉積進行了模擬研究,認為RSM湍流模型及“two-layer”邊界條件的使用,能合理的預(yù)測顆粒物的沉積。模擬結(jié)果表明0.01～50μm范圍的無因次顆粒沉積速率呈現(xiàn)“V”形分布,因此,采用 RSM 模型和“two-layer”模型的邊界處理條件。

1.2 沉積速率及熱泳力

為了與前人的研究結(jié)果進行比較,采用了無因次顆粒沉積速率,通常表示為：

其中Cm為時均顆粒濃度;u*為磨擦速度,可寫為下式：

其中U m為平均風速;f為范寧系數(shù),采用W hite[17]得出的關(guān)聯(lián)式計算 f。

其中k為壁面粗糙度,光滑壁面為0;Re為雷諾數(shù);D為水力直徑。

無因次松弛時間τ+,如下示：

其中μ和ν分別為流體的動力粘度與運動粘度;ρp和dp分別為顆粒密度和直徑;Cc為坎寧漢系數(shù),如下示：

其中λ為氣體分子自由程,在1大氣壓,25℃時為0.065μm。

如果氣流與管壁之間存在溫差,顆粒物會沿著溫度降低的方向運動,即為熱泳現(xiàn)象。顆粒所受熱泳力可表示為如下：

2 計算域

于矩形斷面0.3 m×0.2 m,長度為3m接近實際空調(diào)通風管道尺寸的管道內(nèi)模擬完全發(fā)展湍流顆粒沉積。采用網(wǎng)格生成工具Gambit生成六面體結(jié)構(gòu)網(wǎng)格;對于管道內(nèi)顆粒沉積而言,邊壁處高精度網(wǎng)格的生成十分重要,因此采用Tian和Ahmadi[16]的邊界層劃分方法,即第一個網(wǎng)格距離邊壁處0.05 mm,其后以1.2的倍率增加直至1～1.2個單位長度為止,近壁處網(wǎng)格較為密集,見圖1,整個計算域共劃分150 000個網(wǎng)格單元。

圖1 矩形斷面網(wǎng)格劃分

應(yīng)用拉格朗日隨機軌道模型追蹤離散相顆粒的運動軌跡,顆粒密度為900 kg/m3,個數(shù)為3 000個。為了和以前研究者的實驗或計算結(jié)果進行比較,對于粒徑范圍0.01～50μm的球形顆粒沉積進行了模擬計算,入口風速為5 m/s,溫度15℃,管道出口邊界條件為“outflow”,管壁熱絕緣并認為水力光滑,為“trap”;忽略顆粒間的合并、碰撞及與壁面的反彈。

3 結(jié)果與討論

3.1 流場分析

為了和前人的研究結(jié)果進行比較,對小斷面尺寸(0.075m×0.05m)的通風管道內(nèi)低雷諾數(shù)流動及顆粒沉積也進行了模擬計算,流場見圖2。2種斷面尺寸的管道進口風速均為5 m/s,由于斷面尺寸的減小,相應(yīng)的雷諾數(shù)也減小,分別為 82 135、5 136。由圖2可以看出,雖然進風風速相同,但斷面尺寸的減小導(dǎo)致雷諾數(shù)大幅降低?？梢钥闯龃髷嗝婀艿纼?nèi)流速梯度較小,斷面流速較為均勻;而小斷面管道流速梯度較大,斷面越小,管道流動核心區(qū)流速越大,而近壁區(qū)域流速明顯低于大斷面管道處的流速。而前人如Tian等[16]主要基于小斷面管道內(nèi)顆粒沉積的研究,由此可知這和實際通風管道內(nèi)流體流動情況稍有差距。

3.2 模型驗證

圖3為無因次沉積速率隨無因次松弛時間的變化曲線,塵粒粒徑范圍為0.01～50μm,并將模擬結(jié)果與前人數(shù)據(jù)進行了比較。由于粒徑不同,所受擴散力、慣性力大小也不同,致使顆粒沉積運動呈現(xiàn)3個區(qū)：即擴散區(qū)、過渡區(qū)和慣性區(qū)。對于次微米及納米級顆粒而言,顆粒所受重力很小,慣性力可以忽略,而擴散力則是顆粒沉積的主要機理。對于大粒徑顆粒,重力及慣性力是顆粒沉積的主要作用力。由圖3可以看出,模擬結(jié)果基本符合前人的實驗數(shù)據(jù)及計算結(jié)果,無因次沉積速率隨無因次松弛時間呈現(xiàn)“V”形變化,滿足顆粒沉積的擴散區(qū)、過渡區(qū)及慣性區(qū)分布規(guī)律。但是計算域為矩形斷面(0.3 m×0.2m)的三維計算,更接近于實際空調(diào)管道的斷面尺寸,雷諾數(shù)(82 135)也要遠高于前人的實驗或計算中采用的雷諾數(shù)[16,19](Tian等[16]采用的雷諾數(shù)為6 667)。但雷諾數(shù)的提高不是由風速的提高而引起,而僅由斷面尺寸的增加導(dǎo)致的。一般而言,隨風速的提高,顆粒的沉積速率會相應(yīng)提高,但由斷面尺寸增加而導(dǎo)致雷諾數(shù)提高對顆粒沉積的影響并不顯著。

圖2 不同通風管道斷面流速場

圖3 模型計算結(jié)果的比較驗證

3.3 溫度對顆粒沉積的影響

對于實際空調(diào)系統(tǒng)通風管道而言,送風溫度、濕度因季節(jié)的變化而不同,因主氣流溫度與管壁的溫差而導(dǎo)致的熱泳力將會促進次微米級小微粒的沉積。有關(guān)熱泳力對顆粒沉積的影響,一些研究者曾加以研究,如Yang等[11]認為PM 2.5的沉積效率主要受主氣流與冷壁面的溫差影響,熱泳力是導(dǎo)致PM 2.5沉積的關(guān)鍵因素。He和Ahmadi[20]認為只有顆粒進入邊界層的粘性底層時,布朗擴散才是顆粒沉積的控制機理,熱泳力對小顆粒(如幾個微米級或以下)的傳輸、沉積有重要影響。并且沉積速率及粒徑范圍皆隨著主氣流與壁面溫差的提高而增加。其它研究者[21-22]也都認為溫差導(dǎo)致的熱泳力能促進顆粒的沉積。該研究針對1μm的顆粒探討主氣流與管壁之間溫差對顆粒沉積的影響。圖4為入口空氣溫度24℃時不同壁面溫度所形成主氣流與冷壁面的溫差對顆粒沉積的影響,溫度梯度均為沿著壁面的法線方向,溫度向管壁處降低。由圖4可以看出,當固定入口溫度時,隨著管壁溫度的提高,顆粒沉積速率有下降的趨勢,無論底面、側(cè)壁還是頂壁面都是如此趨勢,表明了主氣流與壁面溫差所產(chǎn)生的熱泳力促進了顆粒的沉積。同時可以看出,側(cè)壁顆粒沉積速率同頂、底面相比,其沉積速率略小,底面的沉積速率最大。然而對于實際空調(diào)系統(tǒng)通風管道而言,都會有保溫措施,管道壁面可視為熱絕緣,主氣流與管壁之間的溫差很小,那么熱泳力對顆粒沉積所起的作用可以被忽略。

圖4 熱泳力對顆粒沉積的影響

然而當入口空氣溫度提高,而管壁保持絕緣時,發(fā)現(xiàn)不同之現(xiàn)象,如圖5示?？梢园l(fā)現(xiàn)隨著入口溫度的提高,底面、側(cè)壁的顆粒沉積速率均略呈下降趨勢,而頂面的顆粒沉積速率則呈現(xiàn)升高趨勢。雖然設(shè)立了管壁絕緣的邊界條件,但管道頂部的氣流溫度略高于管道底部,會產(chǎn)生一個升力。因此,溫差產(chǎn)生的升力作用使較多的顆粒分散于管道空間的上部,所以,管道頂部的顆粒沉積速率有升高趨勢,底部及側(cè)壁的顆粒沉積速率則略呈下降趨勢。另外,隨著溫度的升高,空氣的粘度增加而導(dǎo)致管道內(nèi)氣體流動的雷諾數(shù)有所減小,湍動能有所降低致使顆粒的湍流擴散作用減弱也有關(guān)。

圖5 管道氣流溫度對顆粒沉積的影響

3.4 濕度對顆粒沉積的影響

空調(diào)系統(tǒng)為了滿足室內(nèi)一定的濕度要求,在新風處理單元一般需要加濕過程,經(jīng)過加濕處理后的空氣,其中的顆粒沉積是否有影響尚需研究。在該研究中,于入口空氣16℃時,對通風管道內(nèi)空氣相對濕度40%～80%的環(huán)境條件顆粒沉積進行了模擬研究,見圖6。空氣相對濕度的提高會改變空氣的粘度、密度、熱導(dǎo)率等物性,將空氣相對濕度改變所引起的物性變化計算出來,在模擬計算過程中,認為每個物性都對應(yīng)著一個相對濕度。由圖6可以看出,隨著氣流相對濕度的提高,顆粒沉積速率均相應(yīng)增加。這是由于空氣相對濕度的提高,濕空氣的粘度也會相應(yīng)提高,增加了湍流邊界層的厚度,使得更多的顆粒進入這個低速的粘性邊界層,因此提高了顆粒沉積速率。該研究只是對于絕緣管壁的情況進行的計算,在實際運行中,倘若管壁保溫措施不當,則較高濕度空氣中的水會凝結(jié)于管壁,這就會相應(yīng)提高管壁的粗糙度,按照Lai等[3]的研究結(jié)果可知,粗糙度的提高會促進顆粒的沉積,所以沉積速率會有相應(yīng)的提高,而潮濕的管道內(nèi)環(huán)境,適宜的溫度則為真菌、細菌等提供了繁殖的溫床,對室內(nèi)空氣品質(zhì)會有嚴重的影響。

圖6 管道內(nèi)氣流濕度對顆粒沉積的影響

4 結(jié)論

應(yīng)用拉格朗日隨機軌道模型模擬分析了通風管道內(nèi)溫濕度變化對顆粒物沉積的影響,得出了如下主要結(jié)論：

1)集中式空調(diào)系統(tǒng)通風管道內(nèi)高溫送、回風氣流與低溫管壁形成較大溫差時,熱泳力會加速小顆粒(亞微米、納米級顆粒)的沉積速率,對于實際工程應(yīng)用,若管道具有良好的保溫措施,熱泳力對顆粒沉積的作用則可以被忽略。

2)通風管道內(nèi)送、回風濕度的提高會加速顆粒的沉積速率,對于缺乏良好保溫措施的通風管道,若濕度較大,管壁上凝結(jié)水的出現(xiàn)會增加粗糙度,更會加速顆粒物的沉積。

3)集中式空調(diào)系統(tǒng)通風管道內(nèi)顆粒物的沉積及其衛(wèi)生狀況對室內(nèi)空氣品質(zhì)的影響已逐漸引起重視,若通風管道維護不當,送回風溫濕度的變化會影響顆粒物的沉積,并惡化室內(nèi)空氣品質(zhì)。通過文章的分析研究,有助于空調(diào)運行管理人員提高通風管道內(nèi)顆粒物沉積的認識,促進其對空調(diào)系統(tǒng)通風管道的維護、清潔工作,同時可為評價室內(nèi)人員對顆粒物的接觸、暴露情況提供參考。

[1]CH EONG K W.Deposition of aerosol particles in ductw ork[J].App lied Energy,1997,57(4)：253-261.[2]WU J,ZHAO B.Effect of ventilation ductas a particle filter[J].Building and Environment,2007,42：2523-2529.

[3]LA I A C K,BYRNE M A,GODDARD A J H.Enhanced particle loss in ventilation duct with ribbed surface[J].Building and Environment,2000,35：425-432.

[4]ZHAO B,WU J.Modeling particle deposition from fully developed turbulent flow in ventilation duct[J].A tmospheric Environment,2006,40：457-466.

[5]ZHANG H,AHMADIG.Aerosol particle transport and deposition in vertical and horizontal turbu lent duct flow s[J].Journal of Fluid Mechanics,2000,406：55-80.

[6]BAE G N,LEE C S,PARK S O.M easurements and control o f particle deposition velocity on a horizontal w afer w ith therrnophoretic effect[J].Aerosol Science and Technology,1995,23：321-330.

[7]WIRZBERGER H,LEKH LANAKHER S,SHAPIRO M,et al.Prevention o f particle deposition by means of heating the deposition surface[J].Journal of Aerosol Science,1997,28(S1)：S83-S84.

[8]TSA I R,LIANG L J.Correlation for thermophoretic deposition o f aerosol particles onto cold plates[J].Journal of Aerosol Science,2001,32：473-487.

[9]CHANG Y C,RANADE M B,GENTRY J W.Thermophoretic deposition in flow along an annu lar cross-section：experiment and simu lation[J].Journal of Aerosol Science,1995,26(3)：407-428.

[10]LUO X W,YU S Y.Deposition of particles in turbulent pipe flow[J].China Particuology 2006,4(1)：31-34.

[11]YANG R C,LIU R L,ZHOU T,et al.Kinematic characteristics and thermophoretic deposition of inhalable particles in turbulent duct flow[J].Chinese Journal o f Chem ical Engineering,2008,16(2)：192-197.

[12]DEH BI A.A stochastic Langevin model of turbu lent particle dispersion in the presence of thermophoresis[J].Internationa l Journal of M ultiphase Flow,2009,35：219-226.

[13]ZHANG Z,CHEN Q.Experimentalmeasurements and numerical simulations of particle transport and distribution in ventilated rooms[J].A tmospheric Environment,2006,40：3396-3408.

[14]ZHANG J P,LIA G.Study on particle deposition in vertical square ventilation duct flow s by different models[J].Energy Conversion and Management,2008,49：1008-1018.

[15]LA IA C K,CHEN F Z.M ode ling particle deposition and distribution in a chamber with a tw o-equation Reynolds-averaged Navier–Stokesmodel[J].Journal of Aerosol Science,2006,37：1770-1780.

[16]TIAN L,AHMADIG.Particle deposition in turbu lent duct flow s-comparisons of different model predictions[J].Journal of Aerosol Science,2007,38：377-397.

[17]WH ITE F M.Fluid Mechanics[M].2th ed.New York：M cGraw-H ill,1986.

[18]TALBOT L,CHENG R K,SCHEFER R W,et al.Thermophoresis of particles in a heated boundary layer[J].Journal of Fluid M echanics,1980,101：737-758.[19]FAN F,AHMADIG.A sub layermodel for turbu lent deposition of particles in verticalductswith smooth and rough surfaces[J].Journal of Aerosol Science,1993,24：45-64.

[20]HE C H,AHMADI G,Particle deposition w ith thermophoresis in lam inarand turbulentduct flow s[J].Aerosol Science and Technology,1998,29(6)：525-546.

[21]ROM AY F J,TAKAGAK I S S,PU I D Y H,et al.Thermophoretic deposition of aerosol particles in turbulent pipe flow[J].Journal of Aerosol Science,1998,29(8)：943-959.

[22]LINA JS,TSAIA C J,CHANG C P.Supp ression of particle deposition in tube flow by thermophoresis[J].Journalof Aerosol Science,2004,35：1235-1250.