海流能發(fā)電裝置葉片性能及氣蝕研究

1、海流能發(fā)電裝置葉片性能及氣蝕研究research on the blade performance and cavitation of generate electricity equipment by energy of ocean current摘要：對水平軸螺旋槳式海流能發(fā)電裝置葉片的升阻系數(shù)特性及氣蝕特性進(jìn)行了數(shù)值模擬研究，研究結(jié)果與文獻(xiàn)的試驗結(jié)果較好地吻合，并得出了結(jié)論。關(guān)鍵詞：海流能發(fā)電裝置； 葉片性能；氣蝕；仿真研究abstract: numerical simulation and research on the rising-resistance coefficient ch

2、aracteristic and cavitation characteristic are carried out for blade of generate electricity equipment by energy of ocean current with level-screw type. the research result is accordance with the test result . the conclusion is obtained.key words: generate electricity equipment by energy of ocean cu

3、rrent； blade performance；cavitation ; simulation research1 概述在開發(fā)綠色新能源的過程中，海洋能的開發(fā)利用被認(rèn)為有巨大開發(fā)潛力，作為海洋能的一種海流能，人類可以像風(fēng)力機捕獲風(fēng)能一樣來捕獲海流能1-2。如圖1所示為用于捕獲海(潮)流能的水下風(fēng)車裝置，利用海流以一定的速度和攻角沖擊流線型的槳葉，通過變速機構(gòu)，驅(qū)動發(fā)電機發(fā)電。*基金項目：國家自然科學(xué)基金資助（50505043）；浙江省重大科技攻關(guān)項目（2006c11012）.收稿日期：2007-12-29 杭州市 310027要提高該類型裝置的能量捕獲效率，關(guān)鍵要素之一就是要設(shè)計出高效的葉片

4、結(jié)構(gòu)。葉片在強度滿足水下工作環(huán)境要求的前提下，既要有較高的升力系數(shù)特性，又要有較低的阻力系數(shù)和延遲的失速特性；同時還需防止葉片氣蝕現(xiàn)象的產(chǎn)生，這是一個不同于風(fēng)力機的問題。葉片在水中旋轉(zhuǎn)，易出現(xiàn)氣蝕現(xiàn)象，破裂的氣泡對槳葉的沖擊損壞表面的防腐層，加速槳葉的腐蝕，也會使系統(tǒng)產(chǎn)生振動和噪聲。 本文對葉片的升阻系數(shù)特性及氣蝕特性進(jìn)行數(shù)值模擬研究，旨在對解決上述問題作有益的探討。 圖1“水下風(fēng)車”示意圖2 葉片結(jié)構(gòu)對升阻系數(shù)的影響水平軸螺旋槳式海流能發(fā)電裝置的葉片所處的水動力學(xué)環(huán)境與風(fēng)力機或機翼的空氣動力學(xué)環(huán)境有相似之處，也有其自身的特點，所以一方面要借鑒風(fēng)力機翼型理論，另一方面也要結(jié)合具體的水工作環(huán)境對

5、其進(jìn)行研究。國內(nèi)外關(guān)于翼型繞流以及流體和機翼相互作用問題的研究包括理論研究、數(shù)值模擬和試驗3個方面?；谒嘛L(fēng)車葉片與風(fēng)力機葉片相似的特點均采用大的高寬比結(jié)構(gòu)，它的基本性能可以很好地通過葉素動量（blade element momentum）理論來模擬分析3。葉片表面的三維不可壓縮定常粘性流在柱面坐標(biāo)下的流動控制方程為(1) (2) （3）式中為流體密度；為速度矢量沿軸的3個速度分量；為柱坐標(biāo)3個坐標(biāo)軸；為流體壓力；為動力粘度。對于一個旋轉(zhuǎn)系統(tǒng)，通過對式（1）式（3）積分，可得沿水流方向和垂直于水流方向的動量積分方程5： （4）（5）圖3 升力系數(shù)曲線式中為邊界層動量厚度；為邊界層形狀系數(shù)；為

6、邊界層邊緣速度；為表面摩擦因數(shù)；為葉輪旋轉(zhuǎn)速度；為邊界層位移厚度；為邊界層厚度；為流體沿翼型表面的脫落點。當(dāng)流體流過葉片時會產(chǎn)生分離的現(xiàn)象，從而對葉片產(chǎn)生持續(xù)的作用力；同時在葉片的后部會有漩渦脫落的現(xiàn)象產(chǎn)生，特別當(dāng)葉片后部兩側(cè)交替產(chǎn)生的漩渦周期性脫落會對葉片產(chǎn)生與來流方向垂直的周期性作用力時，持續(xù)和脈動流體力的作用會使葉片產(chǎn)生疲勞、振動和噪聲，嚴(yán)重時會使葉片破壞；反之葉片也會影響流場的分布，實際上對于有限翼展的葉片，其上下表面壓差產(chǎn)生了使葉片旋轉(zhuǎn)的升力，同時由于水流企圖從高壓側(cè)流向低壓側(cè)，所以在葉片的后緣處(見圖2)及兩端形成渦流，渦流造成阻力的增加，即引起一個誘導(dǎo)阻力。圖2翼型渦流升力系數(shù)和

