王龍 李相偉 胡南云 李亮 周毅鈞

摘 要：外場工作風(fēng)力機葉片會受到環(huán)境腐蝕，導(dǎo)致風(fēng)力發(fā)電機氣動性能下降，發(fā)電效率降低。為研究腐蝕翼型的氣動性能變化，建立9種腐蝕模型，計算4種風(fēng)速（30m/s、35m/s、40m/s、45m/s）下腐蝕翼型的氣動性能，分析翼型形狀與腐蝕位置的變化對翼型氣動性能的影響。數(shù)值模擬結(jié)果表明來流風(fēng)速與翼型的升阻力成線性關(guān)系;S803翼型的氣動性能優(yōu)于S802和DU93翼型，DU93翼型發(fā)生腐蝕缺陷，氣動性能表現(xiàn)最差;高速流體流過翼型缺陷位置，在缺陷位置發(fā)生流動分離形成渦旋流場，導(dǎo)致翼型氣動性能降低。研究結(jié)果有利于惡劣環(huán)境下的風(fēng)力機葉片選型及葉片的后期維護。

關(guān)鍵詞：腐蝕翼型;氣動性能;高風(fēng)速;流動分離

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

文章編號：1672-1098（2021）02-0016-07

收稿日期：2020-09-30

基金項目：安徽省自然科學(xué)基金資助項目（1708085QE123）

作者簡介：王龍（1984-），男，安徽蚌埠人，副教授，博士，研究方向：風(fēng)力機設(shè)計及風(fēng)力機振動特性。

Research on Aerodynamic Characteristics of Corroded Airfoil under High Wind Speed

WANG Long1，LI XiangWei1， HU NangYun2， LI Liang3， ZHOU YiJun1

（1.School of Mechanical Engineering， Anhui University of Science and Technology， HuainanAnhui 232001，China;2.SafetyEvaluation DepartmentNanchang Anda Security Technology Consulting Co. Ltd.， NanchangJiangxi 330000， China;3.School of Mechanics and Optoelectronics Physics， Anhui University of Science and Technology， HuainanAnhui 232001，China）

Abstract：In the working field， the blades of wind turbines could be corroded by the environment， which would cause the aerodynamic performance of the wind turbine to decrease. Then the efficiency of power generation would decrease. In order to study the changes in the aerodynamic performance of the corroded airfoil， nine corrosion models were established to calculate the aerodynamic performance of the corroded airfoil under four wind speeds （30m/s， 35m/s， 40m/s， 45m/s）. The influence were analyzed of the change of the shape and position of the corroded airfoil to its aerodynamic performance. Numerical simulation results show that the incoming wind speed has a linear relationship with the lift drag of the airfoil， that among the three types of airfoils， the aerodynamic performance of the Type S803 is better than the Types of S802 and DU93 and the aerodynamic performance of the Type DU93 will be the worst when it has a corroded defect and that when the high-speed fluid flows through the defect location of the airfoil， it will form a vortex flow field resulting in a decrease of the aerodynamic performance. The research results are conducive to the selection of wind turbine blades and the later maintenance in harsh environments.

Key words：corroded airfoil; aerodynamic performance; high wind speed; flow separation

風(fēng)沙、雨水長時間沖刷風(fēng)力機葉片后，會導(dǎo)致葉片表面出現(xiàn)腐蝕現(xiàn)象，改變?nèi)~片表面粗糙度，該現(xiàn)象對風(fēng)力機葉片發(fā)電效率極為不利。一旦該缺陷得不到修復(fù)，有可能會導(dǎo)致葉片出現(xiàn)裂紋擴展的現(xiàn)象，將嚴(yán)重影響風(fēng)力發(fā)電機的發(fā)電穩(wěn)定性和使用壽命。

