摘要：為獲得最優(yōu)的冬季自然通風(fēng)時(shí)段，以蘭州地區(qū)某地板供暖民用住宅為研究對象，考慮圍護(hù)結(jié)構(gòu)對室外溫度波的衰減和延遲效應(yīng)，得到建筑圍護(hù)結(jié)構(gòu)的等熱流邊界條件。在此基礎(chǔ)上，采用FLUENT軟件對不同通風(fēng)時(shí)段的室內(nèi)污染物濃度場進(jìn)行了數(shù)值分析，結(jié)果表明：同一開窗方式時(shí)，不同時(shí)段的通風(fēng)時(shí)長不同，不同進(jìn)風(fēng)溫度對室內(nèi)污染物的時(shí)空分布及遷移特性影響顯著，綜合考慮污染物排除效果及舒適性，發(fā)現(xiàn) 17：00時(shí)開始通風(fēng)為最優(yōu)通風(fēng)時(shí)段。

關(guān)鍵詞：自然通風(fēng)；非穩(wěn)定傳熱；通風(fēng)效率；優(yōu)化

中圖分類號：TU832 文獻(xiàn)標(biāo)志碼：A 文章編號：1674-4764（2015）05-0101-08

Abstract：In order to find out the optimal period of time of natural ventilation in winter， the realizable k-ε turbulent model of CFD is adopted to numerically analyze the pollutant concentration field in a residence in Lanzhou when the unsteady heat transfer through outerwalls is considered. The results indicate that the inlet temperature has significant effect on the indoor pollutants profile and migration， and different ventilation period corresponds to different duration under the same opening patterns. Considering the pollutants elimination efficiency and thermal comfort simultaneously， the ventilation beginning at 17：00 is the optimal ventilation mode.

Key words：natural ventilation； unsteady heat transfer； ventilation efficiency； optimization

隨著經(jīng)濟(jì)的發(fā)展，人們對室內(nèi)空氣品質(zhì)的要求不斷提高。自然通風(fēng)由于具有控制室內(nèi)污染物濃度水平、改善室內(nèi)空氣品質(zhì)、調(diào)節(jié)室內(nèi)熱環(huán)境等功能，并且具有“節(jié)能、環(huán)境效益好”等優(yōu)點(diǎn)，越來越受到人們的重視。文獻(xiàn)[1]對某教學(xué)樓辦公室在秋季利用熱壓和風(fēng)壓共同作用下的自然通風(fēng)與室內(nèi)污染物排除效率關(guān)聯(lián)性進(jìn)行了數(shù)值和實(shí)驗(yàn)研究。文獻(xiàn)[2]對不同污染源散發(fā)條件下的自然置換通風(fēng)進(jìn)行了實(shí)驗(yàn)研究，結(jié)果表明：室內(nèi)外溫差較大時(shí)有利于污染物的排出。文獻(xiàn)[3]實(shí)驗(yàn)研究了重慶某建筑夏季自然通風(fēng)對室內(nèi)污染物排除效果、室內(nèi)溫度、濕度的影響關(guān)系。文獻(xiàn)[4]對德國法蘭克福辦公樓的通風(fēng)監(jiān)測數(shù)據(jù)表明：冬季通風(fēng)時(shí)間超過全天時(shí)間的10%～25%，會使冬季的能量消耗增加。文獻(xiàn)[5]對印度復(fù)雜氣候條件下自然通風(fēng)效果進(jìn)行了仿真模擬，發(fā)現(xiàn)在夏季改變開窗策略可以提高13%的熱舒適度，夜間通風(fēng)顯著改善了室內(nèi)熱舒適性，但冬天只能利用白天通風(fēng)，而且難以保證室內(nèi)的熱舒適性。文獻(xiàn)[6]對意大利4 000名學(xué)生在冬季和夏季待在教室時(shí)的主觀熱感覺進(jìn)行了調(diào)查研究，結(jié)合客觀的自然通風(fēng)條件給出了地中海氣候PMV的預(yù)期因子的建議值。以上研究工作主要是針對夏季或者過渡季展開的。而寒冷地區(qū)冬季氣溫低、晝夜溫度波動大，如何合理利用自然通風(fēng)的優(yōu)越性來達(dá)到既改善室內(nèi)空氣品質(zhì)又不影響居住者的舒適性這一目標(biāo)，是在寒冷地區(qū)實(shí)施自然通風(fēng)的前提。而生活在同一熱工分區(qū)的居住者，由于起居習(xí)慣、生活規(guī)律等存在差異，他們選擇開窗通風(fēng)的時(shí)間段也必然有所不同。目前，關(guān)于寒冷地區(qū)冬季自然通風(fēng)時(shí)段的選擇研究，還未見報(bào)道。本文以蘭州某民用住宅為研究對象，試圖通過自然通風(fēng)數(shù)值試驗(yàn)，探尋適合寒冷地區(qū)氣候特征的冬季自然通風(fēng)最優(yōu)時(shí)段。

