防風抑塵網(wǎng)研究進展

孫昌峰，陳光輝，范軍領，李建隆

（青島科技大學化工學院，山東 青島 266042）

進展與述評

防風抑塵網(wǎng)研究進展

孫昌峰，陳光輝，范軍領，李建隆

（青島科技大學化工學院，山東 青島 266042）

防風抑塵網(wǎng)作為一種能有效控制散堆料場起塵與擴散的手段，近年來得以推廣，并得到學術(shù)界的關注。本文綜述了國內(nèi)外防風抑塵網(wǎng)的研究現(xiàn)狀和進展。著重從作用機理、防風效應的影響因素、數(shù)值模擬進展及防風網(wǎng)的應用等方面對防風網(wǎng)進行了總結(jié)評述；同時介紹了一種新型防風網(wǎng)——導流型防風網(wǎng)，即在常規(guī)平板型防風網(wǎng)網(wǎng)孔上增設導流翅片，使來流風繞過料堆產(chǎn)生上揚以減小其對料堆的直接作用力，并借助數(shù)值模擬手段對新型防風網(wǎng)的抑塵性能進行了研究；最后對防風網(wǎng)目前存在的不足和未來的研究方向進行了探討。

防風抑塵網(wǎng)；防風抑塵效應；數(shù)值模擬；導流翅片

Abstract：Porous fence，which is an effective method to abate dust emission in a bulk material storage yard，has been used in many practical applications and aroused extensive attention in academic fields. In this paper，the current status and progress of porous fence are summarized，including mechanisms，progress of numerical simulation and several influence parameters on the shelter effect. Simultaneously，a porous fence with deflector，which plays a guiding role and directs the air flow up to a certain angle in order to reduce the direct impact on the windward surface of the bulk material pile is introduced，the shelter effect of the new porous fence is investigated numerically. The present imperfection and direction of future research are also discussed.

Key words：porous fence；shelter effect；numerical simulation；deflector

防風網(wǎng)古稱防風柵、防風障[1]，是一種多孔的障礙物，最初是人類用來抵御風沙侵襲的有效手段。廣義的防風網(wǎng)是可以減小風速的任何結(jié)構(gòu)，如木柵、金屬網(wǎng)、樹籬等[2]。而工程上的防風網(wǎng)是指減小風強技術(shù)的工程裝置，主要由具有一定開孔率的金屬網(wǎng)板或者由編制、粘接、擠壓成型的非金屬網(wǎng)片、支撐鋼結(jié)構(gòu)、地下混凝土基礎和相應的輔助噴水裝置以及自控儀表系統(tǒng)組成。

防風網(wǎng)按其移動性能可以分為移動式和固定式[3]。移動式防風網(wǎng)采用電動升降，在使用時將防風網(wǎng)上升到一定的高度，不使用時將防風網(wǎng)降低，不影響其它作業(yè)，其主要應用于移動性較大的現(xiàn)場作業(yè)，如港口堆場煤炭加工等行業(yè)。而固定式防風網(wǎng)一般被加工成型后裝配到框架上，然后連接到固定在地面的鋼支架上。考慮到造價、操作性等因素，目前防風網(wǎng)以固定式為主。按網(wǎng)板形式的不同防風網(wǎng)可分為蝶形、直板形、半圓形等[4]，其中蝶形和直板形較為常見。而根據(jù)材質(zhì)的不同，防風網(wǎng)又可分為鍍鋁鋅網(wǎng)、玻璃鋼網(wǎng)、柔性纖維網(wǎng)等。其中鍍鋁鋅網(wǎng)因具有耐腐蝕、耐濕熱等優(yōu)點，主要應用于港口或近海堆料場。而玻璃鋼網(wǎng)和柔性纖維網(wǎng)則主要用于防風網(wǎng)設計使用年限較短的堆料場。

防風網(wǎng)的基本原理是通過設置多孔透風的屏障來降低來流風速并減弱其湍流強度，達到防風抑塵的目的[5-6]。因具有結(jié)構(gòu)簡單、維護方便、營建迅速等特點，已在農(nóng)、林、工礦等行業(yè)及許多大型港口得到廣泛應用[7-10]。本文作者對近年來國內(nèi)外在防風網(wǎng)方面的研究進行了綜述，并對防風網(wǎng)未來的研究方向進行了展望。

1 抑塵機理

防風網(wǎng)是利用空氣動力學的原理，將網(wǎng)前大尺度、高強度的旋渦梳理成小尺度、弱強度旋渦，從而改變網(wǎng)后的微環(huán)境[11]以達到減弱下游風速[12]和流場湍流度的效果，防止粉塵的飛揚。

