李爭 程立源 高夢海 張文達(dá)

摘 要： 為了提高風(fēng)力機組的整體性能，解決風(fēng)力機在實際運行中受建筑物影響的問題，利用仿真分析軟件Fluent對不同葉片數(shù)的新型鸚鵡螺等角螺線型風(fēng)力機進行氣動性能研究，建立建筑物與風(fēng)力機組排布模型，分析建筑物擾流特性，對比擾流環(huán)境對風(fēng)力機組轉(zhuǎn)矩性能的影響。結(jié)果表明：3個葉片風(fēng)力機的整體性能更優(yōu);建筑物下游出現(xiàn)紊流區(qū)域，切向速度明顯增加，其附近的新型鸚鵡螺等角螺線型風(fēng)力機組轉(zhuǎn)矩性能明顯提升，驗證了建筑物附近安裝鸚鵡螺等角螺線型風(fēng)力機組的可行性。所提風(fēng)力機組排布方式可有效提升風(fēng)力機性能，為風(fēng)力機結(jié)構(gòu)優(yōu)化設(shè)計和建筑物附近風(fēng)力機排布提供參考。

關(guān)鍵詞： 空氣動力學(xué);垂直軸風(fēng)力機;Fluent;轉(zhuǎn)矩;建筑物

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

doi：10.7535/hbkd.2020yx06010

Numerical simulation of nautilus equiangular spiral wind

turbine and analysis of building turbulence

LI Zheng， CHENG Liyuan， GAO Menghai， ZHANG Wenda

（School of Electrical Engineering， Hebei University of Science and Technology， Shijiazhuang， Hebei 050018， China）

Abstract：

In order to improve the overall performance of wind turbines and solve the problem that the wind turbine is affected by the building in actual operation， the simulation analysis software Fluent was used to study the aerodynamic performance of the new Nautilus equiangular spiral wind turbine with different numbers of blades. The layout model of buildings and wind turbines were established， the turbulence characteristics of buildings were analyzed， and the torque performance of wind turbines in the turbulent environment were compared. The results show that the overall performance of the three-blade wind turbine is more superior; there is turbulence area at the lower position of the building， and the tangential speed is significantly increased; the performance of the new Nautilus equiangular spiral wind turbine near the building are significantly improved， which verifies the feasibility of installing the Nautilus equiangular spiral wind turbine near the building. The proposed arrangement of wind turbines can effectively improve the performance of wind turbines， and provide references for the optimal design of wind turbine structures and the arrangement of wind turbines near buildings.

Keywords：aerodynamics; vertical axis wind turbine; Fluent; torque; building

發(fā)量都居世界首位。風(fēng)能是一種可再生的清潔型能源，風(fēng)力發(fā)電可作為可再生能源的主力電源[1-4]。垂直軸風(fēng)力機具有結(jié)構(gòu)簡單、便于維護、占地面積小、噪聲低等優(yōu)勢，吸引了大量的學(xué)者進行研究，進一步促進了垂直軸風(fēng)力機的發(fā)展[5-8]。風(fēng)力機的性能受安裝位置的影響較大，建筑物下游的風(fēng)力機會受到建筑物對來流風(fēng)速大小、方向、湍流強度等的影響。如果能夠在建筑物附近找到風(fēng)力機合適的安裝位置，則可以避免低風(fēng)速區(qū)，利用局部大風(fēng)區(qū)域提高風(fēng)力機的性能。

風(fēng)力機的研究方法主要有實物實驗法、風(fēng)洞實驗法和仿真模擬法。大多數(shù)研究采用數(shù)值模擬仿真的研究方法對風(fēng)力發(fā)電機的空氣動力學(xué)性能進行分析，可以節(jié)約時間和研究成本。文獻[9]對不同水平面上的風(fēng)力機進行分析，兩臺風(fēng)力機相互促進，提升整個風(fēng)力機組的性能。文獻[10—11]分別從風(fēng)力機仿真過程中二維、三維的模型選取進行了綜合分析。二維仿真的計算時間短，仿真結(jié)果優(yōu)于實驗值;三維仿真可以更好地反映出實際情況，但是運算時間長。改變風(fēng)力機的結(jié)構(gòu)和運用新的計算方法可提升風(fēng)力機的性能以及仿真模擬的計算精度。文獻[12]將垂直軸風(fēng)力機與外軸風(fēng)力機相組合，最大程度上減少了功率波動和到達(dá)峰值的時間，提高了風(fēng)力機的效率。文獻[13]對三維建筑物周圍的流場進行了數(shù)值模擬，模擬結(jié)果與風(fēng)洞試驗結(jié)果相似。建筑群間的繞流也比單個建筑物要復(fù)雜得多，建筑物后方的流場存在復(fù)雜漩渦區(qū)，不利于空氣的擴散傳播，嚴(yán)重影響風(fēng)力機的正常運行。文獻[14]根據(jù)超高層建筑的集風(fēng)效應(yīng)，提出在高層建筑上尋找較高風(fēng)速和較低湍流強度的位置安裝風(fēng)力發(fā)電機可以有效降低棄風(fēng)率。這些研究對推動風(fēng)力機的發(fā)展非常有價值，但是關(guān)于哪種建筑物對應(yīng)哪種特定的風(fēng)力機模型的研究相對較少。

