沈國輝　姚旦　樓文娟　邢月龍　潘峰　郭勇

摘要：采用CFD模擬方法研究單山和雙山情況下三維山丘風(fēng)場，研究計算模型表面粗糙度對風(fēng)場的影響，計算不同坡度山體情況下單山的風(fēng)場，進行兩個山體左右排列情況的風(fēng)場計算，分析坡度、風(fēng)向角、間距對雙山風(fēng)場的影響.研究表明：山體的計算模型表面粗糙度增大時，山頂上方半山高度的加速效應(yīng)減弱，山后尾流區(qū)的高度增加；山體橫風(fēng)切面的加速效應(yīng)大于順風(fēng)切面，橫風(fēng)切面內(nèi)半山以上的位置均為風(fēng)速最大值區(qū)域；左右雙山緊貼排列情況下，風(fēng)斜吹時前山山頂?shù)募铀俦却笥诤笊缴巾?，風(fēng)直吹的數(shù)據(jù)在兩者中間.

關(guān)鍵詞：CFD模擬；平均風(fēng)速；風(fēng)場計算；三維山體；峽谷風(fēng)效應(yīng)

中圖分類號：TU312.1 文獻標(biāo)識碼：A

Abstract：CFD numerical simulation method was employed to study the wind field characteristics of threedimensional hills under isolated hill and two adjacent hills condition. The influence of the calculation model surface roughness on the wind field was studied. The wind fields on isolated hill with different slopes were calculated. The wind fields of two adjacent hills in leftright arrangement were also investigated， and the influence of the hill slope， wind azimuth and hill distance on the wind field was studied. The results show that， if the model surface roughness increases， the speedup effect of half hill height above the top of the hill decreases and the height of vortex region increases. The speedup effects on the crosswind plane of an isolated hill are more significant than those on the alongwind plane， and the most unfavorable position ranges from the half hill height to the top of the hill on the crosswind plane. When two adjacent hills are in leftright arrangement with zero distance， the speedup ratios on the front hill are larger than those on the rear hill under yawed azimuth， and the speedup ratio， when wind blows perpendicularly to the two hills， is between these two data.

Key words：CFD simulation； mean wind velocity； wind field calculation； three dimensional hill； wind tunneling effect

風(fēng)流經(jīng)山地時，受到地形的影響，會形成復(fù)雜的山地風(fēng)場，該風(fēng)場與平地風(fēng)場有較大差異.各國規(guī)范[1-4]均給出了山地風(fēng)場的修正方法，中國GB規(guī)范[1]對風(fēng)壓高度變化系數(shù)進行了修正，其他規(guī)范[2-4]對風(fēng)速進行了修正.但多數(shù)規(guī)范[1-3]的修正方法通常針對二維山體情況，沒有給出三維山體情況下的數(shù)據(jù)，同時各國規(guī)范[1-4]對于兩個山體形成的風(fēng)場沒有進行規(guī)定.

山體風(fēng)場特征的研究方法主要有現(xiàn)場實測[5]、風(fēng)洞試驗[6-8]和CFD數(shù)值模擬[6-12].CFD數(shù)值模擬可以按實際尺寸進行計算，可避免風(fēng)洞試驗只能進行縮尺試驗的不足，并具有成本低、速度快、資料完備等優(yōu)點，本文采用該方法進行研究.以往研究中，Kim等[6]對兩個二維山體前后排列的情況進行模擬，分析后山的尾流區(qū)尺寸；Breuer等[7]通過數(shù)值模擬和風(fēng)洞試驗研究了二維周期性山體的風(fēng)場.李朝[8]進行不同山體風(fēng)場的計算，并與風(fēng)洞試驗結(jié)果進行比較；陳平[9]對典型山體的風(fēng)場進行模擬，并分析不同山體組合山地情況；祝志文等[10]對峽谷地貌大跨度橋梁選址處進行數(shù)值模擬，獲得橋址處各位置的風(fēng)場數(shù)據(jù)；肖儀清等[11]針對某復(fù)雜地形進行模擬，并與實測風(fēng)速進行了比較；魏奇科等[12]、李鑫[13]、李正良等[14]對一系列山體的風(fēng)速進行CFD和風(fēng)洞試驗研究，并進行比較.