防風網(wǎng)的抑塵機理與網(wǎng)后顆粒的起塵機理密切相關。Dong等[13]依據(jù)對網(wǎng)后不同距離處顆粒的起動風速的觀測評價了防風網(wǎng)的防風抑塵性能，指出防風網(wǎng)增大了網(wǎng)后顆粒的臨界起動風速，因而能達到抑塵的目的。Lee等[14]研究發(fā)現(xiàn)當來流風達到顆粒的臨界起動風速時，料堆表面的顆粒開始跳躍，且隨著開孔率的降低，網(wǎng)后顆粒的臨界起動速度逐漸增大。叢曉春等測定了不同礦料顆粒的起動風速并考察了揚塵量隨風速的變化關系，結(jié)果表明當來流風速低于顆粒臨界起動風速時，顆粒不起塵；當風速超過臨界起動風速時，起塵量隨著風速的增大呈4次方增大。

目前，研究者認為防風網(wǎng)減小了網(wǎng)后風速，增大了網(wǎng)后顆粒的臨界起動風速從而達到抑塵的目的。林官明等[15]則應用子波分析對防風網(wǎng)后的湍流信號進行了研究，發(fā)現(xiàn)網(wǎng)后旋渦的能量和發(fā)生頻率有顯著地降低，證明了防風網(wǎng)對起塵的抑制除了表現(xiàn)在降低風速以外還表現(xiàn)在對猝發(fā)旋渦的抑制，從而從另一個角度闡述了防風網(wǎng)的抑塵機理。

2 防風抑塵效應的影響因素

影響防風網(wǎng)防風抑塵作用的因素很多，其中防風網(wǎng)的結(jié)構(gòu)對其防風抑塵效應及網(wǎng)后庇護區(qū)的大小起決定作用，了解這些結(jié)構(gòu)因素可以更深入的認識防風網(wǎng)，并指導結(jié)構(gòu)的設計優(yōu)化，提高防風網(wǎng)性能。

2.1 開孔率的影響

開孔率是防風網(wǎng)的開孔透風面積與總面積之比，作為防風網(wǎng)的結(jié)構(gòu)指標[16]，是影響防風網(wǎng)防風抑塵性能最重要的因素[12，17-20]。

Raine等[21]考察了不同開孔率的防風網(wǎng)后平均風速和湍流度的變化。研究發(fā)現(xiàn)，開孔率為20%時，防風網(wǎng)能最有效地減小其背風面的平均風速。與未開孔的防風網(wǎng)相比，開孔的防風網(wǎng)能更好的減小平均風速，而較高開孔率的防風網(wǎng)則能夠提供更好的總體庇護效應。Li Wei等[22]研究發(fā)現(xiàn)開孔率在23%～30%時，近網(wǎng)處具有較小的風速。Lee和Park[23]討論了防風網(wǎng)后料堆表面的壓力變化，發(fā)現(xiàn)開孔率在 40%～50%時，防風網(wǎng)能較大程度的減小料堆表面的壓力，對料堆起塵的抑制作用最強。

防風網(wǎng)的開孔率也是決定網(wǎng)后回流區(qū)變化的重要參數(shù)[24]。Castro[25]和Ranga等[26]均發(fā)現(xiàn)當防風網(wǎng)開孔率大于30%時，網(wǎng)后回流即會消失，而Lee和Kim等[27]則報道了當開孔率大于40%時，由于通過網(wǎng)的強滲流風的影響，網(wǎng)后平均風速增大，回流消失；而對于開孔率為40%的防風網(wǎng)，網(wǎng)后流場具有較好的流動特性，網(wǎng)后流場的湍流強度也較小。

此外，Lee等[28-29]還通過風洞試驗研究了防風網(wǎng)對料堆的庇護效應，研究表明具有最佳庇護效應的防風網(wǎng)其開孔率在30%～40%。而Mercer[30]則報道了防風網(wǎng)的最佳開孔率在25%左右。

綜上所述，開孔率對防風網(wǎng)的防風抑塵性能起著至關重要的作用，但目前國內(nèi)外對開孔率的研究結(jié)果并不一致，這是由于試驗所處的大氣環(huán)境、模擬及實驗采用的模型不同所致。最佳的開孔率一般在20%～50%。

2.2 網(wǎng)高的影響