為了得到較好的風(fēng)力機模型，本文對2，3，4個葉片的風(fēng)力機轉(zhuǎn)矩特性和風(fēng)能利用率進行比較分析，確定風(fēng)力機的葉片數(shù)量，然后分析建筑物附近的繞流特性，并仿真模擬計算建筑物周圍4個特定位置風(fēng)力機組的轉(zhuǎn)矩，觀察風(fēng)力機性能的變化，以期為建筑物附近安裝新型鸚鵡螺等角螺線型風(fēng)力機組的可行性提供理論依據(jù)。

1 風(fēng)力機性能系數(shù)相關(guān)定義

根據(jù)牛頓第二定理[15-16]，

F=ma，

E=12mv2，

式中：F為風(fēng)力;m為空氣質(zhì)量;a為加速度;E為動能。

風(fēng)能能量為

P=12ρSv3，

式中：ρ為空氣密度;S為風(fēng)力機掃風(fēng)面積;v為來流風(fēng)速。

風(fēng)力機功率為

Pw=Tω，

式中：T為風(fēng)力機轉(zhuǎn)矩;ω為風(fēng)力機旋轉(zhuǎn)角速度。

風(fēng)能利用率Cp是衡量垂直軸風(fēng)力機的重要指標(biāo)，其計算式為

Cp=Tω12ρSv3。

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

TSR=Rωv，

式中R為風(fēng)力機半徑。

2 幾何模型及計算域的創(chuàng)建

在對新型鸚鵡螺等角螺線型風(fēng)力機的空氣動力學(xué)性能進行分析時，為了節(jié)約計算資源和減少計算時間，本文只對鸚鵡螺垂直軸風(fēng)力機的風(fēng)輪部分進行數(shù)值模擬。圖1為3種不同葉片數(shù)的風(fēng)輪模型。風(fēng)力機的入風(fēng)口截面為橢圓形并且3個葉片間的螺線角距離相等，風(fēng)力機入口處橢圓的長軸為0.33 m，短軸為0.15 m。鸚鵡螺垂直軸風(fēng)力機葉片的外形是特殊弧形結(jié)構(gòu)， 且風(fēng)輪的整體結(jié)構(gòu)相對圓滑。它的橫截面接近橢圓形，有效減少了風(fēng)力機在旋轉(zhuǎn)時的風(fēng)阻，使得風(fēng)力機可以獲得較大的旋轉(zhuǎn)速度。在風(fēng)輪材料強度允許的情況下，風(fēng)力機能夠以較高速度旋轉(zhuǎn)，產(chǎn)生更多電能。風(fēng)力機在旋轉(zhuǎn)時保證一直存在1個葉片處于迎風(fēng)狀態(tài)，其余葉片的另一側(cè)恰好處于背風(fēng)處。這種阻力型風(fēng)力機在運行過程中有效提升了轉(zhuǎn)矩差，產(chǎn)生更大旋轉(zhuǎn)轉(zhuǎn)矩，有效提高了風(fēng)力機的空氣動力學(xué)性能。 在進行繞流特性分析時選取規(guī)則的長方體墻體代替建筑物墻體。風(fēng)力機和建筑物的參數(shù)如表1所示。

為了盡可能模擬風(fēng)力機實際的運行環(huán)境，準(zhǔn)確計算風(fēng)力機的基本性能，首先構(gòu)建流體仿真中風(fēng)力機的計算域。整個計算域內(nèi)的介質(zhì)為均勻的空氣，計算域的尺寸遠(yuǎn)大于風(fēng)力機模型的大小，如圖2所示，風(fēng)力機的計算域尺寸為20 m×10 m×5 m。

