周翌晨，付海明，楊　會

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

?

植物對街谷流場影響的數(shù)值模擬與評價

周翌晨，付海明，楊會

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

摘要：利用計算流體動力學(xué)(CFD)軟件對三維街谷模型進行數(shù)值模擬，采用附加源項法模擬植物的多孔介質(zhì)特性，對街谷中種植樹木前后的流場變化進行分析，并且從街谷中風(fēng)速分布變化、街區(qū)理想風(fēng)速比例及建筑周圍平均風(fēng)速3個方面進行評價. 結(jié)果表明：在兩種典型風(fēng)向下，街谷中的風(fēng)速在種植樹木后均有不同程度的降低，不合理極端風(fēng)速得到一定改善，可以認為植物對降低風(fēng)速、穩(wěn)定流場、改善行人區(qū)舒適性有一定的意義；同時從降低建筑物周圍風(fēng)速、 減少墻面對流換熱的角度而言，種植樹木對降低建筑能耗也有一定意義.

關(guān)鍵詞：數(shù)值模擬； 街谷流場； 植物； 風(fēng)速； 舒適性

街谷(street canyon)是街區(qū)的最基本單元，其由兩排高大建筑物及中間空曠車道構(gòu)成，由于形如峽谷，故名為街谷.分析街谷的流場特性，對于研究整個街區(qū)甚至是城市區(qū)域都有一定意義[1].綠化布置被認為是改善區(qū)域環(huán)境品質(zhì)的重要手段，植物的顏色和氣味可以大幅度提高人體對室外環(huán)境的滿意度和舒適感[2]，同時植物會對建筑物周圍的流場產(chǎn)生一定的影響，也可以捕集空氣中的塵粒.

計算流體動力學(xué)(CFD)技術(shù)已經(jīng)成為目前模擬科學(xué)研究的重要手段，其能全真模擬還原出現(xiàn)實狀態(tài)中區(qū)域的流場特點.本文采用CFD軟件分別對道路兩旁栽種樹木和無樹木的街谷進行數(shù)值計算，并從不同角度對街谷內(nèi)的流場進行分析評價.

對于街谷中樹木模型的建立，由于還原真實狀態(tài)十分困難，經(jīng)過了很多國內(nèi)外學(xué)者的探索，近年來較多采用的是附加源項法.所謂附加源項法，即在動量和湍流輸運方程中采用附加源/匯項的方法，模擬植物的繞流效應(yīng).文獻[3]通過附加源項法修改動量方程，并模擬了防風(fēng)柵對城市建筑局部風(fēng)環(huán)境的影響.文獻[4]利用附加源項法模擬了樹木對高層建筑局部風(fēng)環(huán)境的影響.文獻[5]采用附加源項法建立了用于城市污染物擴散分析的樹木模型，并對意大利巴里地區(qū)一處典型街區(qū)進行了案例分析.文獻[6]比較了RANS(Reynolds average Navier-Stokes)模型和LES(large eddy simulation)模型對于街谷研究的優(yōu)缺點.文獻[7]評估了LES模型應(yīng)用于含植物的城市街谷中的表現(xiàn)，其中的植物模型均是使用附加源項法表現(xiàn)的.本文基于附加源項法的植物模型，通過不同角度對街谷中存在植物的流場進行評價，為植物改善室外風(fēng)環(huán)境提供理論支持.

1數(shù)值模擬

1.1模型介紹

本文采用兩種模型來分析不同風(fēng)向?qū)止攘鲌龇植嫉挠绊?，如圖1所示. 建筑物長為100 m，寬為40 m，高為20 m，兩棟建筑物間隔為40 m；流域總長為660 m，寬為400 m，高為140 m.本文中植物模型是基于常見的街道樹——香樟，由于其連續(xù)布置于街道兩側(cè)，故可近似為長條形區(qū)域，將該區(qū)域簡化為多孔介質(zhì)區(qū)域，長為100 m，寬為4 m，高為4 m，離地高度為4 m，布置于街區(qū)內(nèi)中軸兩旁.建立坐標系，風(fēng)沿著x軸正方向吹，高度方向為z軸.

