孫麗穎，魏慧嬌

(1．哈爾濱工程大學(xué) 航天與建筑工程學(xué)院，150001哈爾濱，sunliying@hrbeu．edu．cn;2．中國建筑科學(xué)研究院天津分院，300834天津)

孔板送風(fēng)房間內(nèi)污染物分布的模擬與實(shí)驗(yàn)

孫麗穎1，魏慧嬌2

(1．哈爾濱工程大學(xué) 航天與建筑工程學(xué)院，150001哈爾濱，sunliying@hrbeu．edu．cn;2．中國建筑科學(xué)研究院天津分院，300834天津)

為研究孔板送風(fēng)條件下室內(nèi)污染物體積分?jǐn)?shù)的分布規(guī)律，采用實(shí)驗(yàn)測試和數(shù)值模擬兩種方法．建模時(shí)采用基本模型和N點(diǎn)風(fēng)口模型兩種方法進(jìn)行孔板入口邊界條件的處理，模擬結(jié)果表明:兩種處理方法計(jì)算的污染物體積分?jǐn)?shù)分布規(guī)律相近;將實(shí)驗(yàn)數(shù)據(jù)和模擬結(jié)果進(jìn)行對比，其變化趨勢吻合良好，其中N點(diǎn)風(fēng)口模型的模擬結(jié)果更接近實(shí)驗(yàn)情況，從而驗(yàn)證了應(yīng)用該模型研究室內(nèi)污染物分布特性的可行性．

孔板送風(fēng);污染物;數(shù)值模擬

隨著經(jīng)濟(jì)發(fā)展，人們對室內(nèi)裝修的要求不斷提高，各種新型材料在室內(nèi)裝飾裝修中的廣泛采用，導(dǎo)致室內(nèi)空氣品質(zhì)開始惡化，已經(jīng)嚴(yán)重威脅到室內(nèi)人員的身心健康．因此，解決室內(nèi)空氣污染問題日益受到人們的關(guān)注．空調(diào)系統(tǒng)的送風(fēng)對室內(nèi)污染物可以起到去除或稀釋的作用，在送風(fēng)過程中，污染物體積分?jǐn)?shù)的變化規(guī)律與分布情況可以通過實(shí)驗(yàn)與數(shù)值模擬的方法來研究．隨著計(jì)算流體力學(xué)技術(shù)的廣泛應(yīng)用，室內(nèi)空氣品質(zhì)的研究借助于數(shù)值模擬方法，取得了很多成績，李麗等［1］用數(shù)值模擬方法研究了新風(fēng)量對室內(nèi)Total Votatile Organic Compound(簡稱TVOC)體積分?jǐn)?shù)的影響．張明遠(yuǎn)［2］運(yùn)用數(shù)值計(jì)算的方法，分析比較了不同送回風(fēng)方式對室內(nèi)污染物體積分?jǐn)?shù)分布的影響．趙彬等［3］采用三維CFD程序STACH－3對潔凈室的速度場進(jìn)行了數(shù)值模擬，得出的結(jié)論是用“N點(diǎn)風(fēng)口動量模型”描述孔板風(fēng)口可以得到與實(shí)際相符合的模擬結(jié)果．張文勝等［4］對一個(gè)體積為0.45×0.45×0.4(m3)的密閉小室，在送風(fēng)溫度恒定的條件下對小室內(nèi)孔板送風(fēng)氣流分布情況進(jìn)行了實(shí)驗(yàn)研究．李勇［5］利用CFD軟件，對各種通風(fēng)空調(diào)方案下ICU病房內(nèi)溫度場、速度場、體積分?jǐn)?shù)場的數(shù)值模擬結(jié)果進(jìn)行分析比較，得出結(jié)論:百葉送風(fēng)時(shí)病房的排污能力要優(yōu)于孔板送風(fēng)．目前，國內(nèi)外通過數(shù)值模擬分析室內(nèi)氣流組織的研究較多，對孔板送風(fēng)條件下室TVOC體積分?jǐn)?shù)場的CFD模擬分析與實(shí)驗(yàn)驗(yàn)證相對較少，存在的問題是數(shù)值模擬的方法、邊界條件和初始條件設(shè)定的合理性未能很好地用實(shí)驗(yàn)來驗(yàn)證，使得模擬結(jié)果的可信度降低．本文以某一采用孔板送風(fēng)的實(shí)驗(yàn)房間為研究對象，通過實(shí)驗(yàn)測試和數(shù)值模擬的方法，研究孔板送風(fēng)過程中實(shí)驗(yàn)房間內(nèi)污染物TVOC體積分?jǐn)?shù)的分布規(guī)律，在進(jìn)行孔板送風(fēng)口的邊界條件處理時(shí)，采用兩種不同形式的簡化，并將實(shí)驗(yàn)與2種模擬結(jié)果進(jìn)行比較，從而驗(yàn)證所建模型的正確性．

