朱 斌，鄭 堤*，喬 峰，陳俊華

(1.寧波大學機械工程與力學學院，浙江寧波315211;2.浙江大學寧波理工學院機能學院，浙江寧波315100)

0 引言

潮流能的開發(fā)利用是當前一個熱點前沿問題，我國潮流能資源豐富，逐步提高潮流能捕獲效率、高效利用潮流能發(fā)電，對于優(yōu)化我國的能源消耗結(jié)構(gòu)、緩解沿海地區(qū)的能源壓力具有深遠意義。

槳葉是潮流能發(fā)電的關(guān)鍵部件，如將槳葉沿翼展方向分成n 個截面，每個截面都是一個葉素，而這些葉素是具有不同外形特征的二維翼型［1-3］，不同的翼型升阻力系數(shù)和導流性能有很大差別，對潮流能捕獲有很大影響。在水平軸潮流能槳葉的設(shè)計方面，目前普遍采用的是航空翼型和風力發(fā)電機專用翼型，缺少適合潮流流動特點的專用翼型。

國內(nèi)外學者做了大量有關(guān)提高水平軸潮流能捕獲裝置槳葉翼型設(shè)計的研究，并得到了較好的成果。Batten W M J 等［4］曾報道利用實驗驗證的數(shù)值方法設(shè)計水平軸潮流能水輪機，提高了潮流能捕能效率;王儉超等［5-7］基于CFD 軟件針對傳統(tǒng)翼型幾何形狀進行數(shù)值仿真，獲得了不同翼型的升阻力系數(shù)，為翼型選擇提供了參考;羅先武等［8］對翼型的流動特性進行分析，為小型潮流能葉片的設(shè)計提供了參考;余清清等［9］利用變槳距方式提高潮流能捕獲裝置在多變流速下的捕能效率。

本研究在上述研究的基礎(chǔ)上，通過計算仿真分析確定翼型的修改方案，在保證翼型的升力系數(shù)與阻力系數(shù)未發(fā)生較大改變的前提下設(shè)計得到綜合利用壓差升力和沖擊作用力的翼型，提高槳葉整體的能量捕獲效率和自啟動性能。

1 翼型分析

1.1 翼型周圍流場分析

現(xiàn)有的水平軸潮流能槳葉的設(shè)計基本都是基于風力發(fā)電Glanert 設(shè)計模型，將槳葉的流場視為一個流體繞翼型流動產(chǎn)生壓差升力的流場，忽略沖擊作用力。本研究通過分析發(fā)現(xiàn)大質(zhì)量流沖擊一個大面域產(chǎn)生沖擊作用力，在葉輪直徑較小時量級與壓差升力相當，而目前水平軸潮流能發(fā)電所用葉輪直徑受限于近海深度和漲潮落潮影響普遍在5 m 以內(nèi)，遠遠無法達到風力發(fā)電機的葉輪直徑，沖擊作用力不能簡單略去。

經(jīng)分析簡化后的翼型周圍的流場如圖1所示。

圖1 翼型流場

圖1 中，將翼型周圍流場分解為一個繞翼流動產(chǎn)生的壓差升力流場和一個來流直接沖擊翼型產(chǎn)生的沖擊作用力流場。翼型前端水流以合速度V合速度繞流翼型產(chǎn)生壓差升力;翼型下方形成一個射流區(qū)，由于槳葉轉(zhuǎn)動，水流的絕對運動速度方向如V絕對速度所示，相對于槳葉的速度方向如V相對速度所示。

1.2 沖擊捕能分析

1.2.1 沖擊捕能原理

沖擊作用力的產(chǎn)生:對于宏觀低速流體，牛頓碰撞理論有較高精度。流體質(zhì)點流與物面碰撞時，流體質(zhì)點將失去與物面垂直的法向動量，而保持原有的切向動量沿物面向下流去。由于法向動量的變化從而引起流體作用在物體上的力。

1.2.2 沖擊作用力產(chǎn)生條件

