朱貴森, 鄭長川

(1.烏海職業(yè)技術(shù)學(xué)院,內(nèi)蒙古自治區(qū) 烏海 016000; 2.海勃灣發(fā)電廠,內(nèi)蒙古自治區(qū) 烏海 016000)

覆冰對風(fēng)力發(fā)電機葉片氣動性能影響的仿真研究

朱貴森1, 鄭長川2

(1.烏海職業(yè)技術(shù)學(xué)院,內(nèi)蒙古自治區(qū) 烏海 016000; 2.海勃灣發(fā)電廠,內(nèi)蒙古自治區(qū) 烏海 016000)

覆冰對風(fēng)力發(fā)電機葉片會造成一定的傷害，為了研究覆冰對風(fēng)力發(fā)電機葉片氣動性能的影響，使用κ-ε紊流模型對NACA 23012翼型進行了仿真分析，著重分析了在40°攻角情況下翼型覆冰前后的氣動性能，發(fā)現(xiàn)該翼型在仿真條件下升力明顯減小、阻力增大，升阻比由覆冰前的1.594下降到了覆冰后的0.7，機翼附近的壓強云圖顯示翼型上表面已形成了渦流區(qū)域，說明空氣動力性能發(fā)生惡化，翼型產(chǎn)生了失速現(xiàn)象。同時對比了覆冰前后翼型附近的壓強、空氣速度，仿真結(jié)果表明覆冰的存在是翼型周圍的空氣壓強、速度場畸變得更嚴(yán)重。

覆冰； 風(fēng)機葉片； 氣動性能

0 引 言

由于特定季風(fēng)氣候、地形地理條件，我國可利用風(fēng)力資源地區(qū)集中在東北、華北、西北(“三北”)以及東南沿海地帶，具有很大的商業(yè)開發(fā)價值。寒冷的冬春季節(jié)，“三北”地區(qū)的平均氣溫在零度以下，在雨雪天氣中，空氣中的過冷卻水滴在風(fēng)力作用下撞擊在風(fēng)輪機葉片上形成覆冰[1]。覆冰改變了風(fēng)機葉片的線型輪廓，使葉片的空氣動力性能惡化，進而導(dǎo)致風(fēng)輪機力矩不足以帶動風(fēng)力發(fā)電機在額定狀態(tài)在運行，導(dǎo)致電能減產(chǎn)，覆冰嚴(yán)重的情況下將導(dǎo)致葉片變形甚至風(fēng)機損毀，大塊覆冰脫落導(dǎo)致的葉片振動也會對風(fēng)力發(fā)電機葉片產(chǎn)生危害。

葉片覆冰外形與攻角、空氣中過冷卻水的含量、風(fēng)速等因素均有關(guān)系，可以通過數(shù)值模擬的方法選擇符合自然環(huán)境下的控制參數(shù)獲得翼型覆冰外形[2]。

隨著風(fēng)電快速發(fā)展，國內(nèi)研究者對該問題的研究集中在通過選用合適的紊流模型和計算流體力學(xué)方法，對不同風(fēng)機翼型、不同覆冰形態(tài)的空氣動力性能進行仿真，普遍得到覆冰使葉片空氣動力性能惡化的結(jié)論[3-8]。挪威、加拿大、日本等研究者更加注重探討影響葉片覆冰的因素，并嘗試使用三維模型、完善算法結(jié)合風(fēng)洞覆冰試驗結(jié)果對葉片覆冰進行分析[9-13]。

本文分析了覆冰對風(fēng)力發(fā)電機葉片氣動性能的影響，通過仿真軟件使用κ-ε紊流模型分析葉片覆冰前后升力、阻力的變化，分析這種變化產(chǎn)生的后果以及提出相應(yīng)的解決方案。研究覆冰對風(fēng)機葉片氣動性能的影響，對防止覆冰事故有著一定的指導(dǎo)意義。

1 風(fēng)力發(fā)電機葉片的氣動性能

1.1 相關(guān)概念

風(fēng)力發(fā)電機葉片通過捕捉風(fēng)能并將其轉(zhuǎn)換成機械能，通過齒輪箱增速帶動發(fā)電機發(fā)電。風(fēng)能與風(fēng)速的3次方成正比，當(dāng)風(fēng)速較小時不能正常啟動。風(fēng)輪機的轉(zhuǎn)速取決于葉片受力的大小，風(fēng)吹向葉片時，葉片同時受到升力和阻力的作用，來風(fēng)方向與葉片中線的夾角稱為攻角。

