蔣光彪，陽　杰，鐘　海，殷水平，李其勇（中南林業(yè)科技大學，湖南 長沙 410004）

高層建筑風場數(shù)值模擬研究

蔣光彪，陽杰，鐘海，殷水平，李其勇
（中南林業(yè)科技大學，湖南 長沙 410004）

本文以FLUENT軟件和高性能集群作為數(shù)值模擬平臺，通過改變高層建筑物高度進行了并行數(shù)值模擬。并對并行數(shù)值模擬結(jié)果和集群并行計算效率進行了分析，分析結(jié)果表明：隨著建筑物高度的增加，建筑物周圍最大風速和最大風速比增大，建筑物對其周圍空氣流動的影響增強。在并行計算中，當CPU個數(shù)相同時，隨著網(wǎng)格數(shù)量的增加，并行效率呈增加趨勢但增大的幅度越來越??；當網(wǎng)格數(shù)相同時，隨著CPU個數(shù)增加，并行計算效率呈減小趨勢。

FLUENT；高層建筑；數(shù)值模擬；并行計算；風場

高層和超高層建筑的大量涌現(xiàn)引起了的局部大氣邊界層的改變，會對低層氣流和湍流產(chǎn)生顯著的影響，特殊的城市下墊面，具有較大的粗糙度，可以引起更強的機械湍流，改變局部地的風場結(jié)構(gòu)，對日常生活造成影響［1］。風環(huán)境是空氣氣流在建筑內(nèi)外空間的流動狀況及其對建筑物使用的影響，因此，風環(huán)境本身形成了一個系統(tǒng)。若在規(guī)劃和設(shè)計中忽略了風環(huán)境問題，在一般的氣候條件下，它們將直接影響著城市環(huán)境的小氣候和環(huán)境給予人們的舒適感；一旦遇到大風，這種影響往往會形成所謂“城市風害”［2］，影響到建筑物本身的某些使用功能和行人、行車安全等方面的問題；特別對于高層建筑的背風側(cè)會形成較大的回流區(qū)域，如果在設(shè)計之初沒能充分考慮其風環(huán)境，而使污染物或有害排放物能夠進入這一回來區(qū)域，則會造成局部環(huán)境的空氣品質(zhì)惡化，因為這一回流區(qū)域內(nèi)的污染物和有害物容易積聚，很難被稀釋和擴散［3-5］。因此，掌握高層建筑物的風環(huán)境特征，構(gòu)建有利于通風和污染物消散，避免人行高度出現(xiàn)過高風速和過急渦流的風環(huán)境布局，對城市環(huán)境和安全具有重要意義。目前建筑風場的預測方法主要有現(xiàn)場實測、風洞試驗和數(shù)值模擬，其中數(shù)值模擬已成為風環(huán)境預測的最有效手段［6，7］。本文采用CFD數(shù)值模擬的方法，利用標準k-ε湍流模型，模擬單體建筑在不同情況下的的繞流風場，研究其建筑物的風場風速的分布特征。

1　問題的描述及數(shù)值模擬方法

1.1問題的提出

本文設(shè)計了如圖1簡化的方形高層建筑類繞流湍流場數(shù)值模擬實例，通過改變方形建筑物的高度進行k-ε湍流模型模擬，研究了計算流域內(nèi)的速度場和壓力場的分布情況，將不同工況的數(shù)值模擬求解結(jié)果進行對比分析。

本文在研究建筑物高度對建筑物周圍風環(huán)境的影響時，風速取2.5m/s，選擇的計算區(qū)域為：長度×寬度×高度=180m×90m×150m，并且為了讓建筑背后流場能充分發(fā)展，計算模型中建筑物左側(cè)距離入口邊界為60m占總計算流域長度的三分之一，距離出口邊界為120m占總計算流域的三分之二。

1.2數(shù)值方法

1.2.1網(wǎng)格劃分

利用FLUENT軟件包中的GAMBIT模塊生成算例風場網(wǎng)格［8-10］。在GAMBIT軟件中建立高層建筑物的幾何模型，然后在已生成的模型中采用非結(jié)構(gòu)的四面體網(wǎng)格劃分了均勻網(wǎng)格，劃分網(wǎng)格量大約為130萬個。如圖2所示。

圖1　流場幾何模型

圖2　劃分的網(wǎng)格

1.2.2邊界條件

