NPU－WA系列風(fēng)力機翼型設(shè)計與風(fēng)洞實驗

喬志德， 宋文萍， 高永衛(wèi)

（西北工業(yè)大學(xué)，翼型、葉柵空氣動力學(xué)國家重點實驗室，陜西 西安710072）

0 引 言

風(fēng)力機葉片設(shè)計是風(fēng)力發(fā)電機組設(shè)計的一項核心技術(shù)，構(gòu)成葉片的翼型是葉片設(shè)計的基礎(chǔ)，20世紀(jì)90年代以前，風(fēng)力機葉片設(shè)計通常使用已有的傳統(tǒng)翼型，如4位數(shù)字NACA44系列和NACA63或64系列翼型［1－2］，自從1980年代后期以來，西歐和美國進行了專門用于風(fēng)力機的先進翼型設(shè)計研究，瑞典宇航研究院20世紀(jì)90年代設(shè)計了FFA－W3－211、FFA－W3－241和FFA－W3－301三個較厚的翼型，荷蘭Delft大學(xué)在1992年和1993年設(shè)計了相對厚度分別為25%、21%的 DU91－W2－250和 DU93－W2－210翼型，此后設(shè)計了相對厚度分別為18%的DU95－W－180、DU96－W－180 翼 型 和 30% 的 DU97－W－300 翼型，形成了相對厚度18%到40%的DU翼型系列［3］，丹麥RIS?國家實驗室在90年代后期發(fā)展了由RIS?－A1－18，RIS?－A1－21和 RIS?－A1－24三個翼型組成的RIS?風(fēng)力機翼型系列，在2000年之后針對大型風(fēng)力機設(shè)計了具有更高設(shè)計升力的RIS?－B1族翼型［4］。但上述翼型缺乏較高雷諾數(shù)下的實驗驗證，目前還主要用于中、小風(fēng)力機葉片設(shè)計。

1984年美國可再生源國家實驗室（NREL）開展了風(fēng)力機翼型族的設(shè)計研究，到90年代，為各類風(fēng)力機發(fā)展了不同性能的，從根部到葉尖的，能適應(yīng)結(jié)構(gòu)要求的9個翼型族［2，5］，這些新翼型的采用，使風(fēng)電機組年發(fā)電量增加了10%～35%。由于是針對當(dāng)時的中、小型風(fēng)力機設(shè)計的，其設(shè)計升力和最大升力系數(shù)都不很高。

美國Sandia國家實驗室在2004年發(fā)表了大型風(fēng)力機創(chuàng)新葉片的設(shè)計技術(shù)報告［6］，提出了把新翼型作為4項關(guān)鍵技術(shù)之一的大型風(fēng)力機葉片創(chuàng)新設(shè)計新概念，所有4項關(guān)鍵技術(shù)都指向一個主要目標(biāo)：在獲得高空氣動力學(xué)性能的條件下，減少葉片重量。所提出的一個2.4MW，直徑104m的風(fēng)力機設(shè)計方案中，葉片大部分剖面雷諾數(shù)都超過5.0×106。

風(fēng)力機葉片的重量和費用正比于半徑的2.4次方，而發(fā)電量正比于風(fēng)力機半徑的平方，所以隨風(fēng)力機功率增加，風(fēng)力機尺寸將會有更快的增加，更大的尺寸意味著更高的運行雷諾數(shù)、更大的重量、更大的陣風(fēng)風(fēng)載及伴隨的振動和疲勞限制。因此大型風(fēng)力機葉片的主要技術(shù)要求是：減少葉片重量，以減少包括制造費用和運輸成本在內(nèi)的發(fā)電成本，減少慣性載荷、陣風(fēng)載荷以及相應(yīng)的系統(tǒng)載荷；并提高葉片的風(fēng)能捕獲能力。由于大型風(fēng)力機運行工況下葉片主要剖面具有很高的雷諾數(shù)，因此要求翼型在高雷諾數(shù)時具有高的氣動性能，此外，大型風(fēng)力機還要求翼型具有更高的設(shè)計升力，這是因為高設(shè)計升力可以減少實度（減少葉片弦長）以減少葉片面積，從而可以減少葉片重量、節(jié)約制造和運輸成本，并減輕陣風(fēng)載荷和慣性載荷；還有，高設(shè)計升力有利于在低于平均風(fēng)速的使用周期內(nèi)提高風(fēng)能捕獲能力，以增加風(fēng)力機的年發(fā)電量。由上述可見，大型風(fēng)力機翼型應(yīng)該在在更高的雷諾數(shù)和更高的設(shè)計升力下具有更好的空氣動力特性。目前國內(nèi)外風(fēng)力機翼型的設(shè)計升力系數(shù)大都小于1.0，除個別翼型外，風(fēng)力機翼型族的實驗雷諾數(shù)都沒有超過3.0×106，所以研究發(fā)展經(jīng)過實驗驗證的高雷諾數(shù)、高設(shè)計升力風(fēng)力機翼型族，可以導(dǎo)致設(shè)計出具有創(chuàng)新性的高性能大型風(fēng)力機葉片。

