大攻角下開縫位置改變對風(fēng)力機葉片影響研究

王 龍, 孫倫業(yè), 張 瑾

(安徽理工大學(xué) 機械工程學(xué)院,安徽 淮南 232001)

大攻角下開縫位置改變對風(fēng)力機葉片影響研究

王 龍, 孫倫業(yè), 張 瑾

(安徽理工大學(xué) 機械工程學(xué)院,安徽 淮南 232001)

風(fēng)力機葉片在大攻角條件下存在著嚴重的流動分離現(xiàn)象,降低風(fēng)電機組的發(fā)電效率,文章采用計算流體力學(xué)方法研究開縫位置改變對風(fēng)力機葉片性能的影響。研究結(jié)果表明:在一定攻角范圍下,射流會改善流動狀況,縮小渦旋結(jié)構(gòu)影響區(qū)域;開縫位置位于分離點附近時,氣動性能最好,隨著射流位置后移,其控制效果減弱,即使射流位置位于回流區(qū),仍可減弱回流區(qū)影響范圍與強度;射流有利于提高葉片升力系數(shù),改善大攻角下葉片氣動性能及穩(wěn)定性。

射流;計算流體力學(xué);風(fēng)力機葉片;附面層控制

0 引 言

風(fēng)力機葉片作為風(fēng)電機組的重要工作部件,其氣動性能高低直接關(guān)系到發(fā)電功率輸出。在來流大攻角條件下,風(fēng)力機葉片會存在明顯的流動分離現(xiàn)象,對水平軸風(fēng)力機而言,其升力系數(shù)會嚴重下降。因此在非工況設(shè)計條件下,其空氣動力學(xué)問題應(yīng)是研究重點之一[1-2]。

計算流體力學(xué)方法[3-5]可以較高精度地再現(xiàn)風(fēng)力機運行過程中的流場擬序結(jié)構(gòu),可處理其遇到的復(fù)雜流動問題。因此隨著現(xiàn)代計算能力的迅速發(fā)展,計算流體力學(xué)在風(fēng)力機氣動性能分析中占據(jù)的地位和作用愈加重要和明顯。

風(fēng)力機葉片的邊界層流動在大攻角下抗逆壓梯度能力較弱,存在嚴重的流動分離現(xiàn)象,從而對風(fēng)力機升阻比及安全穩(wěn)定性造成嚴重影響。通過流體控制手段改變附面層相干結(jié)構(gòu),較大程度減小渦旋結(jié)構(gòu)影響區(qū)域,提高風(fēng)力機運行安全及穩(wěn)定性,無疑具有重要意義。

文獻[6]對合成射流控制圓柱分離進行實驗研究,研究結(jié)果表明在不同雷諾數(shù)下合成射流均能改善圓柱繞流的分離狀況。

在葉片流動分離機理研究方面,國內(nèi)外學(xué)者均進行相對深入的研究。文獻[7]對NACA0015 翼型進行射流試驗,研究表明在低馬赫數(shù)條件下,射流可以提高升力系數(shù),改善流動狀況。文獻[8]對低雷諾數(shù)下機翼翼型的前緣流動分離機進行研究,發(fā)現(xiàn)前緣橢圓弧靠近葉身位置存在吸力峰,并在強逆壓梯度作用下發(fā)生流動分離;翼型上表面渦系結(jié)構(gòu)包含駐留渦、脫落渦和二次渦,具有強烈的非定常性。

文獻[9]對低雷諾數(shù)下多孔分布式抽吸結(jié)構(gòu)對流動分離的控制進行相關(guān)研究,結(jié)果表明通過在翼型前緣附近合適位置進行附面層抽吸,可以抑制流動分離。

風(fēng)力機葉片表面開縫示意圖如圖1所示。

圖1 風(fēng)力機葉片表面開縫示意圖

本文采用計算流體力學(xué)方法研究開縫位置改變對風(fēng)力機葉片性能的影響,獲取不同攻角下的葉片流場及升力系數(shù),討論了4種不同開縫位置對流場及氣動性能的影響,可為非工況條件下風(fēng)力機葉片設(shè)計提供參考依據(jù)。

1 控制方程及湍流模型

