朱海天， 郝文星， 李 春， 丁勤衛(wèi)， 余 萬(wàn)

(上海理工大學(xué) 能源與動(dòng)力工程學(xué)院， 上海 200093)

目前，環(huán)境污染與能源緊缺問(wèn)題日益嚴(yán)峻，風(fēng)能是一種安全且無(wú)污染的可再生能源[1]，其開(kāi)發(fā)利用受到世界各國(guó)的廣泛關(guān)注。隨著風(fēng)能利用技術(shù)的快速發(fā)展，風(fēng)電迅速發(fā)展成為世界電能產(chǎn)業(yè)愈發(fā)重要的組成部分[2-3]。風(fēng)力機(jī)是一種將風(fēng)的動(dòng)能轉(zhuǎn)化為電能的旋轉(zhuǎn)機(jī)械。目前，風(fēng)力機(jī)組主要安裝于遠(yuǎn)離城市負(fù)荷中心的風(fēng)資源集中區(qū)域，由于當(dāng)?shù)仉娋W(wǎng)接納能力、跨區(qū)輸電能力以及應(yīng)對(duì)風(fēng)速隨機(jī)性和間歇性能力的不足，棄電限電現(xiàn)象較為嚴(yán)重[4]。2014年，中國(guó)平均風(fēng)電棄風(fēng)率約為8%，2015年則高達(dá)18.6%[5]。

為減緩棄電限電現(xiàn)象的發(fā)生，Keirstead等[6]指出，近年來(lái)各種與建筑物結(jié)合的部分風(fēng)力機(jī)被應(yīng)用于城市中，包括適用于多場(chǎng)合的微型風(fēng)力機(jī)、安裝于樓頂?shù)娘L(fēng)力機(jī)以及垂直軸阻力型風(fēng)力機(jī)等，這種風(fēng)能利用形式避免了輸、送電能問(wèn)題，產(chǎn)生的電能可直接用于建筑本身，為解決城市能源緊缺問(wèn)題和棄電限電問(wèn)題提供了技術(shù)支持。因此，充分利用建筑環(huán)境中的風(fēng)能為風(fēng)電發(fā)展提供了一種新思路。較之郊區(qū)和偏遠(yuǎn)地區(qū)，城市來(lái)流風(fēng)具有平均風(fēng)速低和湍流度大等特點(diǎn)，由于建筑物的影響，也存在局部大風(fēng)區(qū)域[7]。

風(fēng)力機(jī)與建筑物結(jié)合的形式基本可分為3類：安裝于建筑頂部的風(fēng)力機(jī)、安裝于建筑物中央的風(fēng)力機(jī)及建筑增強(qiáng)型風(fēng)力機(jī)(BAWT)[8]。1930年， Honnef首先提出了建筑增強(qiáng)型風(fēng)力機(jī)概念設(shè)計(jì)。建筑增強(qiáng)型風(fēng)力機(jī)將風(fēng)力機(jī)與建筑相結(jié)合，利用建筑對(duì)風(fēng)能的強(qiáng)化和集結(jié)作用，從而提高發(fā)電效率[9]。2008年底，巴林世界貿(mào)易中心完工[10]，在雙子塔建筑之間加裝了3個(gè)直徑為29 m的建筑增強(qiáng)型水平軸風(fēng)力機(jī)，這是世界上第一座大型結(jié)合風(fēng)力渦輪的建筑。Li等[11]對(duì)珠江大廈空洞中安裝垂直軸風(fēng)力機(jī)(VAWT)的效益進(jìn)行實(shí)驗(yàn)評(píng)估，結(jié)果表明喇叭口型洞口將增強(qiáng)風(fēng)速，受風(fēng)向和周圍建筑物的影響較為明顯。Heo等[12]對(duì)建筑增強(qiáng)型水平軸風(fēng)力機(jī)進(jìn)行了數(shù)值計(jì)算，結(jié)果表明建筑擴(kuò)散體可大幅提升發(fā)電效率，但該種BAWT受風(fēng)多向性差，當(dāng)偏航角小于30°時(shí)才具有較高的電能輸出。

