喬文佳賀啟濱齊亞茹劉艷華

1西安交通大學人居學院

2深圳市建筑科學研究院有限公司

交通源顆粒物對小區(qū)內(nèi)環(huán)境的影響

喬文佳1賀啟濱2齊亞茹1劉艷華1

1西安交通大學人居學院

2深圳市建筑科學研究院有限公司

利用湍流k-ε模型及離散相模型（DPM）研究了小區(qū)外交通源（線源）排放顆粒物對小區(qū)內(nèi)環(huán)境的影響。對比分析了有、無圍墻情況下，小區(qū)內(nèi)環(huán)境受交通源污染物的影響狀況。結(jié)果表明：圍墻對顆粒物有很好的阻擋作用。因圍墻的存在，顆粒物大部分積聚在小區(qū)外、圍墻附近。雖然部分顆粒物能夠進入小區(qū)，由于其濃度很低，小區(qū)環(huán)境基本不受影響。但是，無圍墻存在時，顆粒物無阻擋地進入小區(qū)，小區(qū)受交通源顆粒物的影響很大，特別是前排建筑物的迎風面和繞流區(qū)域。

交通源顆粒物 小區(qū)內(nèi)環(huán)境 圍墻 離散相模型

0 引言

機動車保有量的持續(xù)增加，不僅帶來了城市中心區(qū)域交通擁堵等嚴重問題，與此同時也帶來了交通源污染的持續(xù)攀升[1～2]，對城市環(huán)境造成嚴重損害。交通源排放的污染物不僅包括氣態(tài)污染物（CO、SO2、NOx等），也包括固態(tài)顆粒物。研究表明，顆粒物可引起呼吸系統(tǒng)疾病、心腦血管疾病等多種疾病。Dockery等在一項對哈佛六城市長期顆粒物污染的研究中，指出：居住在濃度大的城市增加30%心血管疾病死亡的危險[3]；Dominici等分析了204個美國城市的醫(yī)療數(shù)據(jù)，指出PM2.5的濃度每增加10ug/m3，缺血性心臟疾病入院人數(shù)增加0.44%，腦血管疾病0.81%，外周動脈疾病0.86%，心律失常0.57%，心臟衰竭1.28%[4]。

近年來，針對交通源氣態(tài)污染物對城市環(huán)境的影響研究已經(jīng)較為深入，這其中包括現(xiàn)場實測[5～6]、風洞實驗[7]、數(shù)值模擬[8～9]。針對交通源固態(tài)顆粒物污染的研究還較少，并且相關(guān)研究主要集中于街道峽谷中顆粒物的擴散[10～11]，但交通源顆粒物對臨街小區(qū)環(huán)境是否同樣有影響，對此問題卻很少有人研究。故本文建立小區(qū)物理模型研究在有無圍墻兩種情況下，道路交通源顆粒物對小區(qū)環(huán)境的影響問題。

1 物理模型

本文的物理模型如圖1。該小區(qū)由四棟建筑物組成，各建筑物尺寸相同。建筑物長、高、寬皆為10m；各建筑左右、前后間距都是10m；當小區(qū)有圍墻時，圍墻高度為2m，圍墻距建筑物10m。距離圍墻5m有道路。假設(shè)道路上車輛污染源為線污染源，位于地面。污染源的長度為50m，寬度為1m。當小區(qū)無圍墻時，除了去除圍墻外，其參數(shù)相同。計算區(qū)域選擇為：入口距離圍墻150m，出口距離圍墻250m，左右側(cè)距離圍墻皆為100m，高度為50m。采用非結(jié)構(gòu)網(wǎng)格，并在壁面及其附近區(qū)域進行網(wǎng)格加密。加密區(qū)域為距離圍墻外各10m、高度20m的立方體。圖2為結(jié)果分析剖面所在位置的示意圖。

圖1 有圍墻小區(qū)結(jié)構(gòu)圖

圖2 剖面示意圖

2 控制方程與邊界條件