自然界所存在的流體滿足質(zhì)量守恒、動量守恒和能量守恒三大規(guī)律,在忽略徹體力和源項的條件下,其數(shù)學(xué)積分的表現(xiàn)形式為:

(1)

(2)

(3)

其中,Ω為空間內(nèi)任一控制體的體積,由閉合曲面邊界?圍成;ρ、v分別為密度和速度矢量;q為熱流矢量；et為流體的總能量；n為曲面的法矢量；π為應(yīng)力張量的大小,其牛頓流體的本構(gòu)方程為:

(4)

S-A模型為文獻[10]提出的一方程湍流模型,該方程在壓聲速和跨聲速流動領(lǐng)域中的翼型及機翼應(yīng)用廣泛,但由于其只含1個湍流變量,相對于其他2個方程模型來說計算量和耗時均較小。不引入轉(zhuǎn)捩項的守恒型輸運方程可改寫為:

(5)

其中,Cb1=0.135 5;σ=2/3;Cb2=0.622;K=0.41;υ為流體黏度；其余符號表達意思可參見文獻[10]。

2 模型及網(wǎng)格

根據(jù)開縫位置不同(縫所在位置分別為0.5、0.6、0.7、0.8 m處),依次命名為model-1、model-2、model-3、model-4,其中,model-0表示原型模型,即未開縫。本文計算所采用模型如圖2所示,為簡單起見,僅給出了原始模型model-0及開縫位置x=0.5 m的model-1的對比。

圖2 風(fēng)力機葉片模型及開縫模型

計算域劃分圖如圖3所示,計算域的上下邊界為15倍弦長,為保證網(wǎng)格的正交性,計算域前部采用“C”結(jié)構(gòu),流體出口邊界距離原點為30倍弦長。邊界條件由遠場及壁面組成。網(wǎng)格示意圖如圖4所示,包括整體網(wǎng)格及風(fēng)力機葉片前后緣局部網(wǎng)格,近壁面第1層網(wǎng)格距離壁面約為1.0×10-5m,總的網(wǎng)格數(shù)約為8×104個。

圖3 計算區(qū)域

圖4 網(wǎng)格示意圖

3 計算結(jié)果

為保證計算結(jié)果的正確性和合理性,本文首先對有實驗結(jié)果的S809原型葉片進行計算。來流空氣的馬赫數(shù)為0.12,雷諾數(shù)為2×106。其計算所獲取的壁面壓力系數(shù)如圖5所示,其中,實驗結(jié)果為文獻[11]所測量的實驗值;計算結(jié)果為采用課題組自行開發(fā)的軟件[12]所模擬的結(jié)果,可見兩者數(shù)據(jù)吻合較好。

圖5 葉片壁面壓力系數(shù)

本文計算條件如下:① 計算攻角范圍為12°～20°,如圖6所示,計算間隔角度為2°;② 對4個位置(X=0.5、0.6 、0.7 、0.8 m)的開縫及原型模型計算,且縫寬均為2 cm、噴射角度為0°角;③ 風(fēng)速約為40 m/s,遠場壓強為72 465 Pa,縫進口壓力為74 000 Pa。

圖6 計算攻角示意

升力系數(shù)對比如圖7所示,其中,圖7a所示為4種不同開縫位置模型及原始模型葉片在攻角為12°～20°條件下的升力系數(shù)計算結(jié)果,從圖7a中可以看出,開縫后葉片升力系數(shù)均比原型葉片的升力系數(shù)增加較多,其無量綱增加的幅度均在0.2以上;圖7b所示為不同開縫模型升力系數(shù)對比,其中,model-1(X=0.5 m)氣動性能相對最高; model-2性能次之;model-4性能相對最差,與原始模型“model-0”相當(dāng)。從升力系數(shù)隨攻角的變化可以看出,隨著開縫位置后移,其升力系數(shù)逐漸下降,靜態(tài)失速角前移。model-1與model-2在12°～16°之間升力系數(shù)幾乎一致,隨著攻角增加,model-2的升力系數(shù)從1.64降到1.18。

