楊秀峰，鐘 珂，亢燕銘

（1.東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 201620；2.揚(yáng)州大學(xué) 能源與動力工程學(xué)院，江蘇 揚(yáng)州 225127）

熱壓自然通風(fēng)房間的瞬態(tài)污染狀況分析

楊秀峰1，2，鐘 珂1，亢燕銘1

（1.東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 201620；2.揚(yáng)州大學(xué) 能源與動力工程學(xué)院，江蘇 揚(yáng)州 225127）

熱壓自然通風(fēng)向穩(wěn)態(tài)發(fā)展的過程中室內(nèi)空氣污染狀況必然會不斷變化.在熱壓自然通風(fēng)瞬態(tài)模型的基礎(chǔ)上，給出了自然通風(fēng)的污染物輸送模型，理論分析了典型條件下通風(fēng)房間內(nèi)污染物濃度的瞬變過程，討論了主要因素對室內(nèi)瞬態(tài)污染物濃度的影響.研究表明，室內(nèi)熱空氣層的污染物濃度和原有冷空氣層的污染物濃度均是先升高后降低，且前者的最高值大于后者.增大熱源浮升力通量或有效通風(fēng)面積、減小房間面積或房間高度均可以加速室內(nèi)污染物濃度的變化，但是增大熱源浮升力通量或減小房間面積對污染物濃度的峰值沒有影響，而增大有效通風(fēng)面積或減小房間高度卻會使污染物濃度峰值升高.

熱壓自然通風(fēng)；瞬態(tài)；污染物濃度；峰值

室內(nèi)熱源產(chǎn)生的浮升力驅(qū)動的熱壓自然通風(fēng)是自然通風(fēng)的一種重要形式，其基本過程是熱源加熱周圍空氣，形成向上浮升的羽流，羽流浮升力產(chǎn)生的煙囪效應(yīng)使室外空氣由下部通風(fēng)口流入，室內(nèi)上部的熱空氣由上部通風(fēng)口排出，室內(nèi)逐漸形成垂直溫度梯度和空氣密度分層［1－4］.

在室內(nèi)外空氣污染物濃度不同的條件下，自然通風(fēng)過程必然會導(dǎo)致室內(nèi)污染物濃度不斷變化.文獻(xiàn)［5－7］分析了熱壓自然通風(fēng)發(fā)展至穩(wěn)態(tài)后室內(nèi)污染物濃度隨時間的衰減規(guī)律.其中，文獻(xiàn)［5］研究了熱源局部分布條件下自然通風(fēng)氣流對室內(nèi)污染物的沖刷過程，文獻(xiàn)［6］則比較了熱源均勻分布和局部分布情形下通風(fēng)氣流對室內(nèi)污染狀況的影響，二者均假設(shè)室內(nèi)熱空氣層和冷空氣層的污染物均勻分布.文獻(xiàn)［7］根據(jù)自然置換通風(fēng)情形下室內(nèi)水平截面的運(yùn)動特性，分別給出了室內(nèi)熱空氣層和冷空氣層內(nèi)的污染物輸送模型，并與假設(shè)污染物充分混合的理論模型進(jìn)行了比較.然而，對教室和劇院等非長時間使用的公共建筑而言，熱壓自然通風(fēng)向穩(wěn)態(tài)發(fā)展的過程也很重要，這是因為該過程相對于建筑的使用時間而言并不簡短［8－11］.因此，了解熱壓自然通風(fēng)向穩(wěn)態(tài)發(fā)展的過程中室內(nèi)空氣污染狀況的瞬態(tài)變化，對全面評價熱壓自然通風(fēng)對室內(nèi)環(huán)境的影響具有重要的現(xiàn)實意義.

本文旨在分析室內(nèi)外初始溫度相同的條件下，熱壓自然通風(fēng)向穩(wěn)態(tài)發(fā)展的過程中室內(nèi)氣態(tài)污染物濃度的變化規(guī)律.在瞬態(tài)熱壓通風(fēng)理論模型［12］的基礎(chǔ)上，給出了一個熱壓自然通風(fēng)瞬態(tài)發(fā)展過程的氣態(tài)污染物輸送模型，并通過數(shù)值求解模型控制方程，討論分析熱源浮升力通量、有效通風(fēng)面積、房間面積和高度等因素對室內(nèi)瞬態(tài)污染物濃度的影響.