我國過去缺乏以厚翼型為特點的風(fēng)力機翼型設(shè)計與計算方法研究；缺乏可減輕風(fēng)力機葉片陣風(fēng)過載的翼型設(shè)計與計算方法研究，也沒有開展風(fēng)力機翼型族的設(shè)計研究，沒有具有知識產(chǎn)權(quán)的風(fēng)力機翼型族，特別是沒有可供兆瓦級和多兆瓦級大型風(fēng)力機葉片設(shè)計使用的翼型族，嚴(yán)重影響了我國風(fēng)電工業(yè)的發(fā)展。

NPU－WA風(fēng)力機翼型族是針對變距或變轉(zhuǎn)速大型風(fēng)力機葉片，使用計算流體力學(xué)方法設(shè)計的，并在NF－3低速翼型風(fēng)洞中進行了從1.0×106到5.0×106的5個不同雷諾數(shù)的風(fēng)洞實驗，與已有國外風(fēng)力機翼型較低雷諾數(shù)實驗數(shù)據(jù)相比，給出具有更高雷諾數(shù)的、更完整的氣動性能實驗數(shù)據(jù)。

1 NPU－WA翼型設(shè)計方法

在風(fēng)力機翼型設(shè)計中綜合使用了課題組多年來研究發(fā)展的翼型設(shè)計與計算方法，這些方法的詳細(xì)描述見所給出的相應(yīng)參考文獻，這里簡介如下：

（1）反設(shè)計方法

按給定目標(biāo)壓力分布的翼型反設(shè)計方法［13］，用于給定較小迎角（或較低設(shè)計升力系數(shù)）時的目標(biāo)壓力分布設(shè)計翼型；基于亞聲速速勢方程的混合邊界條件翼型設(shè)計方法［7］，用于根據(jù)部分表面給定目標(biāo)壓力分布，部分表面給定幾何外形的設(shè)計要求設(shè)計翼型；N－S方程翼型設(shè)計方法［13］，用于大迎角（或高升力）條件下，給定目標(biāo)壓力分布設(shè)計翼型，適用于各種迎角和雷諾數(shù)。

（2）翼型優(yōu)化設(shè)計方法

通過采用多目標(biāo)優(yōu)化方法基于雷諾平均N－S方程［9－10］、低速線化速勢方程或跨聲速全速勢方程－附面層迭代方法進行優(yōu)化設(shè)計［8］。優(yōu)化目標(biāo)是提高翼型在高升力和高雷諾數(shù)下的升阻比，并將中等升力系數(shù)、較低雷諾數(shù)等非設(shè)計條件以及翼型力矩系數(shù)及翼型相對厚度作為需要滿足的約束條件。

（3）人－機對話修改設(shè)計［13］

通過計算機屏幕直接修改翼型幾何外形，通過對修形后翼型的氣動特性計算檢驗修改效果，由于優(yōu)化方法難以使翼型外形有大的改變，所以允許較大修改量的人－機對話修形是優(yōu)化方法的必要補充，同時也是進行適當(dāng)修形以滿足非設(shè)計條件的補充手段。

（4）校核計算方法［13－15］

對使用上述方法設(shè)計的翼型進行設(shè)計條件和非設(shè)計條件下的校核計算，在小迎角（或低升力）時，主要使用XFOIL計算軟件，在大迎角（或高升力）使用課題組研究發(fā)展的雷諾平均N－S方程計算方法。

圖1給出所設(shè)計的NPU－WA－210翼型在設(shè)計工況下壓力分布與國外同類風(fēng)力機翼型的比較。

圖1 設(shè)計升力和設(shè)計雷諾數(shù)下NPU－WA－210翼型與其他風(fēng)力機翼型壓力分布計算比較Fig.1 Comparison of the design pressure distribution of NPU－WA－210with the other airfoils

2 NPU－WA翼型族翼型的名稱和幾何外形