model-1的最大升力系數(shù)為1.71,其失速角為18°;model-2的最大升力系數(shù)為1.67,在攻角為16°出現(xiàn),較model-1的最大升力系數(shù)小0.04;model-3的最大升力系數(shù)為1.58,在攻角為16°發(fā)生,較model-1的最大升力系數(shù)小0.13;model-4的最大升力系數(shù)為1.40,在攻角為12°出現(xiàn),較model-1的最大升力系數(shù)小0.31,隨后升力系數(shù)隨攻角增加而下降。

圖7 升力系數(shù)對比

18°攻角下不同開縫位置葉片的流場分布如圖8所示。從圖8可以看出,在加入射流以后,model-1、model-2的流場均沒有出現(xiàn)流動分離現(xiàn)象,并且流場分布大體相似;而model-3、model-4在葉背后半部的吸力面上方出現(xiàn)渦旋結(jié)構(gòu)區(qū)域。對比圖8c、圖8d流場分布還可以看出,隨著射流起始位置后移,其改善流場狀態(tài)的能力進一步削弱,渦旋影響區(qū)域顯著增加。在回流區(qū)的影響下,model-3和model-4的升力系數(shù)均開始下降,下降幅度在0.34以上,小于model-1、model-2的升力系數(shù)。

18°攻角下開縫葉片與原型葉片的流場分布對比如圖9所示。從圖9可以看出,2種模型在較大攻角下其流場出現(xiàn)截然不同的流動現(xiàn)象,其中,model-1在射流的作用下不再出現(xiàn)回流區(qū),model-0在葉背后半部吸力面上方出現(xiàn)更大的漩渦區(qū)域,并出現(xiàn)尾緣渦并向下游脫落。開縫與原型葉片回流區(qū)位置流場如圖10所示,對比圖10a、圖10b可以發(fā)現(xiàn)分離點位于射流起始位置附近,附加的高能流體從回流區(qū)上游出發(fā),迫使渦旋結(jié)構(gòu)出現(xiàn)破碎,大渦分解為若干小渦,同時增加附面層低能流體所攜帶的能量,減小附面層厚度,最終達到消除分離回流區(qū)的目的。觀察流場還可以看出,相同攻角下,葉背處漩渦區(qū)域的存在會降低氣動性能。

圖8 攻角18°下不同開縫葉片流線

18°攻角下model-0(原型葉片)與model-1(開縫葉片)壓力系數(shù)對比如圖11所示,從圖11可以看出,在加入射流后,壓力面上的壓力系數(shù)分布差別較小,吸力面上的壓力系數(shù)出現(xiàn)較大不同。原因在于射流出現(xiàn)后,吸力面的流體加速向下游流動,壓強變小,從物理角度來看則是“等效”增加了葉片的彎度與環(huán)量,起到增升效果,即壓力系數(shù)曲線所包含的面積增加,升力系數(shù)變大。

圖9 開縫與原型葉片流場對比

圖10 開縫與原型葉片回流區(qū)位置流場

圖11 開縫與原型葉片壓力系數(shù)對比

通過上述分析可以看出,引入高能氣流在狹縫出口處形成射流,迫使下游附面層低能流體獲取能量,改善附面層流動狀態(tài),對葉片在大攻角下出現(xiàn)的大尺度漩渦及后緣渦起到一定的削弱和清除作用,提高葉片的氣動性能。

4 結(jié) 論

本文采用計算流體力學(xué)方法對4種不同開縫位置風(fēng)力機葉片進行數(shù)值模擬,獲取了不同攻角下的葉片流場及升力系數(shù),并討論了開縫位置對流場及氣動性能的影響,所得結(jié)論如下:

(1) 在來流攻角12°～20°條件下,開縫模型性能均優(yōu)于原始葉片模型,其中,model-1氣動性能相對最高,其最大升力系數(shù)為1.71,較原始模型提高28.6%,失速角為18°。

(2) 大攻角下,射流會抑制流動分離現(xiàn)象發(fā)生,但與射流位置密切相關(guān)。開縫位置位于分離點附近時,氣動性能最好,隨著后移,其控制效果減弱,即使射流位置位于回流區(qū),仍可減弱回流區(qū)影響范圍。

(3) 射流可以增加葉片氣動性能的穩(wěn)定性。高能氣體出口處附近形成高速射流,掃掠附面層低能流體,有利于削弱和清除吸力面上的大尺度渦旋結(jié)構(gòu)及尾緣渦,可以提高風(fēng)力機機械性能。