(a) 風(fēng)縱(豎)向吹過街谷

(b) 風(fēng)橫向吹過街谷

本文使用Ansys-ICEM軟件，用四面體網(wǎng)格劃分物理模型，對于建筑物及樹冠區(qū)域采用局部加密處理.使用Fluent軟件進行數(shù)值計算，湍流模型采用RNGκ-ε模型，壓力-速度耦合方程的解法采用SIMPLE算法.使用Tecplot后處理軟件進行模擬數(shù)據(jù)的采集，使用Origin軟件進行數(shù)據(jù)的統(tǒng)計和分析.

1.2模擬設(shè)置

入口邊界使用UDS(user-defined funtions)編寫的指數(shù)風(fēng)函數(shù)描述速度入口條件，見式(1)，并根據(jù)氣象參數(shù)及文獻[8]，取參考高度z0=10 m時，平均風(fēng)速u(z0)為3.94 m/s.

(1)

式中：u(z)為高度為z處的平均風(fēng)速；z0為參考高度；α為地面粗糙指數(shù)(地面越粗糙，對氣流的滯留作用越明顯)，取α=0.3.

出口邊界采用流量出口邊界條件(outflow)，兩側(cè)和頂部采用對稱邊界條件(symmetry)，地面和建筑墻面采用無滑移邊界條件，粗糙度分別取0.005和0.002 m[9].

1.3控制方程及源項

流域內(nèi)的空氣為不可壓縮理想氣體.由于樹冠區(qū)復(fù)雜的微觀結(jié)構(gòu)，若進行全真建模，對網(wǎng)格劃分以及模擬計算都是很大的考驗.因此，參照文獻[10-12]，本文的研究選擇用源項的形式描述速度及壓力的損失，表現(xiàn)樹冠的多孔介質(zhì)特性.以下給出涉及的控制方程及源項.

連續(xù)性方程為

(2)

式中：ui為速度；xi為空間坐標；i為方向角標，可取1,2,3.

動量方程為

(3)

(4)

湍動能方程為

(5)

湍流耗散率方程為

(6)

式中：p為空氣壓力；ρ為空氣密度；νe為有效黏度；ν為運動黏度；νt為湍流黏度；σk和σε為普朗特數(shù)；Cμ,Cε1,Cε2為無量綱常數(shù)(取值見表1)；S、Sk、Sε分別為動量源項、湍流動能源項和湍流耗散率源項，用于表現(xiàn)氣流流經(jīng)植物的樹葉樹枝后產(chǎn)生的變化，具體公式為

S=-ρCdA|u|ui

(7)

(8)

(9)

表1　各系數(shù)取值

1.4可靠性驗證

對于樹冠模型可靠性驗證，使用源項作為植物模型(樹高為h)，在與文獻[16]相同邊界條件下進行模擬，得到了樹后垂直方向風(fēng)廓線，如圖2所示.由圖2可以看出，使用本文源項進行數(shù)值模擬的結(jié)果趨勢與文獻[16]的研究結(jié)果基本一致，故本文中使用該源項表現(xiàn)植物對風(fēng)的阻力作用是可靠的.

圖2　模型驗證Fig.2　Model validation

2街谷流場描述

對于風(fēng)縱向吹過街谷的流場，如圖3(a)所示，吹向建筑物的氣流受到建筑物的阻擋，其中一部分向下運動，之后遇到地面阻擋后產(chǎn)生與來流方向相反的氣流，形成渦流，于是在建筑物迎風(fēng)面前形成氣流滯留區(qū)；另一部分向上運動，最終在建筑物頂部逃逸.在建筑物的前端拐角處同樣也發(fā)生了分離，這些分離使局部位置獲得了較大的風(fēng)速，尤其是在冬季，會對行人產(chǎn)生極不舒適的感受.同樣，由于街谷的存在，前端拐角的分離氣流向街谷擠壓，使街谷區(qū)域成為高流速區(qū)域，并且當(dāng)街谷越窄時，該現(xiàn)象越明顯.當(dāng)氣流吹過建筑物后，在建筑物背風(fēng)面產(chǎn)生較大范圍的負壓區(qū)，在拐角處氣流發(fā)生分離，形成較大的渦流.