7、阻力系數(shù)是表征葉片性能的重要參數(shù)，文獻(xiàn)6從試驗測量的角度研究了葉片升力系數(shù)(cl)、阻力系數(shù)(cd)及葉片表面的壓力分布特性；文獻(xiàn)4則試驗研究了“水下風(fēng)車”（marine current turbine）在不同流速、不同壓力環(huán)境下的功率系數(shù)、扭矩系數(shù)(ct)與葉輪尖速比(tsr)的關(guān)系。隨著計算流體力學(xué)技術(shù)的成熟，通過數(shù)值模擬方法獲取葉片的各種性能參數(shù)越來越得到應(yīng)用。翼型的雷諾數(shù)計算如下：re式中v為流動的特征速度，這里取4m/s(為了與文獻(xiàn)6的試驗結(jié)果對比)；l為流動的特征長度，這里取葉片弦長0.23m (取該值同樣是為了與文獻(xiàn)6的試驗結(jié)果對比)；為流體的運動粘度，對于水=1.007

8、5;10-6 m2/s。圖3和圖4為二維翼型naca63-215在fluent軟件中求得的升力系數(shù)、阻力系數(shù)圖，與文獻(xiàn)6中試驗結(jié)果對比，二者非常好地實現(xiàn)吻合；但在仿真過程中湍流模型的選擇是個非常重要的問題，在仿真中，小攻角時采用標(biāo)準(zhǔn)的k-e模型，并采用非平衡壁面函數(shù)的近壁面圖4 阻力系數(shù)曲線圖5 不同翼型的升力系數(shù)曲線處理，該模型對于處理貼近壁面附近的流動計算是很有效的；當(dāng)大攻角時，在考慮流體分離及回流現(xiàn)象的基礎(chǔ)上，采用k-模型。仿真結(jié)果表明用數(shù)值模擬的方法也可以獲得葉片翼型合理的升阻特性。下面以naca2415、naca 4415、naca4424、 naca2411為例研究翼型的彎度、厚度

9、對升阻特性的影響，同時通過所得的葉片表面壓力系數(shù)分布研究葉片的氣蝕現(xiàn)象，通過分析力求可以獲得既提高葉片性能又降低葉片氣蝕發(fā)生的合理結(jié)構(gòu)。對圖5和圖6中不同翼型的升力系數(shù)和阻力系數(shù)曲線分析，對相同厚度(15%倍弦長)的naca 4415和naca 2415翼型，naca 2415的升力系數(shù)要小于naca 4415，而其阻力系數(shù)基本相當(dāng)，即同等厚度下，適當(dāng)彎度的增加(naca 2415的彎度要小于naca 4415)可以增大升力系數(shù)而維持阻力系數(shù)基本不變；再比較同等彎度但不同厚度的naca 4415與naca 4424翼型及naca 2415和naca 2411翼型，結(jié)果表明大的厚度在減小升力系數(shù)

10、的同時也增大了阻力系數(shù)，這對葉片的性能是不利的，但隨著彎度的降低，厚度對阻力系數(shù)的影響逐漸減小。葉片的形狀主要根據(jù)動力學(xué)性能和結(jié)構(gòu)設(shè)計來確定。為了盡可能地提高葉片效率，就要增大葉片的升力系數(shù)cl工作范圍和攻角變化范圍，同時結(jié)構(gòu)上也要盡量避免氣蝕的發(fā)生。從效率的角圖6 不同翼型的阻力系數(shù)曲線度看，要減小阻力系數(shù)cd，而從結(jié)構(gòu)強度上考慮，又需要有較厚的截面，尤其在葉片根部，而這又會導(dǎo)致較高的阻力系數(shù)，所以在選擇葉片翼型時，要有合適的彎度，同時在保證葉片結(jié)構(gòu)強度的前提下，應(yīng)選用較小的厚度。3 葉片的氣蝕特性研究不同壓力下，水的汽化溫度是不同的，在常溫下，當(dāng)葉片的局部壓力達(dá)到或低于水的汽化壓力時，葉片

11、就會發(fā)生氣蝕；這是因為葉片工作時，葉元體的迎流面的壓力較高，而葉背處壓力較低，同時不同葉輪半徑處的攻角和入流速度是不同的，迎流速度的提高和攻角的增大使葉背處的壓降更加明顯，所以氣蝕現(xiàn)象的出現(xiàn)與否可以根據(jù)葉片周圍局部壓力分布的大小來判斷。汽化系數(shù)常被用來表征葉片結(jié)構(gòu)是否發(fā)生氣蝕，它被定義為(6)而結(jié)構(gòu)物周圍的壓力系數(shù)定義為式中 p0 參考靜壓力，等于patm+gh，papatm大氣壓力，papv 汽化壓力，pa 水的密度，kg/m3h水的深度，mv 自由來流的速度，m/spl 局部壓力，pa圖8 不同翼型的壓力系數(shù)比較 當(dāng)plpv或cp時，就會出現(xiàn)氣蝕現(xiàn)象。由上式可見，由于1/v2，而v=r，所