2.1 控制方程

本文采用標準k-ε模型模擬風場運動。

連續(xù)方程：

動量方程：

k和ε輸運方程：

式中：vt=Cμk2/ε，Cμ=0.9，Cε1=1.41～1.45，Cε2=1.91～1.92，Ck=0.09～0.11，Cε=0.07～0.09。

本文應(yīng)用拉格朗日方法跟蹤顆粒物的運動。對于顆粒相本文做如下假設(shè)：①假定所有顆粒相都為球形；②由于顆粒相所占體積比率十分小，所以不考慮顆粒的碰撞。通過對顆粒所受的各種外力平衡方程可以得到預(yù)測顆粒的運動軌跡運動方程如下：

式中：FD(u-up)為顆粒的單位質(zhì)量曳力，表達式為

式中：u為流體項速度；up為顆粒速度；μ為流體動力粘性系數(shù)；ρ為流體密度；ρp為顆粒密度（骨架密度）；dp為顆粒直徑；Re為相對雷諾數(shù)。

式中：Fx為顆粒所受的其它力，包含布朗力FB和薩夫曼力FS，即Fx=FB+FS。

布朗力表征的是微小顆粒熱運動所受的力，粒徑越小，布朗力越大，其表達式如下：

S0是強度量，其表達式如下：

薩夫曼剪切提升力是由于流場中的速度梯度引起，其表達式為

式中：ui為流體項速度；upi為顆粒速度；v為流體的運動粘性系數(shù)；ρ為流體密度；ρp為顆粒密度（骨架密度）；dp為顆粒直徑；dij為變形速率張量；kc=2.594為一常數(shù)[12]。

流體渦團特征生存時間表達式如下：

式中：經(jīng)驗常數(shù)CL在k-ε模型中取為0.15；r為服從（0，1）區(qū)間均勻分布的隨機數(shù)。

顆粒穿過流體渦團時間表達式如下：

式中：τ為顆粒松弛時間；Le為渦團特征長度標尺；|u-up|為顆粒與流體速度差。

隨機軌道模型假設(shè)脈動速度服從高斯概率分布。脈動速度表達式如下：

式中：δ為服從正態(tài)分布的隨機數(shù)。

對于本文選取的標準k-ε模型，脈動速度表達式如下：

經(jīng)過流體與顆粒相互作用的時間間隔，顆粒終止在當前渦內(nèi)的計算，并在下一個時間間隔，顆粒進入一個新的渦團，重復(fù)以上步驟獲得顆粒新的速度、軌跡。當顆粒碰到壁面被捕捉或離開計算區(qū)域時，終止軌跡的計算。

2.2 邊界條件

本文計算有無圍墻兩種模型。入口風速風剖面表示為：

式中：zb和uzb分別為參考高度和參考高度上的風速；α為地面粗糙度指數(shù)。這里取zb=10，uzb=2.0，α=0.24。

污染源排放的顆粒物顆粒密度取2000kg/m3，粒徑取1μm，流量取0.003kg/s。

1）入口邊界條件：速度入口；

2）出口邊界條件：質(zhì)量出口邊界條件；

3）壁面邊界條件：地面和墻面無滑移邊界條件；計算域的上邊界、左右邊界取為對稱邊界，相當于滑移邊界條件。

4）顆粒相邊界條件：壁面為reflect邊界條件；出口、入口為escape邊界條件；

5）對控制方程離散格式：二階迎風格式，simple算法非耦合計算。

3 結(jié)果分析

3.1 流場分析

圖3分別是有、無圍墻小區(qū)y=15m剖面速度矢量圖。從中可看出無論有無圍墻都會在建筑物中間區(qū)域及其后排建筑物尾部出現(xiàn)渦旋，并且位置基本相同。不同之處在于對于有圍墻小區(qū)，因圍墻對氣流的阻擋及壓差梯度作用在圖3中位置1、2處會出現(xiàn)渦旋；對于無圍墻小區(qū)因前排建筑物阻擋及壓差作用會在圖3中位置3處出現(xiàn)渦旋。