防風網(wǎng)的高度與庇護范圍密切相關。Torano等[31]考察了露天儲料場顆粒的起塵情況，發(fā)現(xiàn)當網(wǎng)高小于料堆高度時，在防風網(wǎng)至網(wǎng)后2倍料堆高度間的顆粒起塵量最小；而當網(wǎng)高為1倍和1.2倍料堆高度時，在網(wǎng)至網(wǎng)后 3倍料堆高度間的顆粒起塵量最小。陳凱華等[32]則對某鋼鐵廠露天堆料場防風網(wǎng)防風效果進行了數(shù)值模擬，結(jié)果顯示，在風速一定的情況下，持續(xù)增大網(wǎng)的高度并不能達到持續(xù)擴大庇護范圍的效果，網(wǎng)高與受保護料堆的高低和料場的面積有直接的關系，網(wǎng)高取值為受保護料堆高度的 1.5倍較為適宜。Dong等[13]研究發(fā)現(xiàn)防風網(wǎng)的庇護范圍隨網(wǎng)高的增加而增大，當達到一個峰值后，網(wǎng)高再增加，防風網(wǎng)的庇護范圍變化并不明顯。

有研究表明[1]：當防風網(wǎng)的高度為料堆高度的0.6～1.1倍時，網(wǎng)高與抑塵效果成正比；當防風網(wǎng)高度為料堆高度的1.1～1.5倍時，網(wǎng)高對抑塵效果的影響趨于平緩；當防風網(wǎng)高度為料堆高度的 1.5倍以上時，隨著網(wǎng)高的增加抑塵效果無顯著的變化，因此防風網(wǎng)最佳高度為料堆高度的1.1～1.5倍。

2.3 網(wǎng)與料堆距離的影響

與開孔率和網(wǎng)的高度相比，網(wǎng)與料堆距離（網(wǎng)至料堆前堆腳的距離）的影響并不十分顯著。Lee等[23]研究發(fā)現(xiàn)，在開孔率和網(wǎng)高一定的情況下，改變防風網(wǎng)與料堆之間的距離對料堆表面的平均壓力并無大的影響。Li Wei等[22]則報道了當間距大于4倍的網(wǎng)高時，開孔率為50%的防風網(wǎng)的庇護效應接近于未開孔的防風網(wǎng)的庇護效應。而有研究表明[1]：當防風網(wǎng)沿料堆腳放置時，由于網(wǎng)后漩渦的影響，極易導致起塵，同時在堆腳和堆頂處滲流風速的降低并不顯著，而在網(wǎng)后2～3倍料堆高度的距離內(nèi)有一個低風區(qū)，減速效果較好。日本的研究[1]表明防風網(wǎng)與最近料堆的距離可控制在1.0～1.5倍料堆高度之內(nèi)。

此外，Kim 和Lee[33]還考察了防風網(wǎng)網(wǎng)孔大小對網(wǎng)后流場的影響，發(fā)現(xiàn)在同一開孔率下，隨孔徑的減小，防風網(wǎng)對來流風的阻礙作用不斷增大，滲流風的湍流度也不斷增加。Yeh等[34]研究了不同的來流風方向?qū)Ψ里L網(wǎng)性能的影響，結(jié)果表明，防風網(wǎng)的抑塵作用很大程度上取決于來流風的方向，傳統(tǒng)防風網(wǎng)布置方式為在料堆四周呈長方形布置，此時，在來流風向與迎風面呈45°夾角時對網(wǎng)后料堆的抑塵作用較差；而呈八邊形布置時，對于傾斜角度來流風，防風網(wǎng)的抑塵作用有較大改善。

3 防風抑塵網(wǎng)的工程應用

在防風網(wǎng)的工程應用方面，近年來日本、韓國等所做的研究較多，日本從20世紀70年代起，就開始對防風網(wǎng)進行相關研究，并將防風網(wǎng)技術(shù)應用于控制港口露天煤堆場的粉塵污染上[1]。韓國浦項科技大學的Lee和Park[35]對POSCO鋼廠5號煤堆場兩側(cè)設置的防風網(wǎng)的抑塵效果進行了現(xiàn)場測試，結(jié)果發(fā)現(xiàn)網(wǎng)后流場湍流強度減小了50%，總懸浮顆粒（TSP）減少 70%～80%。目前，國外應用防風網(wǎng)的案例有日本電源開發(fā)株式會社下轄的大型發(fā)電廠、中國臺灣臺中火力發(fā)電廠燃煤儲運場、日本東京電廠、荷蘭鹿特丹港務局等，均取得了較好的防塵效果。