文獻[1]采用數(shù)值模擬方法，對FFA-W3-301 和FFA-W3-211翼型的氣動性能進行分析，模擬結(jié)果的升力系數(shù)與阻力系數(shù)與風(fēng)洞實驗[2]的結(jié)果有較好的一致性，同時發(fā)現(xiàn)翼型表面粗糙度對風(fēng)力機風(fēng)能轉(zhuǎn)化率有很大影響。文獻[3-4]等根據(jù)實際工作中風(fēng)力機葉片磨損狀況，建立三種缺陷翼型幾何模型，基于CFD方法，研究風(fēng)力機葉片磨損對其氣動性能的影響，發(fā)現(xiàn)葉片磨損加劇翼型尾緣附近的流動分離，同時磨損程度越大，葉片的升力系數(shù)越小，阻力系數(shù)越大。文獻[5]等人研究海上風(fēng)力機葉片水滴撞擊問題，發(fā)現(xiàn)水滴撞擊會導(dǎo)致葉片前緣出現(xiàn)腐蝕問題，破壞葉片表面完整度，進而降低風(fēng)力機的氣動性能，影響發(fā)電效率。文獻[6]等通過CFD瞬態(tài)模擬方法，研究NACA 64-618翼型腐蝕狀態(tài)下的氣動性能，發(fā)現(xiàn)葉片的升力系數(shù)降低一半，阻力系數(shù)增大三倍，嚴(yán)重影響風(fēng)力發(fā)電效率。文獻[7]等數(shù)值模擬仿真DU91_W2_250型風(fēng)力機葉片的動態(tài)失速特性，發(fā)現(xiàn)折合頻率與振動幅度的增大會加劇葉片的動態(tài)失速，導(dǎo)致氣動性能變差。文獻[8]等在NACA0018翼型表面裝置固定氣動彈片，并分析氣動彈片對翼型的氣動性能影響，研究發(fā)現(xiàn)增加氣動彈片可提高葉片的氣動新能。文獻[9]以NACA4412翼型為基礎(chǔ)，設(shè)計多種傾斜角度葉尖，研究葉尖小翼結(jié)構(gòu)變化對葉片氣動性能的影響，結(jié)果發(fā)現(xiàn)大斜角葉尖小翼會改善葉尖繞流，提高風(fēng)力機輸出功率。文獻[10]在高風(fēng)速條件下研究環(huán)境參數(shù)對翼型的動態(tài)失速的影響，發(fā)現(xiàn)隨著風(fēng)速的增大，葉片升力系數(shù)逐漸增大，最大升力攻角與之成正比。文獻[11]使用N-S方程研究在不同位置開射流口對S809翼型氣動性能的影響，結(jié)果發(fā)現(xiàn)射流位置于3/4弦長處，攻角18°時升阻比可提高一半，極大改善翼型的氣動性能。文獻[12]等采用CFD方法，研究開縫分布位置對風(fēng)力機葉片氣動性能的影響，發(fā)現(xiàn)大攻角下射流會抑制流動分離現(xiàn)象，開縫位置于1/2弦長處對，升力系數(shù)顯著提高。文獻[13]等分析研究同一系列不同翼型厚度對風(fēng)力機氣動性能的影響，結(jié)果發(fā)現(xiàn)NACA0012翼型厚度對風(fēng)力機的各氣動參數(shù)有較大的影響，在同一系列翼型中存在一最佳翼型厚度。文獻[14]研究大氣穩(wěn)定性和表面粗糙度對風(fēng)力發(fā)電機氣動性能的影響，發(fā)現(xiàn)當(dāng)來流不穩(wěn)定時，氣體流過葉片高低不平粗糙度表面時，會在葉片尾緣產(chǎn)生渦旋尾跡，降低葉片升力系數(shù)。文獻[15]通過3D CFD模擬分析NREL VI型空氣動力學(xué)特性，將S809替換S826型機翼后，發(fā)現(xiàn)葉片流動分離現(xiàn)象延遲，其中S826翼型升力系數(shù)增加了23.7%。

總體來看，風(fēng)力機葉片的翼型形狀以及葉片的表面粗糙程度會對實際工況時的風(fēng)力機的發(fā)電效率產(chǎn)生很大影響。目前國內(nèi)外學(xué)者多研究單一翼型發(fā)生腐蝕狀況時的葉片氣動性能，較少關(guān)注同一型號翼型發(fā)生腐蝕和腐蝕位置變化之間的變化狀況。