翼型的選擇對風(fēng)力發(fā)電機的效率有很大影響，使用具有較高升阻比的翼型的風(fēng)力發(fā)電機效率越高。目前，水平軸風(fēng)力發(fā)電機葉片使用的翼型主要有航空翼型及其修正翼型[13]，如NACA系列、NASA LS系列；也有專門的風(fēng)力機翼型，如美國的NERL S系列、荷蘭的DU系列等。本文選用NACA 23012翼型進行覆冰前后的空氣動力性能分析。

風(fēng)力利用系數(shù)是風(fēng)輪機將風(fēng)能轉(zhuǎn)換為機械能的效率，是風(fēng)速、葉片轉(zhuǎn)速、葉片直徑和槳距角有關(guān)的函數(shù)。貝茨理論指出理想風(fēng)輪的最大風(fēng)能利用系數(shù)為0.593。

1.2 風(fēng)力發(fā)電機葉片的氣動性能

翼型的氣動性能是指翼型在與氣體發(fā)生相對運動的過程中，氣體的流動特性對翼型升力、阻力的影響。氣動性能發(fā)生惡化，風(fēng)輪機不能獲得足夠的力矩帶動發(fā)電機正常發(fā)電，嚴(yán)重情況下風(fēng)輪機將出現(xiàn)失速現(xiàn)象產(chǎn)生一系列后果。

失速現(xiàn)象是指當(dāng)攻角超過某一臨界值，翼型附近的其他流動狀況發(fā)生惡化，邊界層受到破壞，在葉片背面尾端出現(xiàn)渦流區(qū)。一般情況下，攻角超過臨界值越多，產(chǎn)生的失速現(xiàn)象越嚴(yán)重，氣體的流動阻力越大。葉片失速導(dǎo)致平衡狀態(tài)破壞，葉片受到的應(yīng)力增加，不及時糾正則可能發(fā)生葉片斷裂事故。

2 覆冰對風(fēng)力發(fā)電機葉片氣動性能影響的仿真研究

2.1 仿真模型

風(fēng)力發(fā)電機在寒冷的雨雪天氣中運行過程中空氣的雷諾系數(shù)大，當(dāng)雷諾系數(shù)大于10 000則可視為完全紊流。標(biāo)準(zhǔn)的κ-ε模型是一個半經(jīng)驗?zāi)Ｐ蚚14]，其有兩個基本方程，一個是湍動能κ的傳輸方程，另一個是擴散率ε的傳輸方程，適用于完全紊流流動。因此本文選用該模型作為葉片覆冰仿真所用的紊流模型。

2.2 覆冰前后葉片升阻力的變化

本文使用標(biāo)準(zhǔn)κ-ε模型對NACA 23012翼型進行覆冰前后的氣動性能仿真，該翼型的設(shè)計升力系數(shù)為0.3，最大彎度位置為翼型的0.15倍弦長處，中弧線為簡單型，相對厚度為12%。分析該翼型在覆冰前后升力、阻力以及翼型附近其他參數(shù)的變化?，F(xiàn)以攻角為40°為例，說明覆冰對葉片空氣動力性能的影響。

本文所有仿真都是在表1所示基礎(chǔ)上完成的。

表1 仿真分析條件

覆冰前后的網(wǎng)格圖如圖1所示

圖1 覆冰前后翼型的網(wǎng)格圖(上圖覆冰前，下圖覆冰后，下同)

對翼型外部繞流問題需要定義一個邊界，邊界離翼型越遠(yuǎn)則邊界對流動的影響越小，計算越精確。本文中，邊界離翼型的距離為10倍翼型弦長。對覆冰前后的翼型使用同樣的方法進行劃分，網(wǎng)格劃分得足夠密集以滿足計算精度的要求。