圖2給出了NPU－WA風(fēng)力機翼型族的幾何外形，根據(jù)設(shè)計技術(shù)要求，該風(fēng)力機翼型族的相對厚度分別為0.15、0.18、0.21、0.25、0.30、0.35和0.40弦長，共7個翼型，分別被命名為：NPU－WA－150、NPUWA－180、NPU－WA－210、NPU－WA－250、NPU－WA－300、NPU－WA－350和 NPU－WA－400。NPU－WA 翼型編號最后三位數(shù)字中的前兩位表示相對厚度，最后一位“零”表示該翼型為初次設(shè)計，若為第一次修改設(shè)計則為1，以此類推。編號中的“NPU”表示該翼型族是由西北工業(yè)大學(xué)研究發(fā)展的，“WA”表示該翼型族是為風(fēng)力機設(shè)計的專用翼型。弦長C為100時，不同相對厚度翼型的后緣厚度見表1。

表1 不同相對厚度翼型的后緣厚度Table 1 Trailing edge thickness of airfoils with different relative thickness

圖2 NPU－WA風(fēng)力機翼型族Fig.2 NPU－WA wind turbine airfoil family

3 NPU－WA風(fēng)力機翼型族的風(fēng)洞實驗

NF－3風(fēng)洞為直流閉口式全鋼結(jié)構(gòu)風(fēng)洞，有三個可更換實驗段，洞體長80m，二元實驗段寬1.6m、高3m、長8m。風(fēng)速為15m／s～120m／s，湍流度0.045%。試驗用的翼型模型均為鋼質(zhì)骨架木質(zhì)結(jié)構(gòu)模型，模型展長1.595m，弦長0.8m，模型在風(fēng)洞中的安裝情況如圖3所示。在翼型模型上下表面開靜壓孔，測量表面的壓力，用于計算翼型的升力和俯仰力矩；在距模型后緣1.2倍弦長處安裝總壓排管，有186根文德利型總壓管和9根靜壓管，測量模型尾跡區(qū)的總壓分布和靜壓，用以計算翼型的阻力。測量寬度范圍為2000mm，可根據(jù)試驗的具體情況進行移動。數(shù)據(jù)采集系統(tǒng)采用美國PSI公司的9816電子掃描閥，共有512個壓力測量通道，采集速度為100Hz／ch，采集精度為±0.05%。該系統(tǒng)用來采集翼型的表面壓力和尾耙的壓力。除了自由轉(zhuǎn)捩實驗外，還進行了固定轉(zhuǎn)捩實驗，固定轉(zhuǎn)捩采用zigzag粗糙帶。

圖3 翼型模型在NF－3風(fēng)洞中的安裝Fig.3 The test model in NF－3wind tunnel

實驗雷諾數(shù)分別為1.0×106、1.5×106、3.0×106、4.0×106和5.0×106。對應(yīng)實驗段風(fēng)速約為20m／s～100m／s。對實驗數(shù)據(jù)進行了包括阻塞修正、流向曲率修正和側(cè)壁干擾修正等風(fēng)洞干擾修正。

4 NPU－WA翼型的氣動特性

4.1 NPU－WA翼型的高升力特性

由圖4可見，主翼型與外側(cè)翼型具有高于1.2的設(shè)計升力，最大升阻比對應(yīng)的升力很接近于最大升力是該翼型族的主要特點之一。

4.2 NPU－WA翼型的高雷諾數(shù)特性

圖5給出雷諾數(shù)對NPU－WA－210翼型最大升阻比和最大升力系數(shù)的影響，由圖可見，翼型的氣動性能隨雷諾數(shù)的變化較為和緩，在直到Re＝5.0×106時仍能保持高的升阻比。

4.3 最大升力系數(shù)對粗糙度的敏感性

最大升力系數(shù)CLmax對粗糙度敏感性可被定義為（CLmaxfr－CLmaxfix）／CLmaxfr，下標(biāo)中的fr表示自由轉(zhuǎn)捩，fix表示固定轉(zhuǎn)捩，一般來說，高設(shè)計升力和大相對厚度都會增加翼型對粗糙度的敏感性，圖6指出，用于葉片外側(cè)剖面的 NPU－WA－150和 NPU－WA－180翼型在所有試驗雷諾數(shù)范圍內(nèi)，最大升力系數(shù)對粗糙度的敏感性都低于10%，分別為0.049～0.076和0.052～0.095，表明其對粗糙度不敏感。主翼型NPU－WA－210在高雷諾數(shù)下，如Re＝5.0×106時敏感性參數(shù)為0.097，Re＝3.0×106時為0.123。

