吳友健 楊 藝,2

（1.廣東海洋大學 機械與動力工程學院；2.南方海洋科學與工程廣東省實驗室）

0 引言

化石能源的過度消耗使得全球氣候變暖問題越來越嚴重，世界各國都在尋找清潔的可再生能源來實現(xiàn)可持續(xù)發(fā)展。風能作為一種環(huán)保能源受到大家的重視[1-2]，風力發(fā)電成為一種新興的發(fā)電技術(shù)，許多國家都積極開展風力發(fā)電機的研究工作[3]。風力機有垂直軸和水平軸兩種[4]，其中升力型垂直軸風力機由于具有無需偏航機構(gòu)可適應(yīng)全向風、造價低、結(jié)構(gòu)簡單[5-6]等優(yōu)點而被受到重視。其工作原理是當來流過升力型葉片時會產(chǎn)生升力，升力沿圓周運動的切線方向分力形成的力矩推動風輪轉(zhuǎn)動從而帶動發(fā)電機工作。葉片是風力發(fā)電機的重要組成部分，葉片翼型氣動性能的優(yōu)劣對風力發(fā)電機的發(fā)電效率有很大的影響，因此對影響翼型氣動特性的因素進行深入研究是十分必要的。

國內(nèi)外學者對此做了相關(guān)研究。呂黎等人采用XFLR5程序來探究NACA4409和NACA0015翼型的氣動性能差異，結(jié)果發(fā)現(xiàn)對稱翼型NACA0015在正負攻角范圍內(nèi)都具有穩(wěn)定的氣動性能[7]。王旱祥等人發(fā)現(xiàn)減小翼型的彎度有利于提升其氣動性能[8]。賈雙林等人對 NACA63A010、NACA0018、NACA23012、NACA4412翼型進行氣動性能分析，結(jié)果發(fā)現(xiàn)在這四種翼型當中，對稱翼型NACA0018的氣動性能最好，最適用于垂直軸風力機[9]。LI Xingxing等人通過增加翼型厚度的方法來獲得一種新翼型，并證明了該新翼型的氣動特性優(yōu)于原翼型[10]。

目前對于影響風力機葉片翼型氣動性能研究主要針對某單一因素開展，本文擬采用數(shù)值計算方法對多個因素進行詳細研究，可以為進一步深入研究翼型氣動特性和翼型設(shè)計優(yōu)化提供系統(tǒng)性參考。

1 翼型氣動性能參數(shù)

圖1為氣流繞流翼型示意圖。翼型受到的合力可分解為升力和阻力，升力方向垂直于流速方向，阻力方向平行于流速方向，來流流過翼型表面時，上翼面速度大，下翼面速度小，根據(jù)伯努利定理，流速大則壓強小，流速小則壓強大，上下翼面存在的壓差即為翼型的升力。來流方向與弦線的夾角α記為攻角，不同攻角下翼型的氣動特性有較大的差異。翼型的升力系數(shù)Cl和阻力系數(shù)Cd是表征翼型氣動特性的重要無量綱參數(shù)，表達式為：

式中，v為來流風速，m/s；c為翼型弦長，m。

升阻比Cl/Cd是指升力系數(shù)與阻力系數(shù)的比值，升阻比是衡量翼型綜合性能的重要參數(shù)，升阻比越高說明翼型的氣動性能越好，風力機發(fā)電效率越高[11]。

圖1 翼型繞流示意圖Fig.1 Airfoil flow diagram

2 翼型幾何參數(shù)

翼型的氣動特性主要取決于其幾何形狀，通常使用翼型的幾何參數(shù)來表征其幾何形狀[12]。翼型的幾何參數(shù)主要包括弦線、厚度、最大厚度、彎度及最大彎度，圖2為翼型幾何參數(shù)示意圖。

圖2 翼型幾何參數(shù)示意圖Fig.2 Airfoil geometry parameter diagram

1）弦線c：翼型前緣至后緣的距離，其長度通常稱為弦長。

2）最大厚度t：翼型上翼面與下翼面的最大距離。

3）最大厚度位置Xt：最大厚度所處位置與前緣的距離。

4）最大彎度f：翼型中弧線與弦線的最大距離。

5）最大彎度位置Xf:最大彎度所處位置與前緣的距離。