我國大陸防風網(wǎng)防塵技術(shù)雖然起步較晚，但近幾年發(fā)展較快。交通部水運科學研究所所做工作較多，對神華天津煤炭碼頭、曹妃甸港、秦皇島煤炭碼頭等防風網(wǎng)工程進行了相關的研究及工程設計工作，現(xiàn)已投入使用，其中于2008年10月竣工的秦皇島港煤三期防風網(wǎng)工程是目前世界上最大的防風網(wǎng)工程，每年可直接減少粉塵排放2000多噸[4]。青島科技大學段振亞等[36]提出了采用防風網(wǎng)和攔沙網(wǎng)的組合形式來抑止二次揚塵的方法，在此基礎上開發(fā)的防風網(wǎng)相繼在某鋼廠 500萬噸球團廠料堆場、西柏坡電廠和山東沾化熱電廠堆煤場投入使用。目前，隨著防風網(wǎng)防塵技術(shù)的日趨成熟，防風網(wǎng)在我國的工程應用范圍正逐漸擴大。

4 數(shù)值模擬研究進展

隨著計算機技術(shù)的快速發(fā)展，越來越多的研究者開始應用數(shù)值模擬方法對防風網(wǎng)后流場進行模擬分析。而應用CFD商業(yè)軟件（如Fluent、CFX）對防風網(wǎng)進行數(shù)值模擬已成為趨勢[37-39]。

4.1 模擬方法

由于防風網(wǎng)后流場表現(xiàn)的是湍流特征，因此其數(shù)值模擬的研究進程很大程度上依賴于湍流模擬的研究進展。目前針對湍流的數(shù)值模擬方法主要是非直接數(shù)值模擬法，即設法對湍流作某種程度的近似和簡化處理，而不直接計算湍流的脈動特性[40]。非直接模擬法主要有兩大類：大渦模擬（LES）法和Renolds平均法（也稱RANS方法）。圖1是湍流數(shù)值模擬方法的分類圖[40]。

大渦模擬法基本思想為用瞬時的N-S方程直接模擬湍流中的大尺度渦，不直接模擬小尺度渦，而小渦對大渦的影響通過近似的模型來考慮。Maruyama[41]利用LES方法對防風網(wǎng)周圍的湍流特性進行了模擬。模擬過程采用亞格子尺度模型使控制方程封閉。由于使用了較粗糙的計算網(wǎng)格，模擬過程的計算量與細網(wǎng)格相比大大減輕，模擬結(jié)果與風洞試驗結(jié)果吻合仍較好。

Renolds平均法是目前使用最廣泛的湍流數(shù)值模擬方法，其核心是不直接求解瞬時的Navier-Stokes方程，而是想法求解時均化的Renolds方程。

Wilson[42]在1985年利用Renolds應力模型對防風網(wǎng)的動力學特性及網(wǎng)周圍的流體運動進行了研究，發(fā)現(xiàn)該模型對近網(wǎng)處湍流運動模擬較準確，而對網(wǎng)上空域內(nèi)高速區(qū)及網(wǎng)后回流區(qū)的模擬結(jié)果不夠穩(wěn)定。在2004年的研究中Wilson[43]分析了RANS模型在模擬網(wǎng)后流場較復雜區(qū)域時出現(xiàn)的不穩(wěn)定性，認為湍流模型的選取對模擬結(jié)果的影響很大。為此，Santiago等[44]比較了3種不同的k-ε湍流模型（標準k-ε、RNG k-ε、realizable k-ε）對防風網(wǎng)后流場模擬的影響。結(jié)果發(fā)現(xiàn)：僅考慮平均風速的變化時，3種模型的模擬結(jié)果均與 Bladely和Mulhearn等[45]的試驗結(jié)果吻合較好，standard k-ε模型對網(wǎng)后回流區(qū)內(nèi)平均風速的模擬較理想，而RNG k-ε和realizable k-ε對網(wǎng)上空域內(nèi)高速區(qū)的模擬較好，但三者的區(qū)別不大。而在對防風網(wǎng)最佳開孔率的模擬研究中，由于考慮了湍流脈動的影響，發(fā)現(xiàn)RNG k-ε和 realizable k-ε模型的模擬結(jié)果要明顯好于 standard k-ε模型。目前，國內(nèi)外對防風網(wǎng)數(shù)值模擬應用較多的湍流模型為 RNG k-ε模型和標準k-ε模型。

由于Renolds平均法避免了直接模擬法計算量大的問題，且模擬效果較好，因此，目前對防風網(wǎng)的數(shù)值模擬所采用的方法大都為Renolds平均法。

圖1 湍流數(shù)值模擬方法及相應的湍流模型

4.2 流場分布

