寒冷地區(qū)冬季自然通風(fēng)時(shí)段選擇優(yōu)化分析

王 燁,王靖文

(蘭州交通大學(xué) a.環(huán)境與市政工程學(xué)院；b.鐵道車輛熱工教育部重點(diǎn)實(shí)驗(yàn)室，蘭州 730070)

寒冷地區(qū)冬季自然通風(fēng)時(shí)段選擇優(yōu)化分析

王 燁a,b,王靖文a

(蘭州交通大學(xué) a.環(huán)境與市政工程學(xué)院；b.鐵道車輛熱工教育部重點(diǎn)實(shí)驗(yàn)室，蘭州 730070)

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

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

隨著經(jīng)濟(jì)的發(fā)展，人們對(duì)室內(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]對(duì)某教學(xué)樓辦公室在秋季利用熱壓和風(fēng)壓共同作用下的自然通風(fēng)與室內(nèi)污染物排除效率關(guān)聯(lián)性進(jìn)行了數(shù)值和實(shí)驗(yàn)研究。文獻(xiàn)[2]對(duì)不同污染源散發(fā)條件下的自然置換通風(fēng)進(jìn)行了實(shí)驗(yàn)研究，結(jié)果表明：室內(nèi)外溫差較大時(shí)有利于污染物的排出。文獻(xiàn)[3]實(shí)驗(yàn)研究了重慶某建筑夏季自然通風(fēng)對(duì)室內(nèi)污染物排除效果、室內(nèi)溫度、濕度的影響關(guān)系。文獻(xiàn)[4]對(duì)德國法蘭克福辦公樓的通風(fēng)監(jiān)測數(shù)據(jù)表明：冬季通風(fēng)時(shí)間超過全天時(shí)間的10%～25%，會(huì)使冬季的能量消耗增加。文獻(xiàn)[5]對(duì)印度復(fù)雜氣候條件下自然通風(fēng)效果進(jìn)行了仿真模擬，發(fā)現(xiàn)在夏季改變開窗策略可以提高13%的熱舒適度，夜間通風(fēng)顯著改善了室內(nèi)熱舒適性，但冬天只能利用白天通風(fēng)，而且難以保證室內(nèi)的熱舒適性。文獻(xiàn)[6]對(duì)意大利4 000名學(xué)生在冬季和夏季待在教室時(shí)的主觀熱感覺進(jìn)行了調(diào)查研究，結(jié)合客觀的自然通風(fēng)條件給出了地中海氣候PMV的預(yù)期因子的建議值。以上研究工作主要是針對(duì)夏季或者過渡季展開的。而寒冷地區(qū)冬季氣溫低、晝夜溫度波動(dòng)大，如何合理利用自然通風(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)道。本文以蘭州某民用住宅為研究對(duì)象，試圖通過自然通風(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 計(jì)算房間平面示意圖Fig. 1 Physical model of the calculating room

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

采用標(biāo)準(zhǔn)k-ε模型、RNGk-ε模型以及Realizablek-ε模型分別對(duì)文獻(xiàn)[7]的通風(fēng)過程進(jìn)行了數(shù)值計(jì)算，并將數(shù)值結(jié)果與文獻(xiàn)[7]的實(shí)驗(yàn)結(jié)果進(jìn)行了比對(duì)。結(jié)果表明：Realizablek-ε模型計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果吻合度最好。所以，選用Realizablek-ε模型進(jìn)行模擬。自然通風(fēng)下室內(nèi)氣流運(yùn)動(dòng)屬于非穩(wěn)態(tài)湍流流動(dòng)，與組分輸運(yùn)方程具有相同的形式[8]，其通用控制方程形式均可表示為

(1)

式中：φ為通用變量，分別表示u、v、w、T、k、ε、τp。方程中各系數(shù)列于表1。

注：μt=cμρk2/ε；cμ為經(jīng)驗(yàn)常數(shù)；Prt為能量方程中的湍流普朗特?cái)?shù)；σk為湍流動(dòng)能方程中的普朗特?cái)?shù)；σε為湍流動(dòng)能耗散率方程中的普朗特?cái)?shù)；στ為組分輸運(yùn)方程中的普朗特?cái)?shù)。

湍流動(dòng)能剪切項(xiàng)：

方程中各系數(shù)取值[9]：cμ=0.09,σk=1.0,σε=1.2,Prt=1.0,στ=1.0,c2=1.9。

S=(2Si,jSi,j)1/2，Si,j=(?ui/?xj+?uj/?xi)/2。

