華鳳皎， 亢燕銘， 鐘 珂

(東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 201620)

建筑體型廓線對建筑物表面粒子沉積特性的影響

華鳳皎， 亢燕銘， 鐘 珂

(東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 201620)

選取典型建筑體型，利用計算流體力學(xué)方法研究了建筑體型對氣溶膠粒子在建筑表面沉積速度和沉積通量的影響．結(jié)果表明：迎風(fēng)面寬度是影響氣溶膠粒子在建筑表面沉積速度和沉積通量的重要因素，且迎風(fēng)面和背風(fēng)面對其的影響規(guī)律不同；建筑表面上氣溶膠粒子的沉積通量隨著其粒徑的變化而變化，且不同建筑表面最小沉積通量對應(yīng)的粒徑不同．

建筑體型； 建筑外表面；慣性沉積； 沉積速度； 沉積通量； 氣溶膠粒子

目前，城市大氣霾污染已成為我國重要的環(huán)境問題．工業(yè)過程、機動車尾氣以及其他人類活動產(chǎn)生的氣溶膠粒子，不僅對人體健康產(chǎn)生威脅，也將對建筑物、構(gòu)筑物以及其他各種市政設(shè)施產(chǎn)生嚴重污染腐蝕和損害，如建筑表面的粒子沉積會影響其圍護結(jié)構(gòu)的熱工性能及美觀要求．近幾年，隨著霾天氣在我國大范圍地區(qū)的頻繁出現(xiàn)，這種氣溶膠污染效應(yīng)也變得越發(fā)嚴重．即使大氣中有害氣溶膠已被有效清除，沉積在建筑物表面的粒子仍然會長時間地影響建筑物附近的空氣質(zhì)量．為了減少氣溶膠粒子造成的危害，須對這些粒子沉積特征有清楚的認識．此外，古建筑的保護也需要了解和預(yù)測建筑各表面的粒子沉積情況，通過對建筑物墻體進行及時清洗，減弱氣溶膠對文物和古跡的腐蝕．

粒子在固體表面的沉積特性研究報道較多.如文獻[1]給出了封閉空間內(nèi)均勻湍流粒子沉積的解析解．文獻[2]對粒子的沉積速度進行了參數(shù)化分析，給出了計算粒子干沉積速度的簡化公式．文獻[3]采用多層箱體模型模擬了包含多種地表形態(tài)的大尺度區(qū)域粒子的沉積通量，但并沒有求解區(qū)域內(nèi)的流動情況．文獻[4]采用離散軌跡模型研究了不同通風(fēng)房間內(nèi)粒子的質(zhì)量濃度分布和沉積特性．文獻[5]研究了城市環(huán)境下墊面、建筑物墻體和屋頂?shù)臍馊苣z粒子沉積速度，指出僅用一個或幾個物理量來表示城市氣溶膠粒子沉積速度的方法過于簡化了對問題的分析．

上述文獻側(cè)重于研究粒子沉積速度與氣流分布的關(guān)系以及室內(nèi)粒子的沉積，均沒有考慮建筑體型對室外粒子在墻面沉積的影響．通常情況下，粒子的沉積效應(yīng)主要受擴散和慣性的共同作用[6]．由于不同建筑物體型對氣流的阻擋效應(yīng)和其周圍流場特征不同，導(dǎo)致了粒子在不同體型建筑周圍隨氣流運動時的動力學(xué)行為的差異．

為考察不同的建筑體型廓線下粒子在建筑物表面的沉積特征，本文將針對城市中常見的點式建筑和條形建筑體型，采用計算流體力學(xué)方法構(gòu)建粒子沉積模型，分析建筑體型對其圍護結(jié)構(gòu)外表面上的粒子沉積特性的影響．

1 數(shù)值計算模型和方法

1．1 數(shù)值計算模型