(1) 覆冰前后翼型升力、阻力的變化。使用Gambit將網(wǎng)格劃分好之后導(dǎo)入Fluent軟件，選擇κ-ε模型以及設(shè)定好其他相關(guān)參數(shù)。進行迭代計算，迭代收斂后可以得到翼型的升力報告和阻力報告。結(jié)果顯示，在攻角為40°風(fēng)速為10 m/s時，翼型覆冰前的升力為93.28 N，阻力為58.5 N；覆冰后，翼型的升力下降為44.26 N，阻力上升為63.19 N；升阻比由覆冰前的1.594下降為覆冰后的0.7。翼型的升力驟然下降、阻力加大，這說明覆冰已經(jīng)使葉片進入了失速狀態(tài)。

(2) 覆冰前后翼型表面壓力分布。覆冰前后翼型上下表面的壓力并沒有明顯的變化，翼型迎流方向上表面空氣流速大，壓強小，由于空氣來流在翼型的前部已經(jīng)發(fā)生分離，翼型上表面中部和尾部壓力受來流影響已不大，因此導(dǎo)致翼型尾部上下表面的壓力趨于一致(見圖2)。

圖2 覆冰前后翼型表面的壓力分布圖

(3) 覆冰前后翼型附近的壓強云圖。攻角為40°時，來流在翼型下表面前端形成高壓區(qū)，同時下表面相對上表面都是高壓區(qū)域，此時機翼附近的壓強云圖如圖3中的上圖。覆冰后，由于覆冰外形的影響，高壓區(qū)出現(xiàn)在覆冰區(qū)域所在的翼型前端，翼型下表面的壓強也明顯小于該區(qū)域的壓強。并且在翼型的上表面區(qū)域出現(xiàn)了渦流區(qū)，致使外部空氣向內(nèi)流動，空氣動力性能發(fā)生惡化，增大了翼型的阻力，翼型產(chǎn)生失速現(xiàn)象。

(4) 覆冰前后翼型附近的速度云圖。覆冰前翼型周圍空氣最大流動速度為17.5 m/s，同樣由于覆冰形狀的影響，附近空氣局部最大速度達(dá)到了19.3 m/s，但這翼型微小區(qū)域的速度增大并沒有改善翼型的氣動性能。由于覆冰的存在，被覆冰擋住來流的覆冰后部區(qū)域的空氣流動變得異?；靵y，翼型上表面尾部空氣流動速度有小范圍減小(見圖4)。

(5) 覆冰前后翼型附近的壓強、氣流速度變化。對比了覆冰前后翼型附近的壓強、速度，其中的壓強數(shù)值是相對于標(biāo)準(zhǔn)大氣壓強的數(shù)值，可以發(fā)現(xiàn)，覆冰后翼型周圍空氣的壓強最大值有顯著增大，壓強最小值也有增大，空氣最大流速也稍有變大。這都是由于覆冰造成的。盡管數(shù)值有所增大，但這些只是翼型附近的某個點的數(shù)值，并不反映整個翼型的受力狀態(tài)。覆冰沒有使翼型的氣動性能變得更優(yōu)，相反，覆冰使翼型周圍的壓強和速度場畸變更加厲害。

圖3 覆冰前后翼型附近的壓強云圖

圖4 覆冰前后翼型附近的速度云圖

表2 覆冰前后翼型附近的壓強、速度比較

3 結(jié) 語

本文使用κ-ε紊流模型研究了覆冰對風(fēng)力發(fā)電機葉片氣動性能的影響，發(fā)現(xiàn)在這種特定的覆冰外形以及40°的來流攻角下，翼型將發(fā)生失速事故。此時需要朝葉片迎風(fēng)面積增大的方向調(diào)節(jié)槳距角以及時恢復(fù)風(fēng)輪機正常運行所需的氣動性能。在相同仿真條件下，對覆冰前后翼型附近的壓強、空氣速度也做出了比較，發(fā)現(xiàn)覆冰的存在是翼型周圍的空氣壓強、速度場畸變得更嚴(yán)重。

由仿真結(jié)論可知，覆冰使翼型的氣動性能發(fā)生惡化，在這種特定的條件下，如果不調(diào)節(jié)槳距角使葉片迎風(fēng)面積增大，風(fēng)輪機將發(fā)生失速事故。若氣動性能已經(jīng)足夠惡化，調(diào)節(jié)槳距角也并不一定能讓風(fēng)輪機回復(fù)正常運行、發(fā)電機達(dá)到額定發(fā)電狀態(tài)，此時風(fēng)輪機將進行停機處理。值得注意的是，覆冰并不必然導(dǎo)致風(fēng)輪機失速，具體是否失速需要對特定的覆冰形狀及其對葉片的氣動性能影響進行分析。