1 熱壓自然通風(fēng)瞬態(tài)發(fā)展的理論模型

由于熱壓自然通風(fēng)瞬態(tài)發(fā)展過程的復(fù)雜性，關(guān)于這方面的理論研究并不多見，其中最具代表性的是文獻(xiàn)［8］給出的理論模型.文獻(xiàn)［8］研究了室內(nèi)外初始溫度相同、通風(fēng)口設(shè)置在房間頂部和底部、熱源位于室內(nèi)地面的情形下，熱壓自然通風(fēng)的瞬態(tài)發(fā)展過程.其將房間內(nèi)的氣流形式看作置換通風(fēng)，并假設(shè)通風(fēng)過程中上部熱空氣層充分混合，溫度均勻分布，給出了通風(fēng)過程中熱空氣層的體積守恒和浮升力守恒方程.事實上，熱空氣層內(nèi)一般存在著溫度差異，且位置越高處空氣溫度越高.文獻(xiàn)［8］也指出該模型關(guān)于熱空氣層充分混合的假設(shè)可能會低估排出空氣的浮升力.

圖1 改進(jìn)模型的空間熱分層示意圖Fig.1 Schematic diagram of the thermal stratification of the modified model

文獻(xiàn)［12］在文獻(xiàn)［8］的基礎(chǔ)上給出了熱壓自然通風(fēng)瞬態(tài)發(fā)展的改進(jìn)理論模型.改進(jìn)模型認(rèn)為室內(nèi)熱空氣層并不是充分混合的，而是由緊靠天花板的熱空氣薄層及其下方的熱空氣層兩部分組成.于是，通風(fēng)發(fā)展過程中室內(nèi)空間可以看作由豎向的3個區(qū)域組成，即下部的冷空氣層、中間的熱空氣層及天花板處的熱空氣薄層，如圖1所示.圖中Ta，T和Tc分別為冷空氣層、熱空氣層和熱空氣薄層的溫度，H，hc和h分別為房間高度和兩個熱分層界面的高度，q為瞬時通風(fēng)量.

其中：B0為熱源的浮升力通量，m4／s3；g′和g′p，h分別為熱空氣層的浮升力和熱分層界面處羽流的浮升力，m／s2；qp，h為熱分層界面處的羽流流量，m3／s；c為常數(shù)，c＝ （6α／5）（18απ2／5）1／3（α＝0.083，為卷吸系數(shù)）.

文獻(xiàn)［12］認(rèn)為，羽流穿過冷空氣層和熱空氣層到達(dá)天花板后沿天花板水平擴(kuò)散過程中與周圍熱空氣（浮升力為g′）相混合，因此，可將擴(kuò)散形成的熱空氣薄層的浮升力近似看作和g′的算術(shù)平均值，即

則通風(fēng)過程中熱空氣層（包括上部熱空氣薄層和中部熱空氣層）的體積守恒方程為

式中：A＊為有效通風(fēng)面積，m2；S為房間面積，m2；t為時間，s.

通風(fēng)過程中熱空氣層的浮升力守恒方程為

定義下列無量綱參數(shù)：

則式（2）和（3）的無量綱形式分別為

初始條件為τ＝0，ζ＝1，δ＝1.

2 污染物輸送模型

通風(fēng)房間的面積S一般遠(yuǎn)大于任意高度的羽流截面積，因此，羽流對周圍流體的卷吸作用可以近似看作水平的［13］.根據(jù)體積守恒，任意高度z處羽流外氣流垂直方向的平均速度v［7］可以表示為

v可能為正、負(fù)或0，其值與截面位置和時間有關(guān).文獻(xiàn)［12］中將熱分層界面到達(dá)最低高度ζover和名義穩(wěn)態(tài)高度ζn的時刻定義為通風(fēng)過程的特征時刻τ2和τ3.在τ＝0的初始時刻，地面處的v為0.在0＜τ＜τ2階段，地面處對應(yīng)的v＞0，熱分層界面高度h處對應(yīng)的v＜0，則熱分層界面和地面之間必然存在v＝0的截面，將其高度記為h0.在τ2時刻，熱分層界面下移至最低位置hover（無量綱高度為ζover），此時q＝qp，h，熱分層界面處的v＝0.在τ2＜τ＜τ3階段，q＞qp，h，熱分層界面處對應(yīng)的v＞0，說明v＝0的截面位于熱分層界面上方.隨著通風(fēng)過程逐漸趨于穩(wěn)態(tài)，熱分層界面和v＝0的截面均緩慢向上趨于穩(wěn)態(tài)熱分層界面處.熱壓自然通風(fēng)發(fā)展過程中，空間內(nèi)典型截面的垂直速度如圖2所示.