上述研究一部分利用大氣邊界層風(fēng)剪切效應(yīng)，充分利用高層建筑的高度優(yōu)勢(shì)捕獲更多的風(fēng)能，另一部分利用建筑周圍的局部大風(fēng)。其中風(fēng)剪切是指風(fēng)速隨高度變化[13]，主要受地面粗糙度和科利奧效應(yīng)影響。筆者提出一種風(fēng)力機(jī)新型安裝方式與安裝位置，將VAWT懸掛固定于建筑擴(kuò)散體間的懸臂梁下端，利用上述2種優(yōu)勢(shì)對(duì)建筑環(huán)境中VAWT進(jìn)行研究。

直線翼垂直軸風(fēng)力機(jī)(SB-VAWT)是最具代表性[14]的垂直軸風(fēng)力機(jī)之一。SB-VAWT在旋轉(zhuǎn)時(shí)具有極復(fù)雜的氣動(dòng)特征[15-16]：一方面，攻角周期性大幅變化導(dǎo)致SB-VAWT氣動(dòng)載荷波動(dòng)較大且直葉片周圍流場(chǎng)結(jié)構(gòu)非常復(fù)雜；另一方面，單制動(dòng)槳盤(pán)假設(shè)已不適用于SB-VAWT，應(yīng)采用雙制動(dòng)槳盤(pán)假設(shè)，因此下游尾跡區(qū)易受到上風(fēng)區(qū)脫落渦軌跡的影響，而且葉片尾跡、主軸尾跡與及葉片支撐桿尾跡之間也發(fā)生相互作用，從而加劇流動(dòng)的復(fù)雜性。因此，研究建筑擴(kuò)散體間SB-VAWT的氣動(dòng)性能十分必要。

筆者采用計(jì)算流體力學(xué)(CFD)方法對(duì)原始SB-VAWT進(jìn)行數(shù)值研究，并與實(shí)驗(yàn)值進(jìn)行對(duì)比，驗(yàn)證了該網(wǎng)格拓?fù)浣Y(jié)構(gòu)及Realizablek-ε湍流模型的適用性。以上海風(fēng)資源特點(diǎn)作為研究條件，研究風(fēng)向?qū)ㄖ鰪?qiáng)型直線翼垂直軸風(fēng)力機(jī)(BASB-VAWT)性能與流場(chǎng)結(jié)構(gòu)的影響，得到其氣動(dòng)特性，以期為不同城市中采用BASB-VAWT提供依據(jù)。

1 氣動(dòng)理論

1.1 建筑環(huán)境中風(fēng)能利用特點(diǎn)

圓形截面建筑擴(kuò)散體周圍風(fēng)場(chǎng)CFD模擬結(jié)果如圖1所示。由于“夾道效應(yīng)”，2座建筑物間的風(fēng)速增大，形成局部強(qiáng)化風(fēng)[13]。

圖1 圓柱體建筑物之間的夾道效應(yīng)

1.2 風(fēng)剪切

風(fēng)剪切效應(yīng)是指風(fēng)速隨海拔高度增加而增大，要利用高層建筑的高度優(yōu)勢(shì)就必須考慮城市中的風(fēng)剪切效應(yīng)，筆者采用以下指數(shù)模型：

V/V0=(H/H0)β

(1)

式中：H為地面高度；H0為地面參考高度；V為H處的風(fēng)速；V0為H0處某參考點(diǎn)的風(fēng)速；β為切變系數(shù)，該系數(shù)與大氣穩(wěn)定度和地面粗糙度有關(guān)，其值如表1所示[17]。

表1 中國(guó)地面粗糙度類別與對(duì)應(yīng)的β值

采用地面粗糙度類別為D的經(jīng)驗(yàn)指數(shù)模型來(lái)描述上海市市區(qū)平均風(fēng)速沿高度的變化規(guī)律。上海市年平均風(fēng)速為3.65 m/s，根據(jù)風(fēng)剪切效應(yīng)可知，200 m高空處的城市來(lái)流風(fēng)速為9 m/s。

1.3 主要?dú)鈩?dòng)參數(shù)

在SB-VAWT中赤道平面攻角α隨相位角的周期性變化關(guān)系如下：

(2)

式中：θ為相位角；λ為葉尖速比。

葉尖速比定義為葉尖切向速度與來(lái)流風(fēng)速的比值：

(3)

式中：ω為角速度；V∞為來(lái)流風(fēng)速；R為旋轉(zhuǎn)半徑。