本文針對單山和雙山情況下三維山丘的風(fēng)場，研究CFD模擬中表面粗糙度對山體風(fēng)場的影響，進行單個山體模型的模擬計算并分析山體坡度的影響，進行兩個山體左右排列情況的計算，分析山體坡度、風(fēng)向角、山體間距對雙山風(fēng)場的影響，研究成果可為復(fù)雜山地的風(fēng)場特性提供參考.需要說明的是，本文只給出山體上順風(fēng)向的平均風(fēng)速，不涉及脈動風(fēng)速，也不涉及其他風(fēng)向的風(fēng)速.

1CFD的建模過程和湍流模型

山體高度H均取100 m，形狀為余弦型，方程為：

z（x，y）=Hcos 2（π（x2+y2）1/2D）（1）

式中：z為高度方向；x，y為水平方向；D為山體底部直徑.對7組坡度山體進行CFD建模計算，山體底部直徑分別為200 m，300 m，400 m，500 m，600 m，700 m和800 m.

計算的流場高度統(tǒng)一取為8H，流場的寬度依據(jù)山體直徑而變化，以典型陡坡山體（D=300 m）為例，流場寬取為6D，上游長度取為3D，下游長度取為4.5D.流場內(nèi)全部采用結(jié)構(gòu)化六面體單元網(wǎng)格進行劃分，網(wǎng)格尺寸在山體表面最密，然后向外擴散.所有工況的阻塞率全部小于3%，網(wǎng)格數(shù)各異，數(shù)量在70萬到340萬之間，典型陡坡網(wǎng)格劃分情況如圖1所示，圖中還給出了來流方向和x，y，z坐標(biāo)定義.

湍流模型選用被廣泛應(yīng)用于模擬大氣邊界層的Realizable kε湍流模型[9-15].計算中地貌取B類（α=0.15），平均風(fēng)速u和湍流度Iu依據(jù)荷載規(guī)范[1]：

u=u10z100.15（2）

Iu=0.14z10-0.15（3）

式中：u10為基本風(fēng)速，計算中取30 m/s.入口處湍流動能k和湍流耗散率ε通過UDF接口編程輸入，表達式如下：

k=32（uIu）2（4）

ε=C3/4μk3/2l（5）

式中：Cμ=0.09，l為湍流積分尺度，參考日本AIJ2004規(guī)范[3]的取值：

l=100z300.530 m

100z<30 m.（6）

式中：zg為梯度風(fēng)高度.入流面設(shè)定為速度入口，出流面設(shè)定為壓力出口，流場兩側(cè)和頂部設(shè)定為對稱面，山體表面和周圍的平地面設(shè)定為壁面，選用非平衡壁面函數(shù)模擬近壁面的流動，數(shù)值離散格式均采用二階迎風(fēng)格式.

2計算模型的山體表面粗糙度選擇

山體表面粗糙度的設(shè)定對山體風(fēng)場的計算有一定的影響，對各種模型表面粗糙度進行試算.以典型陡坡山體（D=300 m，H=100 m）為例，山體表面粗糙高度取7種，分別為0 m，0.05 m，0.1 m，0.2 m，0.5 m，1 m和2 m，地面粗糙度取0 m.圖2給出了0 m，0.1 m和1 m粗糙高度情況下的流線圖，由圖可知：1）隨著粗糙高度增加，山體后部尾流區(qū)的尺寸擴大；2）當(dāng)粗糙高度為0 m時，尾流區(qū)很??；當(dāng)粗糙高度為0.1 m時，尾流區(qū)長度方向延伸到后側(cè)山腳之后D/4的距離，寬度約D/3，高度約0.6H；當(dāng)?shù)孛娲植诟叨葹? m時，尾流區(qū)長度方向延伸到后側(cè)山腳之后D/2的距離，寬度約0.6D，高度約0.8H.

圖7給出了離地面（山體）表面距離為5 m和40 m的風(fēng)速加速比，可以發(fā)現(xiàn)：1）最大加速比出現(xiàn)在山頂位置；2）在離地面（山體）表面距離為5 m處，山側(cè)往迎風(fēng)方向大約15°的位置存在一條加速比極大值線；3）山后尾流區(qū)的加速比很??；4）前山腳附近有較小的減速區(qū)；5）隨著離地面（山體）表面距離增加，加速比趨近于1，但是尾流區(qū)范圍內(nèi)的加速比仍然很小.