為簡化問題，常假設(shè)通風(fēng)口處的氣流為單向流，不考慮流入空氣與室內(nèi)空氣之間的混合，從而將熱壓通風(fēng)房間內(nèi)的氣流形式看作置換通風(fēng)［1，5，8］.于是，通風(fēng)過程中可以將室內(nèi)地面附近的新流入冷空氣層內(nèi)的氣態(tài)污染物看作均勻分布，其濃度與室外空氣中的氣態(tài)污染物濃度ka相等.假設(shè)房間上部的熱空氣層（包括熱空氣薄層和中間熱空氣層）和室內(nèi)原有冷空氣層內(nèi)的氣態(tài)污染物均勻分布，濃度分別記為ku和k1，則在0＜τ≤τ2的通風(fēng)階段，室內(nèi)空間由3個濃度不同的區(qū)域組成.當(dāng)τ＞τ2后，室內(nèi)原有冷空氣層消失，室內(nèi)空間由濃度分別為ku和ka的2個區(qū)域組成.當(dāng)τ＝τ3時，通風(fēng)發(fā)展至穩(wěn)態(tài)，通風(fēng)量和熱分界面高度分別取穩(wěn)態(tài)值qss和hss.

2.1 0＜τ≤τ2的階段

該階段室內(nèi)熱空氣層的污染物質(zhì)量守恒方程為

室內(nèi)原有冷空氣層的污染物質(zhì)量守恒方程為

垂直速度v＝0的截面距離地面的高度h0可由式（12）得到：

將室內(nèi)原有的污染物濃度記為k0，并在此基礎(chǔ)上定義下列無量綱濃度：

則式（10）和（11）的無量綱形式分別為

v＝0的截面的無量綱高度為

2.2 τ＞τ2的階段

如圖2所示，當(dāng)瞬態(tài)通風(fēng)過程發(fā)展至τ2時刻后，室內(nèi)原有冷空氣層消失，熱分層界面由最低高度處緩慢上移，界面處的垂直速度v＞0，而v＝0的截面位于熱分層界面上方.v＝0的截面將室內(nèi)空間分為兩個濃度不同的區(qū)域，截面以上區(qū)域的濃度為ku，截面以下區(qū)域的濃度為ka.上部區(qū)域內(nèi)污染物質(zhì)量的守恒方程可以表示為

式（17）的無量綱形式為

該階段的起始時刻為τ2，初始條件即為τ＝τ2，Ku＝Ku，τ2，其中Ku，τ2為τ2時刻室內(nèi)熱空氣層的污染物濃度.

3 熱壓通風(fēng)房間瞬態(tài)污染狀況分析

聯(lián)立熱壓通風(fēng)瞬態(tài)發(fā)展過程的改進(jìn)模型和上述污染物輸送模型，并利用 Mathematica 7.0軟件數(shù)值求解相關(guān)的微分方程組，數(shù)值方法為顯式RungeKutta法.

圖3所示為a＝0.2條件下熱壓通風(fēng)過程中熱分層界面高度ζ、v＝0截面的高度ζ0以及原有冷空氣層厚度ζ-ζ0的瞬時變化.由圖3可知，通風(fēng)過程開始后熱分層界面由天花板處向下移動，v＝0的截面由地面處向上移動，兩截面在τ2時刻完全重合.這表明通風(fēng)開始后室內(nèi)原有的空氣層越來越薄，直至消失.τ2時刻以后，兩截面均向穩(wěn)態(tài)熱分層界面處緩慢移動，且v＝0的截面稍高于熱分層界面.

圖3 a＝0.2條件下ζ，ζ0和ζ-ζ0的瞬時變化Fig.3 Dimensionless heightζ，ζ0andζ-ζ0，as a function ofτfor a＝0.2