本文采用SSTκ-ω湍流模型，以S802、S803和DU93翼型為基礎(chǔ)，設(shè)計3種不同腐蝕位置，建立9種腐蝕翼型模型，在4種高風(fēng)速下，基于CFD方法數(shù)值分析翼型形狀的變化以及腐蝕位置的變化對葉片氣動性能的影響。所做工作有利于不同環(huán)境下的葉片選型及葉片的后期維護，提高惡劣環(huán)境下的風(fēng)力機發(fā)電效率。

1 腐蝕翼型的等效模型

本文對不同型號翼型以及不同腐蝕位置進行仿真分析，研究形狀與腐蝕位置變化對翼型氣動性能的影響。

模擬腐蝕缺陷模型的時候，是在翼型吸力面的不同位置增加凹槽，凹槽長為0.5mm，寬為1mm，研究不同位置的凹槽對翼型氣動性能的影響，三種缺陷位置分布如圖1所示。

2 數(shù)值計算模型

2.1 網(wǎng)格劃分

本文利用Mesh軟件生成9套網(wǎng)格模型，圖2為整體網(wǎng)格結(jié)構(gòu)。使用三角形非結(jié)構(gòu)化網(wǎng)格，同時對翼型附近流場域進行加密處理保證計算精度的準(zhǔn)確性。

2.2 腐蝕翼型網(wǎng)格

30%弦長腐蝕點記作X=0.3，45%弦長腐蝕點記作X=0.45，60%弦長腐蝕點記作X=0.6。翼型前緣首先受到高速流體沖擊，壓強較大，為保證仿真模擬的精準(zhǔn)度，對腐蝕缺陷位置和翼型前緣、尾緣進行加密處理，各工況下的翼型網(wǎng)格如圖3所示。

3 計算結(jié)果與分析

3.1 不同翼型形狀的影響

其它設(shè)置條件未變，僅改變來流風(fēng)速（30m/s，35m/s，40m/s，45m/s），分析S802、S803和DU93翼型對應(yīng)的腐蝕翼型在風(fēng)速變化工況下的氣動性能，其各自對應(yīng)的升力變化曲線與阻力變化曲線如圖4和圖5所示。各工況下升阻力數(shù)值如表1～表3所示。

隨著風(fēng)速的增大，三種翼型的升力與阻力都逐步變大，且升阻力大小與風(fēng)速大小成線性增長關(guān)系。升力變化：當(dāng)腐蝕缺陷位置為X=0.3時，S802與DU93翼型升力數(shù)值較小，S803翼型的升力最大，氣動性能相對最優(yōu);X=0.45與X=0.46時， S803翼型的升力數(shù)值最大，S802升力小于S803，DU93翼型氣動性能相對較差;S803翼型隨著腐蝕缺陷的變化，升力變化較為平緩，彼此數(shù)值相差不大。氣體沿著翼型前緣開始流動，流動到腐蝕缺陷點位置，會提前產(chǎn)生流動分離現(xiàn)象，降低翼型的氣動性能，其中S803翼型抗腐蝕性能較好，腐蝕位置的變化對其升力變化的影響很小。阻力變化：S802翼型腐蝕缺陷位置為X=0.6時，阻力最小;DU93在X=0.45時，對應(yīng)的阻力最小;S803翼型隨著腐蝕位置的變化，阻力也產(chǎn)生相應(yīng)的變化。

3.2 不同腐蝕位置的影響

翼型不同位置發(fā)生腐蝕缺陷，會改變翼型表面流線結(jié)構(gòu)，流體流過翼型表面時，就會造成氣動性能的改變。為研究S802、S803和DU93翼型抗腐蝕性能，現(xiàn)比較分析三種腐蝕缺陷位置時，翼型在各工況下的氣動性能，其對應(yīng)的升力變化曲線與阻力變化曲線如圖6和圖7所示。