1 物理模型和數(shù)學(xué)模型

1.1 物理模型

所研究住宅物理模型如圖1所示。位于8層建筑的中間樓層，結(jié)構(gòu)尺寸為：X×Y×Z=10.5 m×13.2 m×2.9 m。其中，1#為客廳與餐廳，2#、3#、4#為臥室。C1、C2為北外窗，C3、C4為南外窗。外墻傳熱系數(shù)K=0.46 W/（m2·K）。推拉式窗戶傳熱系數(shù)Kc=2.2 W/（m2·K）。室內(nèi)采用地板供暖。

1.3 邊界條件

1） 入口邊界：綜合IWEC（美國國家氣象數(shù)據(jù)中心）和CSWD（清華大學(xué)和中國氣象局）所提供蘭州冬季氣象參數(shù)，蘭州在冬季以北風(fēng)為主導(dǎo)風(fēng)向，故本文以北外窗C1、C2為計(jì)算區(qū)域的流動入口邊界；11月的月平均風(fēng)速最小，兩種統(tǒng)計(jì)結(jié)果均為0.2 m/s，3月份風(fēng)速最大，兩種統(tǒng)計(jì)結(jié)果均為1.2 m/s。前期研究表明：取1.2 m/s進(jìn)風(fēng)速度不能同時(shí)滿足室內(nèi)采暖溫度和污染物濃度要求，故本文進(jìn)風(fēng)速度取0.4 m/s；進(jìn)風(fēng)溫度取不同通風(fēng)時(shí)段內(nèi)的室外平均溫度。

2）出口邊界：排風(fēng)口C3、C4設(shè)為自由出流邊界條件[10]。

3）壁面邊界：氣固交界面的空氣流動取速度無滑移條件。

4）熱邊界設(shè)置： 采用非穩(wěn)定傳熱方法得到不同通風(fēng)時(shí)段通過圍護(hù)結(jié)構(gòu)的熱流密度平均值，以此為圍護(hù)結(jié)構(gòu)的熱邊界條件。具體數(shù)值見表2。室內(nèi)采暖熱負(fù)荷均由地板提供，設(shè)為定熱流邊界條件，根據(jù)文獻(xiàn)[11]設(shè)定地板的熱流密度為45 W/m2。內(nèi)墻、天花板均設(shè)為絕熱邊界。

5） 污染源條件設(shè)置：來自地板的污染物C7H8假定其沿地板外法線方向向上均勻散發(fā)，根據(jù)文獻(xiàn)[12]對室內(nèi)C7H8允許濃度規(guī)定，結(jié)合散發(fā)時(shí)長以及房間體積，可求得散發(fā)速率為5×10-9 kg/s。選用層流模型模擬了污染物的擴(kuò)散過程（限于篇幅，擴(kuò)散過程另文討論）。結(jié)合人們冬季開關(guān)窗的習(xí)慣，對計(jì)算房間密閉情況下污染物散發(fā)15、17、19、21 h后的室內(nèi)污染物濃度場分別作為4個通風(fēng)時(shí)段自然通風(fēng)數(shù)值模擬的初始條件。