1 數(shù)值模擬基礎(chǔ)及模型建立

1.1 計(jì)算基礎(chǔ)

以哈爾濱工程大學(xué)人工環(huán)境實(shí)驗(yàn)室的空調(diào)房間作為研究對象，房間的物理模型見圖1．房間尺寸為5m×3.5 m×2.8 m，一門一窗位于南墻上．設(shè)實(shí)驗(yàn)房間內(nèi)共有2張桌子和2個(gè)人，根據(jù)實(shí)際情況設(shè)2張桌子是污染物TVOC的主要散發(fā)源．室內(nèi)溫度為18℃，送風(fēng)方式為孔板送風(fēng)，并且整個(gè)頂棚布滿孔板．回風(fēng)口為百葉風(fēng)口，尺寸為0.65 m×0.42 m，位于房間的西墻下方，桌子尺寸為0.7 m×0.4 m×0.6 m．

圖1 房間的物理模型

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

大量實(shí)驗(yàn)表明空調(diào)房間內(nèi)的氣流基本上是湍流，本文采用目前應(yīng)用最廣、接受檢驗(yàn)最多、數(shù)值技術(shù)也最成熟的湍流模型——k－ε二方程模型［6］．采用k－ε二方程模型求解湍流對流換熱問題時(shí)，k－ε方程與動量方程、能量方程、連續(xù)性方程一起構(gòu)成了室內(nèi)空氣流動的控制方程，其通用形式［7］為

式中:φ為通用變量，分別表示動量方程中的速度u，能量方程中的溫度T，體積分?jǐn)?shù)方程中的C，湍流動能方程中的湍流動能K及湍流動能耗散方程中的耗散率ε;Γφ為廣義擴(kuò)散系數(shù);Sφ表示廣義源項(xiàng)．

劃分網(wǎng)格時(shí)在x、y、z方向上網(wǎng)格的長度≤0.11 m，并在人體、桌子、風(fēng)口附近進(jìn)行加密，共生成143 851網(wǎng)格．然后采用控制容積法將控制方程在網(wǎng)格上離散，差分格式使用混合格式，求解算法為SIMPLE算法．

1.3 邊界條件和初始條件

在進(jìn)行數(shù)值計(jì)算時(shí)，計(jì)算所需的邊界條件和初始條件如下:

1)精確的CFD模型可對室內(nèi)污染物體積分?jǐn)?shù)的分布進(jìn)行全面預(yù)測，風(fēng)口模型的設(shè)定是影響計(jì)算精度的主要因素之一．孔板送風(fēng)口的邊界條件簡化采用以下兩種處理方法:模型一采用“N點(diǎn)風(fēng)口模型”［3］，即利用外形面積與原孔板風(fēng)口相等的簡單開口替換復(fù)雜孔板風(fēng)口，以描述其入流邊界條件，并保證入流的質(zhì)量流量和動量流量與實(shí)際一致;模型二將孔板風(fēng)口等效為一個(gè)簡單開口，其面積與孔板風(fēng)口的有效通過面積相等，這樣可確保入流的動量流量和質(zhì)量流量與實(shí)際一致［8］．

2)在滿足工作區(qū)溫濕度和風(fēng)速要求的基礎(chǔ)上，房間送風(fēng)量可取 750 m3/h，送風(fēng)溫度為18℃．

3)因?yàn)閷?shí)驗(yàn)室建在室內(nèi)，可忽略實(shí)驗(yàn)室圍護(hù)結(jié)構(gòu)內(nèi)外的傳熱問題，模型中設(shè)定墻體為沒有厚度的絕熱墻體．

4)室內(nèi)空氣中TVOC的初始體積分?jǐn)?shù)為620×10－9．認(rèn)為室內(nèi)空氣為干空氣，室內(nèi)空氣初始溫度18℃．由于桌子散發(fā)污染物速率很小，短時(shí)間內(nèi)對室內(nèi)環(huán)境的影響較小，所以近似認(rèn)為桌子散發(fā)污染物速率為0．

2 實(shí)驗(yàn)方法與步驟