圖3 y=15m剖面速度矢量圖

圖4分別是有、無圍墻小區(qū)z=1.5m剖面速度矢量圖。進入小區(qū)的氣流在遇到建筑物后分成上下兩股繞流在兩棟建筑物中間和后排建筑物尾部形成渦旋。但由于圍墻的阻擋作用，導致進入小區(qū)的氣流在z=1.5m剖面速度矢量不同，進而影響建筑物中間區(qū)域和后排建筑物尾部的渦旋位置分布。

圖4 z=1.5m剖面速度矢量圖

圖5分別是有、無圍墻小區(qū)z=5m剖面速度矢量圖。來流遇到前排建筑物分成兩股繞流，因壓差作用在建筑物中間區(qū)域和后排建筑物尾部形成渦旋。但由于圍墻的存在影響了小區(qū)流場，渦旋的位置有所變化。

圖5 小區(qū)z=5m剖面速度矢量圖

3.2 濃度分析

為便于分析濃度的相對值大小，圖6～圖8皆采用統(tǒng)一比色卡。圖6分別是有、無圍墻小區(qū)y=15m剖面濃度分布圖。對于有圍墻小區(qū)，圍墻對于顆粒物有很好的阻擋作用，部分顆粒物被阻擋在圍墻外，導致圍墻外邊緣附近濃度極大；只有部分顆粒物爬墻進入小區(qū)，且因為前排建筑物阻擋，使小區(qū)內(nèi)各點的濃度值皆很小。對于無圍墻小區(qū)，顆粒物無阻礙地進入小區(qū)，遇到前排建筑物受到阻擋，部分顆粒物積聚在前排建筑物迎風面附近，導致前排建筑物迎風面附近區(qū)域濃度很大，部分隨流體繞流或是爬墻繼續(xù)進入小區(qū)內(nèi)部，使兩排建筑物中間和后排建筑物尾部區(qū)域濃度偏大。

圖6 小區(qū)y=15m剖面濃度分布

圖7分別是有、無圍墻小區(qū)z=1.5m剖面濃度分布圖。對于有圍墻小區(qū)，因為圍墻的阻擋作用，進入小區(qū)的顆粒物只是部分，所以顆粒物對小區(qū)的環(huán)境影響很小。對于無圍墻小區(qū)，顆粒物無阻擋地進入小區(qū)，遇到前排建筑物隨氣流以繞流方式通過，繞流和壓差影響使前排建筑物迎風面濃度很大，繞流區(qū)域濃度很大。后排建筑物尾部和建筑區(qū)域后部也因同樣原因呈現(xiàn)較高濃度值。

圖7 小區(qū)z=1.5m剖面濃度分布

圖8分別是有、無圍墻小區(qū)z=5m剖面濃度分布圖。對于有圍墻小區(qū)，z=5m剖面濃度值極小。對于無圍墻小區(qū)只有部分區(qū)域因為壓差作用等出現(xiàn)較高濃度值。從中可看出，隨著高度的增加，交通源顆粒物對小區(qū)的影響逐漸減小。

圖8 小區(qū)z=5m剖面濃度分布

4 結(jié)論

利用湍流k-ε模型及離散相模型（DPM）研究了小區(qū)外交通源（線源）排放顆粒物對小區(qū)內(nèi)環(huán)境的影響。分析了小區(qū)外交通源在小區(qū)有、無圍墻兩種情況下對小區(qū)的影響情況，通過模擬計算結(jié)果得到如下結(jié)論：

1）圍墻的存在時，對顆粒物擴散有阻擋作用。部分顆粒物爬過圍墻進入小區(qū)，均勻地向小區(qū)擴散，顆粒物分布范圍大，影響區(qū)域大。但所影響區(qū)域濃度相對值不會很大。