表1給出了典型建筑體型的特征參數(shù)，其中，建筑體型尺寸長為L，寬為W，高度為H．計算區(qū)域的大小對于計算量和模擬結(jié)果的準確性會造成影響，當(dāng)計算區(qū)域長度和高度不小于建筑物高度的5～6倍時，可消除邊界效應(yīng)對計算結(jié)果的影響[7-8]．因此，本研究中計算區(qū)域高度取為6H，入口距建筑物迎風(fēng)面的距離為5H，出口距建筑物背風(fēng)面距離為10H．本文建筑體型case 2 L的計算區(qū)域如圖1所示．

表1 建筑體型尺寸

圖1 模型計算域及網(wǎng)格示意圖Fig．1 Schematic diagram of the computational domain and enlarged mesh

采用Gambit軟件對求解區(qū)域進行非結(jié)構(gòu)化網(wǎng)格劃分，并對建筑外表面附近區(qū)域進行網(wǎng)格局部加密處理，而其他區(qū)域則采用正常尺寸網(wǎng)格來優(yōu)化整體網(wǎng)格數(shù)量并節(jié)約數(shù)值計算成本．模擬計算前，進行了網(wǎng)格獨立性檢驗，即增加網(wǎng)格數(shù)，直至網(wǎng)格數(shù)的增加對計算結(jié)果沒有影響為止[9]．經(jīng)檢驗，最終網(wǎng)格數(shù)量為163～291萬．

1．2 計算方法和邊界條件的設(shè)定

湍流流場可以采用不同的湍流模型進行模擬，為節(jié)約計算成本并獲得簡潔合理的數(shù)值結(jié)果，根據(jù)文獻[10-13]可知，標準κ-ε模型能很好地模擬本文條件下的湍流流場，其中壁面湍流特征用標準壁面函數(shù)描述，采用有限容積法離散控制方程，并對動量、湍流動能κ、能量耗散率ε的輸運方程采用二階迎風(fēng)格式進行差分，以提高計算精度．

污染物輸送質(zhì)量濃度變化的對流擴散方程為

(1)

式中：Ci為污染物質(zhì)量濃度；K為粒子擴散系數(shù)；Si為污染源項；uj為j方向的速度分量；xj為笛卡爾坐標分量．

計算域的進口邊界條件設(shè)為速度入口，風(fēng)速大小采用速度廓線形式．風(fēng)速的垂直特性與地形及地面粗糙度有關(guān)，速度廓線可由指數(shù)律分布[14-16]表示，如式(2)所示．

(2)

式中：z為任一高度， m；zs為參考高度， m；u(zs)為zs處的參考風(fēng)速， m/s， 考慮到城市區(qū)域的穩(wěn)定氣象條件，本文取20 m高度處的風(fēng)速為2 m/s；α為地面粗糙系數(shù)(風(fēng)剖面冪指數(shù))，本文α取0．25[17-18]．速度廓線采用C++語言用戶編程并與Fluent程序提供的UDF(user-defined functions)對接．

計算區(qū)域出口處可視為充分發(fā)展的湍流，采用出流邊界條件(outflow)，即該流出區(qū)域的質(zhì)量流率和流入?yún)^(qū)域的質(zhì)量流率相等．建筑表面及地面采用無滑移條件，其余各面則采用滑移速度邊界．

假設(shè)城市背景氣流中粒子質(zhì)量濃度均勻，因此，設(shè)每種粒徑粒子在計算域進口處空氣中的質(zhì)量濃度為10 μg/m3．城市大氣中，氣溶膠粒子的粒徑范圍很廣，生物有害性粒子尺寸(空氣動力學(xué)直徑)通常在1～10 μm范圍，其中2 μm左右及以下的粒子占絕大多數(shù)[3]，粒子密度約為1 300 kg/m3[19]．故本文選取粒子密度取ρp= 1 500 kg/m3，對粒徑dp=0．5， 1．0， 1．5， 2．5， 3．5， 5．0， 7．0和10．0 μm分別進行計算和分析．