3.2山體坡度對單山風(fēng)場的影響

圖9給出了7種山體坡度下山頂位置的平均風(fēng)速和加速比，山體底部直徑分別為200 m，300 m，400 m，500 m，600 m，700 m和800 m.由圖可知：1）不同坡度山體在山頂位置均出現(xiàn)了較大的加速效應(yīng)；2）離山頂表面約30 m以下，坡度大的山體加速比大；3）離山頂表面約30 m至100 m是一個過渡區(qū)，加速比隨坡度的變化規(guī)律逐漸改變；4）離山頂表面100 m以上，不同坡度山體的平均風(fēng)速較接近.

圖10給出不同山體坡度下山頂位置離山體表面不同高度的最大加速比.由圖可知：1）在離山頂表面5 m高度，隨著標(biāo)準坡度增大，加速比增大且漸趨平緩；2）在離山頂表面10 m，20 m，50 m高度，隨著標(biāo)準坡度增大，加速比先增后減；3）在離山頂表面100 m，200 m高度，隨著標(biāo)準坡度增大，加速比基本上呈略微下降趨勢；4）各離頂高度的加速比曲線在標(biāo)準坡度100%的位置，都已經(jīng)非常平坦.

4陡坡雙山左右緊貼排列（間距為零）的風(fēng)場

4.1典型陡坡雙山的風(fēng)場

計算典型陡坡山體（D=300 m）雙山左右緊貼排列（山體間距0 m）的風(fēng)場，圖11給出了離地面（山體）表面距離5 m高度處的加速比.由圖可知：1）左右排列雙山加速比最大的位置依然是山頂附近；2）前山腳附近有較小的減速區(qū)，且兩減速區(qū)之間的區(qū)域也存在一定的減速；3）左右雙山山后部的尾流區(qū)在山后連成一體，形成一個更大的尾流區(qū).

4.2山體坡度對雙山風(fēng)場的影響

計算兩個山體左右緊貼排列情況7種山體坡度時的風(fēng)場，山體底部直徑分別為200 m，300 m，400 m，500 m，600 m，700 m和800 m，圖12給出7種坡度雙山山坳中的平均風(fēng)速和加速比.由圖可知：1）7種坡度山體山坳位置的風(fēng)速均存在較明顯的加速效應(yīng)；2）離地約200 m以上，隨著坡度減少加速比呈增大趨勢，但增幅并不大；3）離地100 m以下，加速比隨坡度的變化比較復(fù)雜，無明顯規(guī)律.

從圖12（b）還可發(fā)現(xiàn)，離地100 m以下，直徑300 m的加速比是7種坡度山體中的最大值.直徑300 m山體比底部直徑更大山體的加速比大是很好理解的，因為兩山中間的斷面增大了，風(fēng)速自然小了.但怎么理解比底部直徑200 m的加速比大，下面從繞流方式角度進行分析.圖13給出了底部直徑200 m和300 m山體10 m高度的流線，可以發(fā)現(xiàn)對于直徑200 m的山體，有3條流線從外側(cè)繞流；對于直徑300 m山體，只有2條流線從外側(cè)繞流.說明來流流經(jīng)底部直徑200 m山體時經(jīng)峽谷的阻力增大，有較多的來流從雙山的外側(cè)繞流，而不經(jīng)過峽谷，因而峽谷中底部直徑300 m山體的加速比大于直徑200 m山體.

4.3風(fēng)向角對峽谷風(fēng)風(fēng)場的影響

進行6個風(fēng)向角下左右緊貼排列雙山風(fēng)場的計算，風(fēng)向角分別為0°，5°，10°，15°，20°，25°和30°，圖14給出了10°和30°風(fēng)向角下離地面（山體）距離5 m高度處的加速比.由圖可知：1）山頂上的加速比分布同單山情況類似；2）峽谷尾流區(qū)的加速比等值線隨著風(fēng)向角增大發(fā)生偏轉(zhuǎn)，與風(fēng)向大致平行.