力矩系數(shù)Cm與風(fēng)能利用系數(shù)Cp均是衡量VAWT性能的重要指標(biāo)。

(4)

式中：T為轉(zhuǎn)矩；P為風(fēng)力機(jī)功率；ρ為空氣密度；A為掃風(fēng)面積。

實(shí)度σ作為描述SB-VAWT幾何形狀的重要無(wú)因次參數(shù)，其表達(dá)式[1]如下：

σ=Nc/(2R)

(5)

式中：N為葉片數(shù)；c為翼型弦長(zhǎng)。

2 計(jì)算模型

2.1 幾何模型

采用對(duì)低雷諾數(shù)工況有較高風(fēng)能利用系數(shù)的NACA0021對(duì)稱翼型[18]，其弦長(zhǎng)c為85.8 mm，最大厚度為0.21c，最大厚度位置在0.3c處，翼型幾何模型見(jiàn)圖2。

圖2 翼型幾何模型

為與實(shí)驗(yàn)值進(jìn)行比較，采用的幾何模型參數(shù)均與文獻(xiàn)[19]相同。SB-VAWT幾何參數(shù)見(jiàn)表2。

巴林世界貿(mào)易中心的建筑增強(qiáng)型水平軸風(fēng)力機(jī)實(shí)圖與本文BASB-VAWT幾何模型的對(duì)比如圖3所示。其中建筑增強(qiáng)型水平軸風(fēng)力機(jī)放置于懸臂梁前端，而B(niǎo)ASB-VAWT垂直懸掛于懸臂梁下端。

2.2 計(jì)算域結(jié)構(gòu)及網(wǎng)格分布

SB-VAWT與翼型周圍計(jì)算域結(jié)構(gòu)及網(wǎng)格分布如圖4所示。SB-VAWT遠(yuǎn)場(chǎng)計(jì)算域結(jié)構(gòu)由半圓形上游區(qū)與矩形下游尾跡區(qū)組成，如圖4(a)中上游區(qū)特征長(zhǎng)度為50R，下游尾跡區(qū)特征長(zhǎng)度為30R。SB-VAWT近場(chǎng)計(jì)算域結(jié)構(gòu)由Z1和Z22個(gè)網(wǎng)格域構(gòu)成，Z1為靜止結(jié)構(gòu)網(wǎng)格域，Z2為滑移非結(jié)構(gòu)網(wǎng)格域。翼型周圍計(jì)算域結(jié)構(gòu)也由半圓形上游區(qū)與矩形下游尾跡區(qū)組成，如圖4(b)與圖4(c)所示。上游區(qū)特征長(zhǎng)度為0.5R，下游尾跡區(qū)特征長(zhǎng)度為0.3R，翼型周圍網(wǎng)格域共分為6個(gè)域,如圖4(d)所示。域1和域6分別為翼型前緣域和翼型后緣域。定義SB-VAWT計(jì)算域中的半圓和矩形上下邊界ad與bc速度入口邊界條件，來(lái)流風(fēng)速為9 m/s，矩形右邊界cd為壓力出口條件。

表2 SB-VAWT幾何參數(shù)

(a) 巴林世界貿(mào)易中心

(b) BASB-VAWT實(shí)體圖

(c) BASB-VAWT側(cè)視圖

(d) BASB-VAWT俯視圖

(a) 遠(yuǎn)場(chǎng)拓?fù)浣Y(jié)構(gòu)示意圖(b) 二維簡(jiǎn)化圖

(c) 遠(yuǎn)場(chǎng)計(jì)算域網(wǎng)格分布(d) 翼型周圍計(jì)算域結(jié)構(gòu)

(e) 近場(chǎng)計(jì)算域網(wǎng)格分布(f) 翼型周圍網(wǎng)格分布

圖4 計(jì)算域結(jié)構(gòu)及網(wǎng)格分布

Fig.4 Calculation region and mesh distribution of the airfoil

SB-VAWT遠(yuǎn)場(chǎng)計(jì)算域結(jié)構(gòu)網(wǎng)格分布如圖4(c)所示，Z1和Z2近場(chǎng)計(jì)算域混合網(wǎng)格分布如圖4(e)所示。翼型表面布置250個(gè)網(wǎng)格節(jié)點(diǎn)，并對(duì)翼型周圍進(jìn)行加密處理，由翼型表面到遠(yuǎn)場(chǎng)邊界方向布置70個(gè)節(jié)點(diǎn)，如圖4(f)所示?？偩W(wǎng)格數(shù)為421 092。