[1] 李俊峰，施鵬飛，高 虎.中國風(fēng)電發(fā)展報告(2010)[M].?？冢汉Ｄ铣霭嫔?，2010:12-13.

[2] 易 賢,朱國林,王開春,等. 翼型積冰的數(shù)值模擬[J]. 空氣動力學(xué)學(xué)報，2002(4):428-433.

[3] 朱程香,王 瓏,孫志國，等.風(fēng)力機葉片翼型的結(jié)冰數(shù)值模擬研究[J].空氣動力學(xué)學(xué)報,2011,29(4):522-528.

[4] 李聲茂,李 巖.垂直軸風(fēng)力機葉片翼型結(jié)冰的數(shù)值計算[J].動力工程學(xué)報,2011,30(3):214-219.

[5] 王文龍，何 斌，傅 潔，等.風(fēng)力機葉片翼型結(jié)冰的氣動性能數(shù)值模擬[J].可再生能源,2013,31(4):60-70.

[6] 張聘亭,楊 濤,盧緒祥.覆冰對風(fēng)力機專用翼型氣動性能影響的數(shù)值研究[J].動力工程學(xué)報,2011,31(12):955-959.

[7] 鄧曉湖,盧緒祥,李錄平,等. 水平軸風(fēng)力機葉片翼型結(jié)冰的數(shù)值模擬[J]. 電力與能源,2010(5)：266-271.

[8] 陳維建,張大林.瘤狀冰結(jié)冰過程的數(shù)值模擬[J].航空動力學(xué)報,2005,20(3):472-476.

[9] Virk M S, Homola M C, Nicklasson P J. Atmospheric icing on large wind turbine blade[J]. International Journal of Energy and Environment,2012,3(1):1-8.

[10] Yuichiro Imamura,Kazuyuki Toda,Makoto Yamamoto.Numerical simulation of ice accretion on wind turbine blade[J].Computational Mechanics.2008,32(3):16-19.

[11] Villalpando F, Reggio M, Ilinca A. Numerical study of flow around iced wind turbine airfoil[J]. Engineering Applications of Computational Fluid Mechanics, 2012, 6(1):39-45.

[12] Fernando Villalpando,Marcelo Reggio,Adrian Ilinca.Numerical study of flow around iced wind turbine airfoil[J].Engineering Applications of Computational Fluid Mechanics,2012,6(1):39-45.

[13] 李銀然.風(fēng)力機專用翼型的空氣動力特性研究[D].蘭州：蘭州理工大學(xué)，2010:11-15.

[14] 韓占忠.流體工程仿真計算實例與分析[M].北京:北京理工大學(xué)出版社，2009:90-91.

Comparison of Flashover Performance of Typical Post Insulators under Pollution Conditions

ZHU Guisen， ZHENG Changchuan

(1. Wuhai Vocational and Technical College, Wuhai 016000, Inner Mogolia, China； 2. Haibowan Power Plant， Wuhai 016000, Inner Mogolia, China)

The ice will cause some harm to the wind turbine blades. In order to study the effect of icing on the aerodynamic performance of wind turbine blade, using theκ-εturbulence model, NACA 23012 airfoil was simulated and analyzed, the dynamic performances before and after the 40° angle of attack of the airfoil icing gas were emphatically analyzed. It was found that in the simulation of the airfoil under the condition of lift reduces and the resistance increases, the lift drag ratio decreases from 1.594 down to the 0.7. The pressure cloud near wing displays that on the airfoil surface there is a vortex area. It shows that the aerodynamic performance deteriorates, and the airfoil stalls. At the same time, the pressure and air velocity around the airfoil were compared, the simulation results show that the air pressure and velocity fields around the airfoil are more serious.

icing; wind turbine blade; aerodynamic characteristics

2016-11-25

朱貴森(1980-)，男，河北衡水人，碩士，講師，主要從事電廠熱能動力方面的研究。

Tel.：15048349006;E-mail：30210912@qq.com

TP 391

A

1006-7167(2017)08-0129-04