為得出實(shí)際情況下室內(nèi)污染物體積分?jǐn)?shù)分布規(guī)律，對實(shí)驗(yàn)房間內(nèi)的污染物分布進(jìn)行了測試．實(shí)驗(yàn)在圖1所示的房間中距地面1.0 m處布置測點(diǎn)5個(gè)，其位置如圖2所示．圖2中將人和桌子簡化成長方體，圓圈代表測試點(diǎn)．實(shí)驗(yàn)前在實(shí)驗(yàn)室內(nèi)人為布置污染源，污染源采用建筑裝修中常用的粘合劑涂在桌面上，則室內(nèi)的兩張桌子成為污染物TVOC的主要散發(fā)源．為了使室內(nèi)污染物的體積分?jǐn)?shù)符合實(shí)際新裝修建筑室內(nèi)的一般水平，將污染源在室外放置約1 h，再放置室內(nèi)約6 h，以使其污染物充分散發(fā)并均勻分布在室內(nèi)．污染物釋放過程中，實(shí)驗(yàn)室的門窗全部關(guān)閉，以防止人為的擾動．室內(nèi)總揮發(fā)性有機(jī)物(TVOC)的體積分?jǐn)?shù)用PpbRAE VOC檢測儀(PGM－7240X型)檢測，該設(shè)備的測試范圍是 0～9 999×10－9，精度為10%．

主要實(shí)驗(yàn)步驟:

1)按照儀器的使用要求連接儀器，并進(jìn)行使用前的校準(zhǔn)工作．

2)在各測點(diǎn)測試TVOC的初始體積分?jǐn)?shù)，測試時(shí)盡量減少室內(nèi)人員走動，以防止氣流擾動對實(shí)驗(yàn)造成的誤差．

3)打開空調(diào)系統(tǒng)，設(shè)置送風(fēng)溫度為18℃，送風(fēng)相對濕度為36%，風(fēng)量設(shè)置為750 m3/h．開始送風(fēng)時(shí)監(jiān)測測點(diǎn)處的污染物體積分?jǐn)?shù)，當(dāng)體積分?jǐn)?shù)發(fā)生變化時(shí)開始測量各測點(diǎn)的污染物體積分?jǐn)?shù)，本實(shí)驗(yàn)取各測點(diǎn)5 s內(nèi)的污染物體積分?jǐn)?shù)的平均值作為實(shí)驗(yàn)結(jié)果．每組測試時(shí)間間隔為60 s，共計(jì)測試10～15組．

4)當(dāng)室內(nèi)污染物體積分?jǐn)?shù)隨時(shí)間的變化不大時(shí)，停止測量，關(guān)閉實(shí)驗(yàn)儀器，關(guān)閉空調(diào)系統(tǒng)．

圖2 測試點(diǎn)布置(m)

3 模擬結(jié)果分析

3.1 兩種模型不同時(shí)刻模擬結(jié)果比較

圖3為不同時(shí)刻兩種模型模擬結(jié)果比較，可以看出:

1)采用兩種不同模型的模擬結(jié)果，室內(nèi)污染物體積分?jǐn)?shù)分布規(guī)律相近，其中測點(diǎn)1和測點(diǎn)4由于距離污染源較近體積分?jǐn)?shù)較高．采用模型一時(shí)同一時(shí)刻室內(nèi)測點(diǎn)2的體積分?jǐn)?shù)最低，而采用模型二時(shí)測點(diǎn)3的體積分?jǐn)?shù)最低，其主要原因是模型二相對于模型一送風(fēng)口尺寸減小了，使得模型2中的測點(diǎn)3剛好位于送風(fēng)口的正下方，而測點(diǎn)2位于送風(fēng)口的邊緣處．

2)同一時(shí)刻模型二所得室內(nèi)污染物體積分?jǐn)?shù)較低，主要是因?yàn)槟Ｐ投惋L(fēng)口較小，風(fēng)速相對較大．模型二中測點(diǎn)4所處的位置位于體積分?jǐn)?shù)“死角”內(nèi)，所以體積分?jǐn)?shù)偏高．

3)送風(fēng)1 000 s時(shí)，室內(nèi)污染物體積分?jǐn)?shù)變化趨于穩(wěn)定，兩種模型模擬得到的室內(nèi)污染物體積分?jǐn)?shù)基本相同，各測點(diǎn)TVOC體積分?jǐn)?shù)相差1 ～2 ×10－9．

圖3 兩種模型不同時(shí)刻模擬結(jié)果對比

3.2 各測點(diǎn)污染物體積分?jǐn)?shù)逐時(shí)變化情況

