考慮粗糙度敏感位置的鈍尾緣翼型氣動性能研究

張 旭，劉海龍，王格格，李 偉

?

考慮粗糙度敏感位置的鈍尾緣翼型氣動性能研究

張 旭1,2，劉海龍1，王格格1，李 偉2,3

（1. 天津工業大學天津市現代機電裝備技術重點實驗室，天津 300387；2. 建筑安全與環境國家重點實驗室，北京 100013；3. 天津城建大學能源與安全工程學院，天津 300384）

針對考慮粗糙度敏感位置的風力機翼型鈍尾緣改型前后的氣動性能進行研究，揭示鈍尾緣改型對表面粗糙翼型增升效果的影響規律。基于-SST湍流模型，計算表面光滑與粗糙的S822翼型的升、阻力系數，并與試驗結果進行比較；采用坐標旋轉變換與縮放橫縱坐標系數相結合的方法，建立鈍尾緣改型型線數學表達式，分析對稱鈍尾緣改型増升效果得到S822翼型的最佳尾緣厚度；研究吸力面和壓力面布置粗糙度時翼型的氣動性能，獲得上、下翼面的粗糙度敏感位置；對具有粗糙度敏感位置的翼型按最佳尾緣厚度進行鈍尾緣改型，計算改型前后翼型的升、阻力系數和升阻比，并分析尖、鈍尾緣翼型的粗糙度敏感性。結果表明：翼型進行鈍尾緣改型的最佳尾緣厚度為2%弦長；吸力面和壓力面的粗糙度敏感位置分別為距前緣1%弦長和5%弦長處；鈍尾緣改型使升力系數和最大升阻比均明顯升高，顯著改善了表面粗糙翼型的氣動性能，且尖、鈍尾緣翼型的粗糙度敏感性綜合指標值為10.68%和8.15%，降低了翼型對粗糙度位置的敏感性。研究結論可為表面粗糙風力機葉片翼型的設計和優化提供指導。

風力機；翼型；粗糙度；敏感位置；鈍尾緣改型；氣動性能

0 引 言

安裝于高寒、沿海地帶的風力機，經常遇到灰塵、昆蟲以及雨雪等，葉片表面會出現污垢附著現象。污垢會改變葉片氣動外形并增大表面粗糙度，從而導致風能轉化率降低[1-8]。而鈍尾緣改型能增加最大升力和失速攻角，降低最大升力對前緣粗糙的敏感度[9-11]。因此，研究表面粗糙翼型鈍尾緣修改后的氣動性能對提高風力機的風能利用系數具有重要意義。

針對粗糙度對翼型氣動性能的影響，很多學者開展了一系列數值與試驗研究。Khalfallah等[12]實驗分析葉片表面積灰形成的粗糙度對風力機出力的影響，研究表明隨灰塵的積累，升力下降，阻力升高，且積灰降低了風力機的輸出功率。Soltani等[13]通過風洞試驗、Ferrer等[14]利用計算流體動力學（computational fluid dynamics，CFD）方法以及Sagol等[1]采用數值模擬和風洞試驗相結合，均獲得表面粗糙度降低翼型的氣動性能以及翼型前緣對粗糙度最敏感等結論。包能勝等[15]利用粗糙帶貼片試驗，研究壓力面前緣粗糙度對葉型氣動性能的影響，結果表明升力系數在失速之后變化非常大，阻力系數變化很小。陳進等[16]探討了前緣粗糙條件對翼型氣動特性的影響，并優化設計得到一種風力機專用新翼型。李德順等[17]數值研究粗糙度對風力機翼型氣動性能的影響，結果顯示適當增加壓力面后緣粗糙度可提高升力系數，且升阻比始終低于光滑表面的翼型。但吳攀等[18]研究發現該結論并不適用于FFA-W3-211翼型。張駿等[19]采用CFD方法研究風力機積灰和昆蟲引起的粗糙度效應，結果表明前緣粗糙度嚴重影響翼型的氣動性能，壓力面尾緣布置適當粗糙度卻有積極作用。以上文獻并未涉及通過降低粗糙度敏感性提高翼型氣動性能的研究。