為簡化計(jì)算，假定：室內(nèi)空氣不可壓縮，密度的變化采用Boussinesq假設(shè)，忽略流體粘性力作功而引起的耗散熱，室內(nèi)空氣為輻射透明介質(zhì)，不參與輻射換熱。

1.3 邊界條件

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

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

3)壁面邊界：氣固交界面的空氣流動(dòng)取速度無滑移條件。

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è)定地板的熱流密度為45W/m2。內(nèi)墻、天花板均設(shè)為絕熱邊界。

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

表2 圍護(hù)結(jié)構(gòu)熱流密度Table 2 Heat flux of the building envelopes (W·m-2)

2 數(shù)值求解方法

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

采用六面體結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行離散?？紤]進(jìn)、排風(fēng)口處流動(dòng)參數(shù)的高梯度變化，對(duì)這些區(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ù)值方法

采用有限體積法對(duì)控制方程(1)進(jìn)行離散；應(yīng)用SIMPLE算法求解速度/壓力耦合問題[12]；梯度項(xiàng)方程采用Green-Gauss Cell Based格式離散；動(dòng)量方程、能量方程、湍流動(dòng)能方程、湍流動(dòng)能耗散率方程、組分輸運(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-9kmol/m3,并且，地面上方0.1 m到1.1 m間的豎向溫差不大于3.0 ℃，則認(rèn)為該工況是可行的自然通風(fēng)工況。據(jù)此，對(duì)進(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)工況。為對(duì)比分析不同通風(fēng)時(shí)段的室內(nèi)相關(guān)參數(shù)分布，采用這一窗戶開關(guān)方式，增加另外3個(gè)通風(fēng)時(shí)段的起始時(shí)刻分別為13:00、15:00、17:00，從而組成4個(gè)對(duì)比通風(fēng)工況，如表3所示。各工況通風(fēng)時(shí)長的確定過程如圖2所示。

表3 自然通風(fēng)模擬工況Table 3 Simulation cases of natural ventilation

圖2 各工況通風(fēng)時(shí)長確定Fig.2 Ventilation time decision of the 4 cases

2.3.2 收斂準(zhǔn)則 計(jì)算中，同時(shí)滿足以下條件，認(rèn)為計(jì)算已收斂：

1)連續(xù)性方程及動(dòng)量方程殘差設(shè)為10-3；

2)能量方程求解殘差設(shè)為10-6；

3)監(jiān)視位置氣流參數(shù)不再波動(dòng)；

4)進(jìn)出口流體質(zhì)量守恒。

3 計(jì)算結(jié)果及分析

3.1 通風(fēng)效率分析

一般利用工作區(qū)的通風(fēng)效率來評(píng)價(jià)室內(nèi)污染物的排除效果，其定義式如式(2)所示。

(2)

式中:η為通風(fēng)效率；Cg為工作區(qū)平均濃度，kmol/m3；Cp為排風(fēng)口濃度，kmol/m3；Cs為進(jìn)風(fēng)口濃度，kmol/m3。

對(duì)工況1～4，統(tǒng)計(jì)不同通風(fēng)時(shí)長Z=1.1 m平面上的C7H8平均濃度作為工作區(qū)污染物平均濃度，同樣的方法得到排風(fēng)口對(duì)應(yīng)時(shí)長的污染物平均濃度Cp，進(jìn)入室內(nèi)的新風(fēng)中認(rèn)為不含有C7H8，所以，取Cs=0。據(jù)此得到4種工況下的通風(fēng)效率，如圖3(a)所示?？梢钥闯觯?種工況的通風(fēng)效率隨時(shí)間呈現(xiàn)出類似的變化趨勢：通風(fēng)開始的10～15 min內(nèi)通風(fēng)效率迅速上升，接著緩慢下降，最后趨于微弱的水平波動(dòng)，數(shù)值均在0.95～1.0之間波動(dòng)。這是因?yàn)橥L(fēng)開始時(shí)，房間內(nèi)C7H8分布較均勻，隨著通風(fēng)的持續(xù)進(jìn)行，C7H8濃度出現(xiàn)分層現(xiàn)象，Z=1.1 m水平面上的C7H8面平均濃度下降趨勢因此變緩，從而導(dǎo)致排風(fēng)口的C7H8平均濃度下降速率比Z=1.1 m水平面上的C7H8面平均濃度下降速率大。但總體上，工況1的通風(fēng)效率要高于其他各工況的值，通風(fēng)40 min后，工況1與工況4的通風(fēng)效率趨于一致，而工況2和工況3的通風(fēng)效率趨于一致。圖3(b)為各工況通風(fēng)結(jié)束時(shí)的排污效率比較，也印證了這一點(diǎn)。這說明雖然同樣的開關(guān)窗情況對(duì)應(yīng)相同的氣流運(yùn)動(dòng)路徑，但由于各工況進(jìn)風(fēng)溫度不同、室內(nèi)原有的污染物濃度的空間分布不同，導(dǎo)致了冷熱氣體的熱質(zhì)交換過程也必然不同。溫度梯度與濃度梯度的雙重作用，使得自然通風(fēng)中流場、溫度場、污染物濃度場之間的耦合效應(yīng)與排污效率之間的關(guān)聯(lián)性得以體現(xiàn)。