圖4是采用兩種模型時(shí)測點(diǎn)1和測點(diǎn)3污染物體積分?jǐn)?shù)逐時(shí)變化曲線，兩曲線基本重合，體積分?jǐn)?shù)下降趨勢完全相同，在大約第900 s時(shí)體積分?jǐn)?shù)變化曲線的斜率很小，室內(nèi)TVOC體積分?jǐn)?shù)已趨于穩(wěn)定．

3.3 Y方向污染物體積分?jǐn)?shù)分布

孔板送風(fēng)所形成的室內(nèi)氣流組織形式是典型的“活塞流”，送風(fēng)從上部空間逐漸向下推進(jìn)，所以同一時(shí)刻距地面不同高度的測點(diǎn)上污染物體積分?jǐn)?shù)有較大差別，本文取房間內(nèi)測點(diǎn)3位置上距地面不同高度上的5個(gè)測點(diǎn)，分析同一時(shí)刻各測點(diǎn)的體積分?jǐn)?shù)值．

從圖5可以看出:

1)兩種模型的模擬結(jié)果均表明，在送風(fēng)初期，離地面距離越大，測點(diǎn)污染物體積分?jǐn)?shù)越低．采用模型二得出的各測點(diǎn)體積分?jǐn)?shù)隨高度變化較大，60 s時(shí)刻模型一最大與最小體積分?jǐn)?shù)值相差20×10－9，模型二最大與最小體積分?jǐn)?shù)值相差45×10－9．

2)隨著送風(fēng)時(shí)間的延長各測點(diǎn)的體積分?jǐn)?shù)差值逐漸減小，說明室內(nèi)污染物體積分?jǐn)?shù)逐漸趨于均勻狀態(tài)．

3)距地面較近的測點(diǎn)污染物體積分?jǐn)?shù)隨高度的變化并不大，而在距地面1.4 m以上的區(qū)域內(nèi)高度會對污染物體積分?jǐn)?shù)產(chǎn)生很大影響，模型二表現(xiàn)得更加明顯．這說明射流區(qū)污染物體積分?jǐn)?shù)變化較明顯，在工作區(qū)污染物體積分?jǐn)?shù)較均勻．

圖4 測點(diǎn)1和測點(diǎn)3體積分?jǐn)?shù)變化曲線

圖5 y方向體積分?jǐn)?shù)分布對比

3.4 速度場比較

圖6是采用兩種不同模型模擬時(shí)，室內(nèi)(z=0.5 m截面與x=0.25 m截面)的流線分布情況．

圖6 兩種模型模擬的流線分布

采用不同簡化模型進(jìn)行數(shù)值模擬時(shí)，室內(nèi)氣流分布有很大的不同．采用模型一時(shí)室內(nèi)氣流分布均勻，送風(fēng)口下的區(qū)域內(nèi)流線近似平行，較為符合實(shí)際情況下的氣流分布;采用模型二時(shí)，在出風(fēng)口下的區(qū)域內(nèi)流線近似平行，而在出風(fēng)口的兩側(cè)氣流分布出現(xiàn)漩渦，其中z=0.5 m截面左側(cè)氣流因?yàn)槭艿饺说淖钃醵霈F(xiàn)一個(gè)很大的漩渦．正是由于兩種模型模擬的室內(nèi)氣流分布的差異導(dǎo)致室內(nèi)污染物的體積分?jǐn)?shù)分布有所不同．

4 數(shù)值模擬的驗(yàn)證

數(shù)值模擬結(jié)果表明，采用兩種不同的簡化模型進(jìn)行模擬的結(jié)果相差不大，但是哪種模型的模擬結(jié)果更接近實(shí)際情況，還需通過實(shí)驗(yàn)檢驗(yàn)．

圖7是將不同時(shí)刻的模擬值和實(shí)驗(yàn)值繪制在同一坐標(biāo)系下的體積分?jǐn)?shù)分布曲線．

圖7 模擬值和實(shí)驗(yàn)值對比

通過圖7中模擬值與實(shí)驗(yàn)值的對比，可得出下述結(jié)論:

1)模擬結(jié)果與實(shí)驗(yàn)結(jié)果顯示的室內(nèi)空氣中污染物體積分?jǐn)?shù)分布規(guī)律符合較好，對于不同測點(diǎn)，兩者在數(shù)值上存在0.1% ～15%的誤差，測點(diǎn)1與測點(diǎn)4由于靠近污染源，TVOC體積分?jǐn)?shù)比較高，測點(diǎn)2與測點(diǎn)5體積分?jǐn)?shù)較低，而測點(diǎn)3處于平均水平．