設置防風網(wǎng)后，來流風被分為兩部分：一部分沿著防風網(wǎng)的網(wǎng)頁向上流動，在防風網(wǎng)頂部產(chǎn)生邊界層分離，從而形成湍流剪切層[21]；另一部分由透風孔穿過防風網(wǎng)，稱為滲流風，繼續(xù)向下游運動[46]。由于防風網(wǎng)后流場的復雜性，通過風洞試驗深入了解防風網(wǎng)的空氣動力學特性具有一定的局限性，而數(shù)值模擬則可以更直觀的顯示出網(wǎng)后流體的流動特性。

Lee等[47]通過模擬得出了如圖2所示的流線圖，較直觀的顯示出了防風網(wǎng)后的流場。其研究還表明，隨著開孔率的減小，滲流風相應的減少，在網(wǎng)與料堆迎風面之間和料堆背風面處分別有回流出現(xiàn)，且當開孔率為20%時，由于網(wǎng)后風速被大大減小，使得料堆前后表面的壓力值近乎一致。Wang和Takle[48]利用Renolds平均法對防風網(wǎng)前后的流場進行了模擬，考察了流場內(nèi)流體的運動情況。研究發(fā)現(xiàn)：防風網(wǎng)開孔率的大小決定了其后流場的分布，當開孔率在40%～99%時，防風網(wǎng)后流體為非分離流動（unseparated flow），此時，防風網(wǎng)附近流線的彎曲程度隨開孔率的減小而增大，風速呈U型分布，網(wǎng)后無速度回流區(qū)出現(xiàn)。流場被防風網(wǎng)網(wǎng)分為 3個區(qū)域：網(wǎng)前風速減弱區(qū)、網(wǎng)后風速減弱區(qū)和網(wǎng)上部風速加速區(qū)。而最大風速減弱點在網(wǎng)后（4H，1/3H）（H為網(wǎng)高）處，且高度越大，最大風速減弱點越趨近于防風網(wǎng)。當開孔率在6%～30%時，流體為分離流動（separated flow），此時，近網(wǎng)處流線彎曲程度較大，網(wǎng)后形成了近似三角形的速度回流區(qū)，使得風速呈W型分布，且網(wǎng)后存在兩個最大風速減弱點，分別在網(wǎng)后（1H，0.8H）和（3H，0.2H）處。在對湍動能進行分析時發(fā)現(xiàn)，開孔率為 6%、20%、40%、62%的防風網(wǎng)網(wǎng)后均在近地面處存在較弱的湍流區(qū)，而在靠近網(wǎng)高處則存在較強的湍流區(qū)。

圖2 防風網(wǎng)的作用機理

4.3 導流型防風網(wǎng)

來流風經(jīng)過常規(guī)防風網(wǎng)后，滲流風會直接作用于料堆的迎風面，產(chǎn)生對顆粒較強的作用力，迫使顆粒從料堆表面揚起，形成二次起塵。針對常規(guī)防風網(wǎng)的不足，青島科技大學化學工程研究所[49-50]研究開發(fā)了一種新型防風網(wǎng)，通過導流、整流措施，在常規(guī)平板型防風網(wǎng)網(wǎng)孔上增設導流翅片，使來流風在導流翅片的作用下沿料堆向上爬升，減小風直接沖擊料堆迎風面的作用力，從而進一步改善防風網(wǎng)的防風抑塵效果。新型防風網(wǎng)模型如圖3所示，采用長孔或橢圓孔，網(wǎng)孔的一側(cè)被沖壓成翻邊導流翅片，且翅片與防風網(wǎng)平面呈一定角度，因而對來流風起到導流作用。

作者課題組[51]應用CFD模擬軟件Fluent 6.2分別對常規(guī)平板型防風網(wǎng)與新型防風網(wǎng)后的流場進行了模擬。數(shù)值模擬的計算區(qū)域設置為長2000 mm，寬600 mm，高600 mm的長方體空間；防風網(wǎng)厚2 mm，高度H=100 mm，寬度為600 mm，開孔排布近似為菱形分布，開孔率為 38.5%，防風網(wǎng)后設置三角形料堆，其高度為75 mm，網(wǎng)堆距為160 mm。利用網(wǎng)格生成軟件Gambit對模型進行前處理，生成的網(wǎng)格如圖4、圖5所示。

圖3 導流型防風網(wǎng)結(jié)構(gòu)示意圖

圖4 防風網(wǎng)流場模擬的網(wǎng)格劃分

圖5 導流型防風網(wǎng)附近局部網(wǎng)格

圖6 常規(guī)型防風網(wǎng)后流場速度矢量圖

圖7 導流型防風網(wǎng)后流場速度矢量圖