流體質(zhì)點沖擊物面，物面所受的作用力來源于流體速度(大小、方向)的改變，碰撞后質(zhì)點的速度方向逐漸與物面平行，而后質(zhì)點將不在對物面產(chǎn)生作用力。

現(xiàn)有的水平軸潮流能發(fā)電機械翼型采用的多是傳統(tǒng)的航空翼型和風力發(fā)電機專用翼型，這類翼型為得到較大的升力系數(shù)和較小的阻力系數(shù)，翼型下表面中后部位幾乎沒有彎度，且整個翼型的曲率一般都不超過6%，很難捕獲沖擊作用力。

2 優(yōu)化翼型

2.1 優(yōu)化目標

翼型結(jié)構(gòu)優(yōu)化的目標:在保證翼型的升力系數(shù)與阻力系數(shù)不發(fā)生較大的改變的前提下，實現(xiàn)翼型沖擊捕能的最大化?，F(xiàn)有翼型的普遍特點:彎度小，且主要彎曲部分集中在前部，尾部幾乎沒有彎度，不適合沖擊捕能。對于定槳距的潮流能捕獲裝置而言，啟動時，槳葉旋轉(zhuǎn)的切向速度很小，入流攻角很大，一般情況下當攻角大于13°時壓差升力迅速減小，阻力迅速增大，啟動所需轉(zhuǎn)矩大部分來源于沖擊捕能，由于現(xiàn)有翼型不適合沖擊捕能，直接導致自啟動性能不佳。

2.2 結(jié)構(gòu)優(yōu)化方法

本研究通過增加尾部翼型彎度的方法增加沖擊捕能。

2.2.1 確定可以優(yōu)化的部位

本研究借助fluent 對翼型的壓力場進行模擬，觀察翼型上、下表面壓力系數(shù)分布。為了能明顯看到翼型壓力分布的特點，筆者采用高速流體，具體設(shè)置參數(shù)如表1所示。

迭代200 步后收斂的壓力分布云圖和壓力系數(shù)圖如圖2所示。

表1 模擬參數(shù)

圖2 壓力分布

在圖2 中能明顯看出，翼型上表面形成一個負壓，下表面形成一個正壓，整體產(chǎn)生一個托舉力。且壓力主要集中于翼型前部，本研究對圖2(b)中壓力系數(shù)的面域進行計算，翼型后1/5 部分只提供了不到10%的升力，可以對這部分進行優(yōu)化。

2.2.2 翼型結(jié)構(gòu)優(yōu)化

以Naca6412 翼型為原外形，優(yōu)化后外形如圖3所示。圖3 是槳距角為槳距角40°時翼型優(yōu)化的一個例舉。

圖3 翼型結(jié)構(gòu)優(yōu)化

結(jié)構(gòu)設(shè)計約束:

(1)為了不額外增加翼型在水中的阻力，劃了平行于弦線的上、下邊界線，優(yōu)化后翼型必須在該范圍內(nèi)。

(2)尾部增加了彎度，出流速度方向與水平線夾角不得超過90°，超過90°后水流反向流動，既阻礙了前方水流，又大大增加了軸向載荷。

(3)彎度的變化要均勻，彎度突變會使沖擊變成撞擊，產(chǎn)生能量損失。

(4)彎度的增加，不能使流體的邊界層發(fā)生嚴重分離產(chǎn)生紊流，紊流會使阻力大大增加，且會使流體發(fā)生振動。

圖3 中，灰色點為原翼型外形，翼型尾部1/5 增加了彎度，使出流速度與水平線的夾角增加36°。

從宏觀上看，流體流經(jīng)葉輪所產(chǎn)生的轉(zhuǎn)矩可以看作兩個單獨的力所產(chǎn)生的力矩之和，且這兩個力都集中作用在進口和出口處。對不可壓縮的定常流體而言，總轉(zhuǎn)矩為［10］:

其中:V2cosα2=u－v2cosβ2。