此外，國內外學者也在風力機鈍尾緣翼型氣動性能的研究方面取得了一些進展。Baker等[20]用風洞試驗研究鈍尾緣翼型的氣動性能，發現適度增加尾緣厚度可增大升阻比并減小對前緣過渡的敏感性。Standish等[21]采用4種數值解法對鈍尾緣翼型進行數值模擬，研究表明鈍尾緣改型能增加最大升力和升力曲線斜率，降低過早的邊界層轉捩。韓中合等[22]通過改進翼型尾部的Gurney襟翼形成鈍尾緣翼型，分析翼型修改前后風力機的功率輸出特性。劉杰平等[23]利用XFOIL軟件比較不同尾緣加厚方式所得翼型的氣動性能。楊瑞等[24]采用CFD方法模擬薄、鈍尾緣翼型的氣動性能，結果顯示鈍尾緣翼型增大了最大升力系數，并降低前緣污染對升力特性的影響。馬林靜等[25]參數化研究尾緣厚度對風力機翼型氣動性能的影響，發現尾緣厚度在一定范圍內增大時，升、阻力系數增大，升阻比先增后降。徐浩然等[26]數值模擬尾緣對稱加厚翼型的氣動性能，結果表明最大升力系數隨尾緣厚度增大而增大，但厚度應控制在約5%弦長范圍內。然而，以上關于粗糙度和鈍尾緣改型對翼型氣動性能影響的研究多是單獨開展的。針對表面粗糙的翼型進行尾緣對稱加厚，研究考慮粗糙度敏感位置的鈍尾緣改型對翼型氣動性能的影響較少。

本文以美國可再生能源實驗室設計的風力機專用翼型S822為研究對象，建立鈍尾緣改型型線公式并進行改型。考慮鈍尾緣的最佳厚度以及粗糙度敏感位置，采用CFD方法對表面粗糙翼型鈍尾緣修改后的氣動性能進行詳細研究，并分析尖、鈍尾緣翼型的粗糙度敏感性。

1 計算模型及方法

基于CFD軟件Fluent，計算表面光滑與粗糙的S822翼型氣動性能，并與試驗結果進行比較分析。

1.1 幾何模型

S822翼型廣泛應用于風力機葉片的主要功率產生區，具有16%的最大相對厚度，以及最大相對彎度為1.92%。風洞試驗在美國伊利諾伊大學香檳分校的低紊流亞音速風洞進行，雷諾數=5×105，風速24.38 m/s，翼型弦長0.305 m。使用圖1所示Zigzag粗糙帶，分別布置在上、下翼面距前緣2%和5%位置處[27]。粗糙帶厚0.33 mm，沿翼展方向相鄰鋸齒齒尖距離4.98 mm，齒尖角度83°，由此可計算出沿氣流方向相鄰鋸齒齒尖距離為3 mm。

圖1 Zigzag粗糙帶幾何尺寸

為了較好模擬風洞試驗，并與數據進行對比，表面光滑與粗糙的S822翼型幾何模型的弦長同樣選用 0.305 m。采用凸臺描述試驗中的粗糙帶，建立翼型粗糙度模型。凸臺高0.33 mm（即粗糙帶厚），寬3 mm（即沿氣流方向相鄰鋸齒齒尖距離）。

利用Fluent的前處理軟件Gambit生成翼型的幾何模型、計算域和網格，如圖2所示。計算域由一個直徑為50的半圓形和一個長50、寬25的矩形構成，翼型位于半圓中心。表面光滑的S822翼型計算網格采用C型四邊形網格。在光滑翼型上布置320個節點，采用邊界層進行近壁面加密。邊界層首層網格高度10-5m，共劃分15層，壁面率值小于5。表面布置粗糙度的S822翼型幾何形狀較復雜，采用分塊結構化網格處理凸臺區域。沿凸臺側邊引出與半圓邊界相交的線段，將表面粗糙翼型的半圓形計算域分成6部分，每部分均構成四邊形。使用映射網格進行劃分，得到結構化C型網格。在含粗糙度的翼型上布置380個節點，且邊界層同表面光滑的S822翼型，值也小于5。

注：A和B為光滑和粗糙翼型的網格局部放大圖。

1.2 控制方程

風力機翼型繞流為不可壓縮流動[28]，其控制方程為連續性方程

和二維不可壓縮Navier-Stokes方程

（2）

1.3 計算方法

使用Fluent軟件進行翼型的氣動性能數值計算。-SST湍流模型綜合了-與-模型的優勢，適合處理帶有逆壓梯度和分離流動的問題[29]。壓力和速度耦合采用SIMPLE算法，各方程離散格式均采用二階迎風格式。速度和連續性方程的殘差值分別為10-6和10-4，和方程的為10-5。計算域左側及上下兩側進流邊界設定為速度入口，風速由=5×105確定，湍流強度為0.05%。計算域右側出流邊界設定為壓力出口，表壓力給定0 Pa。翼型表面設定為固壁絕熱無滑移邊界條件。

1.4 網格無關性與模型適應性驗證

利用上述計算模型及方法，計算攻角在-5.17o～20o之間變化時表面光滑和粗糙的S822翼型的升、阻力系數，并與文獻[27]的風洞試驗結果進行對比，如圖3所示。表面光滑和粗糙的S822翼型的網格數分別采用74 990、 95 185、112 740和79 214、109 046、128 546。

圖3 光滑和粗糙的S822翼型計算值與試驗值比較