對入流處（inlet）的邊界條件設(shè)定，采用FLUENT中的速度入口邊界條件（velocity inlet），該邊界條件適用于不可壓縮流動；流體在出口處認為是充分發(fā)展的，所以出口邊界定義為自由出口（outflow）。outflow邊界條件用于出流邊界上的壓力或速度都未知的情況，適用于出口處流動時完全發(fā)展的情況；地面和建筑物表面都是實際存在而又固定不動的邊界，在此都采用無滑移的壁面條件（wall）；流域兩側(cè)及頂部定義為對稱邊界條件（symmetry）。

1.2.3控制方程

計算風工程中，鈍體繞流問題的控制方程是粘性不可壓Navier-Stokes方程，考慮到鈍體流場的特性及計算的收斂性和計算效率，本文選用“標準” k-ε二方程模型作為湍流模型。

其中，Gk是由于平均速度梯度引起的湍流動能k的產(chǎn)生項：Gb是由于浮力引起的湍流動能k的產(chǎn)生項：。

2　單體建筑高度對周圍風環(huán)境影響的模擬分析

2.1建筑模型工況設(shè)置

本節(jié)研究了單體建筑的高度變化對高層建筑物室外風場的影響，在本例的數(shù)值模擬計算中，建筑物模型的幾何尺寸高度分別取30m、60m、90m 和120m的四種情況，如表1所示。

表1　建筑模型工況設(shè)置

2.2模擬結(jié)果與對比分析

圖3是高度不同（表1工況）的建筑物建筑中心平面內(nèi)的風速等值線圖。從圖中可以看到，隨著建筑高度增加，建筑背后風影區(qū)高度增加。高度為30m的建筑，在計算區(qū)域內(nèi)，最大風速約為3.50m/s，最大風速比為1.38；高度為60m的建筑，在計算區(qū)域內(nèi)的最大風速約為4.0m/s，最大風速比為1.60；高度為90m的建筑，在計算區(qū)域內(nèi)的最大風速為4.27m/s，最大風速比為1.71。高度為120m的建筑，在計算區(qū)域內(nèi)，最大風速約為4.66m/s，最大風速比為 1.86，可見隨著建筑物高度的增加，建筑周圍的最大風速、最大風速比在增大，所以當建筑的高度增加時，建筑物對周圍空氣流動的影響增強。

圖3　Y=0m豎直剖面風速等值線圖

圖4是高度不同（表1工況）的建筑物的豎直平面流線圖，從圖中可以看到，隨著建筑高度的增加，建筑背后的漩渦隨建筑物的高度的增加而升高，并且漩渦越來越不明顯。

圖4　Y=0m豎直剖面風速流線圖

3　集群并行計算效率分析

3.1不同網(wǎng)格數(shù)下的并行計算效率比較分析

為了比較不同的網(wǎng)格數(shù)量對集群計算效率的影響，在CPU核數(shù)相同的情況下，網(wǎng)格數(shù)量方面選取10萬個、50萬個、100萬個、150萬個和200萬個5種工況，通過PBS作業(yè)提交系統(tǒng)，向集群提交并行計算作業(yè)。分別得到了這5組工況下在集群上計算所需要的總時間。通過計算得到了它們各自的并行效率。具體的計算結(jié)果如表2所示。

表2　CPU個數(shù)為12時集群并行計算工況

圖5是計算執(zhí)行時間與并行效率和網(wǎng)格數(shù)量之間關(guān)系的曲線圖，從圖中可以看出在CPU核數(shù)相同的條件下，任務的計算執(zhí)行時間隨任務的網(wǎng)

圖5　計算時間和并行效率與網(wǎng)格數(shù)量關(guān)系圖

格數(shù)量的增多而變長；集群效率隨任務網(wǎng)格數(shù)量的增多而變大，但隨網(wǎng)格數(shù)量的不斷增大，集群效率增大的趨勢越來越緩。

3.2不同計算節(jié)點數(shù)下的并行計算效率比較分析

為了比較CPU個數(shù)對計算效率影響的分析，設(shè)定工況B和工況C中的網(wǎng)格數(shù)分別都為80萬個和100萬個，8組不同的CPU個數(shù)工況；通過向集群提交作業(yè)得到了它們并行計算的結(jié)果，具體計算結(jié)果分別如表3和表4所示。

表3　網(wǎng)格數(shù)為80萬個時集群并行計算工況

表4　網(wǎng)格數(shù)為100萬時集群并行計算工況

結(jié)合表3和表4的數(shù)值結(jié)果，將工況B和工況C中計算執(zhí)行時間和CPU核數(shù)在同一坐標軸中繪圖描點得到了計算運行時間和CPU核數(shù)量之間關(guān)系的曲線圖，如圖6所示。從圖中可以看出任務在80萬網(wǎng)格和100萬網(wǎng)格時，它們運行計算時間隨CPU核數(shù)的增多具有相同的變化趨勢。很明顯的看出當網(wǎng)格數(shù)量不變的時候，任務的計算執(zhí)行時間隨CPU核數(shù)的增多先是快速減小而后有緩慢增多的趨勢。所以說一個任務不是所分配的CPU核數(shù)越多，它的計算執(zhí)行時間就越短。