2)第1 000 s時(shí)，室內(nèi)TVOC體積分?jǐn)?shù)基本達(dá)到穩(wěn)定，實(shí)驗(yàn)所測得的平均值(142×10－9)與模擬結(jié)果(137×10－9)有一定差別，但誤差較小，進(jìn)一步證實(shí)了用數(shù)值模擬的方法研究室內(nèi)污染物分布特性的可行性．

3)實(shí)驗(yàn)測得數(shù)據(jù)與采用模型一所得的計(jì)算結(jié)果更為接近，因?yàn)樵趯?shí)際情況中送風(fēng)小孔均勻布滿整個(gè)天花板，“N點(diǎn)風(fēng)口模型”利用外形面積與原孔板風(fēng)口相等的簡單開口替換復(fù)雜孔板風(fēng)口，更符合實(shí)際情況，更適合進(jìn)行孔板送風(fēng)條件下污染物分布的數(shù)值研究．

5 結(jié)論

1)孔板送風(fēng)對室內(nèi)污染物有明顯稀釋作用．對于本文研究的空調(diào)房間，送風(fēng)量為750 m3/h時(shí)，房間開始送風(fēng)120 s后人坐姿呼吸層上的污染物體積分?jǐn)?shù)已經(jīng)降到了國家規(guī)定值500×10－9之下，送風(fēng)900 s后，污染物體積分?jǐn)?shù)趨于穩(wěn)定．

2)模擬和實(shí)驗(yàn)結(jié)果表明，室內(nèi)污染物體積分?jǐn)?shù)分布和送回風(fēng)口、污染源的位置有很大關(guān)系，距離送風(fēng)口較近的區(qū)域內(nèi)污染物體積分?jǐn)?shù)較低，污染源及其周邊的污染物體積分?jǐn)?shù)較高．

3)孔板送風(fēng)采用兩種入流邊界條件處理方法的模擬結(jié)果均和實(shí)驗(yàn)結(jié)果吻合較好，模型一的模擬結(jié)果更接近實(shí)驗(yàn)情況，因此建議采用模型一的“N點(diǎn)風(fēng)口模型”進(jìn)行孔板送風(fēng)條件下污染物分布的數(shù)值研究．

［1］李麗，徐文華．新風(fēng)量對室內(nèi)VOCs濃度分布影響的模擬分析［J］．潔凈與空調(diào)技術(shù)，2003(3):16－20．

［2］張明遠(yuǎn)．室內(nèi)高大空間污染物分布數(shù)值模擬［D］．北京:華北電力大學(xué)，2003．

［3］趙彬，曹莉，李先庭．潔凈室孔板型風(fēng)口入流邊界條件的處理方法［J］．清華大學(xué)學(xué)報(bào):自然科學(xué)版， 2003，43(5):690 －692．

［4］張文勝，安大偉，董瑞，等．密閉恒溫小室流場均勻性的研究［J］．山東建筑大學(xué)學(xué)報(bào)， 2006，21(5):430－433．

［5］李勇．ICU病房空調(diào)通風(fēng)系統(tǒng)數(shù)值模擬［D］．武漢:華中科技大學(xué)，2006．

［6］鹿世化，黃虎，李奇賀．?dāng)?shù)值模擬用于室內(nèi)空氣質(zhì)量控制的研究進(jìn)展［J］．暖通空調(diào)， 2007，37(5):40 －46．

［7］陶文銓．?dāng)?shù)值傳熱學(xué)［M］．2版．西安:西安交通大學(xué)出版社，2001．

［8］NIELSEN P V．Description of supply openings in numerical models for room air distribution ［J］．ASHRAE Transactions， 1992，98(1):963 －971．

Numerical simulation and experimental investigation of contamination distribution in a chamber with perforated ceiling air supply

SUN Li-ying1，WEI Hui-jiao2

(1．College of Aerospace and Civil Engineering，Harbin Engineering University，150001 Harbin，China，sunliying@hrbeu．edu．cn;2．China Academy of Building Research(Tianjin Branch)，300834 Tianjin，China)

To study the contamination distribution in a chamber with perforated ceiling air supply，simulative analysis and experimental tests are carried out in this paper．In modeling，the Basic model and N-spot air model are both used for defining the orifice entrance boundary condition．The simulation results show that the distribution of contamination concentration agrees well with the test result．Especially，the simulation result of the N-spot air model is closer to the experimental one，which proves the feasibility of this model in analyzing indoor contamination distribution．

perforated ceiling air supply;contamination;numerical simulation

TU831.3

A

0367－6234(2011)08－0135－05

2010－03－10．

中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(HEUCF100206)．

孫麗穎(1973—)，女，博士，副教授．

(編輯 趙麗瑩)