NACA翼型是美國宇航局開發(fā)的系列翼型，通常使用記號NACAXYZZ來表示，其中X表示翼型的最大彎度占弦長的百分比，Y表示翼型的最大彎度位置與前緣的距離占弦長的百分比；ZZ表示翼型的最大厚度占弦長的百分比[13]。以NACA2415為例，2代表翼型的最大相對彎度為弦長的2%，4代表翼型最大彎度位置與前緣距離為弦長的4%，15代表翼型的最大厚度為弦長的15%。通常我們將翼型的厚度與弦長的比值（t/c）稱為相對厚度，彎度與弦長的比值（f/c）稱為相對彎度。

3 數(shù)值模型

3.1 流場數(shù)值模型

風力機在低速下流動屬于定常流動，空氣可視為不可壓縮氣體，由于不考慮熱量的傳遞故忽略能量方程，連續(xù)性方程和二維不可壓縮N-S方程為[14]：

連續(xù)性方程：

N-S方程：

S-A（Spalart-Allmaras）是一個適用于低雷諾數(shù)流動單方程模型，通常被用于求解受到壓力梯度作用的流動問題[15]，其輸運方程為[16]：

3.2 數(shù)值求解

利用ICEM建立幾何模型和進行結(jié)構(gòu)化網(wǎng)格劃分，網(wǎng)格數(shù)量為41 406。由于翼型近壁面處流場的參數(shù)變化梯度變化比較大，故在進行網(wǎng)格劃分時對翼型近壁面處網(wǎng)格進行加密[17]。如圖3所示，左側(cè)為翼型整體網(wǎng)格分布圖，右側(cè)為翼型近壁面區(qū)域網(wǎng)格分布圖。進口設(shè)置為速度進口，出口為壓力出口，翼型壁面為無滑移壁面邊界。建立基于S-A模型的數(shù)值計算，離散格式為二階迎風格式，壓力和速度耦合采用SIMPLE算法。當所有殘差的最大值均小于10-5時，可以認為計算收斂。

圖3 翼型整體網(wǎng)格分布Fig.3 Airfoil overall mesh distribution

4 結(jié)果分析與討論

4.1 數(shù)值方法驗證

采用本文數(shù)值模擬方法模擬風洞測試的NACA0012翼型[18]氣動性能，并將其結(jié)果與風洞試驗值進行對比，驗證數(shù)值方法的正確性。文獻報導的風洞試驗，研究雷諾數(shù) Re=106，攻角范圍為-18°-18°時NACA0012翼型的升力系數(shù)和阻力系數(shù)。選取其中部分攻角下的氣動性能參數(shù)進行驗證。圖4和圖5為翼型升力系數(shù)、阻力系數(shù)隨攻角變化曲線，結(jié)果表明數(shù)值方法得到的升力系數(shù)與實驗值趨勢相同且誤差很小，阻力系數(shù)有一定的差異，但誤差在允許范圍內(nèi)，說明數(shù)值模擬方法可以準確有效模擬翼型氣動性能。

圖4 NACA0012升力系數(shù)隨攻角變化Fig.4 NACA0012 lift coefficient changes with angle of attack

圖5 NACA0012阻力系數(shù)隨攻角變化Fig.5 NACA0012 drag coefficient changes with angle of attack

4.2 影響翼型氣動特性的因素

翼型的氣動特性隨著攻角的變化而改變，故攻角是影響其氣動特性的關(guān)鍵因素。除此之外，翼型的相對厚度、相對彎度及雷諾數(shù)都是影響翼型氣動特性的因素。

4.2.1 相對厚度（t/c）對翼型氣動特性的影響

選取NACA中的三種對稱翼型進行研究，分別為NACA0014、NACA0016、NACA0018，其中，00代表翼型的相對彎度為零，即中弧線與弦線重合；翼型的最大厚度分別為弦長的14%，16%，18%，圖6為三種翼型的幾何輪廓示意圖。

圖6 不同相對厚度翼型幾何輪廓對比Fig.6 Comparison of geometric profiles of different relative thickness airfoils