模擬采用Renold時均方法，利用standard k-ε湍流模型使方程組封閉。圖6、圖7為防風網(wǎng)后流場速度矢量圖。可明顯看出：導流板具有良好的導向作用。與常規(guī)防風網(wǎng)比較，導流型防風網(wǎng)后料堆表面速度較小，滲流風速減小40%左右；頂部附近風速減小19%，壓力波動最大可減小60%；湍流強度變化較小，防風抑塵效果好。

5 結(jié) 語

防風網(wǎng)作為一種能控制起塵和擴散的有效手段，目的是減小其下風向的風速。目前對防風網(wǎng)防風抑塵效果影響因素的考察主要集中在開孔率、網(wǎng)高、網(wǎng)與料堆的距離等方面；開孔率是影響防風網(wǎng)作用的關鍵因素，但國內(nèi)外對開孔率的研究得出的結(jié)果并不一致，一般認為，最佳開孔率均在20%～50%；網(wǎng)高是影響防風網(wǎng)庇護范圍的重要因素，隨網(wǎng)高的增加防風網(wǎng)的庇護范圍增大，但持續(xù)增大網(wǎng)的高度并不能持續(xù)擴大庇護范圍；與開孔率和網(wǎng)高相比，網(wǎng)與料堆距離的影響并不顯著，但在網(wǎng)后2～3倍料堆高度的范圍內(nèi)庇護作用最好。

目前研究者們對防風網(wǎng)的研究取得了長足的進步，研究方向主要集中在防風網(wǎng)防風抑塵方面，如防風網(wǎng)形式、防風網(wǎng)布置、防風網(wǎng)防風抑塵性能的影響因素等方面[52-56]。防風網(wǎng)抑塵技術(shù)逐漸成熟，其工程應用已越來越廣泛。而隨著計算機技術(shù)的不斷發(fā)展，利用CFD軟件對防風網(wǎng)進行模擬已成為研究熱點，其研究結(jié)果對人們深入認識防風網(wǎng)并指導裝置的設計優(yōu)化意義重大。但這些研究大多集中在流體動力學方面，關于防風網(wǎng)結(jié)構(gòu)設計方法、結(jié)構(gòu)穩(wěn)定性、抗風設計研究[57]和承載安全度分析方法以及新型防風網(wǎng)的開發(fā)等的研究報道還較少，應是今后研究的重點。此外，在理論上應建立更確切的流場數(shù)學模型、網(wǎng)后揚塵量的計算模型以及顆粒揚塵的判據(jù)以對防風網(wǎng)的抑塵機理進行更深入的研究。

[1] 趙海珍，梁學功，馬愛進，等.防風網(wǎng)防塵技術(shù)及其在我國大型煤炭港口的應用與發(fā)展對策[J].環(huán)境科學研究，2007，20（2）： 68-71.

[2] Ucar T，Hall F R.Windbreaks as a pesticide drift mitigation strategy：a review[J].Pest Management Science，2001，57（8）：663-675.

[3] 陳建華，詹水芬.港口散貨堆場防風網(wǎng)防塵技術(shù)研究和應用[J].珠江水運，2008，（3）：44-46.

[4] 馬蘭.港口散貨堆場防風網(wǎng)板工藝及建設方案優(yōu)化研究[D]. 天津：天津大學，2009.

[5] 張光玉，陳立，王奇志，等. 秦皇島港煤堆場防風網(wǎng)風洞試驗研究[J].交通環(huán)保，2003，24（1）：4-6.

[6] 陳凱，朱鳳榮，鈕珍南，等.防風網(wǎng)作用效果的的風洞試驗評估[J].北京大學學報，2006，45（2）：636-640.

[7] 祁祥，程復，武烽東，等.防風網(wǎng)防風效能初步研究[J].水土保持研究，2007，14（1）：191-195.

[8] Wang H，Takle E S.On three-dimensionality of shelterbelt structure and its influences on shelter effects[J].Boundary-Layer Meteorology，1996，79（1）：83-105.

[9] Judd M J，Raupach M R，F(xiàn)innigan J J.A wind tunnel study of turbulent flow around single and multiple windbreaks，Part I：Velocity fields[J].Boundary-Layer Meteorology，1996，80（1）：127-165.

[10] Takahashi S，Du M Y，Wu P M，et al.Three dimensional numerical simulation of the flow over complex terrian with windbreak hedge[J].Environmental Modelling and Software，1998，13（3-4）：257-265.

[11] Wang H，Takle E S.Numerical simulations of shelterbelt effects on wind direction[J].Journal of Applied Meteorology，1995，34（10）：2206-2219.