2 數(shù)值求解方法

2.1 網(wǎng)格劃分及獨(dú)立性驗(yàn)證

采用六面體結(jié)構(gòu)網(wǎng)格對計(jì)算區(qū)域進(jìn)行離散?？紤]進(jìn)、排風(fēng)口處流動參數(shù)的高梯度變化，對這些區(qū)域的網(wǎng)格進(jìn)行了局部加密。分別采用3套網(wǎng)格（網(wǎng)格數(shù)分別為520 590、672 138和868 434）進(jìn)行了數(shù)值試驗(yàn)，所得（X=8.98 m，Y=3 m）處的速度和溫度計(jì)算結(jié)果均吻合得很好，說明本文所得解是網(wǎng)格獨(dú)立性的?？紤]計(jì)算的經(jīng)濟(jì)性，決定選用672 138作為后續(xù)計(jì)算的網(wǎng)格數(shù)。

2.2 數(shù)值方法

采用有限體積法對控制方程（1）進(jìn)行離散；應(yīng)用SIMPLE算法求解速度/壓力耦合問題[12]；梯度項(xiàng)方程采用Green-Gauss Cell Based格式離散；動量方程、能量方程、湍流動能方程、湍流動能耗散率方程、組分輸運(yùn)方程均采用二階迎風(fēng)格式離散。

2.3 模擬工況及收斂準(zhǔn)則

2.3.1 工況確定

為了確定合理的通風(fēng)工況，規(guī)定同時(shí)滿足：通風(fēng)1 h后距地面1.1 m高度水平面（規(guī)定為工作區(qū)高度）上溫度的平均值不低于16.0 ℃和該平面上污染物濃度不高于文獻(xiàn)[13]的要求值2.17×10-9 kmol/m3，并且，地面上方0.1 m到1.1 m間的豎向溫差不大于3.0 ℃，則認(rèn)為該工況是可行的自然通風(fēng)工況。據(jù)此，對進(jìn)風(fēng)溫度為0 ℃，進(jìn)風(fēng)速度為0.4 m/s，與不同進(jìn)、排風(fēng)口開啟度組合，形成的多工況進(jìn)行了數(shù)值試驗(yàn)，結(jié)果表明：窗戶C1、C2、C3的開啟度均為0.1 m寬（關(guān)閉窗戶C4）時(shí)為最佳通風(fēng)工況。為對比分析不同通風(fēng)時(shí)段的室內(nèi)相關(guān)參數(shù)分布，采用這一窗戶開關(guān)方式，增加另外3個通風(fēng)時(shí)段的起始時(shí)刻分別為13：00、15：00、17：00，從而組成4個對比通風(fēng)工況，如表3所示。各工況通風(fēng)時(shí)長的確定過程如圖2所示。