12、以忽略了工作環(huán)境的溫度變化引起的pv變化后，與葉輪的轉(zhuǎn)速、直徑及葉片的浸水深度h有關(guān)。葉輪在旋轉(zhuǎn)一周過程中葉片的浸水深度是不一樣的(在大型的長葉片結(jié)構(gòu)中尤為明顯)，不同深度對應(yīng)的壓力也不同，所以葉片上的同一位置點的氣蝕數(shù)在旋轉(zhuǎn)過程中也是不一樣的，而且這種現(xiàn)象在葉尖處更為突出。圖9 汽化水體積分布圖圖7 naca63215翼型不同攻角時的壓力系數(shù)從理論上講，通過增加葉輪的工作深度及降低葉輪的轉(zhuǎn)速有助于避免氣蝕的發(fā)生，但受工作條件及安裝難度的限制，葉輪深度不可以太深；另一方面，對于大型水下風(fēng)車而言，如果僅僅為了避免氣蝕而將葉輪轉(zhuǎn)速控制在很低時，那么功率一定時，葉片將需要承受很大的扭矩；所以有必要研

13、究葉輪在高性能工作狀態(tài)下可以防止葉片氣蝕的合理結(jié)構(gòu)。下面仍以前述翼型為例，研究葉片結(jié)構(gòu)氣蝕的規(guī)律性。根據(jù)前面所述條件，并設(shè)定式（6）中，pat=101320 pa，為1025kg/m3，pv3175 pa，可求得10m水深處的氣蝕數(shù)為24.5。圖7為翼型為naca63215處于 10m水深處，攻角分別是0°、5°、10°和15°時的壓力系數(shù)曲線，可以看出隨著攻角的增大，系數(shù)cp逐漸增大，但仍低于氣蝕數(shù)，所以不會發(fā)生氣蝕。圖8分別為naca2415、naca4415、aca4424、 naca2411翼型的壓力系數(shù)曲線，從圖中可見，同樣攻角的條件下，隨著葉

14、片翼型厚度及彎度的減小，壓力系數(shù)逐漸增大，即氣蝕的可能性越來越大。圖9是naca2415翼型在攻角為15，流速為8m/s時的氣蝕圖，它展示的是汽化水在翼型上表面前緣后的體積分布。對于水下風(fēng)車的葉片結(jié)構(gòu)，要改善其氣蝕特性，一方面需要改善葉片翼型結(jié)構(gòu)，比如適當(dāng)增加葉片的厚度、增大彎度等，尋求探索新的翼型結(jié)構(gòu)；另一方面，采用抗剝蝕材料或抗剝蝕涂料，目前已知抗剝蝕性能好的材料是鈦合金，俄羅斯已有用鈦合金制作螺旋槳的例子，也可采用表面涂鈦的方法，國內(nèi)外都有可用于螺旋槳抗剝蝕的涂料7。4 結(jié)論通常設(shè)計元件的一般過程為設(shè)計、樣機性能試驗、制造。本文采用cfd方法通過計算機對葉片性能進(jìn)行了仿真研究，其結(jié)果與相

15、關(guān)文獻(xiàn)6的試驗結(jié)果相比較，認(rèn)為該方法能夠很好地在圖紙設(shè)計階段預(yù)測葉片的性能及其周圍流動產(chǎn)生的漩渦、氣穴、邊界層分離、尾流等不良現(xiàn)象，將可能發(fā)生故障的隱患消滅在圖紙設(shè)計階段。通過仿真表明，在水下風(fēng)車葉輪的工作環(huán)境中，大厚度葉片翼型不利于獲得高性能的升阻系數(shù)特性，但對防止氣蝕有較好的效果；彎度的減小可以減弱厚度對阻力系數(shù)的影響，但卻不利于升力系數(shù)的提高，也不利于改善葉片氣蝕特性；所以在葉片的設(shè)計過程中，要均衡二者利弊，實現(xiàn)二者最優(yōu)組合。參 考 文 獻(xiàn)1 j. vanzwieten, f.r. driscoll, a. leonessa and g. deane. design of a proto

16、type ocean current turbine part i-mathematical modeling and dynamics simulation. ocean engineering, in press, corrected proof, available online 19 january 2006.2 fraenkel, p. l. power from marine turbines. proc. instn mech. engrs, part a, 2002.3 barnsley mj, wellicome jf. wind tunnel investigation o

17、f stall aerodynamics for a 1.0m horizontal axis rotor. j wind eng ind aerodyn, 1992(39).4 a.s. bahaj, a.f. molland, j.r. chaplin, w.m.j. batten. power and thrust measurements of marine current turbines under various hydrodynamic flow conditions in a cavitation tunnel and a towing tank. renewable energy, 2007, 32(3):407-426.5 danmei hu, ouyang hua, zhaohui du. a study on stall-delay for horizontal axis wind turbine. renewable energy, 2006, 31(6)：821-836.6 molland, a.f., bahaj, a.s., chaplin, j.r., et al. measurements and predictions of