由圖3可知，表面光滑的S822翼型的升、阻力系數在74 990、95 185和112 740網格數下均相差很小，表面粗糙的S822翼型在79 214、109 046和28 546網格數下亦如此。考慮到節省計算資源，光滑和粗糙的翼型的網格數分別采用95 185和109 046比較合適，均可得到精確的計算結果。

由圖3a可知，光滑翼型的升、阻力系數分別在8.19°和10.19°攻角之前與試驗值吻合良好，在8.19°和10.19°攻角之后高于、低于試驗值。由圖3b可知，粗糙翼型的升力系數在攻角小于4.06°時與試驗值比較接近，在4.06°攻角之后略高于試驗值；阻力系數在11.19°攻角之前與試驗值吻合很好。總體來看，表面光滑和粗糙的翼型的氣動性能計算值均與試驗值基本相同，且總體變化趨勢一致，故可采用本文數值模型和方法進行模擬計算。

2 鈍尾緣改型公式及最佳尾緣厚度

建立鈍尾緣改型型線的數學表達式，數值計算S822翼型的不同尾緣厚度鈍尾緣改型的氣動性能。分析尾緣厚度對翼型的升、阻力系數和升阻比的影響，獲得鈍尾緣改型時最佳的尾緣厚度。

2.1 鈍尾緣改型公式

原始S822翼型的型線如圖4a所示。在不改變翼型的最大相對厚度及其位置、彎度和弦長的情況下，利用坐標旋轉變換以及縮放橫縱坐標系數建立圖4a所示鈍尾緣改型型線的數學表達式。

a. S822翼型鈍尾緣改型示意圖

a. Curve design of S822 airfoil and its blunt trailing-edge modification

注：：虛線為經過坐標旋轉變換得到的翼型型線；r為的長度，m；α為與軸的夾角，rad；β和φ為逆、順時針旋轉角度，rad；h為尾緣厚度，m；、、為尖、鈍尾緣和坐標旋轉后的翼型控制點坐標；S822_1、S822_2、S822_3、S822_4為尾緣厚度為1%c、2%c、3%c和4%c的對稱鈍尾緣翼型；c為翼型弦長，m。

（5）

將式（3）帶入式（4）和式（5），可得

（7）

經過坐標旋轉變換，翼型的弦長變短。此時，可通過上、下翼面型線的橫坐標分別乘以因式和保證改型后翼型弦長不變，其中和為上、下翼面旋轉后尾緣點的橫坐標。因此，鈍尾緣改型上、下翼面橫坐標的表達式為

為使鈍尾緣改型與原翼型具有相同的最大相對厚度及其位置，基于有限元方法中形函數的思想，將坐標旋轉后翼型上、下翼面的縱坐標分別減去和加上，從而得到鈍尾緣改型上、下翼面縱坐標的表達式

（9）

采用式（8）和式（9）對S822翼型進行尾緣厚度為1%、2%、3%和4%的對稱鈍尾緣改型，翼型型線如圖4b所示。并用S822_1、S822_2、S822_3和S822_4分別表示上述4種尾緣厚度情況。

2.2 最佳尾緣厚度

圖5為S822翼型及其鈍尾緣改型的升、阻力系數和升阻比隨攻角變化的曲線。由圖5a可知，4種改型的升力系數在0.03°攻角之前比較接近，0.03°攻角之后隨尾緣厚度增加呈先增大后減小的趨勢，S822_1、S822_2和S822_3改型的升力系數均高于原型和S822_4改型，且S822_4改型的升力系數在5.09°攻角之后低于S822翼型。由圖5b可知，阻力系數在-5.17o～20o攻角范圍內十分接近。

圖5 不同尾緣厚度翼型的升、阻力系數和升阻比

由圖5c可知，4種改型的升阻比在1°攻角之前很接近，且均高于S822翼型；在1°～8.19°攻角范圍內，升阻比隨尾緣厚度增加而先增后減，S822_2改型的最大。在8.19°攻角之后，S822_1、S822_2和S822_3改型的升阻比與原型相差不大，且S822_3改型在9.17°～12.22°攻角之間高于原型及其他改型。S822_4改型的升阻比在3.07°攻角之后低于其他翼型。S822_2改型的最大升阻比高于原型及其他改型。綜上所述，鈍尾緣改型時，并不是尾緣厚度越大翼型的氣動性能越好。S822翼型對稱鈍尾緣改型的最佳尾緣厚度為2%。

3 粗糙度位置對翼型氣動性能的影響

為了研究粗糙度位置對翼型氣動性能的影響，分別在翼型吸力面與壓力面的距前緣1%、2%、5%、10%、20%、30%、40%、50%、60%、90%、100%位置處布置一個高0.33 mm、寬3 mm的凸臺，進行翼型的升、阻力系數和升阻比計算。

3.1 吸力面布置粗糙度翼型的氣動性能