圖4所示為a為0.2、室外污染物濃度Ka＝0條件下，熱壓通風(fēng)過程中室內(nèi)污染物濃度的瞬時變化.圖4表明，室內(nèi)熱空氣層的污染物濃度Ku和室內(nèi)原有冷空氣層的污染物濃度K1均先升高后降低，且K1先于Ku到達(dá)最高值，其最高值大于Ku的最高值.對于室內(nèi)原有冷空氣層而言，雖然通風(fēng)過程中進(jìn)入該層的羽流中不含污染物（Ka＝0），而離開該層的羽流中污染物濃度為K1，使得該層的污染物總量越來越少，但由于通風(fēng)起始階段該層的厚度ζ-ζ0迅速減?。ㄒ妶D3），使該層的污染物濃度顯著升高.隨著通風(fēng)的不斷進(jìn)行，室內(nèi)原有熱空氣層的污染物總量持續(xù)減小，該層的厚度變化漸趨平緩，使得該層的污染物濃度K1開始下降.對室內(nèi)熱空氣層而言，通風(fēng)起始階段進(jìn)入該層的羽流流量明顯大于離開該層的風(fēng)量，且羽流的污染物濃度（與K1相同）不小于離開該層空氣的污染物濃度（與Ku相同），從而使室內(nèi)熱空氣層內(nèi)的污染物總量顯著增加，污染物濃度上升.隨著通風(fēng)過程的進(jìn)行，雖然進(jìn)入熱空氣層的羽流流量仍大于通風(fēng)量，但羽流中的污染物濃度已開始低于排風(fēng)氣流的污染物濃度，從而使熱空氣層的污染物濃度Ku逐漸降低.當(dāng)Ku取最大值時，dKudτ＝0，則由式（7）和（14）可知Ku＝K1，這意味著圖4中兩曲線的交點即為Ku曲線的最高點.

圖4 a＝0.2、Ka＝0條件下室內(nèi)污染物濃度的瞬時變化Fig.4 Dimensionless pollutant concentration as a function of τfor a＝ 0.2and Ka＝0

圖5所示為a為0.2、室外污染物濃度Ka分別為0，0.2和0.5的條件下，熱壓通風(fēng)過程中室內(nèi)無量綱污染物濃度Ku和K1的變化過程.由圖5可知，室外污染物濃度Ka不同時，室內(nèi)瞬時污染物濃度的變化趨勢相似，且Ka越大，室內(nèi)瞬時污染物濃度越大，Ku和K1到達(dá)最高值的無量綱時間τ越長.圖5還表明，Ka越大，τ2時刻室內(nèi)原有冷空氣層的污染物濃度K1與Ka之間的偏離越大，這說明當(dāng)室外污染物濃度較大時，即使通風(fēng)過程進(jìn)行到室內(nèi)原有冷空氣層即將消失時，該層的污染物濃度也不能降低至室外污染物的濃度水平.

圖6所示為室外污染物濃度Ka為0、無量綱有效通風(fēng)面積a分別為0.2，0.5和1.0條件下熱壓通風(fēng)過程中室內(nèi)無量綱污染物濃度Ku和K1的變化過程.由圖6可知，a越大，Ku和K1的最高值越大，且由初始濃度上升到最高值以及由最高值下降到接近室外濃度的過程所需的無量綱時間越短.說明a越大，室內(nèi)污染物濃度Ku和K1隨無量綱時間τ的變化越快.

了解室內(nèi)污染物濃度隨時間t的演變過程更具有實際意義.由上述污染物輸送模型可知，室內(nèi)污染物濃度Ku和K1隨無量綱時間τ的變化取決于a和Ka以及瞬態(tài)通風(fēng)參數(shù)ζ，ζc，δ和δc.而由瞬態(tài)通風(fēng)過程的理論模型可知，瞬態(tài)通風(fēng)參數(shù)僅取決于無量綱有效通風(fēng)面積a，即在室外污染物濃度Ka一定的條件下，Ku和K1隨τ的變化僅與a有關(guān).因此，通風(fēng)過程中室內(nèi)污染物濃度Ku和K1隨時間t的變化與B0，S，A＊和H有關(guān).本文在典型建筑參數(shù)［5］的基礎(chǔ)上，利用瞬態(tài)熱壓通風(fēng)的理論模型［12］和污染物輸送模型分析B0，S，A＊和H對室內(nèi)瞬態(tài)污染狀況的影響.為簡單起見，以下分析過程中假設(shè)Ka＝0.