3.2 污染物濃度場分析

下面給出2.3.1節(jié)所確定的工況1～4通風(fēng)結(jié)束時(shí)的模擬結(jié)果。

圖4為不同工況時(shí)客廳進(jìn)風(fēng)窗口中軸線(X=1.15 m)縱截面上C7H8濃度分布?？梢钥闯觯?個(gè)工況表現(xiàn)出了相似的污染物分布特征：自地面朝天花板方向，C7H8濃度由低到高呈層狀分布，地面附近區(qū)域C7H8濃度梯度較大。這是因?yàn)闇囟容^低的室外新風(fēng)進(jìn)入房間后，沿窗臺(tái)下內(nèi)墻壁形成下降氣流，與地板附近的熱空氣相遇后，便朝遠(yuǎn)離窗口的方向運(yùn)動(dòng)，沿地面形成“冷氣湖”，其厚度自北向南逐漸變薄，地面附近的污染物同時(shí)被稀釋。在人員活動(dòng)區(qū)高度范圍內(nèi)，工況1和工況4的污染物濃度均未超過文獻(xiàn)[13]規(guī)定的上限值，工況2和工況3有局部區(qū)域污染物超標(biāo)現(xiàn)象，而且，在天花板的北部區(qū)域有大面積污染物滯留現(xiàn)象。這是因?yàn)樵诿芏炔钚纬傻母∩ψ饔孟拢廴疚镫S著自然對(duì)流邊界層的運(yùn)動(dòng)，在北外墻上部和天花板所成的角落處形成了聚集，導(dǎo)致該區(qū)域污染物濃度偏高。

圖3 通風(fēng)效率變化曲線Fig.3 Ventilation efficiency profile

圖4 X=1.15 m截面C7H8濃度場Fig.4 C7H8 concentration profile in X=1.15 m section

圖5為不同工況時(shí)Y=3 m橫截面(包含了3#房間、4#房間以及客廳部分區(qū)域)上的C7H8濃度分布?？梢钥闯?，1#房間西墻和北墻形成的角落處有污染物滯留現(xiàn)象，這是因?yàn)檫M(jìn)入室內(nèi)的冷氣流自北朝南的運(yùn)動(dòng)導(dǎo)致該角落處產(chǎn)生了回流。但各工況下1#房間內(nèi)污染物濃度均未超標(biāo)。2#房間內(nèi)的污染物濃度最低，4#房間次之，但各房間內(nèi)污染物濃度均未超標(biāo)。3#房間內(nèi)的污染物濃度梯度較1#和2#房間復(fù)雜，工況2和工況3時(shí)，在3#房間上部區(qū)域有較大范圍的污染物集聚現(xiàn)象，這是由漩渦區(qū)以及熱滯留所致?？傮w上，該截面上人員活動(dòng)區(qū)的污染物濃度在各工況下均未超標(biāo)。

圖5 Y=3 m截面C7H8濃度場 Fig.5 C7H8 concentration profile in Y=3 m section

圖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)動(dòng)，有一部分污染物便被“壓制”在了回流區(qū)及客廳西墻附近區(qū)域[14]。客廳絕大部分區(qū)域空氣質(zhì)量符合文獻(xiàn)[13]要求，工況4時(shí)客廳的空氣質(zhì)量明顯要優(yōu)于其他3個(gè)工況。各工況下，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)變化，從而對(duì)污染物遷移產(chǎn)生影響，但這一影響遠(yuǎn)不及溫度場對(duì)污染物濃度場影響。3#房間作為各工況通風(fēng)的唯一出流通道，對(duì)于污染物的時(shí)空分配起著重要作用，其污染物濃度梯度也最為復(fù)雜。在北墻內(nèi)側(cè)不同范圍的污染物集聚可能與溫度場以及該區(qū)域的回流對(duì)污染物擴(kuò)散過程影響有關(guān)，但各工況下4#房間內(nèi)污染物濃度均未超標(biāo)。