(b) 橫向吹過街谷

對于橫向吹過街谷的流場，如圖3(b)所示，迎風(fēng)建筑與縱向情況中的建筑相似，前端有部分氣流滯留區(qū)，拐角處流速增大.街谷中的氣流較為雜亂，主要來源是迎風(fēng)建筑背風(fēng)拐角處的渦流和從街谷頂部產(chǎn)生的繞流，街谷中的流場情況主要與街谷的寬度有關(guān).后部建筑物背風(fēng)側(cè)有較大的渦流區(qū)，該區(qū)域風(fēng)速偏低，易形成氣流死區(qū)，在夏季會造成悶熱等不舒適感.

圖4(a)和4(b)分別為離地1.5 m高度風(fēng)縱向吹過街道兩旁不含樹木和種植樹木的街谷形成的流場，左端為速度入口方向. 由圖4(a)和4(b)可知，街谷中栽種樹木后，行人高度1.5 m處風(fēng)速普遍低于不含樹木的街谷，局部區(qū)域的風(fēng)速有較大改善.圖4(c)和4(d)分別為離地1.5 m高度風(fēng)橫向吹過街道兩旁不含樹木和種植樹木的街谷形成的流場，產(chǎn)生的現(xiàn)象和風(fēng)縱向吹過街谷類似.不同之處在于當(dāng)風(fēng)縱向吹過街谷時，樹木對街谷中風(fēng)速的削弱作用強于風(fēng)橫向吹過街谷.

(a) 無樹，縱向，z=1.5 m

(b) 有樹，縱向，z=1.5 m

(c) 無樹，橫向，z=1.5 m

(d) 有樹，橫向，z=1.5 m

3對街谷周圍流場的評價方法

(1) 沿街谷中軸速度分布及街谷中心垂直方向速度分布.取街谷中軸1.5 m高度水平方向，觀察種樹前后風(fēng)速分布變化；取街谷中心處垂直方向，觀察種植樹木前后風(fēng)速分布變化.該數(shù)據(jù)可以直接反映樹木對于來流強風(fēng)的阻擋效果.

(2) 流速特性分布指標. 與空氣分布特性指標(ADPI)類似，該指標是所有樣本點中理想風(fēng)速樣本點所占的比例.小于5 m/s的風(fēng)速被認為是舒適風(fēng)速[17]，而理想風(fēng)速應(yīng)比舒適風(fēng)速要求更高，由于夏季室外風(fēng)速過小會給人造成悶熱的感覺，故本文中的理想風(fēng)速取1～4 m/s，即小于1 m/s和大于4 m/s 的風(fēng)速被認為是不理想風(fēng)速.該指標可以反映種植樹木前后流場對行人舒適度的影響.

(3) 墻面平均風(fēng)速. 取距離墻面1 m位置的切面作為參考面，計算該面上所有樣本點風(fēng)速的平均值.該參數(shù)主要反映建筑物墻面上的風(fēng)速大小，這對于建筑能耗有很大意義，尤其是在冬季，建筑物周圍的風(fēng)速直接影響墻面與空氣的對流換熱，從而影響建筑物室內(nèi)的熱負荷.因此，降低建筑物周圍的平均風(fēng)速可以間接降低建筑能耗.

4模擬結(jié)果與評價分析

4.1風(fēng)速分布對比

4.1.1街谷中軸1.5 m高度水平方向風(fēng)速分布