1．3 粒子輸運理論基礎(chǔ)和計算方法

粒子在空氣中的輸運模擬采用歐拉法，選用混合模型求解第二相粒子的動量和體積分數(shù)．大氣環(huán)境中粒子可視為稀薄相，只考慮流場對顆粒物的單向耦合作用．壁面處的粒子沉積通量J可用式(3)計算．

(3)

式中：CVOF-dA為靠近壁面控制體中粒子相體積分數(shù)；vd-dA為某一點粒子的沉積速度； dA為控制體壁面的面積；Awall為建筑物表面面積．

建筑物不同外表面上的粒子沉積速度vd為垂直面上的沉積速度vdv

(4)

朝上的水平面沉積速度vdu

(5)

式中：vs為粒子的終端沉降速度， m/s；u*為摩擦速度， m/s； 參量I的具體表達式見文獻[20]；u*為通過流體的湍流動能k， m2/s，由式(6)計算得到．

(6)

1．4 數(shù)值模型的合理性驗證

為了保證數(shù)值模擬方法分析建筑流場分布特征的可靠性，需要對所采用數(shù)學(xué)模型進行驗證．本文利用漢堡大學(xué)氣象研究所提供的BLASIUS風(fēng)洞試驗數(shù)據(jù)庫[21]對速度場進行驗證．模擬參數(shù)與試驗保持一致，模型采用邊長為125 mm的立方體，入口風(fēng)速采用速度廓線，參考高度500 mm處的參考風(fēng)速為6 m/s，風(fēng)速廓線指數(shù)為0．22．采用文獻[22]的試驗結(jié)果驗證濃度場，具體如下：風(fēng)洞中的模型尺寸為60 mm×60 mm，其長度與風(fēng)洞長度相等，充分發(fā)展區(qū)的風(fēng)速廓線指數(shù)為0．20，參考高度650 mm處的參考風(fēng)速為3 m/s．污染物質(zhì)量濃度采用無因次濃度C*來表示，如式(7)所示．

(7)

式中：Ci為i點處的污染物質(zhì)量濃度， kg/m3；q0為污染物散發(fā)強度， kg/s；Ls為污染源長度，取0．9 m．

建筑迎風(fēng)區(qū)和背風(fēng)區(qū)的數(shù)值計算結(jié)果與風(fēng)洞試驗結(jié)果如圖2所示，其中u(z)為沿?zé)o量綱高度Z(Z=z/H)方向的氣流速度分量．由圖2(a)可以看出，本文模擬結(jié)果在建筑物高度范圍內(nèi)與風(fēng)洞試驗結(jié)果吻合得比較好．圖2(b)給出的污染物質(zhì)量濃度分布驗證說明，迎風(fēng)區(qū)質(zhì)量濃度計算值與實測值符合得較好，背風(fēng)區(qū)幾乎完全吻合．因此，本文采用的數(shù)值模型可以很好地預(yù)測城市中單體建筑物附近的氣流流動狀況和粒子質(zhì)量濃度分布情況．

(a) 氣流速度

(b) 質(zhì)量濃度

2 模擬結(jié)果與分析

2．1 建筑體型對建筑附近氣流速度的影響

對表1給出的各建筑體型逐一進行了模擬．在case 1L和case 3L兩種典型情況下，在z=5 m高度處建筑周圍的氣流流線分布如圖3所示．

圖3 典型建筑z=5 m平面上的氣流流線分布Fig．3 Streamline distributions on typical planes of the buildings at z=5 m

由圖3可以看出，由于建筑物對來流的阻擋，導(dǎo)致在建筑側(cè)面及背風(fēng)面區(qū)域均會形成渦流．由于case 3L的迎風(fēng)面寬度較大，其背風(fēng)面形成的回漩渦流區(qū)域也明顯大于case 1L． case 3L側(cè)面近壁區(qū)完全處于渦流狀態(tài)，而case 1L側(cè)面較長，渦流主要分布在擾流區(qū)前端．二次流的分布不同可能會導(dǎo)致粒子沉積到建筑物表面的概率不同，從而造成建筑表面粒子沉積通量的差異．圖3還表明，不同體型的建筑物附近的氣流流動速度不同，可能會導(dǎo)致建筑物不同表面的粒子沉積速度有所不同．