采用Realizablek-ε湍流模型，在模擬旋轉(zhuǎn)流動(dòng)方面該湍流模型優(yōu)于Standardk-ε湍流模型，也能較精確地捕捉到流動(dòng)分離點(diǎn)[20]。y+值接近30～100才能使壁面函數(shù)較好地捕捉到邊界層的流動(dòng)[21]。故本文網(wǎng)格模型選取的y+在30～50內(nèi)。

2.3 網(wǎng)格無(wú)關(guān)性驗(yàn)證

BASB-VAWT網(wǎng)格無(wú)關(guān)性驗(yàn)證結(jié)果見(jiàn)表3。以大網(wǎng)格數(shù)約66萬(wàn)計(jì)算結(jié)果為基準(zhǔn)，結(jié)果表明：在約42萬(wàn)總網(wǎng)格數(shù)時(shí)的計(jì)算值與約66萬(wàn)總網(wǎng)格數(shù)時(shí)的計(jì)算值偏差較小，故本網(wǎng)格模型是可信的。因此，選取約為42萬(wàn)網(wǎng)格數(shù)的模型作為最終計(jì)算模型。

表3 網(wǎng)格無(wú)關(guān)性驗(yàn)證結(jié)果

2.4 實(shí)驗(yàn)驗(yàn)證

原始SB-VAWT二維風(fēng)能利用系數(shù)Cp計(jì)算值與實(shí)驗(yàn)值[19]的對(duì)比如圖5所示。由圖5可知，在低葉尖速比時(shí)，風(fēng)能利用系數(shù)計(jì)算值與實(shí)驗(yàn)值的差值略大，其主要原因有：(1) 文獻(xiàn)[19]中的實(shí)驗(yàn)值未考慮風(fēng)洞阻塞效應(yīng)且未排除轉(zhuǎn)軸與力矩傳感器間的摩擦阻力力矩；(2) 筆者僅考慮二維SB-VAWT，即設(shè)為無(wú)限長(zhǎng)葉片，然而Li等[22]指出，二維SB-VAWT計(jì)算值相比三維SB-VAWT計(jì)算值與實(shí)驗(yàn)值具有較大差值，這是由于二維SB-VAWT葉片未考慮沿展向葉素表面壓差逐漸降低的因素；(3) 二維SB-VAWT未考慮葉片支撐桿對(duì)SB-VAWT的影響。

圖5 二維風(fēng)能利用系數(shù)計(jì)算值與實(shí)驗(yàn)值的對(duì)比圖

結(jié)合上述3種降低風(fēng)能利用系數(shù)的實(shí)際因素，故本文二維CFD計(jì)算值與實(shí)驗(yàn)值存在誤差，但本文CFD計(jì)算值仍可捕捉到最大風(fēng)能利用系數(shù)對(duì)應(yīng)的最佳葉尖速比，其值為2.62。因此，二維CFD計(jì)算值仍可定性且定量地反映風(fēng)向?qū)Σ煌珺ASB-VAWT的影響。

2.5 風(fēng)玫瑰圖

以上海市近3年天氣報(bào)告作為風(fēng)向研究條件，如圖6所示共8個(gè)風(fēng)向，包括北風(fēng)(N)、東北風(fēng)(NE)、東風(fēng)(E)、東南風(fēng)(SE)、南風(fēng)(S)、西南風(fēng)(SW)、西風(fēng)(W)及西北風(fēng)(NW)。其中，縱軸代表出現(xiàn)此風(fēng)向的頻率。

(a) 2014年

(b) 2015年

(c) 2016年

由圖6可知，上海市東南風(fēng)、東風(fēng)及東北風(fēng)工況頻數(shù)較多，近3年平均風(fēng)速為3.65 m/s。

2.6 建筑擴(kuò)散體型