當(dāng)A＊＝0.139 2m2，S＝7.14m2，H＝3.0m，B0＝0.001 4，0.002 8，0.004 2m4／s3時，通風(fēng)過程中室內(nèi)熱空氣層的污染物濃度Ku和原有冷空氣層的污染物K1隨時間t的變化過程如圖7所示.圖7表明，B0變化對通風(fēng)過程中Ku和K1的最高值沒用影響，但B0越大，Ku和K1由初始濃度上升至最高值的過程和由最高值下降的過程均越短，即室內(nèi)污染物濃度的變化越迅速.這是因為當(dāng)B0增大而A＊，H和S不變時，無量綱通風(fēng)面積a不變，則無量綱濃度Ku和K1隨無量綱時間τ的變化過程不變.而由式（4）可知，B0越大，某無量綱時刻τ對應(yīng)的時間t越小，這表明雖然室內(nèi)污染物濃度變化到某值所需的無量綱時間τ不變，但所需的時間t越短.

當(dāng)A＊＝0.139 2m2，B0＝0.002 8m4／s3，H＝3.0m，S＝3.57，7.14，10.7m2時，熱壓通風(fēng)導(dǎo)致的室內(nèi)污染物濃度Ku和K1隨時間t的變化曲線如圖8所示.由圖8可知，S減小對室內(nèi)瞬態(tài)污染物濃度的影響與B0增大的情形類似，即S越小，室內(nèi)污染物濃度變化到某一水平所經(jīng)歷的時間t越短，但S變化不會對Ku和K1的峰值產(chǎn)生影響.由熱壓通風(fēng)過程的理論模型和污染物輸送模型可以得知，當(dāng)S變小而B0，A＊和H不變時，某個無量綱時刻τ對應(yīng)的熱壓通風(fēng)特性參數(shù)和無量綱污染物濃度均保持不變，而式（4）表明該無量綱時刻τ對應(yīng)的時間t減小.因此，S越小，熱壓通風(fēng)氣流對室內(nèi)污染物的沖刷過程越迅速有效.

當(dāng)B0＝0.002 8m4／s3，S＝7.14m2，H＝3m，A＊＝0.046 4，0.092 8，0.139 2m2時，熱 壓 通 風(fēng) 房 間 的瞬態(tài)污染物濃度曲線如圖9所示.圖9表明，A＊越大，室內(nèi)熱空氣層和原有冷空氣層的污染物濃度隨時間的變化率越大，Ku和K1的上升階段和下降階段均顯著縮短，這與圖7中B0增大及圖8中S減小的情形類似.但當(dāng)A＊增大而其他條件不變時，室內(nèi)污染物濃度Ku和K1的峰值也相應(yīng)增大，該現(xiàn)象與圖7和8的現(xiàn)象并不相同.

當(dāng)B0＝0.002 8m4／s3，S＝7.14m2，A＊＝0.139 2m2，H＝2.5，3.0，4.0，5.0m時，熱壓通風(fēng)過程中室內(nèi)的瞬態(tài)污染物濃度如圖10所示.由圖10可知，當(dāng)空間高度H減小而其他因素不變時，熱壓通風(fēng)可以使室內(nèi)污染物濃度在相同時間內(nèi)的變化幅度明顯增大.此外，圖10還表明，H越小，通風(fēng)過程中室內(nèi)污染物濃度的最高值越大，這與圖9的現(xiàn)象類似.但不同H對應(yīng)的Ku和K1變化曲線均分別經(jīng)過除初始點之外的共同點，該現(xiàn)象與圖9明顯不同.

4 結(jié) 語

本文在熱壓通風(fēng)瞬態(tài)發(fā)展改進(jìn)模型的基礎(chǔ)上，給出了通風(fēng)房間的污染物輸送模型，并通過聯(lián)立求解兩模型的相關(guān)方程，討論分析了典型條件下熱壓通風(fēng)發(fā)展過程中室內(nèi)污染物濃度的變化規(guī)律.本文的主要結(jié)論如下：

（1）通風(fēng)過程中室內(nèi)熱空氣層的污染物濃度Ku和室內(nèi)原有冷空氣層的污染物濃度K1均是先升高后降低，且K1先于Ku到達(dá)最高值，其最高值大于Ku的最高值，當(dāng)Ku取其最高值時，K1與Ku相等；

（2）室外污染物濃度Ka越大，室內(nèi)瞬時污染物濃度越大，Ku和K1到達(dá)最高值的無量綱時間τ越長；

（3）無量綱有效通風(fēng)面積a越大，Ku和K1的最高值均越大，而Ku和K1的上升過程和下降過程所需的無量綱時間越短；