圖6為吸力面布置粗糙度的翼型的升、阻力系數和升阻比隨粗糙度位置變化的曲線。由圖6a和圖6b可知，升、阻力系數分別在7.16°和10.19°攻角之前非常接近。S822_1%_s、S822_2%_s和S822_5%_s翼型的升力系數在7.16°攻角之后隨粗糙度位置靠近前緣而減小并低于光滑翼型，而阻力系數在10.19°攻角之后呈遞增趨勢且高于光滑翼型。吸力面其他位置粗糙的翼型的升力系數在7.16°～16°攻角之間和阻力系數均與光滑翼型十分接近，且升力系數在攻角大于16°之后明顯高于光滑翼型。

注：S822_1%c_s、S822_2%c_s、S822_5%c_s、S822_10%c_s、S822_20%c_s、S822_30%c_s、S822_40%c_s、S822_50%c_s、S822_60%c_s為吸力面1%c、2%c、5%c、10%c、20%c、30%c、40%c、50%c、60%c位置處布置粗糙度的翼型。

由圖6c可知，S822_1%_s、S822_2%_s和S822_ 5%_s翼型的升阻比在2°攻角之前相差不大，且均高于 光滑翼型，2°攻角之后隨粗糙度位置靠近前緣而減小，并低于光滑翼型。吸力面其他位置粗糙的翼型的升阻比在6.17°攻角之前和12.22°攻角之后均與光滑翼型比較接近，6.17°～12.22°攻角之間低于光滑翼型。最大升阻比隨吸力面粗糙度位置靠近前緣呈遞減趨勢，且低于光滑翼型。

3.2 壓力面布置粗糙度翼型的氣動性能

圖7為壓力面布置粗糙度的翼型的升、阻力系數和升阻比隨粗糙度位置變化的曲線。由圖7a和圖7c可知，S822_1%_p、S822_2%_p、S822_5%_p和S822_100%_ p翼型的升力系數和升阻比均在3.07°攻角之前非常接近，且高于光滑翼型；壓力面其他位置粗糙的翼型的升力系數和升阻比在攻角小于18°時與光滑翼型十分接近，18°攻角之后明顯高于光滑翼型。S822_100%_p翼型的升力系數在3.07°～6.17°和10.19°～16°攻角范圍內以及升阻比在3.07°～6.17°攻角之間高于壓力面其他位置粗糙和光滑的翼型。S822_5%_p翼型的升力系數在9.17°～18°攻角之間以及升阻比在9.17°～13.3°攻角范圍內均最小。壓力面粗糙的翼型的升阻比在7.16°～9.17°和13.3°～18°攻角范圍內分別低于和非常接近光滑翼型，最大升阻比相差不大且低于光滑翼型。由圖7b可知，阻力系數在18°攻角之前十分接近，18°攻角之后低于光滑翼型。

注：S822_1%c_p、S822_2%c_p、S822_5%c_p、S822_10%c_p、S822_20%c_p、S822_30%c_p、S822_40%c_p、S822_50%c_p、S822_60%c_p、S822_90%c_p、S822_100%c_p為壓力面1%c、2%c、5%c、10%c、20%c、30%c、40%c、50%c、60%c、90%c、100%c位置處布置粗糙度的翼型。

圖7a和圖7c描述壓力面粗糙翼型的升力系數和升阻比隨攻角變化的規律，并未清晰顯示同一攻角下不同粗糙度位置對翼型氣動性能的影響。因此，以9.17°～12.22°的4個攻角為例，分析升力系數和升阻比隨粗糙度位置變化的情況，如圖7d和圖7e所示。由圖7d和圖7e可知，升力系數隨攻角增加呈遞增趨勢，升阻比卻呈遞減趨勢。升力系數和升阻比在粗糙度距前緣小于1%時快速減小，1%～5%范圍內基本呈先增大后減小的趨勢，5%位置處達到最小值。升力系數在粗糙度距前緣10%～50%之間時相差不大，50%～100%范圍內呈先減小后增大的趨勢。升阻比在粗糙度距前緣10%～50%范圍內基本呈先增大后減小的趨勢，60%～100%之間時比較接近，且低于光滑翼型。

分析可知，吸力面距前緣1%位置處布置粗糙度后，升、阻力系數和升阻比變化最劇烈，因而吸力面粗糙度敏感位置為距前緣1%處。壓力面距前緣5%和100%位置處布置粗糙度后，升、阻力系數和升阻比變化都比較明顯。但實際上翼型的前緣較易污染，研究前緣粗糙對翼型氣動性能的影響更具有意義。因此，壓力面粗糙度的敏感位置為距前緣5%處。

4 表面粗糙翼型改型后氣動性能計算及粗糙度敏感性分析

利用建立的鈍尾緣改型公式對表面粗糙的S822翼型進行尾緣厚度為2%的對稱鈍尾緣改型，研究考慮粗糙度敏感位置的鈍尾緣改型對翼型的升、阻力系數和升阻比的影響規律，并進行尖、鈍尾緣翼型的粗糙度敏感性分析。

