龔玉祥，周曉亮，程明哲

(浙江運(yùn)達(dá)風(fēng)電股份有限公司，風(fēng)力發(fā)電系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 310012)

風(fēng)力機(jī)葉片的大型化趨勢(shì)使得葉片長(zhǎng)度越來(lái)越長(zhǎng)，為了保證葉片根部的強(qiáng)度，葉片的結(jié)構(gòu)設(shè)計(jì)上大多選用大厚度的翼型，這就導(dǎo)致了在大攻角下，大厚度的翼型更加容易發(fā)生氣動(dòng)分離現(xiàn)象，導(dǎo)致葉片吸收風(fēng)能的能力降低，從而降低了風(fēng)輪效率[1]；因此，控制流動(dòng)分離和減小葉片阻力成為風(fēng)力機(jī)氣動(dòng)研究領(lǐng)域的熱門問(wèn)題之一。

1947年，美國(guó)聯(lián)合飛機(jī)公司的Bmynes和Taylor首次提出了將渦流發(fā)生器(VGs)應(yīng)用于推遲飛機(jī)機(jī)翼的邊界層分離[2]。VGs 安裝在葉片上表面邊界層內(nèi)，并與當(dāng)?shù)貋?lái)流保持一定的側(cè)向夾角。VGs產(chǎn)生的尾渦擾動(dòng)分離區(qū)內(nèi)的氣流，使邊界層上部高能氣流與近壁低能氣流混合，增加近壁流體的動(dòng)量和能量，延緩分離，同時(shí)VGs 技術(shù)安裝工藝簡(jiǎn)單，經(jīng)濟(jì)性高；因此，被應(yīng)用于大型風(fēng)力機(jī)上[3-4]。Lin通過(guò)風(fēng)洞試驗(yàn)和數(shù)值模擬，得到了渦流發(fā)生器安裝位置應(yīng)該離分離點(diǎn)較近的結(jié)論[5]。Johansen等在數(shù)值計(jì)算中將渦流發(fā)生器兩側(cè)的邊界層條件設(shè)為對(duì)稱邊界及周期性邊界，研究了VGs 對(duì)風(fēng)力機(jī)葉片氣動(dòng)性能的影響規(guī)律[6]。張磊等采用CFD方法分析了VGs對(duì)風(fēng)力機(jī)專用翼型DU 91-W2-250氣動(dòng)性能的影響規(guī)律，對(duì)比了14°攻角下光滑翼型和VGs 翼型的流動(dòng)特性，揭示了VGs對(duì)翼型的流動(dòng)控制機(jī)制[7]。趙振宙等考慮了轉(zhuǎn)捩效應(yīng)的影響，采用SST全湍流模型和Gamma-Theta轉(zhuǎn)捩模型對(duì)翼型DU 91-W2-250進(jìn)行了數(shù)值模擬，研究結(jié)果表明轉(zhuǎn)捩模型計(jì)算結(jié)果較試驗(yàn)值更加吻合[8-9]。

在翼型上安裝渦流發(fā)生器的研究中使用的翼型弦長(zhǎng)大都為0.6 m，而實(shí)際風(fēng)力機(jī)葉片的厚翼型截面對(duì)應(yīng)的弦長(zhǎng)可達(dá)4 m[10]；因此，本文在不改變渦流發(fā)生器尺寸和安裝位置的情況下，研究了弦長(zhǎng)變化對(duì)翼型氣動(dòng)特性的影響規(guī)律。

1 計(jì)算模型

1.1 幾何模型和網(wǎng)格劃分

翼型選用DU系列風(fēng)力機(jī)專用翼型DU97-W2-300為研究對(duì)象，此翼型試驗(yàn)?zāi)Ｐ拖议L(zhǎng)為0.6 m，裝有VGs的試驗(yàn)數(shù)據(jù)由荷蘭Delft大學(xué)通過(guò)風(fēng)洞實(shí)驗(yàn)得到[11]，可驗(yàn)證數(shù)值計(jì)算方法的可靠性，模型如圖1所示。

圖1 幾何模型