2.1 網(wǎng)格劃分

網(wǎng)格劃分對計算結(jié)果的準(zhǔn)確性起著至關(guān)重要的作用。鸚鵡螺垂直軸風(fēng)力機葉片外形是特殊的曲面結(jié)構(gòu)，在網(wǎng)格劃分時采用六面體及四面體結(jié)構(gòu)會降低計算精度，存在劃分時間長等缺點。本文采用非結(jié)構(gòu)性的網(wǎng)格劃分方法，以有效降低網(wǎng)格的出錯率并提高網(wǎng)格劃分的效率[17-18]，最終確定網(wǎng)格數(shù)量為25萬個。網(wǎng)格劃分時增加了旋轉(zhuǎn)域和風(fēng)輪葉片的網(wǎng)格密度，該部分網(wǎng)格數(shù)量有14萬個，后處理的過程中可以清晰觀察出整個計算域中云圖分布特點。網(wǎng)格劃分如圖3所示。

2.2 計算方法

在實際求解過程中選擇合適的湍流模型進行仿真計算。仿真模擬時合適的湍流模型不僅可以增加計算效率，還可以將計算精度控制在合理的范圍內(nèi)，使計算結(jié)果更加與實際情況相符。與RNG k-ε模型相比，Standard k-ε模型適用范圍廣，計算量合適，有較多的數(shù)據(jù)積累和適當(dāng)?shù)挠嬎憔?，因此選擇Standard k-ε模型進行計算[19-20]。入口速度v為6 m/s，出口壓力為0 Pa，計算域4面為墻體，流體材質(zhì)為氣體。在進行垂直軸風(fēng)力機的數(shù)值仿真計算時，當(dāng)風(fēng)輪的切向力系數(shù)呈現(xiàn)出明顯的周期性變化時，視為計算結(jié)果收斂。

3 結(jié)果分析

3.1 風(fēng)力機分析

假設(shè)來流風(fēng)速固定為6 m/s，風(fēng)力機旋轉(zhuǎn)速度為5°/s，且2，3，4葉片的風(fēng)力機分別旋轉(zhuǎn)1周，繪制出風(fēng)力機的轉(zhuǎn)矩波形圖，如圖4所示，3種風(fēng)力機的風(fēng)能利用率如圖5所示。

從圖4可知，風(fēng)力機產(chǎn)生轉(zhuǎn)矩的大小并不隨著葉片數(shù)的增加而增大，風(fēng)力機轉(zhuǎn)矩的平穩(wěn)性也不與葉片數(shù)有明顯的線性關(guān)系。當(dāng)風(fēng)力機為2個葉片時，風(fēng)力機產(chǎn)生的轉(zhuǎn)矩波動最大，工作時至少有半個周期風(fēng)力機處于背風(fēng)狀態(tài)，轉(zhuǎn)矩值下降。風(fēng)力機為3個葉片時，風(fēng)力機旋轉(zhuǎn)1周產(chǎn)生的平均轉(zhuǎn)矩位于其他2種風(fēng)力機之上，3個葉片的風(fēng)力機的轉(zhuǎn)矩峰值也較大。風(fēng)力機為4個葉片時，轉(zhuǎn)矩波動與3個葉片的風(fēng)力機相比較小，轉(zhuǎn)矩峰值與轉(zhuǎn)矩平均值均小于3個葉片的風(fēng)力機。

從圖5可知，葉片數(shù)量為2，3，4時，風(fēng)能利用率的曲線逐漸向低尖速比的方向移動。葉片數(shù)為2時，風(fēng)力機運行的尖速比范圍與其他二者相比較大，但其提供的啟動力矩小，自啟動的性能較差;葉片數(shù)為3，4時，隨著葉片數(shù)量的增加，風(fēng)能利用系數(shù)對應(yīng)的尖速比范圍變窄，最大的風(fēng)能利用率Cp也變小。在葉片增加的過程中，最大的風(fēng)能利用率Cp對應(yīng)的尖速比也在減小。由于鸚鵡螺垂直軸風(fēng)力機葉片的特殊結(jié)構(gòu)，實度也會隨著葉片數(shù)量的增多而增加，在一定程度上增加風(fēng)輪的阻力。到達(dá)葉片上的風(fēng)速變小，進一步降低葉片產(chǎn)生的旋轉(zhuǎn)力矩，風(fēng)能利用率Cp下降，并且向低尖速比的方向移動;隨著葉片數(shù)的增加，低速旋轉(zhuǎn)時風(fēng)力機更容易產(chǎn)生較大的旋轉(zhuǎn)轉(zhuǎn)矩，也更容易啟動。綜合考慮風(fēng)輪的轉(zhuǎn)矩特性、壓力特征、在較寬范圍內(nèi)的尖速比下的風(fēng)能利用系數(shù)[21]以及結(jié)合材料和成本問題，認(rèn)為鸚鵡螺垂直軸風(fēng)力機的葉片數(shù)量為3時最佳。