內(nèi)側(cè)翼型最大升力系數(shù)對粗糙度敏感性通常是不很重要的，可以允許較高的敏感性。

4.4 厚翼型的氣動特性

圖7給出了用于葉片內(nèi)側(cè)剖面的NPU－WA－300翼型空氣動力學(xué)特性隨雷諾數(shù)的變化。

5 結(jié) 論

風(fēng)洞實驗表明，NPU－WA翼型族達到了在高雷諾數(shù)、高升力條件下實現(xiàn)高升阻比和外側(cè)翼型對粗糙度不敏感的主要設(shè)計要求。為我國自主研發(fā)大型風(fēng)力機提供了可供實際使用的翼型幾何數(shù)據(jù)和雷諾數(shù)從1.0×106～5.0×106的風(fēng)洞實驗數(shù)據(jù)。

［1］ABBOTT F T，Von DOENHOFF A E.Theory of wing sections［M］.Dover Publications，Inc.，New York，1959.

［2］BERTAGNOLIO F，SORENSEN N，JOHANSEN J，et al.Wind turbine airfoil catalogue［R］.Ris?National Laboratory，Roskide，Denmark，August 2001.

［3］TIMMER W A，Van ROOIJ R P J O M.Summary of the delft university wind turbine dedicated airfoils［R］.AIAA－2003－0352，2003.

［4］FUGLSANG P，BAK C，GAUNAA M，et al.Design and verification of the RIS?－B1Airfoil Family for Wind Turbines［R］.AIAA－2004－668，2004.

［5］Dan SOMERS M.The S830，831，and S832Airfoils［R］.NREL／SR－500－36339，2005.

［6］WINDPACT Blade System Design Studies.Innovative design approaches for large wind turbine blades，final report［R］.Sandia National Laboratory，SAND－2004－0074，2004.

［7］QIAO Z D.Subsonic airfoil design with mixed boundary conditions［A］.Numerical and Applied Mathematics［M］.Scientific Publishing Co.，IMACS，1989：209－212.

［8］ZHONG B W，QIAO Z D.Multiobjective optimization design of transonic airfoils［R］.ICAS－94－2.1.1，1994.

［9］錢瑞戰(zhàn)，喬志德，宋文萍.基于N－S方程的跨聲速翼型多目標(biāo)多約束優(yōu)化設(shè)計［J］.空氣動力學(xué)學(xué)報，2000，18（3）：350－365.（QIAN R Z，QIAO Z D，SONG W P.Multiobjective optimization design of transonic airfoils with constraints based on navier－stokes equations［J］.ACTAAerodynamicaSinica，2000，18（3）：350－365.）

［10］QIAO Z D，QIN X L，YANG X D.Wing design by solving adjoint equations［R］.AIAA－2002－0263，2002.

［11］QIAO Z D，YANG X D，QIN X L，et al.Numerical optimization design of wings by solving adjoint equations［R］.ICAS2002－1104.1，2002.

［12］楊旭東，喬志德.基于共軛方程方法的跨聲速機翼氣動優(yōu)化設(shè)計［J］.航空學(xué)報，2003，24（1）：1－5.（YANG X D，QIAO Z D.Optimum aerodynamic design of transonic wing based on adjoint equation－ns method［J］.ACTAAeronauticaetAstronauticaSinica，2003，24（1）：1－5.）

［13］喬志德.先進翼型的CFD設(shè)計及應(yīng)用［A］.現(xiàn)代力學(xué)與科技進步［M］.莊逢甘主編.北京：清華大學(xué)出版社，1997：984－987.（QIAO Z D.CFD Design and Application of Advanced Airfoil［A］.Modern Mechanics and Scientific and Technological Progress［M］.ZHUANG F G，Beijing：Qinghua University Press，1997：984－987.）

［14］喬志德.Computational aerodynamics of wing design［A］.計算流體力學(xué)新進展［M］.DUBOIS F，鄔華謨編.高等教育出版社，2001：273－294.（QIAO Z D.Computational aerodynamics of wing design［A］.Recent Progress in Computational Flui－d Dynamics［C］.DUBOIS F，WU H M，Higher Edu－cation Press，2004：273－294.）

［15］DRELA M.XFOIL：an analysis and design system for low reynolds number airfiols［A］.Conference on Low Reynolds Number Airfoil Aerodynamic［C］，University of Notre Dame，June 1989.

［16］DRELA M.XFOIL 6.94［Z］.MIT Aero ＆ Astro，USA，2001.