[12] Fang F M，Wang D Y.On the flow around a vertical porous fence[J].Journal of Wind Engineering and Industrial Aerodynamics，1997，67：415-424.

[13] Dong Zhibao，Qian Guangqiang，Luo Wanyin，et al.Threshold velocity for wind erosion： The effects of porous fences [J].Environmental Geology，2006，51（3）：471-475.

[14] Lee S J，Park K C，Park C W.Wind tunnel observations about the shelter effect of porous fences on the sand particle movements [J].Atmospheric Environment，2002，36（9）：1453- 1463.

[15] 林官明，葉文虎.防風網(wǎng)泄流區(qū)湍流的子波分析[J].北京大學學報：自然科學版，39（5）：732-735.

[16] Loeffler A E，Gordon A M，Gillespie T J.Optical porosity and windspeed reduction by coniferous windbreaks in Southern Ontario[J].Agroforestry Systems，1992，17（2）：119-133.

[17] Perera M.Shelter behind two-dimensional solid and porous fences[J].Journal of Wind Engineering and Industrial Aerodynamics，1981，8（1-2）：93-104.

[18] Heisler G M，Dewalle D R.Effects of windbreak structure on wind flow[J].Agriculture，Ecosystems and Environment，1988，22：41-69.

[19] Wilson J D.On the choice of a windbreak porosity profile [J].Boundary-Layer Meteorology，1987，38（1）：37-49.

[20] Wilson J D.Numerical studies of flow through a wind break [J].Journal of Wind Engineering and Industrial Aerodynamics，1985，21：119-154.

[21] Boldes U，Colman J，Maranon Di Leo J.About the penetration of a horizontal axis cylindrical vortex into the nearby downwind region of a vertical porous fence[J].Journal of Wind Engineering and Industrial Aerodynamics，2003，91（7）：859-872.

[22] Li Wei，Wang Fang，Bell Simon.Simulating the sheltering effects of windbreaks in urban outdoor open space[J].Journal of Wind Engineering and Industrial Aerodynamics，2007，95（7）：533-549.

[23] Lee S J，Park C W.Surface-pressure variations on a triangular prism by porous fences in a simulated atmospheric boundary layer [J]. Journal of Wind Engineering and Industrial Aerodynamics，1998，73（1）：45-58.

[24] Seginer IDO.Windbreak drag calculated from the horizontal velocity field[J].Boundary -Layer Meteorology，1972，3（1）：87-97.

[25] Castro I P.Wake characteristics of two-dimensional perforated plates normal to an air-stream[J].Journal of Fluid Mechanics，1971，46（3）：599-609.

[26] Ranga K G，Garde R J，Singh S K，et al.Experimental study on characteristics of flow past porous fences[J].Journal of Wind Engineering and Industrial Aerodynamics，1988，29（1-3）：155-163.

[27] Kim H B，Lee S J.The structure of turbulent shear flow around a two-dimensional porous fence having a bottom gap[J].Journal of Fluids and Structures，2002，16（3）：317-329.

[28] Park C W，Lee S J. Experimental study on surface pressure and flow structure around a triangular prism located behind a porousfence[J].Journal of Wind Engineering and Industrial Aerodynamics，2003，91（1-2）：165-184.

[29] Lee S J，Kim H B.Velocity field measurements of flow around a triangular prism behind a porous fence[J].Journal of Wind Engineering and Industrial Aerodynamics，1998，77：521-530.

[30] Mercer G N.Modelling to determine the optimal porosity of shelterbelts for the capture of agricultural spray drift [J].Environmental Modelling and Software，2009，24（11）：1349-1352.

[31] Torano J，Torno S，Diego I，et al.Dust emission calculations in open storage piles protected by means of barriers，CFD and experimental tests[J].Environmental Fluid Mechanics，2009，9（5）：493-507.

[32] 陳凱華，宋存義，李強，等.鋼鐵廠露天料堆場防風抑塵墻效果的數(shù)值模擬[J].環(huán)境工程學報，2008，2（3）：404 – 407.

[33] Kim H B，Lee S J.Hole diameter effect on flow characteristics of wake behind porous fences having the same porosity[J].Fluid Dynamics Research，2001，28（6）：449-464.

[34] Yeh C P，Tsai C H，Yang R J.An investigation into the sheltering performance of porous windbreaks under various wind directions[J].Journal of Wind Engineering and Industrial Aerodynamics，2010.

[35] Lee S J，Park C W.The shelter effect of porous fences on coal piles in POSCO open storage yard[J].Journal of Wind Engineering and Industrial Aerodynamics，2000，84（1）：101-118.