3.2 建筑物擾流分析

城市環(huán)境中的建筑物嚴(yán)重影響風(fēng)速的流動方向，建筑物附近會產(chǎn)生部分的漩渦紊流區(qū)域，而建筑物周圍出現(xiàn)的切向風(fēng)速會嚴(yán)重影響風(fēng)力機的運行狀態(tài)和性能。本文利用Fluent對建筑物進行數(shù)值仿真模擬，通過觀察建筑物2個特殊位置的風(fēng)速矢量圖，分析建筑物對來流風(fēng)速的干擾范圍和速度變化波動，驗證建筑物附近安裝新型鸚鵡螺垂直軸風(fēng)力機的可行性，增加風(fēng)力機的經(jīng)濟效益。建筑物附近的風(fēng)速矢量圖如圖6所示。

根據(jù)圖6的速度矢量圖，可以觀察到建筑物對來流風(fēng)速在水平和垂直方向的影響非常大。當(dāng)上游的風(fēng)遇到建筑物時，在建筑物前緣出現(xiàn)流動分離現(xiàn)象。風(fēng)沿著建筑物周圍向下游擴散時，會出現(xiàn)環(huán)繞現(xiàn)象。建筑物兩側(cè)和頂端向下的風(fēng)速流線相互交織，一起向建筑物的下游傳播擴散。建筑物后方的湍流尾流變窄，風(fēng)向、風(fēng)速等變化十分紊亂，在下游流場形成漩渦。建筑物的正后方風(fēng)速會出現(xiàn)較大的漩渦，隨著距離的增加漩渦逐漸變小，另外該區(qū)域風(fēng)速低、風(fēng)速變化頻繁。在整個計算域內(nèi)最高風(fēng)速為9.408 m/s，是來流風(fēng)速的1.568倍。同時可以發(fā)現(xiàn)，風(fēng)在建筑物頂端和兩側(cè)都出現(xiàn)了明顯的分離擾流現(xiàn)象，因此建筑物邊緣附近的切向風(fēng)速會大于來流風(fēng)速。在距建筑物頂端較低位置時風(fēng)速會出現(xiàn)下降，也就是距離建筑物頂部的較低范圍內(nèi)的風(fēng)速會小于來流風(fēng)速。

3.3 風(fēng)力機組的數(shù)值模擬

通過對建筑物的速度流場進行數(shù)值分析，可以間接證明在建筑物附近安裝風(fēng)力機能在一定程度上提高風(fēng)力機性能。為進一步驗證建筑物附近安裝風(fēng)力機組的可行性，對風(fēng)力機組在特定位置的空氣動力學(xué)性能進行數(shù)值模擬，圖7是建筑物周圍的風(fēng)力機組分布圖。建筑物前方4R處為1號風(fēng)力機，左右兩側(cè)5R處分別為2號和4號風(fēng)力機，正上方4R處為3號風(fēng)力機。

圖8為風(fēng)力機的壓力云圖，建筑物與1號風(fēng)力機之間的距離相對較近，1號風(fēng)力機受到建筑物的反壓影響，由于建筑物的阻礙，來風(fēng)來不及擴散，其整個風(fēng)輪的葉片受到較大壓力。整個風(fēng)力機壓力增大會增加風(fēng)力機的振動頻率，對風(fēng)力機的旋轉(zhuǎn)產(chǎn)生消極影響，降低風(fēng)力機的穩(wěn)定性。由于建筑物附近壓強擴散的影響，2號風(fēng)力機受到的壓強要稍稍大于4號風(fēng)力機受到的壓強。3號風(fēng)力機與建筑物的距離相對較大，建筑物后方的負(fù)壓區(qū)域和正壓區(qū)域形成的壓力差會產(chǎn)生大風(fēng)區(qū)域， 3號風(fēng)力機能夠獲得更多的風(fēng)能，從而降低棄風(fēng)率。