2）無圍墻存在時，顆粒物因無阻擋全部進入小區(qū)。遇到前排建筑物受阻擋，所以前排建筑物受影響較大。部分顆粒物受阻擋后跟隨流體發(fā)生繞流，繞流區(qū)域在建筑物附近，故造成小區(qū)左右通道中靠近建筑區(qū)域濃度較大。顆粒物跟隨流體通過建筑區(qū)域后均勻擴散，在建筑物后側(cè)因有尾渦存在，濃度也會偏大。

3）無論是否有圍墻，隨著高度增加，顆粒物濃度都將減小。

[1]殷杉,蔡靜萍,陳麗萍.交通綠化帶植物配置對空氣顆粒物的凈化效益[J].生態(tài)學報,2007,27(11):4590-4596

[2]杜剛,張見昕.遼寧省典型城市大氣中顆粒物污染狀況及分布特征研究[J].現(xiàn)代科學儀器,2006,(2):69-70

[3]Dockery DW,Pope III CA,Xu X,et al.An association between air pollution and mortality in six U.S.cities[J].N Engl J Med, 1993,329:1753-1759

[4]Dominici F,Peng RD,Bell ML,et al.Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases[J].JAMA,2006,295:1127-1134

[5]Xie S D,Zhang Y H,Li Q,et a1.Spatial distribution of traffic-rel -ated pollutant concentrations in street canyons[J].Atmospheric Environment,2003,37(23):3213-3214

[6]Niachou K,Livada I,Santamouris M.Experimental study of tem -perature and airflow distribution Inside an urban street canyon during hot summer weather conditions-Part II:Airflow analysis [J].Building and Environment,2008,43(8):1393-1403

[7]ChristofG,Ruck B.Influence of trees Oil the dispersion of pollutants in an urban street canyon-experimental investigation of the flow and concentration field[J].Atmospheric Environment,2007, 41(16):3287-3302

[8]Huang Y D,Hu X N,Zeng N B.Impact of wedge-shaped roofs on airflow and pollutant dispersion inside urban street canyons[J]. Building and Environment,2009,44(12):2335-2347

[9]黃遠東,吳文權(quán),范炳全,等.不同建筑物寬高比的街道峽谷內(nèi)部氣流場數(shù)值模擬[J].上海理工大學學報,2005,27(3):203-206

[10]羅昔聯(lián),顧兆林.基于DPM模擬的街道峽谷顆粒物擴散特性研究[J].中國科學院研究生院學報,2007,24(5):578-583

[11]付志民,孫在.街道峽谷內(nèi)顆粒物擴散特性的DPM數(shù)值模擬[J].中國計量學院學報,2011,22(4):322-326

[12]李先庭,趙彬.室內(nèi)空氣流動數(shù)值模擬[M].北京:機械工業(yè)出版社,2009

Im pa c t of Tra ffic Exha us ting Pa rtic le s on Com m unity Environm e nt

QIAO Wen-jia1,HE Qi-bin2,QI Ya-ru1,LIU Yan-hua1
1 School of Human Settlements and Civil Engineering,Xi'an Jiaotong University
2 Shenzhen Institute of Building Research Co.Ltd.

In this paper,the impact of traffic exhausting particles on community environment was simulated using the standard-turbulence model and Dispersed Phased Model(DPM).The community near street which with and without fences are simulated.The results show that fences have a barrier function to particles:most of the particles accumulate nearby the area of fences outside the community for fences exist.Minority particles enter into the community,while the impact on the community environment is little.The particles will enter into the community completely for no fences exist.The community environment will be greatly impact by the traffic exhausting particles,especially the front buildings windward and bilateral area.

traffic exhausting particles,community environment,fences,DPM

1003-0344（2015）02-037-5

2013-10-12

喬文佳（1988～），男，碩士研究生；西安市咸寧西路28號西安交通大學興慶校區(qū)人居學院（710049）；E-mail:798656402@qq.com

"十二五"國家科技支撐計劃（2012BAJ09B01）；"十二五"國家科技支撐計劃（2012BAJ06B03）