圖6 Z=1.1 m截面C7H8濃度場 Fig.6 C7H8 concentration profile in Z=1.1 m section

3.3 舒適性評(píng)價(jià)

舒適性是居住者對(duì)客觀環(huán)境的主觀反映，是室內(nèi)流場、溫度場、相對(duì)濕度以及居住者自身?xiàng)l件綜合作用的結(jié)果。圖7給出了工況1～4沿房間高度不同水平面上的舒適性評(píng)價(jià)結(jié)果。

圖7 PMV平均值沿房間高度變化曲線 Fig. 7 Average value of PMV in vertical direction

可以看出，4種工況的舒適性表現(xiàn)出了相似的變化規(guī)律：PMV值沿房間高度方向逐漸增大，這主要是室內(nèi)溫度分層所致。根據(jù)PMV-PPD的熱舒適評(píng)價(jià)指標(biāo)及ISO7730標(biāo)準(zhǔn)規(guī)定：以PPD≤10%為設(shè)計(jì)依據(jù)[15]，即90%以上的人感到滿意的熱環(huán)境為熱舒適環(huán)境，此時(shí)，對(duì)應(yīng)的PMV= -0.5～+0.5。但考慮中國的實(shí)際經(jīng)濟(jì)狀況，一般計(jì)算時(shí)都將該值放寬至PPD≤20%[16]，此時(shí)，對(duì)應(yīng)的PMV= -0.75～+0.75。所以，根據(jù)中國對(duì)舒適性要求的PMV值以及文獻(xiàn)[15]中關(guān)于熱感覺的7級(jí)標(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對(duì)應(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能滿足國際上對(duì)舒適性要求的PMV值，但房間整體偏涼。對(duì)于寒冷地區(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)于排污效率和舒適性評(píng)價(jià)的討論得到了這樣一個(gè)結(jié)論：感覺舒適的環(huán)境未必對(duì)健康有利，反之亦然。從圖3(a)可以看出，工況1整個(gè)通風(fēng)過程的排污效率均是最高的，工況4次之，工況3最低。從圖3(b)可以看出，通風(fēng)結(jié)束時(shí)，工況1、4的排污效率相當(dāng)，均高于其它工況。綜合通風(fēng)過程的排污效率變化以及室內(nèi)舒適性評(píng)價(jià)結(jié)果，工況4為最佳通風(fēng)時(shí)段。居住者如果能通過改變衣著來適應(yīng)室內(nèi)偏涼的環(huán)境，工況1也是比較理想的通風(fēng)時(shí)段。

4 結(jié) 論

以蘭州某民用住宅為研究對(duì)象，采用Realizablek-ε模型對(duì)冬季自然通風(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)出了相似的變化趨勢，但對(duì)污染物的排除效果存在差異。相同的開窗方式，進(jìn)風(fēng)溫度對(duì)室內(nèi)污染物的時(shí)空分布有顯著影響。

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

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

[1] 周懷改. 通風(fēng)對(duì)建筑物室內(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)下門窗開啟對(duì)室內(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,?akcF 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ì)算流體動(dòng)力學(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ù)對(duì)室外氣象參數(shù)瞬時(shí)變化的動(dòng)態(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)

(編輯 王秀玲)

Optimization analysis of the period of time of natural ventilation in winter in cold zone

WangYea,b,WangJingwena

(a. School of Environmental and Municipal Engineering;b. Key Laboratory of Railway Vehicle Thermal Engineering,Ministry of Education of China,Lanzhou Jiaotong University,Lanzhou 730070,P.R. China)

In order to find out the optimal period of time of natural ventilation in winter,the realizablek-ε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.

natural ventilation; unsteady heat transfer; ventilation efficiency; optimization

10.11835/j.issn.1674-4764.2015.05.015

2015-06-29 基金項(xiàng)目：國家自然科學(xué)基金(51266004、51476073)；甘肅省自然科學(xué)基金(1308RJZA151)

王燁(1972-),男,教授,博士,主要從事暖通空調(diào)、強(qiáng)化傳熱研究，(E-mail)wangye@mail.lzjtu.cn。

Foundation item：National Natural Science Foundation of China(No.51266004,51476073);National Natural Science Foundation of Gansu Province of China(No.1308RJZA151)

TU832

A

1674-4764(2015)05-0101-08

Received:2015-06-29

Author brief：Wang Ye (1972-),PhD,professor,main research interests:heating ventlation and air conditioning,heat transfer enhancement，(E-mail)wangye@mail.lzjtu.cn.