為定量分析建筑各表面附近氣流速度的差異，圖4給出了距離建筑物背風(fēng)面、迎風(fēng)面和側(cè)面0．25 m處投影平面的垂直中心線上氣流速度沿?zé)o量綱高度Z的變化情況．

(a) 迎風(fēng)面

(b) 背風(fēng)面

(c) 側(cè)面

由圖4(a)可見，在迎風(fēng)面附近，高度相同、體型不同的建筑物附近氣流速度變化曲線相似，僅在Z<0．3的近地區(qū)域略有不同．高度不同的建筑物附近氣流速度變化趨勢也基本一致，這表明建筑體型和高度對迎風(fēng)面中心附近氣流速度的影響很?。畯膱D4(b)和4(c)中可看出，建筑物背風(fēng)面和側(cè)面的氣流速度受建筑體型和迎風(fēng)面寬度影響很大．不同體型的建筑物背風(fēng)面在Z<0．3的近地區(qū)域氣流速度基本相等，上部流速隨著迎風(fēng)面寬度增大而增大，且建筑物越高，規(guī)律越明顯．另外，較高建筑的背風(fēng)區(qū)氣流速度也高于相同平面結(jié)構(gòu)的較低建筑．

由圖4(c)還可以看到，在case 2L/case 2H和case 3L/case 3H兩種情況下，建筑側(cè)面近地區(qū)域的氣流速度很接近，且遠大于case 1L/case 1H的氣流速度．這是因為前者建筑側(cè)面尺寸較小，建筑側(cè)面完全位于渦流區(qū)內(nèi)，流速受渦流影響而明顯增大，而后者側(cè)面部分的前端區(qū)域位于渦流區(qū)內(nèi)，而所考察的速度樣本點位置已在渦流區(qū)外，故氣流速度較?。@一差別將會對粒子在建筑物表面的沉積特性產(chǎn)生影響．

2．2 建筑各表面上的粒子沉積通量

由數(shù)值模擬結(jié)果和式(3)沉積通量的計算可以得到表1中6種情形的迎風(fēng)面、側(cè)面、背風(fēng)面和屋頂上粒子沉積通量J與粒徑dp的關(guān)系，結(jié)果如圖5所示．

由圖5(a)可知，不論建筑高低，均表現(xiàn)為迎風(fēng)面越寬，迎風(fēng)面上的沉積通量越小，此外，粒子粒徑越小，沉積通量受建筑體型影響越小.這是因為小粒子的沉積主要受擴散機理控制，而大粒子的沉積主要受慣性力控制．由于建筑體型是通過影響周圍流場來影響粒子沉積效果的，因此只能影響到受慣性力控制的大粒子的沉積效果．

(a) 迎風(fēng)面和側(cè)面

(b) 背風(fēng)面

(c) 屋頂

(d) 建筑物總表面

另外，比較不同建筑體型中迎風(fēng)面最小沉積通量對應(yīng)的粒子直徑可以看到，除了迎風(fēng)面面積較小的case 1H迎風(fēng)面上的最小沉積通量對應(yīng)粒徑為2．5 μm外，其他體型建筑對應(yīng)粒徑均為3．5 μm．這是因為case 1H對氣流的阻擋作用小，迎風(fēng)面附近氣流平均速度大于其他建筑體型的情況，粒子受到的慣性作用大，故導(dǎo)致最小沉積通量對應(yīng)的粒徑小于其他建筑體型．由圖5(a)還可以看到，雖然側(cè)面粒子沉積通量略小于迎風(fēng)面，但沉積通量隨粒徑的變化曲線與迎風(fēng)面非常相似，最小沉積通量對應(yīng)的粒徑也是3．5 μm．