4.1 鈍尾緣改型對吸力面粗糙翼型氣動性能的影響

在吸力面1%位置處布置粗糙度，翼型鈍尾緣改型前后的升、阻力系數和升阻比隨攻角變化的曲線如圖8a、8b、8c所示。由圖8a可知，相較于尖尾緣S822_1%_s翼型，S822_1%_s_2改型的升力系數明顯升高。由圖8b可知，改型的阻力系數在9.17°攻角之前與S822_1%_s翼型非常接近，9.17°攻角之后高于S822_1%_s翼型。由圖8c可知，S822_1%_s_2改型的升阻比在攻角小于11.19°時高于S822_1%_s翼型，11.19°攻角之后低于S822_1%_s翼型。吸力面粗糙翼型的最大升阻比在鈍尾緣改型后升高14.6%。

4.2 鈍尾緣改型對壓力面粗糙翼型氣動性能的影響

在壓力面5%位置處布置粗糙度，翼型鈍尾緣改型前后的升、阻力系數和升阻比隨攻角變化的曲線如圖8d、8e、8f所示。由圖8d和圖8e可知，S822_5%_p_2改型的升、阻力系數分別在0.03°和10.19°攻角之前與S822_5%_p翼型非常接近，0.03°和10.19°攻角之后高于S822_5%_p翼型，而20°攻角時改型的升力系數低于S822_5%_p翼型。

由圖8f可知，攻角小于1°以及在13.23°～16°之間變化時，S822_5%_p_2改型和S822_5%_p翼型的升阻比十分接近；在1°～13.23°攻角范圍內和16°攻角之后，改型的升阻比分別高于、低于S822_5%_p翼型。S822_5%_ p_2改型的最大升阻比明顯高于S822_5%_p翼型。

4.3 鈍尾緣改型對表面粗糙翼型氣動性能的影響

在吸力面1%和壓力面5%兩個位置處布置粗糙度，翼型鈍尾緣改型前后的升、阻力系數和升阻比隨攻角變化的曲線如圖9所示。

圖8 吸力面或壓力面粗糙翼型鈍尾緣（尾緣厚度為2%c）改型前后的升、阻力系數和升阻比

圖9 粗糙翼型鈍尾緣改型前后的升、阻力系數和升阻比

由圖9a可知，S822_1%_s & 5%_p_2改型的升力系數高于S822_1%_s & 5%_p翼型。由圖9b可知，改型的阻力系數在10.19°攻角之前與S822_1%_s & 5%_p翼型非常接近，10.19°攻角之后高于S822_1%_s & 5%_p翼型。由圖9c可知，S822_1%_s & 5%_p_2改型的升阻比在11.19°攻角之前高于S822_1%_s & 5%_p翼型，11.19°攻角之后低于S822_1%_s & 5%_p翼型。鈍尾緣改型后表面粗糙翼型的最大升阻比提升16.9%。

4.4 鈍尾緣翼型粗糙度敏感性分析

為了定量地評價鈍尾緣改型對翼型粗糙度敏感性的影響，需要找出一個評價指標。根據風力機葉片翼型的特點，取升力系數下降率和升阻比下降率的加權平均值作為粗糙度敏感性評價指標[30]。

定義升力系數下降率

定義升阻比下降率

根據文獻[30]可知，升力系數下降率和升阻比下降率的權重系數分別為0.49和0.51，即綜合指標。根據表面光滑和粗糙的尖、鈍尾緣翼型的升力系數與升阻比計算結果，計算升力系數下降率、升阻比下降率以及綜合指標值，如表1所示。由表1可以看出，鈍尾緣翼型綜合指標值為8.15%，低于尖尾緣翼型的10.68%，說明鈍尾緣改型可使翼型具有更低的粗糙度敏感性。

表1 尖、鈍尾緣翼型粗糙度敏感性綜合指標值

注： CC′、(C/C)max、(C′/C′)max為光滑和粗糙翼型的升力系數、最大升阻比，δC、(C/C)為升力系數和升阻比的下降率。

Note:C, C′, (C/C)maxand (C′/C′)maxare the lift coefficients and the maximum lift-drag ratios of airfoils with smooth and rough surfaces;δC,(C/C) are the decline rates of lift coefficient and lift-drag ratio.

5 結 論

通過數值模擬研究，關于考慮粗糙度敏感位置的鈍尾緣改型的增升效果可得到如下結論：

1）通過坐標旋轉變換與縮放坐標系數，建立鈍尾緣改型型線表達式；鈍尾緣翼型的升力系數和升阻比隨尾緣厚度增加而先增后減，且最大升阻比在尾緣厚度為2%（為翼型弦長）時最高。