由圖7可得，當攻角小于12°時，三種翼型的升力系數(shù)隨著相對厚度的增加而減小，在12°攻角時三種翼型的升力系數(shù)達到最大值；隨著攻角的繼續(xù)增大，升力系數(shù)開始下降，在16°攻角時達到最小值0.65。由圖8可得，在0°-12°攻角范圍內(nèi)三種翼型的阻力系數(shù)差異不大，阻力系數(shù)曲線接近重合。當攻角大于12°時，阻力系數(shù)開始快速上升，且翼型的相對厚度越小阻力系數(shù)越大。當攻角超過某一臨界角后，翼型的升力系數(shù)下降、阻力系數(shù)上升的現(xiàn)象稱為失速，此時三種翼型失速臨界攻角近似為12°。

圖7 不同相對厚度翼型升力系數(shù)隨攻角變化曲線Fig.7 Different relative thickness airfoil lift coefficient changes with angle of attack

圖8 不同相對厚度翼型阻力系數(shù)隨攻角變化曲線Fig.8 Different relative thickness airfoil drag coefficient changes with angle of attack

由圖9可得，在0°～20°攻角范圍內(nèi)，升阻比呈現(xiàn)出先增大后減小的變化趨勢。三種翼型的最大升阻比均在8°攻角時出現(xiàn)，且相對厚度越小，最大升阻比越大。隨著攻角的繼續(xù)增加，翼型的升阻比開始下降，相對厚度對翼型升阻比的影響逐漸減小。

圖9 不同相對厚度翼型升阻比隨攻角變化曲線Fig.9 Different relative thickness airfoil lift-to-drag ratio changes with angle of attack

圖10和圖11分別為三種翼型在6°和16°攻角時的流線圖。由圖可知，在6°攻角時，翼型沒有出現(xiàn)旋渦，而在16°攻角時，翼型已經(jīng)出現(xiàn)較大的分離渦，氣流分離現(xiàn)象十分嚴重，翼型進入嚴重的失速狀態(tài)。在16°攻角時，三種翼型都出現(xiàn)較大的旋渦，其中NACA0014的分離點與前緣距離0.15c（c為弦長），NACA0016的分離點與前緣距離0.20c，NACA0018的分離點與前緣距離0.28c。由此可得，當翼型處于失速狀態(tài)時，相對厚度越小，渦的范圍越大，與前緣的距離越近，翼型的失速越嚴重。當攻角超過失速角后，翼型的邊界層會發(fā)生分離，其尾部形成旋渦，導致壓差阻力大大增加，翼型的氣動性能急劇下降進入失速狀態(tài)。這也解釋了圖7和圖8中當攻角超過12°后，為何翼型的升力系數(shù)下降、阻力系數(shù)顯著上升。

綜上所述，在攻角小于失速攻角時，相對厚度小的翼型可以獲得更好的氣動性能效果；當攻角大于失速角后，相對厚度大的翼型的氣動性能更佳。

圖10 6°攻角時翼型速度流線圖Fig.10 Airfoil speed streamline diagram at 6°angle of attack

圖11 16°攻角時翼型速度流線圖Fig.11 Airfoil speed streamline diagram at 16°angle of attack

4.2.2 相對彎度（f/c）對翼型氣動特性的影響

選取NACA0018、NACA2418、NACA4418三種翼型來分析相對彎度對翼型氣動特性的影響，三種翼型最大彎度占弦長的百分比分別為0%、2%、4%，圖12為三種翼型的幾何輪廓示意圖。

圖12 不同相對彎度翼型幾何輪廓對比Fig.12 Comparison of geometric profiles of different relative camber airfoils

由圖13可知，在0°～20°攻角范圍內(nèi)，翼型的升力系數(shù)隨著相對彎度的增加而增大。在12°攻角時，三種翼型的升力系數(shù)均達到最大值，隨后升力系數(shù)開始下降。由圖14可知，在0°～12°攻角范圍內(nèi)，三種翼型的阻力系數(shù)差異較小；當攻角大于12°時，三種翼型的阻力系數(shù)開始快速上升，但差異仍然不大。

圖13 不同相對彎度翼型升力系數(shù)隨攻角變化Fig.13 Different relative camber airfoil lift coefficient changes with angle of attack

圖14 不同相對彎度翼型阻力系數(shù)隨攻角變化Fig.14 Different relative camber airfoil drag coefficient changes with angle of attack