式中:u—槳葉旋轉(zhuǎn)切向速度，v2—出口處流體相對翼型速度，β2—出口處翼型下底面與旋轉(zhuǎn)平面夾角，ρ—流體密度，Q—流量，r1—入口處流體等效半徑，r2—出口處流體等效半徑，V1—進口絕對速度，V2—出口絕對速度，α1—進口絕對速度與旋轉(zhuǎn)方向夾角，α2—出口絕對速度與旋轉(zhuǎn)方向夾角。

翼型尾部彎度增加后，主要減小了出口處翼型下底面與旋轉(zhuǎn)平面夾角，使cosβ2增大，控制轉(zhuǎn)速即u 不變，V2cosα2則隨之減小，由式(1)知轉(zhuǎn)矩T 增大。

2.2.3 優(yōu)化后翼型的升阻力系數(shù)模擬

本研究在網(wǎng)格劃分中對翼型附近圓形區(qū)域做了加密處理。

結(jié)構(gòu)優(yōu)化必須保證升力系數(shù)未發(fā)生較大的減小，阻力系數(shù)未有較大的增加。Naca6412 翼型和優(yōu)化翼型在相同雷諾數(shù)下模擬得到的升阻力系數(shù)如圖4所示。在10°攻角時達到最大升力系數(shù)，升力系數(shù)達到1.37，與原Naca6412 翼型的最大升力系數(shù)1.375 相近，阻力系數(shù)為0.192，略高于Naca6412 翼型的0.172。升阻力系數(shù)均未出現(xiàn)較大變化，符合要求。

圖4 升阻力系數(shù)

翼型彎度增加后有可能發(fā)生嚴重的邊界層分離，大大增加阻力，并發(fā)生振動，筆者繪制的翼型上、下壁面切應力的x 分量如圖5所示。

圖5 切應力x 方向分布

逆向流動可以通過壁面切應力的x 分量是否為負值來判斷，由圖5 可見切應力均大于零，未發(fā)生邊界層分離產(chǎn)生紊流，符合設(shè)計要求。

3 實驗驗證

3.1 槳葉設(shè)計

槳葉的基本設(shè)計參數(shù)如表2所示。

表2 槳葉的基本設(shè)計參數(shù)

本研究利用Matlab 迭代求解式(2 ～5)得槳葉弦長、扭角［11］，槳葉外形參數(shù)如圖6所示。

圖6 槳葉外形參數(shù)

槳葉UG 三維造型如圖7所示。槳葉可手動變槳距角。

圖7 UG 三維造型

式中:a—軸向誘導系數(shù)，b—切向誘導系數(shù)，CL—升力系數(shù)，r—半徑，φ—入流角，β—扭角，l—弦長。

3.2 實驗裝置

潮流能發(fā)電裝置安裝于大型水槽試驗平臺，安裝槳葉完畢后的照片如圖8所示。

圖8 實驗裝置

水槽長100 m，寬4 m，深3.5 m，平臺與鏈條連接，由伺服電機做速度控制。潮流能發(fā)電裝置由槳葉、一級鏈輪、行星增速器、200 W三相永磁同步發(fā)電機、負載控制、數(shù)據(jù)采集模塊組成。行星增速器選用速比為5∶1，電氣系統(tǒng)由整流器、逆變器、大功率電阻等組成，西門子PLC 為主控系統(tǒng)，筆者利用西門子A/D 轉(zhuǎn)換模塊將三相電量采集卡采集到的模擬信號轉(zhuǎn)換為數(shù)字信號，借助西門子OPC 和上位機LabVIEW 編寫的監(jiān)控控制程序通信，可以在電腦上直觀顯示以及記錄運行時的各個參數(shù)。

3.3 測試結(jié)果

在實驗過程中，本研究記錄捕能裝置的三相輸出電壓與電流，并轉(zhuǎn)化為功率數(shù)據(jù)記錄如圖9所示。

圖9 最大有用功率