當(dāng)腐蝕缺陷位置為X=0.3時，S803翼型相較于S802 與DU93，相同工況下升力數(shù)值最大，氣動性能最好，其中S802氣動升力略小于S803，DU93相對氣動升力最差。當(dāng)腐蝕缺陷位置為X=0.45與X=0.6時，S803翼型氣動升力最大，氣動性能最優(yōu)，相較于X=0.3的工況S802與S803氣動升力數(shù)值相差很小。這表明，當(dāng)翼型表面發(fā)生意外腐蝕，造成翼型結(jié)構(gòu)缺陷時，隨著風(fēng)速的增大，三種翼型的氣動升力業(yè)隨之增大，其中S803翼型抗腐蝕干擾最強，氣動升力性能最好。高風(fēng)速流體沿著翼型前緣快速向后緣流動，當(dāng)流體經(jīng)過X=0.3時，流體處于快速上升過程，DU93翼型結(jié)構(gòu)受到腐蝕缺陷影響較大，氣動升力最小，氣動性能相對最差;當(dāng)流體經(jīng)過X=0.45與X=0.6處時，S803與S802氣動性能相差較小，氣動性能均優(yōu)于DU93翼型。

當(dāng)高速空氣流體流經(jīng)三種腐蝕缺陷位置時，會于缺陷處產(chǎn)生渦旋增加氣動阻力，降低了翼型原有的氣動性能。DU93翼型產(chǎn)生的阻力最大，其氣動性能低于S802與S803翼型。X=0.3時，S802翼型產(chǎn)生的氣動阻力最小;X=0.45與X=0.6時，S803翼型產(chǎn)生的氣動阻力最小，S802氣動阻力略大于S803。

3.3 流場分析

由上可知，S803翼型的氣動升阻力性能優(yōu)于S802與DU93模型，為了進一步分析S802與DU93翼型結(jié)構(gòu)發(fā)生腐蝕缺陷時氣動性能變差原因，選取腐蝕缺陷位置為X=0.3與X=0.45，來流速度為30m/s與45m/s，其對應(yīng)的流場與壓強分布如圖8與圖9所示。

高速流體首先撞擊到翼型前緣，對翼型前緣造成較大的正壓強，然后沿著翼型的上表面和下表面向后流動。因下表面結(jié)構(gòu)完整沒有缺陷，流體沿著光滑的下表面向后流動時，流線較為平緩，對壁面產(chǎn)生的壓強較小且呈均勻分布。流體從翼型上表面向后流動，經(jīng)過腐蝕缺陷處時，在缺陷位置產(chǎn)生較強的渦旋，致使翼型腐蝕位置產(chǎn)生較大的負壓強，導(dǎo)致翼型腐蝕位置出現(xiàn)局部負壓場，嚴(yán)重影響翼型的升阻力，極大降低了翼型的氣動性能。

DU93翼型，X=0.3、v=30m/s時壓強變化范圍為-629～529Pa;X=0.3、v=45m/s時壓強變化范圍為-1 435～1 193Pa;X=0.45、v=30m/s時壓強變化范圍為-670～580Pa。S802翼型，X=0.3、v=30m/s時壓強變化范圍為-471～516Pa;X=0.3、v=45m/s時最大壓強變化范圍為-1 092～1 163Pa; X=0.45、 v=30m/s時壓強變化范圍為-452～515Pa。 來流速度越大， 流體對翼型造成的正負壓強越大， 相應(yīng)的升阻力也越大。 當(dāng)腐蝕位置從X=0.3后移到X=0.45，DU93翼型的正負壓強變大，升阻力數(shù)值波動變大，而S802翼型，流場流線較為平穩(wěn)，壓強變化不大，升阻力數(shù)值波動平緩。結(jié)果表明當(dāng)翼型發(fā)生腐蝕缺陷時，S802翼型的氣動性能整體優(yōu)于DU93翼型。

4 結(jié)論