將建筑擴(kuò)散體截面分為4種：正圓形、橢圓形、矩形及菱形。幾何模型如圖7所示，由于4種建筑擴(kuò)散體均為中心點(diǎn)對(duì)稱模型，故僅需研究3種不同風(fēng)向下的工況：與建筑擴(kuò)散體對(duì)稱軸平行的風(fēng)向、與建筑擴(kuò)散體對(duì)稱軸之間角度呈45°的風(fēng)向以及與建筑擴(kuò)散體對(duì)稱軸垂直的風(fēng)向，因而工作量大大降低。因需要引入氣象資料，故將這3種風(fēng)向分別稱為西風(fēng)、西北風(fēng)及北風(fēng)，如圖7所示的方向位置。

正圓形擴(kuò)散體SB-VAWT周圍局部流場(chǎng)網(wǎng)格分布如圖8所示。其他BASB-VAWT的周圍流場(chǎng)網(wǎng)格分布不再贅述。

(a) 正圓形擴(kuò)散體(b) 橢圓形擴(kuò)散體

(c) 矩形擴(kuò)散體(d) 菱形擴(kuò)散體

圖7 4種建筑擴(kuò)散體

Fig.7 Four models of the building diffuser

圖8 正圓形擴(kuò)散體BASB-VAWT網(wǎng)格分布局部放大圖

3 結(jié)果與分析

3.1 BASB-VAWT周圍流場(chǎng)分析

圖9和圖10為各風(fēng)向下建筑擴(kuò)散體的壓力云圖、流線圖和BASB-VAWT速度云圖，其中來(lái)流風(fēng)速為9 m/s。

由圖9可知，建筑前后壓差造成局部大風(fēng)區(qū)域，建筑外廓使一部分流體集聚于兩建筑之間。在西北風(fēng)工況下，矩形擴(kuò)散體間建筑產(chǎn)生的大流動(dòng)分離發(fā)展至大部分流道中，而其他截面擴(kuò)散體間建筑產(chǎn)生的渦流并沒(méi)有大幅發(fā)展至流道中。在北風(fēng)工況下，流場(chǎng)結(jié)構(gòu)極為復(fù)雜，建筑背風(fēng)面存在大分離渦。

根據(jù)最佳捕獲風(fēng)能的位置安裝BASB-VAWT。在西風(fēng)工況下，建筑擴(kuò)散體間的風(fēng)速明顯大于入口風(fēng)速，風(fēng)力機(jī)可捕獲更多的風(fēng)能。在西北風(fēng)工況下，由于建筑擴(kuò)散體產(chǎn)生的流動(dòng)分離部分發(fā)展至風(fēng)力機(jī)旋轉(zhuǎn)區(qū)域內(nèi)，對(duì)風(fēng)力機(jī)產(chǎn)生一定影響。在北風(fēng)工況下，由于建筑擴(kuò)散體間的風(fēng)速較低，且建筑背風(fēng)面存在大分離渦，風(fēng)力機(jī)將承受更劇烈的載荷波動(dòng)。

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

3.2 BASB-VAWT氣動(dòng)性能分析

考慮BASB-VAWT在不同建筑擴(kuò)散體下不同風(fēng)向工況時(shí)的性能，對(duì)西風(fēng)、西北風(fēng)及北風(fēng)工況時(shí)BASB-VAWT一個(gè)旋轉(zhuǎn)周期內(nèi)的周圍流場(chǎng)進(jìn)行數(shù)值模擬，其λ為2.62。BASB-VAWT風(fēng)能利用系數(shù)計(jì)算結(jié)果如圖11所示。

由圖11可知，在西風(fēng)工況下，正圓形BASB-VAWT風(fēng)能利用系數(shù)較高，其值為0.783 30，風(fēng)能利用系數(shù)提升至原始SB-VAWT的2.47倍；橢圓形BASB-VAWT風(fēng)能利用系數(shù)較低，為0.678 73；在西北風(fēng)工況下，正圓形BASB-VAWT風(fēng)能利用系數(shù)較高，其值為0.389 20；菱形BASB-VAWT風(fēng)能利用系數(shù)較低，為0.286 96；在北風(fēng)工況下，由于建筑擴(kuò)散體之間風(fēng)速過(guò)低，風(fēng)力機(jī)效率極低。