由圖15可知，翼型的升阻比呈現(xiàn)出先增大后減小的趨勢，且翼型的相對彎度越大，不同翼型的最大升阻比越大，其最大升阻比所對應(yīng)的攻角越小。當攻角超過失速角后，翼型升阻比受相對彎度的影響開始逐漸減小。

圖15 不同相對彎度翼型升阻比隨攻角變化Fig.15 Different relative camber airfoil lift-to-drag ratio changes with angle of attack

圖16 6°攻角時翼型速度流線圖Fig.16 Airfoil speed streamline diagram at 6°angle of attack

圖17 16°攻角時翼型速度流線圖Fig.17 Airfoil speed streamline diagram at 16°angle of attack

圖16和圖17分別為翼型在6°和16°攻角時的速度流線圖，在6°攻角時相對彎度越大，流動分離出現(xiàn)得越早，渦的形成也越早。在16°攻角時，三種翼型都出現(xiàn)嚴重的氣流分離，其中NACA0018的分離點與前緣距離0.26c（c為弦長），NACA2418的分離點與前緣距離0.36c，NACA4418的分離點與前緣距離0.48c。由此可得，當翼型未失速時翼型的彎度越大其尾部渦的出現(xiàn)越早；當翼型處于失速狀態(tài)時，彎度越小渦的范圍越大，與前緣的距離越近，翼型的失速越為嚴重。

綜上所述，在0°～20°攻角范圍內(nèi)，相對彎度越大翼型的氣動性能越佳。

4.2.3 雷諾數(shù)對翼型氣動特性的影響

雷諾數(shù)是一個可以用于判斷流體流態(tài)的無量綱參數(shù)，其實質(zhì)是指流體的慣性力與粘性力的比值，通常用Re表示。雷諾數(shù)對葉片翼型的氣動特性有較大的影響[19]，因此研究不同雷諾數(shù)下翼型的氣動參數(shù)十分有必要，圖18～圖20為翼型在不同雷諾數(shù)下的NACA0018翼型氣動特性參數(shù)隨攻角變化曲線。

由圖18可得，當攻角小于5°時，不同雷諾數(shù)下的升力系數(shù)沒有太大的差異，曲線幾乎接近重合。隨著攻角的繼續(xù)增大，升力系數(shù)曲線開始分離，且雷諾數(shù)越大升力系數(shù)越大。由圖19可知，翼型的阻力系數(shù)隨著雷諾數(shù)的增加而減小，并且翼型所對應(yīng)的失速攻角越大。從整體上看，翼型的升阻比隨著雷諾數(shù)增加而增大，且最大升阻比所對應(yīng)攻角越大。當攻角超過失速角后，雷諾數(shù)對升阻比的影響開始逐漸減小，升阻比曲線有重合的趨勢。

圖18 不同雷諾數(shù)時NACA0018翼型升力系數(shù)隨攻角變化Fig.18 NACA0018 airfoil lift coefficient changes with angle of attack with different Reynolds numbers

圖19 不同雷諾數(shù)時NACA0018翼型阻力系數(shù)隨攻角變化Fig.19 NACA0018 airfoil drag coefficient changes with angle of attack with different Reynolds numbers

圖20 不同雷諾數(shù)時NACA0018翼型升阻比隨攻角變化Fig.20 NACA0018 airfoil lift-to-drag ratio changes with angle of attack with different Reynolds numbers

5 結(jié)論

本文系統(tǒng)研究影響翼型氣動性能的因素，主要得到如下結(jié)論：

1）基于FLUENT軟件建立S-A模型對攻角范圍為-18°-18°時NACA0012進行數(shù)值計算，并與文獻風洞試驗得到的升力系數(shù)和阻力系數(shù)數(shù)據(jù)進行比較，結(jié)果表明計算值與試驗值吻合，驗證了數(shù)值方法的正確性。

2）當攻角小于失速角12°時，翼型的升阻比隨著相對厚度的增大而減小；隨著相對彎度的增大而增大。當翼型處于失速狀態(tài)時，翼型的相對厚度和相對彎度越小，其尾部出現(xiàn)的渦越大，與前緣的距離越近，流動分離的速度越快。

3）比較不同雷諾數(shù)NACA0018翼型的氣動特性參數(shù)，結(jié)果表明雷諾數(shù)越大，翼型的升力系數(shù)、升阻比越大，翼型的氣動性能越好。