2）吸力面1%、2%和5%位置處粗糙的翼型的升力系數隨粗糙度位置靠近前緣而減小，且與升阻比均低于光滑翼型，阻力系數增大并高于光滑翼型；其他位置粗糙翼型的升、阻力系數與光滑翼型很接近。壓力面5%位置處粗糙的翼型的升力系數在9.17°～18°攻角之間低于其他位置粗糙和光滑的翼型；壓力面粗糙翼型的最大升阻比相差不大且低于光滑翼型。吸、壓力面的粗糙度敏感位置分別為1%和5%處。

3）具有粗糙度敏感位置的翼型鈍尾緣改型后，升力系數和最大升阻比均明顯升高；尖、鈍尾緣翼型的粗糙度敏感性綜合指標值分別為10.68%和8.15%，鈍尾緣改型使翼型對粗糙度位置的敏感性降低。

[1] Sagol E, Reggio M, Ilinca A. Issues concerning roughness on wind turbine blades[J]. Renewable and Sustainable Energy Reviews, 2013, 23(4): 514－525.

[2] Corten G P, Veldkamp H F. Insects cause double stall[C]//European Wind Energy Conference. Copenhagen, Denmark: 1999.

[3] Corten G P, Veldkamp H F. Insects can halve wind-turbine power[J]. Nature, 2001, 412(6842): 41－42.

[4] Morgan C, Bossanyi E, Seifert H. Assessment of safety risks arising from wind turbine icing [C]//BOREAS IV-Wind Energy Production in Cold Climate. Hetta, Finland: 1998: 113－121.

[5] Fu P, Farzaneh M. A CFD approach for modeling the rime-ice accretion process on a horizontal-axis wind turbine[J]. Journal of Wing Engineering and Industrial Aerodynamics, 2010, 98(4/5): 181－188.

[6] Lamraoui F, Fortin G, Benoit R, et al. Atmospheric icing impact on wind turbine production[J]. Cold Regions Science and Technology, 2014, 100(4): 36－49.

[7] Shu L C, Liang J, Hu Q, et al.Study on small wind turbine icing and its performance[J]. Cold Regions Science and Technology, 2017, 134(1): 11－19.

[8] 焦靈燕，汪建文，賀玲麗. 粗糙度對風力機翼型氣動性能影響的模擬研究[J]. 可再生能源，2014，32(12)：1816－1820. Jiao Lingyan, Wang Jianwen, He Lingli. Simulation study on effect of surface roughness on aerodynamic performance of wind turbine airfoil[J]. Renewable Energy Resources, 2014, 32(12): 1816－1820. (in Chinese with English abstract)

[9] Jackson K J, Zuteck M D, van Dam C P, et al. Innovative design approaches for large wind turbine blades[J]. Office of Scientific & Technical Information Technical Reports, 2003, 8(2): 141－171.

[10] Zhang X, Li W, Liu H L. Numerical simulation of the effect of relative thickness on aerodynamic performance improvement of asymmetrical blunt trailing-edge modification[J]. Renewable Energy, 2015, 80: 489－497.

[11] Yoo H S, Lee J C. Numerical analysis of NACA64-418 airfoil with blunt trailing edge[J]. International Journal of Aeronautical and Space Sciences, 2015, 16(4): 493－499.

[12] Khalfallah M G, Koliub A M. Effect of dust on the performance of wind turbines[J]. Desalination, 2007, 209(1): 209－220.

[13] Soltani M R, Birjandi A H, Moorani M S. Effect of surface contamination on the performance of a section of a wind turbine blade[J]. Scientia Iranica, 2011, 18(3): 349－357.

[14] Ferrer E, Munduate X. CFD predictions of transition and distributed roughness over a wind turbine airfoil[C]//47thAIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition. Orlando, Florida: 2009, 1－13.

[15] 包能勝，倪維斗. 風力機翼型前緣表面粗糙度對氣動性能影響[J]. 太陽能學報，2008，29(12)：1465－1470. Bao Nengsheng, Ni Weidou. Influence of additional rough strap of wind turbine airfoil leading edge surface on aerodynamic performance[J]. Acta Energiae Solaris Sinica, 2008, 29(12): 1465－1470. (in Chinese with English abstract)

[16] 陳進，張石強，王旭東，等. 基于粗糙度敏感性研究的風力機專用翼型設計[J]. 空氣動力學學報，2011，29(2)： 142－149. Chen Jin, Zhang Shiqiang, Wang Xudong, et al. Dedicated wind turbine airfoil design based on the roughness sensitivity considerations[J]. Acta Aerodynamica Sinica, 2011, 29(2): 142－149. (in Chinese with English abstract)

[17] 李德順, 李仁年, 楊從新, 等. 粗糙度對風力機翼型氣動性能影響的數值預測[J]. 農業機械學報, 2011, 42(5): 111－115. Li Deshun, Li Rennian, Yang Congxin, et al. Numerical prediction of the effect of surface roughness on aerodynamic performance of a wind turbine airfoil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(5): 111－115. (in Chinese with English abstract)