為進(jìn)一步得到風(fēng)向?qū)ASB-VAWT年風(fēng)能利用系數(shù)的影響，結(jié)合上海市近3年風(fēng)玫瑰圖，通過(guò)加權(quán)平均得出BASB-VAWT年風(fēng)能利用系數(shù)，如圖12所示。其中，在北風(fēng)工況下，4種BASB-VAWT年風(fēng)能利用系數(shù)均較低，且葉片周圍流場(chǎng)極為復(fù)雜，故年風(fēng)能利用系數(shù)并未考慮北風(fēng)工況。同理，在西北風(fēng)工況下，矩形擴(kuò)散體自身的強(qiáng)渦流直接影響風(fēng)力機(jī)的葉片，故在研究年風(fēng)能利用系數(shù)中并未考慮此工況。

由圖12可知，在上海市市區(qū)，正圓形BASB-VAWT具有較高的年風(fēng)能利用系數(shù)，近3年其均值為0.385 92；而矩形BASB-VAWT年風(fēng)能利用系數(shù)較低，近3年其均值為0.140 48。

為進(jìn)一步研究BASB-VAWT在西風(fēng)及西北風(fēng)工況下的載荷波動(dòng)，對(duì)4種BASB-VAWT進(jìn)行了對(duì)比分析，由于北風(fēng)工況下的流場(chǎng)結(jié)構(gòu)復(fù)雜，故對(duì)其力矩系數(shù)進(jìn)行研究無(wú)意義。BASB-VAWT力矩系數(shù)與相位角θ的關(guān)系見(jiàn)圖13。

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

西風(fēng)西北風(fēng)

北風(fēng)

圖11 BASB-VAWT風(fēng)能利用系數(shù)

由圖13可知，在西風(fēng)工況下，BASB-VAWT將承受比原始SB-VAWT更大幅的載荷波動(dòng)。在西北風(fēng)工況下，矩形BASB-VAWT的載荷波動(dòng)較原始SB-VAWT與其他擴(kuò)散體BASB-VAWT的載荷波動(dòng)更為劇烈，正圓形、橢圓形BASB-VAWT與原始SB-VAWT曲線相近。而菱形BASB-VAWT的性能較差，這是由于在菱形擴(kuò)散體間形成局部大風(fēng)區(qū)域較為滯后。

圖12 BASB-VAWT年風(fēng)能利用系數(shù)

4 結(jié) 論

(1) BASB-VAWT氣動(dòng)性能對(duì)風(fēng)向敏感，在西風(fēng)工況下，其性能大幅提升；在西北風(fēng)工況下，性能稍高于原始SB-VAWT；在北風(fēng)工況下，流場(chǎng)結(jié)構(gòu)非常復(fù)雜且建筑產(chǎn)生的脫落渦發(fā)展至風(fēng)力機(jī)旋轉(zhuǎn)半徑內(nèi)，其性能遠(yuǎn)低于原始SB-VAWT。

(a) 西風(fēng)工況下

(b) 西北風(fēng)工況下

(2) BASB-VAWT氣動(dòng)性能對(duì)建筑擴(kuò)散體排布及建筑平面外廓敏感，較之橢圓形和正圓形建筑，矩形建筑的風(fēng)力機(jī)性能較劣，而正圓形建筑優(yōu)于長(zhǎng)半軸平行于風(fēng)向的橢圓形建筑。

參考文獻(xiàn)：

[1] 李春, 葉舟, 高偉, 等. 現(xiàn)代大型風(fēng)力機(jī)設(shè)計(jì)原理[M]. 上海: 上?？茖W(xué)技術(shù)出版社, 2013.

[2] 王東華, 葉舟, 郝文星, 等. 海上漂浮式風(fēng)力機(jī)Spar平臺(tái)波頻與慢漂響應(yīng)性能分析[J].動(dòng)力工程學(xué)報(bào), 2016, 36(11): 907-913.

WANG Donghua, YE Zhou, HAO Wenxing, et al. Analysis on wave frequency and slow drift response of a Spar platform for offshore floating wind turbines[J].JournalofChineseSocietyofPowerEngineering, 2016, 36(11): 907-913.

[3] SCHEURICH F, FLETCHER T M, BROWN R E. Simulating the aerodynamic performance and wake dynamics of a vertical-axis wind turbine[J].WindEnergy, 2011, 14(2): 159-177.

[4] 國(guó)家能源局. 關(guān)于做好2013年風(fēng)電并網(wǎng)和消納相關(guān)工作的通知[R]. 北京：國(guó)家能源局，2013.