由圖5(b)可見，背風(fēng)面上沉積通量隨粒徑的變化規(guī)律為，迎風(fēng)面越寬，背風(fēng)面粒子沉積通量越大．此外，對比圖5(a)和5(b)還可以看到，背風(fēng)面上的沉積通量明顯小于迎風(fēng)面，同時最小沉積通量對應(yīng)粒徑較大，約為5．0～7．0 μm．這是因為背風(fēng)面附近氣流速度遠低于迎風(fēng)面，故粒子在這里受到的慣性力較小，于是最小沉積通量對應(yīng)的dp必然增大．由圖5(c)可看出，在重力沉降的作用下，屋頂上的沉積通量隨著粒徑的增加單調(diào)上升．由于屋頂表面附近的局地風(fēng)速受建筑體型的影響較小，故不同建筑體型的沉積通量差別很小．由圖5(d)可知，由于屋頂水平面上的沉積通量遠大于其他垂直表面的沉積通量，建筑物表面總的沉積通量主要取決于屋頂上的沉積通量，因此，其與屋頂上呈現(xiàn)出相同的沉積規(guī)律．

2．3 建筑體型對沉積速度的影響

根據(jù)氣溶膠動力學(xué)原理[23-24]，粒子的沉積速度與建筑物附近氣流速度密切相關(guān)．為了對建筑物表面沉積情況進行總體評價，在距離建筑物表面0．25 m的投影平面上均勻取點，每個點間距1 m，取值點總數(shù)為442～2 502，提取這些取樣點處的數(shù)值結(jié)果來討論氣溶膠粒子在外墻面的沉積規(guī)律． 根據(jù)流場的數(shù)值模擬結(jié)果和式(4)和(5)，在不同建筑體型下，dp=0．5和10．0 μm的粒子在建筑迎風(fēng)面、背風(fēng)面、側(cè)面和屋頂上的沉積速度統(tǒng)計結(jié)果如圖6所示．

(a) 迎風(fēng)面

(b) 背風(fēng)面

(c) 側(cè)面

(d) 屋頂

由圖6(a)可見，在兩種建筑高度下，均表現(xiàn)為迎風(fēng)面越窄，粒子沉積速度越大的特征．由圖6(b)可知，背風(fēng)面與迎風(fēng)面表現(xiàn)出的規(guī)律則恰好相反．這一規(guī)律在高建筑(即H=40 m)上表現(xiàn)得更為明顯．這可能是由于迎風(fēng)面尺度越大，背風(fēng)區(qū)的湍動能亦越大，故式(6)中的摩擦速度將增加，從而使得沉積速度增大．

圖6(c)中表明，粒子在建筑側(cè)面上沉積速度與迎風(fēng)面沉積速度在量級上相當(dāng)，但與建筑體型和迎風(fēng)面寬度的關(guān)系并沒有非常明確的相關(guān)性．這很可能是因為，不同體型的側(cè)面寬度不同，而建筑側(cè)面空氣流動情況受到迎風(fēng)面和側(cè)面尺寸的共同影響．由圖6(d)可見，dp=10．0 μm的粒子在屋頂上的水平沉積速度明顯高于dp=0．5 μm粒子的沉積速度.這是因為大粒子還會受到重力沉積的作用，從而強化了大粒子的沉積．此外，屋頂上的沉積速度與建筑體型和迎風(fēng)面寬度的相關(guān)性與背風(fēng)面相同，即沉積速度隨著迎風(fēng)面寬度的增大而增大．綜合以上分析，可以認為迎風(fēng)面寬度W與沉積速度vdv存在明確的相關(guān)性．

為評價建筑物迎風(fēng)面寬度W對粒子沉積效果的影響，現(xiàn)以迎風(fēng)面寬度最小的case 1L和case 1H為基準，定義其他建筑模型表面不同dp粒子的沉積速度和沉積通量變化率η，如式(8)所示．