圖6　執(zhí)行時間和CPU個數(shù)之間關(guān)系圖

圖7　并行效率和CPU個數(shù)之間的關(guān)系圖

結(jié)合表3和表4的數(shù)值結(jié)果，將工況B和工況C中的并行效率和CPU核數(shù)在同一坐標軸中繪圖描點得到了并行效率和CPU核數(shù)量之間關(guān)系的曲線圖，如圖7所示。從圖中可以看出任務在80萬網(wǎng)格和100萬網(wǎng)格時，它們的并行效率隨CPU個數(shù)的增多而下降，并且具有相同程度的下降趨勢。

4　結(jié)論

本文利用FLUENT軟件和高性能集群對高層建筑風場大規(guī)模并行數(shù)值模擬，并對計算結(jié)果的分析，總結(jié)如下：

（1）在相同風速下，建筑的長度和寬度相同時，隨著建筑高度的增加，建筑對其周圍氣流的影響增強，建筑周圍最大風速和最大風速比增大，建筑背后風影區(qū)面積增大。

（2）增加網(wǎng)格數(shù)量會增加計算時間，但是可以提高集群計算效率。

（3）隨著結(jié)點數(shù)目CPU的增加，總執(zhí)行時間先是逐漸快速減少，然后慢慢增加；而并行效率呈下降趨勢。

［1］胡曉峰，周孝清，卜增文等.基于室外風環(huán)境CFD模擬的建筑規(guī)劃設(shè)計［J］.工程建設(shè)與設(shè)計，2007（04）：14-18.

［2］龔晨.住宅小區(qū)人行高度風環(huán)境數(shù)值模擬研究［D］：廣東工業(yè)大學，2012.

［3］周允基，郜冶.北京某國際金融大廈風荷載與火災流場的數(shù)值模擬［J］.空氣動力學學報，2002（03）：351-360.

［4］梁向麗，符永正，周傳輝.空氣動力陰影區(qū)對高層建筑周圍污染物濃度分布的影響［J］.工業(yè)安全與環(huán)保，2003（12）：21-23.

［5］郜冶，周允基.風對高層建筑之間煙霧和污染物質(zhì)運動及散布的影響［J］.哈爾濱工程大學學報，2002（01）：113-119.

［6］劉申會.高層建筑的風荷載特性數(shù)值模擬研究［D］.湖南大學，2013.

［7］Murakami S.Current status and future trends in computational wind engineering［J］.Journal of Wind Engineering&Industrial Aerodynamics，1997：67.

［8］李進良，李承曦，胡仁喜.精通FLUENT6.3流場分析［M］.北京：化學工業(yè)出版社，2009.

［9］王瑞金，張凱，王剛.FLUENT技術(shù)基礎(chǔ)與應用實例［M］.北京：清華大學出版社，2007.

［10］溫正，石良辰，任毅如.FLUENT流體計算應用教程［M］.北京：清華大學出版社，2009.

（責任編輯：張洪瑋）

Numerical Simulation on the High-rise Building Wind Field Based on FLUENT and Cluster

JIANG Guang-biao，YANG Jie，ZHONG Hai，YIN Shui-pin，LI Qi-yong
（Central South University of Forestry and Technology，Changsha Hunan 410004）

In this paper，the actual high building is simulated based on FLUENT software and High Performance Cluster platform，and parallel numerical simulations were done by changing the height of high-rise buildings.The parallel numerical simulation results and efficiency of parallel computing clusters were analyzed，which we can get the following conclusions：With the increase of building height，the speed of wind and the speed of wind ratio increased at maximum，and effect of air flow become stronger around the building.When the number of CPU is the same，with the increase in the number grid，efficiency show an increasing tendency，but the increase of the magnitude become smaller；When the grid number is the same，with the increase of CPU number the efficiency of parallel computing show decreasing tendency in parallel computing.

FLUENT；high-rise building；High Performance Cluster；numerical simulation；parallel computing

TU974

A

10.3969/j.issn.1672-7304.2015.03.018

1672-7304（2015）03-0045-04

蔣光彪（1975-），男，湖南永州人，副教授，博士，研究方向：計算流體力學。

湖南省自科基金項目（14JJ2013），中南林業(yè)科技大學人才引進項目（104-0096）。