圖9為建筑物周圍4臺風(fēng)力機的轉(zhuǎn)矩圖，1號風(fēng)力機與建筑物之間的壓力差較小，此范圍內(nèi)的風(fēng)速變化復(fù)雜，風(fēng)能相互抵消，風(fēng)力機產(chǎn)生的轉(zhuǎn)矩較小。建筑物兩側(cè)的速度流場對稱，2號和4號風(fēng)力機轉(zhuǎn)矩相差較小，與1號風(fēng)力機相比轉(zhuǎn)矩大，能夠產(chǎn)生更多的經(jīng)濟效益。建筑物上方的3號風(fēng)力機的轉(zhuǎn)矩是4臺風(fēng)力機中最大的。綜合分析，除了1號風(fēng)力機的轉(zhuǎn)矩偏小，其他3臺風(fēng)力機的轉(zhuǎn)矩都大于1號位置的轉(zhuǎn)矩，說明建筑物附近的風(fēng)力機應(yīng)盡量安裝在建筑物兩側(cè)和上方的位置。另外，也證實：在建筑物附近安裝新型鸚鵡螺垂直軸風(fēng)力機組是可行的，而且在合適的位置可以增加風(fēng)力機的動力性能。

4 結(jié) 語

本文對2，3，4葉片的鸚鵡螺風(fēng)力機進行數(shù)值仿真，對比風(fēng)力機的轉(zhuǎn)矩性能和風(fēng)能利用率，得出3個葉片的鸚鵡螺等角螺線型風(fēng)力機性能最優(yōu)。在建筑物工況下模擬，發(fā)現(xiàn)建筑物擾流對下游的風(fēng)速和流向產(chǎn)生影響，其中最大風(fēng)速約是來流風(fēng)速的1.5倍。對建筑物附近風(fēng)力機組的壓力以及轉(zhuǎn)矩特性分析發(fā)現(xiàn)，1號風(fēng)力機位置風(fēng)速變化紊亂，風(fēng)力機使用壽命會降低，2，3，4號風(fēng)力機轉(zhuǎn)矩特性明顯提升。本文中的排布方式可以有效提高風(fēng)力機組的性能，使風(fēng)能得到更有效的利用。本文僅針對建筑物特定位置風(fēng)力機組進行了研究，未來將會對多種排布方式展開研究，進一步提升城市風(fēng)能利用效率。

參考文獻/References：

[1]賀婷婷，趙曉麗.風(fēng)力發(fā)電的成本效益分析——以張家口為例[J].中外能源，2019，24（5）：21-25.

HE Tingting，ZHAO Xiaoli.Cost-benefit analysis of wind power generation： Taking Zhangjiakou as an example[J].Sino-Global Energy，2019，24（5）：21-25.

[2]GRECO L， TESTA C. Wind turbine unsteady aerodynamics and performance by a free-wake panel method[J]. Renewable Energy，2020.doi： 10.1016/j.renene.2020.08.002.

[3]JUNG C，SCHINDLER D. Introducing a new approach for wind energy potential assessment under climate change at the wind turbine scale[J]. Energy Conversion and Management，2020，225： 113425.

[4]常德生，袁超.張家口-承德地區(qū)風(fēng)能和太陽能資源開發(fā)策略研究[J].中國高新科技，2020（10）：34-37.

CHANG Desheng，YUAN Chao.Research on the development wind strategy and solar energy resources in Zhangjiakou-Chengde area[J].China High-Tech，2020（10）：34-37.

[5]LI Zheng，ZHANG Wenda，HAN Ruihua，et al. Design and performance analysis of distributed equal angle spiral vertical axis wind turbine [J]. Recent Patents on Engineering，2020，14（1）：120-132.

[6]張少波，張二兵.論垂直軸風(fēng)力發(fā)電機在基站上的應(yīng)用前景[J].信息通信，2018，31（9）：275-277.

[7]潘登宇，李琦芬，胡丹梅，等.垂直軸風(fēng)力機模塊化研究進展及其應(yīng)用展望[J].熱能動力工程，2017，32（sup1）：6-13.

PAN Dengyu，LI Qifen，HU Danmei，et al.Progress and application of vertical axis wind turbine modulation investigation[J].Journal of Engineering for Thermal Energy and Power，2017，32（sup1）：6-13.