(8)

式中：Y代表沉積速度vdv或沉積通量J；x=2， 3，分別為表1所示的case 2和case 3類型的建筑；i=L， H，分別表示低建筑和高建筑．

根據(jù)以上定義，圖7給出了不同建筑體型迎風(fēng)面上的粒子沉積通量變化率η與其粒徑的關(guān)系．

圖7 不同建筑迎風(fēng)面上不同粒子沉積通量 變化率隨其粒徑的變化Fig．7 Changed ratios of deposition flux with dpon windward of different buildings

由圖7可以看到，建筑迎風(fēng)面寬度的增大將使得所有粒徑的粒子在迎風(fēng)面的沉積通量均產(chǎn)生不同程度的減小，這與上述迎風(fēng)面寬度對粒子沉積速度的影響的結(jié)論相吻合．由圖7還可以看到，迎風(fēng)面寬度對粒子沉積通量的影響程度隨著粒徑和建筑高度的增大而增加．

為了比較建筑體型對粒子沉積速度和沉積通量的影響程度，圖8分別給出了dp=0．5和10．0 μm的粒子在建筑迎風(fēng)面和背風(fēng)面上沉積速度與沉積通量變化率η．

(a) 迎風(fēng)面

(b) 背風(fēng)面

由圖8(a)可以看到，迎風(fēng)面上的粒子沉積速度與沉積通量的變化率η的變化規(guī)律保持一致，隨著迎風(fēng)面寬度和建筑高度的增加，迎風(fēng)面上的粒子沉積速度和沉積通量都減小，且粒徑越大越明顯．由圖8(b)可以看到：對于較高的建筑，迎風(fēng)面寬度較大的體型使得背風(fēng)面上大粒子和小粒子沉積速度和沉積通量都增大，但沉積通量的增大幅度明顯低于沉積速度；對于較低建筑，背風(fēng)面沉積速度隨著迎風(fēng)面寬度增大而增加，但沉積通量卻出現(xiàn)了減小現(xiàn)象．這是因為沉積通量由沉積速度和粒子質(zhì)量濃度共同決定，沉積通量沒有完全隨著沉積速度的增大而增加，這可能是建筑表面粒子質(zhì)量濃度隨迎風(fēng)面寬度變化的規(guī)律與沉積速度不一致造成的．為此，圖9給出了不同建筑體型下建筑中心剖面上dp=0．5 μm粒子的質(zhì)量濃度分布．

圖9 不同建筑體型中心剖面上的粒子質(zhì)量濃度分布Fig．9 Particle concentration distributions on centre section of different buildings

由圖9可以看到，由于沉積效應(yīng)對粒子的清除作用，隨著迎風(fēng)面寬度的增大，背風(fēng)面附近的粒子質(zhì)量濃度下降，因此，使得粒子沉積通量隨迎風(fēng)面寬度增加的幅度減小(對于較高建筑)．盡管隨著建筑迎風(fēng)面寬度增加，粒子沉積速度增大，但當(dāng)沉積速度的增大量不足以抵消粒子質(zhì)量濃度的減小量時，沉積通量隨迎風(fēng)面寬度增大的變化趨勢就不再與沉積速度相同．

3 結(jié) 語

建筑物表面的粒子沉積不僅嚴重影響建筑外觀，而且在雨雪沖刷和太陽輻射的共同作用下，還會因外墻材料的腐蝕而改變圍護結(jié)構(gòu)的熱工特性和室內(nèi)熱環(huán)境．由于建筑物的體型對其周圍流場有重要的畸變作用，進而影響到隨背景空氣輸送的大氣氣溶膠在建筑物表面上的沉積特性，但相關(guān)的研究還很缺乏．因此，本文應(yīng)用粒子沉積模型和計算流體力學(xué)方法，對幾種典型體型廓線的建筑表面上粒子沉積特性進行了數(shù)值計算和分析，討論了不同表面上發(fā)生粒子沉積的特征．主要結(jié)論如下：