[18] 吳攀, 李春, 葉舟, 等. 粗糙度對風力機專用翼型氣動性能影響[J]. 流體機械, 2014, 42(1): 17－21, 62. Wu Pan, Li Chun, Ye Zhou, et al. Influence of roughness on aerodynamic performance of dedicated wind turbine airfoil [J]. Fluid Machinery, 2014, 42(1): 17－21, 62. (in Chinese with English abstract)

[19] 張駿, 袁奇, 吳聰, 等. 大型風力機葉片表面粗糙度效應數值研究[J]. 中國電機工程學報, 2014, 34(20): 3384－3391. Zhang Jun, Yuan Qi, Wu Cong, et al. Numerical simulation on the effect of surface roughness for large wind turbine blades[J]. Proceedings of the CSEE, 2014, 34(20): 3384－3391. (in Chinese with English abstract)

[20] Baker J P, Mayda E A, van Dam C P. Experimental analysis of thick blunt trailing-edge wind turbine airfoils [J]. Journal of Solar Energy Engineering, 2006, 128(4): 422－431.

[21] Standish K J, van Dam C P. Aerodynamic analysis of blunt trailing edge airfoils [J]. Journal of Solar Energy Engineering, 2003, 125(4): 479－487.

[22] 韓中合，焦紅瑞. 加裝鈍尾緣改善風力機槳葉氣動性能的研究[J]. 動力工程，2009，29(11)：1073－1077.Han Zhonghe, Jiao Hongrui. Study on aerodynamic performance of airfoil by mounting a blunt trailing edge to wind blade[J]. Journal of Power Engineering, 2009, 29(11): 1073－1077. (in Chinese with English abstract)

[23] 劉杰平，陳培，張衛民. 后緣加厚方式對典型風力機翼型氣動性能的影響[J]. 太陽能學報，2009，30(8)：1092－1096. Liu Jieping, Chen Pei, Zhang Weimin. Influence of enlarging trailing edge thickness on the representative wind turbine airfoil’s aerodynamic performances[J]. Acta Energiae Solaris Sinica, 2009, 30(8): 1092－1096. (in Chinese with English abstract)

[24] 楊瑞，李仁年，張士昂，等. 鈍尾緣風力機翼型氣動性能計算分析[J]. 機械工程學報，2010，46(2)：106－110. Yang Rui, Li Rennian, Zhang Shiang, et al. Computational analyses on aerodynamic characteristics of flatback wind turbine airfoils [J]. Journal of Mechanical Engineering, 2010, 46(2): 106－110. (in Chinese with English abstract)

[25] 馬林靜，陳江，杜剛，等. 尾緣厚度對風力機翼型氣動特性影響參數化研究[J]. 太陽能學報，2010，31(8)：1060－1067. Ma Linjing, Chen Jiang, Du Gang, et al. Parametric research on influence of trailing edge’s thickness to aerodynamic performance for wind turbine airfoils[J]. Acta Energiae Solaris Sinica, 2010, 31(8): 1060－1067. (in Chinese with English abstract)