[5] 李俊峰, 蔡豐波, 喬黎明, 等. 2015年中國(guó)風(fēng)電發(fā)展報(bào)告[R]. 北京: 中國(guó)循環(huán)經(jīng)濟(jì)協(xié)會(huì)可再生能源專業(yè)委員會(huì), 2015.

[6] KEIRSTEAD J, JENNINGS M, SIVAKUMAR A. A review of urban energy system models: approaches, challenges and opportunities[J].RenewableandSustainableEnergyReviews, 2012, 16(6): 3847-3866.

[7] 姜瑜君, 桑建國(guó), 張伯寅. 高層建筑的風(fēng)環(huán)境評(píng)估[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 42(1): 68-73.

JIANG Yujun, SANG Jianguo, ZHANG Boyin. An evaluation on wind environment around skyscrapers[J].ActaScientiarumNaturaliumUniversitatisPekinensis, 2006, 42(1): 68-73.

[8] MULLER G, JENTSCH M F, STODDART E. Vertical axis resistance type wind turbines for use in buildings[J].RenewableEnergy, 2009, 34: 1407-1412.

[9] HEYMANN M. Signs of hubris: the shaping of wind technology styles in Germany, Denmark, and the United States, 1940-1990[J].TechnologyandCulture, 1998, 39(4): 641-670.

[10] 袁行飛, 張玉. 建筑環(huán)境中的風(fēng)能利用研究進(jìn)展[J].自然資源學(xué)報(bào), 2011, 26(5): 891-898.

YUAN Xingfei, ZHANG Yu. Progress in the study of wind energy utilization in the built environment[J].JournalofNaturalResources, 2011, 26(5): 891-898.

[11] LI Q S, SHU Z R, CHEN F B. Performance assessment of tall building-integrated wind turbines for power generation[J].AppliedEnergy, 2016, 165: 777-788.

[12] HEO Y G, CHOI N J, CHOI K H, et al. CFD study on aerodynamic power output of a 110 kW building augmented wind turbine[J].EnergyandBuildings, 2016, 129: 162-173.

[13] GUALTIERI G. Atmospheric stability varying wind shear coefficients to improve wind resource extrapolation: a temporal analysis[J].RenewableEnergy, 2016, 87: 376-390.

[14] LI Qing'an, MAEDA T. Study on flow around straight-bladed vertical axis wind turbine under low tip speed ratio[J].JournalofFluidScienceandTechnology, 2014, 9(3): JFST0051.

[15] MOHAMED M H. Impacts of solidity and hybrid system in small wind turbines performance[J].Energy, 2013, 57: 495-504.

[16] XIAO Qing, LIU Wendi, INCECIK A. Flow control for VATT by fixed and oscillating flap[J].RenewableEnergy, 2013, 51: 141-152.

[17] 中華人民共和國(guó)建設(shè)部. 建筑結(jié)構(gòu)荷載規(guī)范(2006版): GB 50009—2001[S]. 北京: 中國(guó)建筑工業(yè)出版社, 2002.

[18] MOHAMED M H. Performance investigation of H-rotor Darrieus turbine with new airfoil shapes[J].Energy, 2012, 47: 522-530.

[19] CASTELLI M R, ARDIZZON G, BATTISTI L, et al. Modeling strategy and numerical validation for a Darrieus vertical axis micro-wind turbine[C]//Proceedingsof2010InternationalMechanicalEngineeringCongressandExposition. Vancouver, Canada: British Columbia, 2010.

[20] LI Qing'an, MAEDA T, KAMADA Y, et al. Wind tunnel and numerical study of a straight-bladed vertical axis wind turbine in three-dimensional analysis (part Ⅱ: for predicting flow field and performance)[J].Energy, 2016, 104: 395-307.

[21] STOLARSKI T, NAKASONE Y, YOSHIMOTO S. Engineering analysis with ANSYS software[M]. Amsterdam, Holland: Butterworth-Heinemann, 2007.

[22] LI Qing'an, MAEDA T, KAMADA Y, et al. Wind tunnel and numerical study of a straight-bladed vertical axis wind turbine in three-dimensional analysis (part I: for predicting aerodynamic loads and performance)[J].Energy, 2016, 106: 443-452.