[36] 段振亞，石文梅，鄭文娟，等.防風網(wǎng)抑塵機理研究及工程應用進展[J].石油化工設備，2010，39（3），40-44.

[37] Diego I，Pelegry A，Trono S，et al.Simultaneous CFD evaluation of wind flow and dust emission in open storage piles[J].Applied Mathematical Modelling，2009，33（7）：3197-3207.

[38] Alhajraf S.Computational fluid dynamic modeling of drifting particles at porous fences[J].Environmental Modelling and Software，2004，19（2）：163-170.

[39] Bitog J P，Lee I B，Shin M H，et al.Numerical simulation of an array of fences in Saemangeum reclaimed land[J]. Atmospheric Environment，2009，43（30）：4612-4621.

[40] 王福軍.計算流體動力學分析-CFD軟件原理與應用[M]. 北京：清華大學出版社，2004：116-118.

[41] Maruyama T.Large eddy simulation of turbulent flow around a windbreak[J]. Journal of Wind Engineering and Industrial Aerodynamics，2008，96（10-11）：1998-2006.

[42] Wilson J D.Numerical studies of flow through a windbreak [J]. Journal of Wind Engineering and Industrial Aerodynamics，1985，21：119-154.

[43] Wilson J D.Oblique，stratified winds about a shelter fence.Part Ⅱ：Comparision of measurements with numerical models [J].Journal of Applied Meteorology，2004，43（8）：1149-1167.

[44] Santiago J L，Martin F，Cuerva A，et al.Experimental and numerical study of wind flow behind windbreaks[J].Atmospheric Environment，2007，41（30）：6406-6420.

[45] Bladely E F，Mulhearn P J.Development of velocity and shear stress distributions in the wake of a porous shelter fence[J].Journal of Wind Engineering and Industrial Aerodynamics，1983，14：145-156.

[46] Raine J K，Stevenson D C.Wind protection by model fences in a simulated atmospheric boundary layer[J]. Journal of Wind Engineering and Industrial Aerodynamics，1977，2（2）：159-180.

[47] Lee S J，Lim H C.A numerical study on flow around a triangular prism located behind a porous fence[J].Fluid Dynamics Research，2001，28（3）：209-221.

[48] Wang H，Takle E S.A numerical simulation of boundary layer flows near shelterbelts[J].Boundary-Layer Meteorology，1995，75（1）：141-173.

[49] 陳光輝，李建隆，王偉文，等.露天散狀物料堆放場的防風抑塵方 法 及 其 防 風 抑 塵 網(wǎng) 和 用 途 ： 中 國 ，200810249627.5[P].2009-7-22.

[50] 范軍領，王偉文，張攀，等.防風抑塵網(wǎng)：中國，200820232839[P].2009-12-16.

[51] 董繼鵬.強風流過散堆料場的流場模擬與抑塵研究[D].青島：青島科技大學，2009.

[52] 陳愛英.防風網(wǎng)的數(shù)值模擬[D].上海：同濟大學應用數(shù)學系，2007.

[53] Frank C，Ruck B.Double-arranged mound-mounted shelterbelts：Influence of porosity on wind reduction in the intermediate field[J].Environmental Fluid Mechanics，2005，5：267-292.

[54] 李建隆，董繼鵬，陳光輝，等.防風抑塵網(wǎng)開孔形式對流場的影響[J].環(huán)境工程學報，2009，3（9）：1725-1728.

[55] Cleugh H A.Effects of windbreaks on airflow，microclimates and crop yields[J].Agroforestry System，1998，41（1）：55-84.

[56] Dong Zhibao，Luo Wanyin，Qian Guangqiang，et al.A wind tunnel simulation of the turbulencefields behind upright porous wind fences[J].Journal of Arid Environments，2010，74（2）：193-207.

[57] 徐洪濤，何勇，廖海黎，等.防風網(wǎng)氣動特性參數(shù)的試驗研究[J].安全與環(huán)境學報，2010，10（1）：70-74.

Progress of porous fence

SUN Changfeng，CHEN Guanghui，F(xiàn)AN Junling，LI Jianlong
（Department of Chemical Engineering，Qingdao University of Science and Teconology，Qingdao 266042，Shandong，China）

X 513

A

1000–6613（2011）04–0871–07

2010-12-08；修改稿日期：2011-01-16。

孫昌峰（1985—），男，碩士研究生。聯(lián)系人：李建隆，教授，博士生導師，主要從事多相流體流動與分離。E-mail ljlong@qust.edu .cn。