圖15給出了山頂位置離山頂表面5 m，10 m和20 m高度的加速比，圖中風(fēng)向角為正值表示前山，負值表示后山.由圖15可知：1）對于這3個高度處的加速比，前山（風(fēng)向角為正）的數(shù)據(jù)大于0°（左右排列）情況，0°（左右排列）的數(shù)據(jù)大于后山；2）這3個高度處的加速比均在15°風(fēng)向角時達到最大.

5山體間距對左右排列雙山風(fēng)場的影響

計算典型陡坡山體（D=300 m）左右排列雙山在5種間距情況的風(fēng)場計算，間距分別為0 m，50 m，100 m，200 m和300 m.圖16給出0 m，50 m和200 m間距下過雙山軸線橫風(fēng)切面內(nèi)的風(fēng)速加速比，由圖可知：1）在不同間距下山體上的加速比數(shù)據(jù)較接近；2）兩山正中間（y=0）處的風(fēng)速隨著山腳間距增大而減小.

圖17給出了不同山體間距下內(nèi)側(cè)山腳位置（例：山腳間距為200 m時，內(nèi)側(cè)山腳即x=0 m，y=±100 m的位置）的風(fēng)速加速比.由圖可知：1）5種山體間距下，內(nèi)側(cè)山腳位置的風(fēng)速均存在明顯的加速效應(yīng)；2）山腳間距為0 m時，加速效應(yīng)最強，隨著山腳間距的增大，加速效應(yīng)減弱；3）山腳間距為200 m和300 m時內(nèi)側(cè)山腳的加速效應(yīng)很接近.

6結(jié)論

本文采用CFD數(shù)值方法研究單山和雙山情況下的三維風(fēng)場特性，得到以下結(jié)論：

1）CFD數(shù)值模擬中，表面粗糙度的設(shè)定對計算結(jié)果的計算有一定的影響，當(dāng)山體表面的粗糙度增大時，山頂上方一半山體高度范圍內(nèi)的風(fēng)速加速效應(yīng)減弱，山后尾流區(qū)的高度增加.

2）單個山體風(fēng)場的計算結(jié)果表明，山體橫風(fēng)切面的加速效應(yīng)大于順風(fēng)切面，橫風(fēng)切面內(nèi)半山以上的位置均出現(xiàn)較大的風(fēng)速增大現(xiàn)象，該現(xiàn)象即為“孤峰繞流效應(yīng)”；加速比在貼近山體位置最大，沿著高度方向遞減，距離山體越遠加速比越趨近于1；山后尾流區(qū)的加速比很小，前山腳附近有較小的減速區(qū)；離山頂表面約30 m以下，坡度大的山體加速比大.

3）雙山左右緊貼排列的計算結(jié)果表明，加速比最大的位置在山頂附近；山坳中離地200 m以上，隨著坡度減少加速比呈增大趨勢，但增幅并不大，離地100 m以下，直徑300 m的加速比是7種坡度山體中的最大值；峽谷尾流區(qū)的加速比等值線隨著風(fēng)向角增大發(fā)生偏轉(zhuǎn)，與風(fēng)向大致平行；在有風(fēng)向角情況下，前山山頂?shù)募铀俦却笥诤笊缴巾敚?°風(fēng)向角的數(shù)據(jù)在兩者中間.

4）不同山體間距下，山體上的加速比數(shù)據(jù)較接近，但在山體中間存在差異；兩山正中間處的風(fēng)速隨著山腳間距增大而減小；山腳間距為0 m時，內(nèi)側(cè)山腳位置的風(fēng)速加速效應(yīng)最強，隨著山腳間距的增大，加速效應(yīng)減弱，當(dāng)距離為200 m和300 m時內(nèi)側(cè)山腳的加速效應(yīng)很接近.

參考文獻

[1]GB 50009-2012 建筑結(jié)構(gòu)荷載規(guī)范[S]. 北京：中國建筑工業(yè)出版社，2012：32-33.

GB 50009-2012 Load code for the design of building structures [S]. Beijing： China Architecture & Building Press， 2012：32-33. （In Chinese）

[2]AS/NZS 1170.0 Australian/New Zealand standard[S]. Sydney： AS/NZS， 2002：23-25.

[3]AIJ 2004 Recommendations for loads on buildings[S]. Tokyo： Architectural Institute of Japan， 2004：611-613.