[8]ABUGAZIA M，EL DAMATTY A A，DAI K S，et al. Numerical model for analysis of wind turbines under tornadoes[J]. Engineering Structures，2020，223： 111157.

[9]李爭，齊偉強，于瀟雪.分布式能源用垂直軸風(fēng)力機的結(jié)構(gòu)優(yōu)化設(shè)計[J].機械設(shè)計與制造，2020（4）：36-40.

LI Zheng，QI Weiqiang，YU Xiaoxue.Optimal structure design of vertical axis wind turbine for distributed energy applications[J].Machinery Design & Manufacture，2020（4）：36-40.

[10]馬璐，陳明陽，王曉東，等.垂直軸風(fēng)力機三維效應(yīng)的數(shù)值模擬研究[J].可再生能源，2020，38（3）：326-332.

MA Lu，CHEN Mingyang，WANG Xiaodong，et al.Numerical simulation of vertical axis wind turbine in 3D effect[J].Renewable Energy Resources，2020，38（3）：326-332.

[11]劉陳，運洪祿，呂續(xù)艦.基于CFD的二維垂直軸風(fēng)力機性能計算[J].太陽能學(xué)報，2020，41（2）：144-151.

LIU Chen，YUN Honglu，LYU Xujian.CFD based study on two-dimension vertical axis wind turbines[J].Acta Energiae Solaris Sinica，2020，41（2）：144-151.

[12]FERDOUES M S， EBRAHIMI S， VIJAYARAGHAVAN K. Multi-objective optimization of the design and operating point of a new external axis wind turbine[J]. Energy， 2017， 125：643-653.

[13]王遠(yuǎn)成，吳文權(quán)，黃遠(yuǎn)東，等.建筑物周圍繞流流場的三維數(shù)值模擬[J].工程熱物理學(xué)報，2004，25（2）：235-237.

WANG Yuancheng，WU Wenquan，HUANG Yuandong，et al.3D-numerical modeling of flow around buildings[J].Journal of Engineering Thermophysics，2004，25（2）：235-237.

[14]李雪琪. 超高層建筑集風(fēng)效果研究——以武漢市超高層為例[D]. 武漢： 華中科技大學(xué)，2019.

LI Xueqi. Study on the Effect of Collecting Wind in Skyscraper： Taking Skyscraper in Wuhan as an Example[D]. Wuhan： Huazhong University of Science and Technology，2019.

[15]CHEN T Y，CHEN Y Y. Developing a vortical stator assembly to improve the performance of drag-type vertical axis wind turbines[J]. Journal of Mechanics，2015，31（6）：693-699.

[16]POUREH A，NOBAKHTI A. Scheduling of H ∞ controller in horizontal axis wind turbines[J]. Control Engineering Practice，2020，38（7）：102-104.

[17]YAGAWA G， YOSHIMURA S， NAKAO K. Automatic Mesh Generation of Complex Geometries Based on Fuzzy Knowledge Processing and Computational Geometry[M]. [S.l.]： IOS Press， 1995.

[18]楊耀文，劉正興.三維接觸問題中的網(wǎng)格自動劃分和單元生成[J].上海交通大學(xué)學(xué)報，1996，30（5）：131-135.

YANG Yaowen，LIU Zhengxing.Automatic grid and element generation for three-dimensional cantact problem[J].Journal of Shanghai Jiao Tong University，1996，30（5）：131-135.

[19]孫欣，彭劍峰，趙凱，等.基于Fluent軟件的泵前進氣管部流場模擬[J].能源與環(huán)保，2019，41（9）：111-116.

SUN Xin，PENG Jianfeng，ZHAO Kai，et al.Flow field simulation of pump front intake manifold based on Fluent software[J].China Energy and Environmental Protection，2019，41（9）：111-116.

[20]ZHAO Xu，ZHOU Ping，LIANG Xiao，et al. The aerodynamic coupling design and wind tunnel test of contra-rotating wind turbines[J]. Renewable Energy，2020，38（7）：146-148.

[21]謝典，顧煜炯，楊宏宇，等.螺旋型垂直軸風(fēng)力機的氣動性能研究及結(jié)構(gòu)參數(shù)優(yōu)化[J].農(nóng)業(yè)工程學(xué)報，2017，33（22）：262-268.

XIE Dian，GU Yujiong，YANG Hongyu，et al.Aerodynamic performance research and structural parameter optimization of helical vertical axis wind turbine[J].Transactions of the Chinese Society of Agricultural Engineering，2017，33（22）：262-268.