[1] 中國工程院.中國能源中長期 (2030-2050) 發(fā)展戰(zhàn)略研究 (可再生能源卷) [M].北京:科學(xué)出版社,2011.

[2] 賀德馨.中國風(fēng)能發(fā)展戰(zhàn)略研究[J].中國工程科學(xué),2011,13(6):95-100.

[3] 鐘偉,王同光.SST湍流模型參數(shù)校正對風(fēng)力機CFD模擬的改進[J].2014,35(9):1743-1748.

[4] 周正貴.壓氣機/風(fēng)扇葉片自動優(yōu)化設(shè)計的研究現(xiàn)狀和關(guān)鍵技術(shù)[J].航空學(xué)報,2008,29(2):257-266.

[5] OZGEN S,CANIBEK M.Ice accretion on multi-element airfoils using extended messinger model[J].Heat Mass Transfer,2009,45:305-322.[6] 葉建,鄒正平,陸利蓬,等.低雷諾數(shù)下翼型前緣流動分離機制的研究[J].北京航空航天大學(xué)學(xué)報,2004,30(8):693-697.

[7] 張旺龍,譚俊杰,陳志華,等.低雷諾數(shù)下翼型分離流動抽吸控制特性研究[J].空氣動力學(xué)報,2015,33(1):113-119.

[8] 王晉軍,馮立好,徐超軍.合成射流控制圓柱分離及繞流結(jié)構(gòu)的實驗研究[J].中國科學(xué)E輯：技術(shù)科學(xué),2007,37(7):944-951.

[9] SEIFERT A,DARABI A,WYGNANSKI I.Delay of airfoil stall by periodic excitation[J].Journal of Aircraft,1996,33(4):691-698.

[10] SPALART P,ALLMARAS S.A One-equation turbulence model for aerodynamic flows[C]//30th aerospace Sciences Meeting and Exhibit. [S.l.：s.n.],1992:1.

[11] CHANG Y L,YANG S L,ARICI O.Flow field computation of the NREL S809 airfoil using various turbulence models[C]//Energy Week `96:American Society of Mechanical Engineers and American Petroleum Institute Energy Week Conference and Exhibition,Houston,TX (United States),21 Jan-2 Feb 1996.New York:American Society of Mechanical Engineerings,1996:372133.

[12] WANG Long,SUN Lunye,WU Guang,et al.Research on algorithm of blade vibration for general wind turbine [C]//Proc SPIE9903,Seventh International Symposium on Precision Mechanical Measurements,2016:990325-8.

(責(zé)任編輯 張 镅)

Study of the impact of the change of slotted position at high angle of attack on wind turbine blades

WANG Long, SUN Lunye, ZHANG Jin

(School of Mechanical Engineering, Anhui University of Science and Technology, Huainan 232001, China)

There is serious flow separation phenomenon for the wind turbine blade in the condition of high angle of attack that reduces the power generation efficiency of wind turbine. The impact of the change of the slotted position on wind turbine blades was studied by using the computational fluid dynamics. The result indicated that the jet improved the flow situation, reducing the affected area of the vortex structure, under a certain range of the angle of attack. The aerodynamic performance was the best when the slotted position was located near the separation point. The control effect was weakened with the jet position moving backward. Even if the jet was located in the recirculation zone, the scope and intensity of the recirculation zone could be reduced. The jet could improve the lift coefficient of blades and improve the aerodynamic performance and stability of blades at high angle of attack.

jet; computational fluid dynamics; wind turbine blade; boundary layer control

2016-02-29；

2016-07-04

國家自然科學(xué)基金青年科學(xué)基金資助項目(51505003);安徽省自然科學(xué)基金資助項目(1708085QE123;1708085QA17;1508085QE94);安徽理工大學(xué)國家自然(社會)科學(xué)基金預(yù)研資助項目(10029);安徽理工大學(xué)青年教師科學(xué)研究基金資助項目(12664)和安徽理工大學(xué)引進人才科研啟動基金資助項目(ZY041)

王 龍(1984-),男，安徽蚌埠人,博士,安徽理工大學(xué)講師.

10.3969/j.issn.1003-5060.2017.08.006

TK83

A

1003-5060(2017)08-1037-05