圖6為不同工況時(shí)Z=1.1 m水平面上的C7H8濃度分布?？梢钥闯觯鞴r下客廳C7H8濃度分布表現(xiàn)出了較一致的特征：靠近客廳西墻附近區(qū)域C7H8濃度較高，沿著墻體呈帶狀分布。這是因?yàn)槭彝獾牡蜏貧怏w進(jìn)入室內(nèi)后，絕大部分很快下沉并與地面附近的污染物進(jìn)行熱質(zhì)交換，只有少部分氣體在客廳西北角落處形成了回流區(qū)，與主流區(qū)新鮮空氣的熱質(zhì)交換遠(yuǎn)不及其他區(qū)域充分。所以，隨著低溫氣體自北朝南方向運(yùn)動，有一部分污染物便被“壓制”在了回流區(qū)及客廳西墻附近區(qū)域[14]?？蛷d絕大部分區(qū)域空氣質(zhì)量符合文獻(xiàn)[13]要求，工況4時(shí)客廳的空氣質(zhì)量明顯要優(yōu)于其他3個工況。各工況下，2#房間內(nèi)污染物的分布差異顯著：工況1時(shí)2#房間內(nèi)污染物分布均勻，均未超標(biāo)，只是在該房間的東南角處有小范圍超標(biāo)；工況2和3時(shí)，2#房間有近一半面積的污染物集聚并超標(biāo)；工況4時(shí)，在2#房間靠近門有小范圍污染物集聚并超標(biāo)。這是因?yàn)楣r1和工況4的進(jìn)風(fēng)溫度要低于工況2和工況3的進(jìn)風(fēng)溫度，2#房間的門使過流斷面突然縮小引起流場結(jié)構(gòu)變化，從而對污染物遷移產(chǎn)生影響，但這一影響遠(yuǎn)不及溫度場對污染物濃度場影響。3#房間作為各工況通風(fēng)的唯一出流通道，對于污染物的時(shí)空分配起著重要作用，其污染物濃度梯度也最為復(fù)雜。在北墻內(nèi)側(cè)不同范圍的污染物集聚可能與溫度場以及該區(qū)域的回流對污染物擴(kuò)散過程影響有關(guān)，但各工況下4#房間內(nèi)污染物濃度均未超標(biāo)。

可以看出，4種工況的舒適性表現(xiàn)出了相似的變化規(guī)律：PMV值沿房間高度方向逐漸增大，這主要是室內(nèi)溫度分層所致。根據(jù)PMV-PPD的熱舒適評價(jià)指標(biāo)及ISO7730標(biāo)準(zhǔn)規(guī)定：以PPD≤10%為設(shè)計(jì)依據(jù)[15]，即90%以上的人感到滿意的熱環(huán)境為熱舒適環(huán)境，此時(shí)，對應(yīng)的PMV= -0.5～+0.5。但考慮中國的實(shí)際經(jīng)濟(jì)狀況，一般計(jì)算時(shí)都將該值放寬至PPD≤20%[16]，此時(shí)，對應(yīng)的PMV= -0.75～+0.75。所以，根據(jù)中國對舒適性要求的PMV值以及文獻(xiàn)[15]中關(guān)于熱感覺的7級標(biāo)度（-3冷，-2涼，-1稍涼，0不冷不熱，+l稍暖，+2暖，+3熱）之規(guī)定，工況1通風(fēng)過程的PMV值自0.1 m至1.7 m高度處均小于0，室內(nèi)從頭到腳整體偏涼；工況2，在0.1～0.92 m高度范圍內(nèi)是舒適的；工況3，在0.1～0.58 m高度范圍內(nèi)是舒適的；工況4，在0.1～1.7 m高度范圍內(nèi)是舒適的，就是說，按工況4通風(fēng)時(shí)，居住者從頭到腳都是舒適的。工況1～3對應(yīng)的舒適區(qū)高度均低于一般建筑室內(nèi)確定PMV值的水平面高度1.1 m。在1.1 m水平面上，工況1的PMV值為-0.3，介于“熱中性”與“稍涼”之間；工況2的PMV值為0.83，“稍暖”；工況3的PMV值為1.1，“偏暖”；工況4的PMV值為0.62，介于“熱中性”與“稍暖”之間。所以，考察Z=1.1 m水平面上的舒適性時(shí)只有工況1能滿足國際上對舒適性要求的PMV值，但房間整體偏涼。對于寒冷地區(qū)的冬季，居住者卻更喜歡待在工況2、3、4這樣“偏暖”的室內(nèi)。實(shí)際上，冬季室內(nèi)溫度偏低，反而有利于抑制污染物的釋放[17]。

