曲寶麟 桂小紅 程魯康 高日鵬 施永博 李舜堯 宋增光

（中國礦業(yè)大學(xué)（北京），北京 100083）

0 引言

近年來，能源問題愈加被重視，對風(fēng)能的開發(fā)利用變得迫在眉睫。目前，研究風(fēng)力機(jī)氣動性能的方法有數(shù)值模擬計算、風(fēng)洞試驗和風(fēng)場測試等，其中，風(fēng)洞實驗與風(fēng)場測試是研究風(fēng)力機(jī)氣動性能最直接明顯的方法?？紤]到風(fēng)洞實驗與風(fēng)場測試受試驗條件和技術(shù)水平限制，成本高、周期長。因此，理論計算在風(fēng)力機(jī)初設(shè)階段的性能分析中發(fā)揮了重要作用[1-4]。

本文使用商用計算流體動力學(xué)(CFD)軟件FLUENT對Savonius風(fēng)力機(jī)(圖1)空氣流域進(jìn)行計算分析，得到壓力場、速度場等云圖，直觀分析其氣動性能，從而為Savonius風(fēng)機(jī)的設(shè)計提供一定的理論依據(jù)。

圖1 Savonius風(fēng)力機(jī)外形

1 計算模型的參數(shù)和建模過程

經(jīng)考慮實際情況，假設(shè)風(fēng)機(jī)額定輸出功率為80W，設(shè)計風(fēng)速為10m/s，轉(zhuǎn)速為356rad/min，葉片總高度為1000mm，葉片直徑為250mm，葉片厚度為10mm，轉(zhuǎn)軸直徑為60mm的Savonius風(fēng)力機(jī)二維模型。設(shè)來流速度為10m/s。模型如圖2所示。

圖2 葉片二維模型

通過CAD對葉片建模，導(dǎo)入ICEM對半圓形葉片進(jìn)行網(wǎng)格劃分，并定義邊界和密度。依據(jù)葉片大小，在葉片的表面生成一個風(fēng)洞，風(fēng)洞的形狀為長方形。

將工作段直接從葉片處延伸到遠(yuǎn)處，這樣具備較好的計算精度，因為風(fēng)動限制了葉片周圍氣體的運動，對模擬精度有一定程度的影響。設(shè)葉片中心與來流入口距離為1000mm，與來流出口距離為4000mm。

先確定葉片和風(fēng)洞各點的位置，將點連接成線，在風(fēng)洞表面生成面，將風(fēng)洞與葉片之間的面減去，拓?fù)錂z查幾何形狀是否丟失，在風(fēng)洞表面生成一個Block，再對Block進(jìn)行切分，切分出葉片的形狀并刪除相應(yīng)的Block，對風(fēng)洞、葉片邊界與Block進(jìn)行點和線的關(guān)聯(lián)，對風(fēng)洞進(jìn)行結(jié)構(gòu)網(wǎng)格單元劃分，考慮到計算機(jī)性能，將節(jié)點間距調(diào)整為較小，在保證網(wǎng)格數(shù)量的基礎(chǔ)上，提高計算精度，對網(wǎng)格進(jìn)行檢查調(diào)整，優(yōu)化網(wǎng)格質(zhì)量，使之滿足計算要求，整個風(fēng)洞網(wǎng)格劃分結(jié)果如圖3所示，葉片周圍的網(wǎng)格劃分如圖4所示。劃分網(wǎng)格后定義邊條件,將風(fēng)力機(jī)的葉片和轉(zhuǎn)軸定義為wall,風(fēng)洞表面定義為FLUID。最后保存為mesh文件備用。

圖4 葉片周圍的網(wǎng)格

2 流動控制方程

本研究中數(shù)值模擬對象為低速流動下葉片周圍的非定常流場，可以不計算能量方程[3]。

忽略熱傳遞的影響，且假設(shè)風(fēng)機(jī)葉片周圍氣體不可壓縮。數(shù)值模擬的基本方程為不可壓縮的非定常二維N-S方程[4]:

湍流模型采用κ-ε低雷諾數(shù)雙方程湍流模型[4]：

其中，κ為湍動能，ε為湍流耗散量。

3 數(shù)值模擬計算

選擇二維單精度處理。同時對重疊、不規(guī)則排列的網(wǎng)格進(jìn)行交換處理。定義湍流模型為具有非平衡壁面函數(shù)的標(biāo)準(zhǔn)雙方程模型。

對葉片進(jìn)行設(shè)置時將其選定為旋轉(zhuǎn)壁面形式，運動方式為相對原點運動，以37rad/s的角速度旋轉(zhuǎn)，將求解方式設(shè)置為二階迎風(fēng)格式，計算采用基于有限體積法的simple系列算法。將殘差收斂級數(shù)設(shè)為0.000001進(jìn)行迭代計算。迭代收斂后的結(jié)果見圖5和圖6。

圖5 殘差收斂圖示

圖6 扭矩收斂圖示

4 結(jié)果分析

經(jīng)過約3000次的迭代運算，殘差和扭矩的誤差都低于設(shè)定值。依據(jù)后處理軟件，可得流場及葉片周圍的速度云圖和壓力云圖。選擇整個流場區(qū)域display壓力云圖和速度云圖，結(jié)果如圖7和圖8所示。

圖7 壓力分布云圖

圖8 速度分布云圖

從圖7可以看出葉片處的壓力波動較大。左側(cè)葉片等壓面要大于右側(cè)葉片的等壓面，左側(cè)葉片迎風(fēng)側(cè)壓力較大，最大值為1.89e+03Pa；葉片背風(fēng)側(cè)較小，最小值為-1.80e+03Pa。故而形成壓力差，使其繞中心轉(zhuǎn)軸逆時針轉(zhuǎn)動。由圖8可得Savonius風(fēng)力機(jī)葉片所在的區(qū)域速度較小，而速度梯度較大，即葉片左右兩側(cè)向轉(zhuǎn)軸中心流速遞減。

圖9、圖10為速度矢量圖，通過分析可知：流體在抵達(dá)葉片之前，流速穩(wěn)定，在抵達(dá)葉片后，流速發(fā)生顯著變化。速度矢量分布圖相對穩(wěn)定，具有較小的扭矩系數(shù)，故該風(fēng)機(jī)葉片的風(fēng)能利用率較低。

圖9 速度矢量圖

圖10 葉片周圍速度矢量圖

5 結(jié)語

設(shè)計了一種適合低風(fēng)速啟動的小型垂直軸風(fēng)力發(fā)電機(jī)葉片，針對Savonius阻力型風(fēng)力擾流流動特點，運用ICEM進(jìn)行了二維建模及網(wǎng)格劃分，Savonius風(fēng)力機(jī)葉片的流動情況進(jìn)行數(shù)值模擬，得到其流場和壓力場，通過對其氣動性能進(jìn)行分析后得到以下結(jié)論:

1）風(fēng)流在抵達(dá)葉片之前速度急劇降低，在葉片的迎風(fēng)面處的風(fēng)速為零，表明Savonius風(fēng)力機(jī)具有較大的阻力系數(shù)。

2）左側(cè)葉片迎、背風(fēng)側(cè)的壓力差大于右側(cè)葉片，這樣就形成了左右葉片的壓力差，從而推動風(fēng)力機(jī)葉片繞中心轉(zhuǎn)軸逆時針轉(zhuǎn)動，揭示了Savonius風(fēng)力機(jī)的工作原理。

3）右側(cè)葉片迎風(fēng)側(cè)風(fēng)流一部分轉(zhuǎn)化為動力，一部分轉(zhuǎn)化為阻力，降低了風(fēng)能的利用率，產(chǎn)生了較大的能量損失，導(dǎo)致風(fēng)能利用率較低，表明該風(fēng)力機(jī)具有較大的改進(jìn)提升空間。