街谷中種植樹木和不含樹木兩種情形下繪制水平方向風(fēng)速分布曲線如圖5所示.由圖5(a)可知，由于街口的構(gòu)造使風(fēng)速產(chǎn)生驟升，但隨后逐漸降低.不同的是，不種植樹木時，風(fēng)速隨后的趨勢是先減小后增大，結(jié)合圖4(a)可知，風(fēng)速增大的原因可能是拐角的高速氣流在街區(qū)后半段產(chǎn)生的影響.而種植樹木時，這些高速氣流由于樹木阻擋，風(fēng)速逐漸減小.將圖5(a)中兩曲線對比可知，種植樹木時風(fēng)速均小于不種植樹木，可見樹木對于氣流的阻擋作用明顯，整個街區(qū)的后部氣流平穩(wěn)舒適. 由圖5(b)可知，氣流在街口位置發(fā)生分離，在建筑物后的街區(qū)內(nèi)產(chǎn)生了渦流，由于兩邊都有街口，氣流進入街區(qū)后會發(fā)生碰撞和混合，加之湍流運動本身具有很大的波動性和不穩(wěn)定性，因此氣流橫向吹過街谷的情況下流場比較復(fù)雜和隨機，并沒有很好的對稱性，這在之前的流場描述中已有提到.但是對比種植和不種植樹木兩種情況，依然能夠看出在種植樹木的情況下，相同位置處風(fēng)速要低于不種植樹木的情況.

(a) z=1.5 m，y=200 m，縱向

(b) z=1.5 m，x=210 m，橫向

4.1.2街谷中心垂直方向風(fēng)速分布

街谷中種植樹木和不種樹木兩種情形下繪制垂直方向風(fēng)速分布曲線如圖6所示.從圖6(a)和6(b)可知，不管是縱向還是橫向，種植樹木情況下的風(fēng)速均小于不種植樹木的情況.不同的是，氣流縱向吹過街谷時，風(fēng)速的趨勢比較單一，基本是逐漸增加的.植物存在時增加速度較慢，氣流較為穩(wěn)定，可以明顯看出樹木對于氣流的阻擋作用.而對于氣流橫向吹過街谷的情況，風(fēng)速的變化隨著高度變化較為復(fù)雜，這與街谷內(nèi)存在較多渦流有關(guān)，并且在高度較大的區(qū)域，但植物的存在依然可以對風(fēng)速產(chǎn)生影響.

(a) x=210 m，y=200 m，縱向

(b) x=210 m，y=200 m，橫向

4.2流速特性分布指標

由于植物的存在不僅僅只作用于街區(qū)內(nèi)部，還對上下游建筑物周圍流場有影響.因此將統(tǒng)計區(qū)域擴大，選取街區(qū)及建筑物周圍10 m范圍作為統(tǒng)計區(qū)，水平高度依然是行人高度1.5 m(即z=1.5 m平面上街谷及離建筑物10 m內(nèi)的所有樣本點，要求選取的樣本點分布均勻).統(tǒng)計整個平面內(nèi)的所有樣本點，整理得到表2.由表2可知，當(dāng)氣流縱向吹過街谷且無種植樹木時，共有16 882個模擬樣本點被統(tǒng)計，平均風(fēng)速為2.82 m/s，理想風(fēng)速的模擬樣本點有14 560個，占總體的86.25%，即有86.25%的樣本點風(fēng)速較為理想.統(tǒng)計植物存在時的模擬結(jié)果，模擬樣本點還是16 882個，平均為2.38 m/s，理想風(fēng)速樣本點的個數(shù)有15 242，占了90.29%，比不種植樹木時提高了4.1%.由此可見，種植樹木有改善不合理風(fēng)速的效果.

表2　理想風(fēng)速統(tǒng)計

在氣流橫向吹過街谷時，效果則不太明顯.平均風(fēng)速在種植樹木后減少了0.21 m/s，理想風(fēng)速的比例只從原先的81.46%提高到了81.90%，改善的幅度很有限.這是由于風(fēng)并沒有直接吹過樹木，大部分流過樹木的氣流是建筑物兩側(cè)和頂部形成的渦流，因此植物對于氣流的改善作用就會小得多.

4.3墻面平均風(fēng)速對比

取距建筑物墻面(包括屋面)1 m位置作為參考面(即墻面向外平移1 m后得到的面)，分別編號1,

2,…，20，共20個面.統(tǒng)計參考面上所有點的風(fēng)速值，取該參考面風(fēng)速的平均值作為對應(yīng)墻面的平均風(fēng)速值.縱向和橫向模型中參考面編號與墻面的對應(yīng)關(guān)系分別如圖7所示.