計(jì)算段翼型的展向長(zhǎng)度與試驗(yàn)?zāi)Ｐ捅３忠恢?，長(zhǎng)度為0.175 m，計(jì)算弦長(zhǎng)分別為0.6 、1和1.5 m，在翼型弦長(zhǎng)20%處共安裝5組VGs，具體幾何尺寸如圖2所示。VGs的形狀為三角形，其弦向有效長(zhǎng)度為17 mm，高度5 mm，攻角為16.4°，每組內(nèi)的兩個(gè)三角形葉尖距離為10 mm，相鄰的兩組VGs的排列周期為35 mm。

采用ANSYS ICEM劃分結(jié)構(gòu)化網(wǎng)格，對(duì)VGs處進(jìn)行網(wǎng)格加密處理，第1層網(wǎng)格高度為0.01 mm，網(wǎng)格增長(zhǎng)率為1.1，保證y+≤1.5，網(wǎng)格如圖3所示。

圖2 渦流發(fā)生器尺寸參數(shù)

圖3 網(wǎng)格劃分

1.2 數(shù)值模擬方法

3種不同弦長(zhǎng)的翼型計(jì)算雷諾數(shù)相同，計(jì)算雷諾數(shù)為Re=2.0×106，攻角范圍為0°到20°，所采用的數(shù)值計(jì)算方法與文獻(xiàn)[8]相同，即使用轉(zhuǎn)捩模型進(jìn)行數(shù)值模擬計(jì)算，定義三維翼型兩側(cè)的斷面為周期性邊界條件。

2 計(jì)算結(jié)果分析

2.1 模型驗(yàn)證

圖4為轉(zhuǎn)捩模型計(jì)算的升力系數(shù)曲線與試驗(yàn)數(shù)據(jù)的對(duì)比圖。從試驗(yàn)數(shù)據(jù)可以看出，VGs使翼型失速現(xiàn)象明顯延遲，從無(wú)VGs翼型的臨界攻角12.37°延遲到有VGs 翼型的16.51°，將最大升力系數(shù)由1.55提升到了1.97。從計(jì)算結(jié)果可以看出，轉(zhuǎn)捩模型計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合良好。在攻角范圍為0°≤α≤18°時(shí)，轉(zhuǎn)捩模型計(jì)算誤差僅為5.28%，在α>18°時(shí)轉(zhuǎn)捩模型計(jì)算誤差僅為6.51%。由此可見(jiàn)，轉(zhuǎn)捩模型對(duì)于渦流發(fā)生器的數(shù)值計(jì)算精度高，其模擬結(jié)果能真實(shí)地反映VGs翼型的氣動(dòng)特性。

圖4 升力系數(shù)曲線

2.2 VGs對(duì)流場(chǎng)和表面壓力分布的影響

從升力系數(shù)的變化趨勢(shì)可以看出，α=18°位于VGs翼型的升力系數(shù)曲線的臨界攻角附近，此時(shí)翼型表面可能會(huì)出現(xiàn)氣流分離現(xiàn)象，截取如圖5所示的位置1和位置2對(duì)應(yīng)的截面，截面1為每組VGs的中心線，截面2為相鄰2組VGs的中心線。

圖5 位置1和位置2

圖6為截面1和截面2的流線圖。從兩個(gè)截面的流線圖可以發(fā)現(xiàn)兩個(gè)截面的流場(chǎng)是不同的，截面1尾緣處未出現(xiàn)明顯分離區(qū)域，截面2尾緣處出現(xiàn)明顯的分離區(qū)域，這表明翼型段上VGs的下游氣流在沿展向分布上是不均勻的。

圖7為截面1和截面2與翼型交線的壓力分布圖?？梢钥闯鰞蓚€(gè)截面的壓力分布只在安裝VGs附近是不同的，交線1上壓力分布出現(xiàn)負(fù)壓峰值，而交線2上出現(xiàn)負(fù)壓谷值。流線圖和壓力分布揭示了流動(dòng)產(chǎn)生渦流的根本原因，對(duì)于單個(gè)渦流發(fā)生器而言，當(dāng)氣流經(jīng)過(guò)時(shí)，其兩側(cè)變成了壓力面和吸力面，造成兩側(cè)壓力的不平衡，形成高能量的渦流，注入到近壁低能量氣流后，增加了近壁流體的動(dòng)量和能量，從而抑制邊界層的分離，推遲翼型失速。