從圖9 中看出，優(yōu)化后的槳葉(6 葉片)在平臺移動速度達到1.2 m/s 附近時功率上升趨勢變緩，槳葉(3 葉片)在1.1 m/s 附近功率上升趨勢變緩，即壓差升力達到最大，當大于臨界速度時，翼型開始進入失速狀態(tài)，而沖擊作用力隨平臺速度增加而增加，整體表現(xiàn)為功率緩慢增加。將優(yōu)化后槳葉與原始槳葉對比發(fā)現(xiàn)，在自啟動性能上優(yōu)化后槳葉優(yōu)于原始槳葉;在臨界速度前兩者功率相差不大，之后發(fā)電功率逐漸拉大，可見在撞擊捕能上優(yōu)化后槳葉占優(yōu)勢。

從6 葉片與3 葉片槳葉的對比也可以發(fā)現(xiàn)，在葉輪直徑較小的情況下，適當增加密實度大大增加捕能效果，在超過臨界速度以后更加顯著，這也說明沖擊作用力在小葉輪中的作用很大。

4 結(jié)束語

通過理論分析、模擬仿真以及槳葉的制作與實驗驗證，本研究可以得出以下結(jié)論:

(1)槳葉轉(zhuǎn)矩是流體的繞翼型流動產(chǎn)生的壓差升力和沖擊作用力共同作用的結(jié)果。葉輪半徑越小沖擊的影響越不能忽略。

(2)翼型的設(shè)計時綜合考慮壓差升力和沖擊捕能效果是有必要的。

(3)仿真結(jié)果表明優(yōu)化后的翼型在升阻力系數(shù)上未發(fā)生較大變化，理論分析表明在沖擊捕能上有一定提高。

(4)實驗結(jié)果表明，翼型達到設(shè)計要求，定槳距槳葉無論在自啟動性能上，還是在穩(wěn)定運行時的捕能效率上都有一定提高。但限于實驗測試設(shè)備，還未對動扭矩、推力等數(shù)據(jù)進行測量。

致 謝

實驗裝備的加工是浙大寧波理工學院工程訓練中心完成的，盛剛、李賢義工程師承擔了大量加工工作，對他們謹致謝意。

實驗工作是在浙大寧波理工學院機能學院海洋裝備研發(fā)中心完成的，謹致謝意。

［1］徐學濤.水平軸潮流能發(fā)電機械葉片設(shè)計及特性分析［D］.杭州:浙江大學機械學院，2012.

［2］鐵 庚，祁文軍.水平軸風力葉片機翼型的氣動特性數(shù)值模擬［J］.流體機械，2013，41(3):29-33.

［3］姜海波，曹樹良，李艷茹.水平軸風力機葉片扭角和弦長的理想分布［J］.太陽能學報，2013，34(1):1-6.

［4］BATTEN W M J，BAHAJ A S，MOLLAND A F，et al.The prediction of the hydrodynamic performance of marine current turbines［J］.Renewable Energy，2008，33(5):1085-1096.

［5］王儉超.水平軸潮流水輪機葉片設(shè)計和模型試驗研究［D］.青島:中國海洋大學海洋科學學院，2011.

［6］KIM B，KIM W，BAE S，et al.Aerodynamic design and performance analysis of multi-mw class wind turbine blade［J］.Journal of Mechanical Science and Technology，2011，25(8):1995-2002.

［7］HARTWANGER D，HORVAT A.3D Modelling of a Wind Turbine using CFD［C］// NAFEMS UK Conference，Cheltenham:［s.n.］，2008:1-4.

［8］羅先武，王 鑫.小型潮流水輪機葉片對稱翼型的流動特性分析［J］.水力發(fā)電學報，2010，29(2):239-243.

［9］余清清，孫 勇，周 浩.基于轉(zhuǎn)矩補償?shù)娘L電變槳加載系統(tǒng)的研究及仿真［J］.機電工程，2014，31(11):1457-1460.

［10］FINNEMORE E J，F(xiàn)RANZINI J B.Fluid mechanics With engineering applications［M］.America:McGraw-Hill，2005.

［11］馬 舜，李 偉，劉宏偉，等.水平軸潮流能發(fā)電系統(tǒng)能量捕獲機構(gòu)研究［J］.機械工程學報，2010，46(18):150-157.