(1) 由于不同體型建筑對來流的阻擋作用不同，形成了不同的流場特征．本文模擬結(jié)果表明，建筑迎風(fēng)面尺寸對流場的影響作用大于建筑高度的影響效果．因此，建筑高度增大雖然增大了迎風(fēng)面面積，但對粒子沉積速度的影響沒有迎風(fēng)面寬度變化的作用明顯．

(2) 隨著建筑物迎風(fēng)面寬度的增大，建筑迎風(fēng)面上粒子沉積速度將會下降，但在背風(fēng)面上的沉積速度則會升高．建筑物迎風(fēng)面上的粒子沉積通量隨迎風(fēng)面寬度的變化規(guī)律與沉積速度一致，但對于背風(fēng)面，由于墻面附近粒子質(zhì)量濃度受迎風(fēng)面寬度的影響趨勢與沉積速度不同，沉積通量與沉積速度變化規(guī)律不完全一致．

(3) 由于不同建筑表面附近氣流慣性作用不同，故不同表面上最小沉積通量對應(yīng)的粒徑也不同．由于背風(fēng)面處在建筑尾流區(qū)，氣流平均速度低，故粒子的慣性和湍流擴散效應(yīng)均較弱，因此，背風(fēng)面最小沉積通量對應(yīng)的粒徑將大于迎風(fēng)面和側(cè)面．

[1] CORNER J， PENDLEBURY E D． The coagulation and deposition of a stirred aerosol[J]． Proceeding of the Physical Society， 1951，64(8)：645-654．

[2] ZHANG L， GONG S， PADRO J， et al． A size-segregated particle dry deposition scheme for an atmospheric aerosol module[J]． Atmospheric Environment， 2001，35(3)：549-560．

[3] MONFORTI F， BELLASIO R， CLAI G， et al． An evaluation of particle deposition fluxes to cultural heritage sites in Florence， Italy[J]． Science of the Total Environment， 2004，334/335(1)：61-72．

[4] ZHAO B， ZHANG Y， LI X T， et al． Comparison of indoor aerosol particle concentration and deposition in different ventilated rooms by numerical method[J]． Building and Environment， 2004，39(1)：1-8．

[5] JONSSON L， KARLSSON E， JONSSON P． Aspects of particulate dry deposition in the urban environment[J]． Journal of Hazardous Materials， 2008，153(1)：229-243．

[6] GREENFIELD S M． Rain scavenging of radioactive particulate matter from the atmosphere[J]． Journal of Atmospheric Sciences， 1957，14(2)：115-125．

[7] TOMINAGA Y， MOCHIDA A， YOSHIE R， et al． AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings[J]． Journal of Wind Engineering and Industrial Aerodynamics， 2008，96(10)：1749-1761．

[8] NAKAMURA Y，OKE T R． Wind temperature and stability conditions in an east-west oriented urban canyon[J]． Atmospheric Environment， 1988，22(12)：2691-2700．

[9] LEI Y G， HE Y L． Effects of baffle inclination angle on flow and heat transfer of a heat exchanger with helical baffles[J]． Chemical Engineering and Processing： Process Intensification， 2008，47(12)：2336-2345．

[10] LAUNDER B E， SPALDING D B． The numerical computation of turbulent flows[J]． Computer Methods in Applied Mechanics and Engineering， 1974，3(2)：269-289．

[11] NAZRIDOUST K， AHMADI G． Airflow and pollutant transport in street canyons[J]． Journal of Wind Engineering and Industrial Aerodynamics， 2006，94(6)：491-522．

[12] DI SABATINO S， BUCCOLIERI R， PULVIRENTI B， et al． Flow and pollutant dispersion in street canyons using FLUENT and ADMS-Urban[J]． Environ Mental Modeling & Assessment， 2007，13(3)：369-381．