（4）增大B0、減小S、增大A＊或減小H均可以加速熱壓通風(fēng)房間內(nèi)污染物濃度的變化，但是增大B0或減小S對污染物濃度的峰值沒有影響，而增大A＊或減小H卻會使污染物濃度峰值升高.

參 考 文 獻(xiàn)

［1］LINDEN P F，LANE－SERFF G F，SMEED D A.Emptying filing boxes：The fluid mechanics of natural ventilation ［J］.Journal of Fluid Mechanics，1990，212：309－335.

［2］ANDERSEN K T. Theoretical considerations on natural ventilation by thermal buoyancy［J］.ASHRAE Transactions，1995，101（2）：1103－1117.

［3］LI Y，DELSANTE A，SYMONS J.Simulation tools for analysing natural ventilation of large openings ［J］.AIRAH Journal，1997，51（11）：21－28.

［4］LI Y.Buoyancy－driven natural ventilation in a thermally stratified one－zone building ［J］.Building and Environment，2000，35（3）：207－214.

［5］HUNT G R，KAYE N B.Pollutant flushing with natural displacement ventilation［J］.Building and Environment，2006，41（9）：1190－1197.

［6］KAYE N B，HUNT G R.Heat source modeling and natural ventilation efficiency［J］.Building and Environment，2007，42（4）：1624－1631.

［7］BOLSTER D T，LINDEN P F.Contaminants in ventilated filling boxes［J］.Journal of Fluid Mechanics，2007，591：97－116.

［8］KAYE N B，HUNT G R.Time－dependent flows in an emptying filling box ［J］.Journal of Fluid Mechanics，2004，520：135－156.

［9］FITZGERALD S D， WOODS A W. Transient natural ventilation of a room with a distributed heat source［J］.Journal of Fluid Mechanics，2007，591：21－42.

［10］BOLSTER D，MAILLARD A，LINDEN P.The response of natural displacement ventilation to time－varying heat sources［J］.Energy and Buildings，2008，40（12）：2099－2110.

［11］FITZGERALD S D， WOODS A W. Transient natural ventilation of a space with localised heating ［J］.Building and Environment，2010，45（12）：2778－2789.

［12］YANG X F，WANG G X，ZHONG K，et al.Transient pollutant flushing of buoyancy－driven natural ventilation ［J］.Building Simulation，2012，5（2）：147－155.

［13］BAINES W D，TURNER J S.Turbulent buoyant convection from a source in a confined region ［J］.Journal of Fluid Mechanics，1969，37：51－80.

Analysis on Transient Air Pollution in Buoyancy－Induced Naturally Ventilated Enclosures

YANGXiu－feng1，2，ZHONGKe1，KANGYan－ming1
（1.School of Environmental Science and Engineering，Donghua University，Shanghai 201620，China；2.School of Energy and Power Engineering，Yangzhou University，Yangzhou Jiangsu 225127，China）

Indoor air pollution is influenced continuously by the airflow during the development of thermal natural ventilation to its steady state.The theoretical model of pollutant transport during natural ventilation is presented on the basic of the transient model of thermal natural ventilation to analyze the time evolution of indoor pollutant concentration and to explore the effects of main factors on the transient pollutant concentration during ventilation.The results indicate that the pollutant concentrations of the warm air layer and the original air layer both increase earlier and drop later.The maximal concentration of the warm air layer is shown to be greater than that of the original air layer.The variation of indoor pollutant concentration with time can be accelerated by increasing the source buoyancy flux or effective vent area，or decreasing the enclosure floor area or enclosure height.The increase of the source buoyancy flux or the decrease of the enclosure floor area would not affect the peak value of the transient pollutant concentration，which may be elevated by increasing of effective vent area or decreasing of the enclosure height.

thermal natural ventilation；transient；pollutant concentration；peak value

TU 834.1

A

1671－0444（2012）06－0750－08

2012－06－01

江蘇省高校自然科學(xué)基金資助項目（09KJB560004）；國家自然科學(xué)基金資助項目（51278094）

楊秀峰（1976—），男，江蘇泰興人，講師，博士研究生，研究方向為室內(nèi)空氣品質(zhì).E－mail：xfyang＠yzu.edu.cn

鐘 珂（聯(lián)系人），女，教授，E－mail：zhongkeyx＠dhu.edu.cn