[4]ANSI/ASCE 7-10 Minimum design loads for buildings and other structures[S]. New York： ASCE， 2010：315-316.

[5]WALMSLEY J L， TAYLOR P A. Boundarylayer flow over topography： impacts of the askervein study[J]. Boundarylayer Meteorological， 1996， 78：291-320.

[6]KIM H G， LEE C M， LIM H C， et al. An experimental and numerical study on the flow over twodimensional hills[J]. Journal of Wind Engineering and Industrial Aerodynamics， 1997， 66（1）： 17-33.

[7]BREUER M， PELLER N， RAPP C， et al. Flow over periodic hillsnumerical and experimental study in a wide range of reynolds numbers[J]. Computers & Fluids， 2009， 38（2）： 433-457.

[8]李朝. 近地湍流風(fēng)場的CFD模擬研究[D]. 深圳： 哈爾濱工業(yè)大學(xué)深圳研究生院， 2010：42-53.

LI Chao. Study on CFD simulation of turbulence wind field near ground[D]. Shenzhen： Shenzhen Graduate School， Harbin Institute of Technology， 2010：42-53. （In Chinese）

[9]陳平. 地形對山地丘陵地區(qū)風(fēng)場影響的數(shù)值研究[D]. 杭州： 浙江大學(xué)建筑工程學(xué)院， 2007：23-35.

CHEN Ping. Numerical study of terrain influence on the airflow over hilly land[D]. Hangzhou： College of Civil Engineering and Architecture， Zhejiang University， 2007：23-35. （In Chinese）

[10]祝志文， 張士寧， 劉震卿， 等. 橋址峽谷地貌風(fēng)場特性的CFD模擬[J]. 湖南大學(xué)學(xué)報：自然科學(xué)報， 2011， 38（10）：13-17.

ZHU Zhiwen， ZHANG Shining， LIU Zhenqing， et al. CFD simulation of wind field at bridge site on gorge terrain[J]. Journal of Hunan University：Natural Sciences， 2011， 38（10）：13-17. （In Chinese）

[11]肖儀清， 李朝， 歐進萍， 等. 復(fù)雜地形風(fēng)能評估的CFD方法[J]. 華南理工大學(xué)學(xué)報：自然科學(xué)版， 2009， 37（9）： 30-35.

XIAO Yiqing， LI Chao， OU Jinping， et al. CFD approach ot evaluation of wind energy in complex terrain[J]. Journal of South China University of Technology：Natural Science Edition， 2009， 37（9）： 30-35. （In Chinese）

[12]魏奇科， 李正良， 孫毅. 山地風(fēng)加速效應(yīng)的計算模型[J]. 華南理工大學(xué)學(xué)報：自然科學(xué)版， 2010， 38（11）： 54-58.

WEI Qike， LI Zhengliang， SUN Yi. Calculation model of speedup effect of wind velocity in hilly terrain[J]. Journal of South China University of Technology：Natural Science Edition， 2010， 38（11）： 54-58. （In Chinese）

[13]李鑫. 山地地形的近地風(fēng)場特性研究[D]. 重慶：重慶大學(xué)土木工程學(xué)院， 2010：21-40.

LI Xin. Study on characteristics of wind field in ground layer on hilly terrain[D]. Chongqing： College of Civil Engineering， Chongqing University， 2010：21-40. （In Chinese）

[14]李正良， 魏奇科， 孫毅. 復(fù)雜山地風(fēng)場幅值特性試驗研究[J]. 工程力學(xué)， 2012， 29（3）：184-191.

LI Zhengliang，WEI Qike， SUN Yi. Experimental research on amplitude characteristics of complex hilly terrain wind field[J]. Engineering Mechanics， 2012， 29（3）：184-191. （In Chinese）

[15]姚旦. 山丘地形風(fēng)場特性及對輸電塔的風(fēng)荷載作用研究[D].杭州： 浙江大學(xué)建筑工程學(xué)院， 2014：9-39.

YAO Dan. Research on characteristics of wind field on hilly terrain and its wind load effect on lattice transmission towers[D]. Hangzhou： College of Civil Engineering and Architecture， Zhejiang University， 2014：9-39. （In Chinese）