(a) 縱向　(b) 橫向

縱向情況墻面平均風(fēng)速值數(shù)據(jù)如表3所示. 由表3可知，當(dāng)氣流縱向吹過街谷時，種植樹木后，墻面3和6上的平均風(fēng)速均有0.5 m/s左右的降低，其余墻面沒有較大影響.同樣地，橫向情況墻面平均風(fēng)速值數(shù)據(jù)如表4所示. 由表4可知，當(dāng)風(fēng)橫向吹過街谷時，在墻面14、17、20處種樹木后風(fēng)速均有0.4 m/s左右的降低，其余墻面速度變化可以忽略不計.

表3　縱向情況墻面平均風(fēng)速值

表4　橫向情況墻面平均風(fēng)速值

結(jié)合表3和4以及圖7可以看出，平均風(fēng)速有所降低的墻面都是與植物最接近的墻面，墻面20是由于受到從街谷中向上逃逸氣流的影響，因此也有所降低.因此，可以認為若在每個墻面附近均種植樹木，墻面的平均風(fēng)速均會降低，而墻體表面風(fēng)速的降低對于冬季降低建筑物墻體的對流換熱有很大的意義，一定程度上可以降低建筑能耗.

5結(jié)語

本文提出了3個針對大空間流場中行人區(qū)舒適性的評價角度，并運用CFD數(shù)值模擬技術(shù)進行分析評價，得到以下結(jié)論：

(1) 不管是氣流縱向吹過街谷還是橫向吹過街谷，從水平及垂直速度分布情況均可以看出，植物對氣流有一定阻擋作用，可以在一定程度上降低風(fēng)速.縱向情況下樹木對氣流的阻擋作用更明顯，對街谷中的環(huán)境改善也較為明顯，并且這一改善作用會持續(xù)到街谷尾部之后的區(qū)域.

(2) 從理想風(fēng)速比例角度，縱向情況下植物對于建筑周圍氣流環(huán)境的改善作用更為明顯，理想風(fēng)速比例由86.25%提高到90.29%，橫向情況下提高的幅度相對較弱，這與街谷內(nèi)部復(fù)雜的氣流狀況有關(guān).

(3) 值得注意的是，對于建筑物墻面上的平均風(fēng)速，不管縱向還是橫向情況下，種植樹木后均有所降低，最高降幅可達20%左右，這對于冬季降低建筑能耗有一定意義.當(dāng)然，根據(jù)種植樹木的種類及種植密度的不同，降低的程度和改善的效果會有所不同.

綜上所述，不管風(fēng)向如何，在街谷中種植樹木對于改善行人區(qū)舒適性有一定的意義.本文旨在證明植物對于街谷流場的影響，也許某些風(fēng)向情況下效果并不明顯，但是依然可以通過增加綠化率或是種植葉面積密度較大的植物來加強植物對流場的影響.

需要指出的是，雖然本文研究的是樹木改變流場以達到改善人體舒適性的目的，但是樹木對人體舒適性的改善主要來源還是其對太陽輻射的吸收和對熱環(huán)境的改善.另一方面，本文中的街谷模型只是眾多街谷樣式中的一種，實際情況下，由于街谷兩旁的建筑物高度、兩邊的對稱程度、上游建筑的布局以及樹木的布置位置等影響，都會形成更加復(fù)雜的流場.另外，本文中的部分情況有一定的特殊性，例如樹木四季的變化引起葉面積密度的變化，冬季時可能并沒有那么理想的效果，更深入的分析還有待之后的課題研究.

參考文獻

[1] VOS P E J,MAIHEU B, VANKERKOM J, et al. Improving local air quality in cities: To tree or not to tree?[J].Environmental Pollution, 2012, 183:113-122.

[2] 冷寒冰,胡永紅,周鑫，等.室外植物對人體舒適度及環(huán)境滿意率的影響[J].中南大學(xué)學(xué)報(自然科學(xué)版),2012,43(11):4574-4580.