本文在定常工況下對高風(fēng)速下腐蝕翼型的氣動性能進行研究，分析對比三種翼型在不同風(fēng)速、不同腐蝕位置下的氣動性能變化。發(fā)現(xiàn)翼型形狀的變化與腐蝕位置的改變都會造成翼型氣動性能的變化。

（1）翼型的升阻力大小與來流風(fēng)速成線性增長關(guān)系，流速越大，翼型產(chǎn)生的升力越大。

（2）隨著腐蝕缺陷位置的改變，S803翼型升阻力變化很小，氣動性能力變化較為平緩。S803翼型整體氣動性能最優(yōu)，而當(dāng)腐蝕位置為X=0.45與X=0.6時，S803升力比S802高1.2%左右，二者氣動性能差距很小。

（3）翼型發(fā)生腐蝕缺陷，高速流體流過翼型缺陷位置時，發(fā)生流場分離，在翼型缺陷位置形成渦旋流場，產(chǎn)生極強的的負壓流場域，最終導(dǎo)致翼型氣動性能降低。

參考文獻：

[1] 任年鑫，歐進萍.大型風(fēng)力機二維翼型氣動性能數(shù)值模擬[J].太陽能學(xué)報，2009，30（8）：1 087-1 091.

[2] BERTAGNOLIO F，RASMUSSEN F.A stochastic model for the simulation of wind turbine blades in static stall[J].Wind Energy， 2010， 13（4）：323-338.

[3] 李德順，李仁年，楊從新，等.粗糙度對風(fēng)力機翼型氣動性能影響的數(shù)值預(yù)測[J].農(nóng)業(yè)機械學(xué)報，2011，42（5）：111-115.

[4] 李德順，王成澤，李銀然，等.葉片前緣磨損形貌特征對風(fēng)力機翼型氣動性能的影響[J].農(nóng)業(yè)工程學(xué)報，2017，33（22）：269-275.

[5] SLOT H M ， GELINK E R M ， RENTROP C ， et al. Leading edge erosion of coated wind turbine blades： Review of coating life models[J]. Renewable Energy， 2015， 80：837-848.

[6] HAN W ， KIM J ， KIM B . Effects of contamination and erosion at the leading edge of blade tip airfoils on the annual energy production of wind turbines[J]. Renewable Energy， 2017， 115：817-823.

[7] 李爽，張磊，宋娟娟，等.風(fēng)力機翼型動態(tài)失速過程邊界層分離規(guī)律研究[J].工程熱物理學(xué)報，2019，40（9）：2 043-2 050.

[8] 郝文星，李春，陳福東，等.翼型附加氣動彈片流動控制研究[J].太陽能學(xué)報，2020，41（9）：302-309.

[9] 蔡新，林世發(fā)，胡莉，等.不同傾斜角葉尖小翼水平軸風(fēng)力機氣動性能[J].河海大學(xué)學(xué)報（自然科學(xué)版），2020，48（5）：460-466.

[10] WANG L ， WANG C ， SUN L . Study on Influence of Environmental Parameters on Dynamic Stall Characteristics of Wind Turbine Blades[J]. Journal of The Institution of Engineers （India） Series C， 2020，101：441-450.

[11] 蘇子昂，孫曉晶.S809翼型環(huán)量控制的射流參數(shù)研究[J].太陽能學(xué)報，2020，41（7）：359-366.

[12] 王龍，李亮，孫倫業(yè)，等.射流參數(shù)對風(fēng)力機葉片氣動性能的影響[J].流體機械，2017，45（7）：28-33.

[13] 廖書學(xué)，李春，聶佳斌，等.不同翼型對垂直軸風(fēng)力機性能的影響[J].機械設(shè)計與研究，2011，27（3）：108-111.

[14] CHURCHFIELD M J，LEE S，MICHALAKES J，et al. A numerical study of the effects of atmospheric and wake turbulence on wind turbine dynamics[J]. Journal of Turbulence， 2012， 13（14）：1-32.

[15] ZHONG J ， LI J . Aerodynamic performance prediction of NREL phase VI blade adopting biplane airfoil[J]. Energy，2020（206）：118182.

（責(zé)任編輯：李 麗）