圖6 攻角為18°下不同截面的流線圖

圖7 攻角為18°下不同截面壓力分布圖

2.3 翼型弦長(zhǎng)變化對(duì)流場(chǎng)的影響

圖8為安裝相同VGs，翼型不同弦長(zhǎng)的升力系數(shù)變化規(guī)律圖?？梢钥闯觯惭b了相同的渦流發(fā)生器，3種不同弦長(zhǎng)翼型的升力系數(shù)較光滑翼型均得到了提升，對(duì)翼型弦長(zhǎng)為0.6 m的升力系數(shù)的提升最為有效。在攻角范圍為0°≤α≤14°時(shí)，3種不同弦長(zhǎng)翼型的升阻力系數(shù)變化趨勢(shì)基本相同，在α>14°時(shí)，相同攻角下，隨著翼型弦長(zhǎng)增大，翼型升力系數(shù)減小。當(dāng)攻角為16°時(shí)，弦長(zhǎng)1 m的翼型升力系數(shù)較0.6 m時(shí)減小了4.35%，弦長(zhǎng)1.5 m的翼型升力系數(shù)較0.6 m時(shí)減小了10.82%。

圖8 弦長(zhǎng)對(duì)升力系數(shù)的影響

基于以上分析，截取得到3種弦長(zhǎng)下攻角分別為14°、16°和18°時(shí)截面2的流線圖，如圖9所示。可以看出：當(dāng)攻角為14°時(shí)，3種弦長(zhǎng)的翼型尾緣處均未有明顯的分離區(qū)域；當(dāng)攻角為16°和18°時(shí)，隨著翼型弦長(zhǎng)的增大，翼型尾緣的分離區(qū)域均逐漸變大。這是因?yàn)樵诶字Z數(shù)相同的情況下，翼型弦長(zhǎng)的增大意味著來(lái)流速度的減小，當(dāng)流體沿著翼型表面流動(dòng)并掠過(guò)VGs時(shí)，形成的渦流的動(dòng)量和能量減小，與近壁低能流體混合后的流體能量將會(huì)減小， 同時(shí)混合后的流體所經(jīng)過(guò)的距離增大，這就導(dǎo)致了翼型尾緣分離區(qū)域的增大?；诖?，預(yù)測(cè)當(dāng)翼型的弦長(zhǎng)增大到4 m時(shí)，安裝本文中的渦流發(fā)生器于弦長(zhǎng)20%處，其性能將進(jìn)一步降低甚至可能失去推遲翼型失速的作用；因此，在設(shè)計(jì)適合于風(fēng)場(chǎng)中運(yùn)行的風(fēng)力機(jī)葉片的渦流發(fā)生器時(shí)，應(yīng)當(dāng)考慮葉片截面的實(shí)際弦長(zhǎng)。

(a)α=14°

(b) α=16°

(c) α=18°

3 結(jié)論

1) 轉(zhuǎn)捩模型計(jì)算的升力系數(shù)與試驗(yàn)值吻合良好，計(jì)算精度高，可以準(zhǔn)確地反映VGs翼型的氣動(dòng)特性。VGs翼型下游流動(dòng)的不均勻性形成的渦流使得近壁流體的能量增強(qiáng)，從而抑制邊界層的分離，推遲失速。

2)在攻角范圍為0°≤α≤14°時(shí)，3種不同弦長(zhǎng)翼型的升阻力系數(shù)變化趨勢(shì)基本一致；在攻角大于14°時(shí)，相同攻角和雷諾數(shù)時(shí)，翼型弦長(zhǎng)增大，升力系數(shù)減小，翼型的尾緣分離區(qū)域增大。本文的研究結(jié)果為設(shè)計(jì)適合于實(shí)際風(fēng)場(chǎng)中的風(fēng)力機(jī)葉片的渦流發(fā)生器提供了指導(dǎo)。