[26] 徐浩然，楊華，劉超. 尾緣加厚的DU系列翼型氣動性能數值分析[J]. 農業工程學報，2014，30(17)：101－108.Xu Haoran, Yang Hua, Liu Chao. Numerical value analysis on aerodynamic performance of DU series airfoils with thickened trailing edge[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(17): 101－108. (in Chinese with English abstract)

[27] Selig M S, Mcgranahan B D. Wind tunnel aerodynamic tests of six airfoils for use on small wind turbines [R]. Golden, Colorado: National Renewable Energy Laboratory, 2003.

[28] 楊從新，巫發明，張玉良. 基于滑移網格的垂直軸風力機非定常數值模擬[J]. 農業機械學報，2009，40(6)：98－102. Yang Congxin, Wu Faming, Zhang Yuliang. Numerical simulation on unsteady rotated flow of a vertical axis wind turbine based on moving meshes[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(6): 98－102. (in Chinese with English abstract)

[29] Menter F R. Two-equation eddy-viscosity turbulence models for engineering applications [J]. AIAA Journal, 1994, 32(8): 1598－1605.

[30] 黃宸武，楊科，劉強，等. 風力機葉片21%相對厚度翼型粗糙敏感性研究[J]. 工程熱物理學報，2012，33(6)： 953－956. Huang Chenwu, Yang Ke, Liu Qiang, et al. Investigation on roughness sensitivity for 21% relative thichness airfoil of wind turbine blade[J]. Journal of Engineering Thermophysics, 2012, 33(6): 953－956. (in Chinese with English abstract)

Aerodynamic performance of blunt trailing-edge airfoil considering roughness sensitivity position

Zhang Xu1,2, Liu Hailong1, Wang Gege1, Li Wei2,3

(1.300387,; 2.100013,; 3.300384,)

Wind turbine is often exposed to dramatically different operational conditions, from icy environments to deserts with sand storms, and there are contaminants in these environments, like dust, dirt, ice, and even insects. These contaminants change the aerodynamic shape of blade and increase the surface roughness, which results in the lower utilization rate of wind energy.The aerodynamic performance of wind turbine blade can be improved through the airfoil modification, so the blunt trailing-edge structure is adopted during the design of an airfoil.Compared with the original airfoil, the blunt trailing-edge modification with bigger trailing-edge thickness and cross-section area not only has a great improvement in the maximum lift coefficient and the stall angle of attack, but also makes the maximum lift less sensitive to the leading-edge roughness.Therefore, it is of great significance to study the aerodynamic performance of blunt trailing-edge modification of the airfoil with rough surface for the improvement of the power utilization coefficient of wind turbine. The aerodynamic performance of the airfoil with rough surface and the blunt trailing-edge modification have been numerically and experimentally investigated in recent years. However, these 2 problems have been discussed separately, and the effects of the blunt trailing-edge modificationon the aerodynamic performance improvement of wind turbine airfoil have been less investigated considering the roughness sensitivity position.In the present study, the aerodynamic performance of wind turbine airfoil and its blunt trailing-edge modification considering the roughness sensitivity position was numerically investigated to reveal the effect of the blunt trailing-edge modification on the lift enhancement of airfoil with rough surface. The dedicated wind turbine airfoil S822 from National Renewable Energy Laboratory (NREL) was used for the simulation. The lift and drag coefficients of S822 airfoil with smooth or rough surfaces were calculated by theSST turbulence model, and were compared with the aerodynamic data from wind tunnel tests, which offered a good opportunity to examine the capability of CFD (computational fluid dynamics) simulation. The mathematical expression of the blunt trailing-edge airfoil profile was established using the coordinates’ rotation combined with the zoom coefficient of coordinates, and the airfoil S822 was modified to be symmetrical blunt trailing-edge airfoil. The lift enhancement of modified airfoils was analyzed to get the best trailing-edge thickness. In order to obtain the roughness sensitivity position of suction and pressure surfaces, the aerodynamic performance of the airfoil with rough surface was studied. The lift and drag coefficients and the lift-drag ratio were calculated for the airfoils with the roughness sensitivity position and their symmetrical modifications with the best trailing-edge thickness. Andthe roughness sensitivity of sharp and blunt trailing-edge airfoils was also analyzed. The results indicated that the best trailing-edge thickness was 2% of chord length for symmetrical blunt trailing-edge airfoil. The roughness sensitivity positions of suction and pressure surfaces were 1% and 5% of chord length away from the leading-edge, respectively. After the blunt trailing-edge modification, the lift coefficient and the maximum lift-drag ratio of the airfoil with the roughness sensitivity position significantly increased. The lift-drag ratio of the blunt trailing-edge airfoil was higher than that of the original airfoilfor the angle of attack less than 11.19° when the suction surface of airfoil is rough, and so does the airfoil with rough suction and pressure surfaces. It is the same change ruler as abovefor the airfoil with rough pressure surface at different angles of attack ranging from 1° to 13.23°. The blunt trailing-edge modification makes the lift coefficient and the maximum lift-drag ratio significantly increase, which remarkably improves the aerodynamic performance of rough airfoil. The compositive index of the roughness sensitivity was 10.68% and 8.15% for sharp and blunt trailing-edge airfoils, respectively. The modification reduces the airfoil’s sensitivity to the roughness position. The research provides a significant guidance for designing and optimizing the wind turbine airfoil under rough blade surface conditions.

wind turbines; airfoils; roughness; sensitivity position; blunt trailing-edge modification; aerodynamic performance

10.11975/j.issn.1002-6819.2017.08.011

TK83

A

1002-6819(2017)-08-0082-08

2016-08-10

2017-04-17

國家自然科學基金（11402168）；建筑安全與環境國家重點實驗室開放基金項目（BSBE2015-03，BSBE2014-08）；天津市自然科學基金面上項目（17JCYBJC20800，15JCYBJC48600）

張 旭，女，河北安平人，副教授，博士后，主要從事風力機空氣動力學和結構力學研究。天津 天津工業大學天津市現代機電裝備技術重點實驗室，300387。Email：zhangxu@tjpu.edu.cn

張 旭，劉海龍，王格格，李 偉.考慮粗糙度敏感位置的鈍尾緣翼型氣動性能研究[J]. 農業工程學報，2017，33(8)：82－89. doi：10.11975/j.issn.1002-6819.2017.08.011 http://www.tcsae.org

Zhang Xu, Liu Hailong, Wang Gege, Li Wei. Aerodynamic performance of blunt trailing-edge airfoil considering roughness sensitivity position[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(8): 82－89. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.08.011 http://www.tcsae.org