[13] NG W Y， CHAU C K． A modeling investigation of the impact of street and building configurations on personal air pollutant exposure in isolated deep urban canyons[J]． Science of the Total Environment， 2014，468/469：429-448．

[14] CHEUNG J O P， LIU C H． CFD simulations of natural ventilation behaviour in high-rise buildings in regular and staggered arrangements at various spacings[J]． Energy and Buildings， 2011，43(5)：1149-1158．

[15] LEUNG K K， LIU C H， WONG C C C， et al． On the study of ventilation and pollutant removel over idealized two-dimensional urban street canyons[J]． Building Simulation， 2012，5(4)：359-369．

[16] RAJAPAKSHA I， NAGAI H， OKUMIYA M． A ventilated courtyard as a passive cooling strategy in the warm humid tropics[J]． Renewable Energy， 2003，28(11)：1755-1778．

[17] MOHAMED F Y．A wind tunnel study on the effect of thermal stability on flow and dispersion of rooftop stack emissions in the near wake of a building[J]． Atmospheric Environment， 2013，65：89-100．

[18] WANG H J， CHEN Q Y． A new empirical model for predicting single-sided， wind-driven natural ventilation in buildings[J]． Energy and Buildings， 2012，54：386-394．

[19] TISA L S， KOSHIKAWA T， GERHARDT P．Wet and dry bacterial spore densities determined by buoyant sedimentation[J]． Applied and Environmental Microbiology， 1982，43(6)：1307-1310．

[20] LAI A C K， NAZAROFF W W． Modeling indoor particle deposition from turbulent flow onto smooth surfaces[J]． Journal of Aerosol Science， 2000，31(4)：463-476．

[21] SCHATZMANN M． Compilation of experimental data for validation of microscale dispersion models[DB]． (1999-05-15)[2015-03-20]． http：//www．mi．zmaw．de/index．php?id=628．

[22] MERONEY R N， PAVAGEAU M， RAFADALIS S， et al． Study of line source characteristics for 2- D physical modelling of pollutant dispersion in street canyon[J]． Journal of Wind Engineering and Industrial Aerodynamics， 1996，62(1)：37-56．

[23] FRIEDLANDER S K， JOHNSTONE H F． Deposition of suspended particles from turbulent gas streams[J]． Journal of Industrial and Engineering Chemistry， 1957，49(7)：1151-1156．

[24] FUCHS N A． The mechanics of aerosols[M]． New York： Dover Publications Inc，1964：159-170．

(責(zé)任編輯：劉園園)

Effects of Building Configuration on the Particle Deposition Features on Building Surfaces

HUAFengjiao，KANGYanming，ZHONGKe

(School of Environmental Science and Engineering， Donghua University， Shanghai 201620， China)

The influences of building configuration on the deposition velocity and deposition flux of aerosol particles on the building surface were investigated by using computational fluid dynamics for typical building configurations． The numerical results show that the width of windward side of a building is a key factor to influence deposition velocity and deposition flux． Deposition velocities and deposition fluxes of aerosol particles are different on windward and leeward sides of the building． Deposition fluxes of aerosol particles on the building surfaces vary with particle diameters． Particle diameters corresponding to minimum deposition flux are also different for different building surfaces．

building configuration； building surface； inertia deposition； deposition velocity； deposition flux; aerosol particles

1671-0444(2017)01-0115-08

2016-03-23

國家自然科學(xué)基金資助項目(40975093)；上海市教委科研創(chuàng)新重點資助項目(14ZZ073)

華鳳皎(1989—)，女，山東煙臺人，博士研究生，研究方向為城市大氣環(huán)境與室內(nèi)空氣品質(zhì). E-mail: 2111173@mail.dhu.edu.cn 鐘 珂(聯(lián)系人)，女，教授，E-mail:zhongkeyx@dhu.edu.cn

X 513； P 426．6

A