綜上所述，工況1、2、4均能滿足符合我國實(shí)際的舒適性要求，工況4的舒適度最好。

3.4 最佳通風(fēng)時(shí)段確定

自然通風(fēng)是將當(dāng)?shù)貧夂驐l件和居住者的生活習(xí)慣緊密結(jié)合的一種改善室內(nèi)空氣品質(zhì)的措施。上文關(guān)于排污效率和舒適性評價(jià)的討論得到了這樣一個結(jié)論：感覺舒適的環(huán)境未必對健康有利，反之亦然。從圖3（a）可以看出，工況1整個通風(fēng)過程的排污效率均是最高的，工況4次之，工況3最低。從圖3（b）可以看出，通風(fēng)結(jié)束時(shí)，工況1、4的排污效率相當(dāng)，均高于其它工況。綜合通風(fēng)過程的排污效率變化以及室內(nèi)舒適性評價(jià)結(jié)果，工況4為最佳通風(fēng)時(shí)段。居住者如果能通過改變衣著來適應(yīng)室內(nèi)偏涼的環(huán)境，工況1也是比較理想的通風(fēng)時(shí)段。

4 結(jié) 論

以蘭州某民用住宅為研究對象，采用Realizable k-ε 模型對冬季自然通風(fēng)時(shí)室內(nèi)污染物濃度場進(jìn)行了數(shù)值分析，得到了如下主要結(jié)論：

1）寒冷地區(qū)冬季自然通風(fēng)可以在不同時(shí)段進(jìn)行，但每一時(shí)段的通風(fēng)時(shí)長存在差異。

2） 不同通風(fēng)時(shí)段的排污效率隨時(shí)間呈現(xiàn)出了相似的變化趨勢，但對污染物的排除效果存在差異。相同的開窗方式，進(jìn)風(fēng)溫度對室內(nèi)污染物的時(shí)空分布有顯著影響。

3）所研究的4個通風(fēng)時(shí)段，工況1、2、4均能滿足符合我國實(shí)際的舒適性要求，工況3稍偏暖，工況4的舒適度最好。結(jié)合排污效率指標(biāo)，17：00時(shí)開窗通風(fēng)，即工況4為最優(yōu)通風(fēng)工況。

4）本文是針對特定結(jié)構(gòu)建筑在特定氣候條件下尋求不同通風(fēng)時(shí)段的通風(fēng)模式、實(shí)現(xiàn)室內(nèi)空氣質(zhì)量動態(tài)化改善的初步探索，后續(xù)研究將會考慮更多因素，使得所得通風(fēng)模式對建筑類型、氣候條件、供暖方式等因素具有更好的適應(yīng)性。

參考文獻(xiàn)：

[1] 周懷改. 通風(fēng)對建筑物室內(nèi)污染物濃度分布的影響[D]. 重慶： 重慶大學(xué)， 2007.

Zhou H G. Effect of ventilation on indoor pollutant concentration distribution [D].Chongqing： Chongqing University， 2007.（in Chinese）

[2] 楊秀峰. 自然置換通風(fēng)條件下室內(nèi)空氣污染的演化規(guī)律研究[D]. 上海： 東華大學(xué)， 2013.

Yang X F. Investigations on the evolution of indoor air pollution in natural displacement ventilation room [D].Shanghai：Donghua University， 2013.（in Chinese）

[3] 劉慶. 自然通風(fēng)下門窗開啟對室內(nèi)環(huán)境的影響研究[D]. 重慶： 重慶大學(xué)， 2014.

Liu Q. Research of indoor environment affected by opening door and window in the condition of natural ventilation [D].Chongqing： Chongqing University， 2014.（in Chinese）

[4] Schakib-Ekbatan K， akc F C， Schweiker M， et al. Does the occupant behavior match the energy concept of the building analysis of a German naturally ventilated office building [J]. Building and Environment， 2015， 84： 142-150.