[3] LI W, WANG F, BELL S. Simulating the sheltering effects of windbreaks in urban outdoor open space[J].Journal of Wind Engineering and Industrial Aerodynamics,2007,95(7):533-549.

[4] 楊易，顧明，金新陽.風(fēng)環(huán)境數(shù)值模擬中模擬植被的數(shù)學(xué)模型與應(yīng)用[J].同濟大學(xué)學(xué)報(自然科學(xué)版),2010,38(9): 1266-1270.

[5] BUCCOLIERI R， SALIM S M， LEO L S. Analysis of local scale tree-atmosphere interaction on pollutant concentration in idealized street canyons and application to a real urban junction[J]. Atmospheric Environment, 2011,45(9):1702-1713.

[6] SALIM S M, CHEAH S C, CHAN A. Numerical simulation of dispersion in urban street canyons with avenue-like tree plantings:Comparison between RANS and LES[J].Building and Environment, 2011, 46:1735-1746.

[7] MOONEN P, GROMKE C. Performance assessment of large eddy simulation (LES) for modeling dispersion in an urban street canyon with tree planting[J]. Atmospheric Environment,2013,75:66-76.

[8] SALIMA S M, BUCCOLIERI R, CHANA A, et al.Large eddy simulation of the aerodynamic effects of trees on pollutant concentrations in street canyons[J]. Procedia Environmental Sciences, 2011(4):17-24.

[9] NG W Y, CHAU C.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.

[10] ENDALEW A M, HERTOG M, DELELE M A. CFD modeling and wind tunnel validation of airflow through plant canopies using 3D canopy architecture[J]. International Journal of Heat and Fluid Flow, 2009, 30(2):356-368.

[11] ENDALEW A M, HERTOG M, GEBREHIWOT M G, et al. Modelling airflow within model plant canopies using an integrated approach[J].Computers and Electronics in Agriculture, 2009, 66(1):9-24.

[12] AMORIM J H, RODRIGUES V, TAVARES R, et al. CFD modelling of the aerodynamic effect of trees on urban air pollution dispersion[J].Science of the Total Environment, 2013, 461/462:541-551

[13] 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.

[14] SANZ C. A note onκ-εmodelling of vegetation canopy air-flows[J]. Boundary Layer Meteorol, 2003, 108(1): 191-197.

[15] KATUL G G, MAHRT L, POGGI D, et al.One-and two-equation models for canopy turbulence[J]. Boundary-Layer Meteorol, 2004, 113(1):81-109.

[16] ZENG P T, HIDENORI T. A first-order closure model for the wind flow within and above vegetation canopies[J].Agricultural and Forest Meteorology, 2000, 103(3): 301-313.

[17] 綠色建筑評價標準:GB/T 50378—2014[S].

文章編號：1671-0444(2016)03-0419-07

收稿日期：2015-05-25

基金項目：國家自然科學(xué)基金資助項目(41371445)

作者簡介：周翌晨(1990—)，男，浙江紹興人，碩士研究生，研究方向為通風(fēng)空調(diào)與氣體凈化.E-mail：zyc427@126.com 付海明(聯(lián)系人)，男，高級工程師，E-mail：fhm@dhu.edu.cn

中圖分類號：S 731.2

文獻標志碼：A

Numerical Simulation and Evaluation of Flow Field of the Street Canyon with Plants

ZHOUYi-chen,FUHai-ming,YANGHui

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

Abstract:Based on the CFD(computational fluid dynamics)method, using additional source term,3D numerical calculation was used to simulate air flow in street canyon to study the influence of plants in the street canyon.The flow field of street with plants was analyzed and evaluated on three aspects, including wind speed in street, the ideal wind velocity distribution performance index and the average wind speed around the buildings.The results indicate that street with plants can reduce wind speed and improve pedestrian comfort in different wind directions cases. And from the perspective of reducing the wind speed around the buildings and convective heat transfer of the wall, planting trees is significant for reducing building energy consumption in the winter.

Key words:numerical simulation; flow field of street canyon; plants; wind velocity; comfort