[5] Rajasekar E， Anupama U，Venkateswaran R. Thermal comfort beyond building design. An investigation in naturally ventilated residential apartments in a hot-dry climate [J]. Advances in Building Energy Research， 2014， 8（2）： 196-215.

[6] d’Ambrosio A F R， Ianniello E， Palella B I. PMV-PPD and acceptability in naturally ventilated schools [J]. Building and Environment， 2013， 67： 129-137.

[7] Blay D， Mergui S， Niculae C. Confined turbulent mixed convection in the presence of a horizontal buoyant wall jet[J]. Fundamentals of Mixed Convection， 1992， 213：65-72.

[8] 李先庭， 江億. 用計(jì)算流體動力學(xué)方法求解通風(fēng)房間的空氣年齡[J]. 清華大學(xué)學(xué)報(bào)， 1998， 38（5）： 28-31.

Li X T，Jiang Y. Calculating air age in a ventilated room with CFD method [J]. Journal of Tsinghua University， 1998， 38（5）： 28-31.（in Chinese）

[9] 陶文銓. 數(shù)值傳熱學(xué)[M].2版. 西安： 西安交通大學(xué)出版社， 2001.

Tao W Q. Numerical heat transfer [M].Second Edition. Xi'an： Xi'an Jiaotong University Press， 2001.（in Chinese）

[10] 張卓鵬. 廣州地區(qū)圍合式住區(qū)室內(nèi)自然通風(fēng)研究[D]. 廣州： 華南理工大學(xué)， 2013.

Zhang Z P. Research on indoor naturalventilation of enclosed residential districts inGuangzhou[D]. Guangzhou： South China University of Technology， 2013.（in Chinese）

[11] 陸耀慶. 實(shí)用供熱空調(diào)設(shè)計(jì)手冊[M].2版. 北京： 中國建筑工業(yè)出版社， 2007.

Lu Y Q.Practical heating air conditioning design manual [M]. Second Edition. Beijing： China Building Industry Press， 2007.（in Chinese）

[12] Patankar S V. Numerical heat transfer and fluid flow [M]. New York： Hemisphere， 1980： 330-351.

[13] 中國國家質(zhì)量檢驗(yàn)檢疫總局.GB/T 18883—2002 室內(nèi)空氣質(zhì)量標(biāo)準(zhǔn)[S].北京：中國標(biāo)準(zhǔn)出版社，2003.

General Administration of Quality Supervision， Inspection and Quarantine ofChina. GB/T 18883—2002 Standards of indoor air quality [S]. Beijing： Standards Press of China， 2003.（in Chinese）

[14] 王燁， 張文霞， 胡文婷. 室內(nèi)環(huán)境參數(shù)對室外氣象參數(shù)瞬時(shí)變化的動態(tài)響應(yīng)研究[J].重慶大學(xué)學(xué)報(bào)， 2015， 38（3） ： 8-14.

Wang Y， Zhang W X， Hu W T. Dynamic response study of the indoor environmental parameters for variable outdoor meteorological conditions [J]. Journal of Chongqing University， 2015， 38（3）： 8-14.（in Chinese）

[15] ISO International Standard 7730. Moderate thermal environment determination of the PMV and PPD indices and specification of the conditions for thermal comfort[S]. Geneva： International Standard Organization， 1984.

[16] 王昭俊.嚴(yán)寒地區(qū)居室熱環(huán)境與熱舒適性研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2002：27-57.

Wang Z J. Study on indoor thermal environment and thermal comfort of residence in cold area [D]. Harbin：Harbin University of Technology， 2002：27-57.（in Chinese）

[17] 徐東群. 居住環(huán)境空氣污染與健康[M]. 北京：化學(xué)工業(yè)出版社， 2005.

Xu D Q. Air contamination and health of residential environment [M]. Beijing： Chemical Industry Press， 2005.（in Chinese）

（編輯 王秀玲）