廖 微,王兵振,張 巍,段云琪
(國家海洋技術(shù)中心,天津 300112)
水平軸潮流能水輪機(jī)數(shù)值模擬方法對比分析
廖微,王兵振,張巍,段云琪
(國家海洋技術(shù)中心,天津300112)
葉素動量理論和CFD方法是水平軸潮流能水輪機(jī)性能分析中運(yùn)用較為廣泛的數(shù)值模擬方法,文中結(jié)合小尺寸水輪機(jī)模型試驗,對比分析了葉素動量理論和CFD方法在水輪機(jī)性能分析中的準(zhǔn)確性和適用性。驗證結(jié)果表明:葉素動量理論和CFD方法均能對水輪機(jī)的性能進(jìn)行預(yù)報,且CFD精度高于動量葉素理論;大尖速比時,動量葉素理論偏離較高,不再適合性能預(yù)報;在小尖速比下,建議采用RNG k-ε模型的CFD方法進(jìn)行分析計算;動量葉素理論適合設(shè)計初期設(shè)計方案的對比分析,而CFD方法適合對設(shè)計結(jié)果的驗證校核和詳細(xì)分析。
水平軸水輪機(jī);數(shù)值模擬;葉素動量理論;CFD方法;水池試驗
隨著經(jīng)濟(jì)發(fā)展與能源短缺、環(huán)境污染等矛盾的日益突出,潮流能作為一種能量密度高、儲備能量大、可再生的清潔能源受到各國政府青睞。據(jù)統(tǒng)計,我國的潮流能儲量達(dá)到13 948.5 MW[1],是世界上潮流資源較為豐富的國家之一,潮流能發(fā)電技術(shù)在我國有較好的應(yīng)用前景。
水輪機(jī)是潮流能系統(tǒng)中能量轉(zhuǎn)換的核心構(gòu)件,對水輪機(jī)性能的準(zhǔn)確預(yù)報是水輪機(jī)設(shè)計的關(guān)鍵。葉素動量理論和CFD方法是水輪機(jī)水動力分析中運(yùn)用最為廣泛的兩種方法。Bahaj[2]等采用葉素動量理論對水平軸水輪機(jī)進(jìn)行數(shù)值模擬及試驗驗證;王樹杰[3]等采用CFD對水輪機(jī)進(jìn)行數(shù)值模擬,對比分析了不同湍流模型的應(yīng)用;Ju Hyun[4]等將CFD計算結(jié)果與試驗值、VLM方法、BEM方法進(jìn)行了比較;盛其虎[5]等總結(jié)了潮流能水輪機(jī)水動力分析的常用方法。本文結(jié)合小尺寸水輪機(jī)模型試驗,對比分析了葉素動量理論和CFD方法在水輪機(jī)性能分析上的準(zhǔn)確性和適用性。
本文進(jìn)行分析驗證的水平軸潮流能水輪機(jī)模型的直徑為0.8 m,水輪機(jī)葉片翼型采用NACA63-424,設(shè)計尖速比為6,理論最佳功率系數(shù)為0.4,葉片數(shù)為2。葉片的設(shè)計結(jié)果如圖1所示。
圖1 模型葉片設(shè)計結(jié)果
葉素動量理論是水輪機(jī)設(shè)計與性能分析中應(yīng)用最為廣泛的方法,它結(jié)合葉素理論和動量理論,將葉片沿翼展方向分割為若干個葉素,并假定相鄰葉素間沒有干擾,葉素上的流動為二維流動,因此每個葉素所受的力與力矩只與該處葉素剖面翼型升阻力相關(guān),受力大小可采用葉素掃略圓環(huán)時發(fā)生的流體動量變化進(jìn)行計算。
采用葉素動量理論對水輪機(jī)性能進(jìn)行計算是一個對葉素處流體流動的軸向誘導(dǎo)系數(shù)和周向誘導(dǎo)系數(shù)不斷迭代求解的過程。單個葉素的流體速度和受力情況如圖2所示。假定葉素位于半徑為r的位置,葉素徑向?qū)挾葹閐r。根據(jù)葉素動量理論以及Prandtl修正[6],該葉素處流體流動的軸向誘導(dǎo)速度和周向誘導(dǎo)速度為:
圖2 葉素流體速度分布和受力示意圖
式中:a為軸向誘導(dǎo)系數(shù);b為周向誘導(dǎo)系數(shù);Ω為葉輪旋轉(zhuǎn)角速度;U為無窮遠(yuǎn)處來流速度。CL,CD分別為葉素翼形的升力和阻力系數(shù);c為葉素的弦長;Ftip為葉尖損失系數(shù);Fhub為輪轂損失系數(shù);Rhub為輪轂半徑;N為葉片數(shù);φ為水流速度和葉輪旋轉(zhuǎn)平面之間的角度即葉素入流角。
當(dāng)軸向誘導(dǎo)系數(shù)大于0.5時,葉輪工作在湍流狀態(tài),動量理論不再適用,需要對a進(jìn)行修正。本文采用文獻(xiàn)[7]提出的方法,通過推力系數(shù)CT的經(jīng)驗公式對a進(jìn)行修正:
當(dāng)求解得a,b后,即可求得葉素所受力、力矩與功率系數(shù),對所有葉素所求量進(jìn)行積分即可求解整個葉輪的受力與功率系數(shù)。
隨著計算機(jī)硬件水平的提高和計算流體力學(xué)的發(fā)展(CFD),采用基于 FLuent,CFX[8],STAR CCM+[9]等商業(yè)軟件的CFD方法得到了廣泛運(yùn)用。根據(jù)水輪機(jī)的工作環(huán)境,CFD求解模式通常為雷諾時均的不可壓縮粘性流體的定常流動,并引入湍流模型對流體運(yùn)動方程組進(jìn)行封閉求解。其中,在對潮流能水輪機(jī)的性能分析中,運(yùn)用最為廣泛的為k-ωSST湍流模型和RNG k-ε湍流模型。k-ωSST湍流模型被證明與試驗值具有較好的吻合度[9-10],其結(jié)合了標(biāo)準(zhǔn)k-ω湍流模型和標(biāo)準(zhǔn)k-ε湍流模型的優(yōu)勢,在處理近壁面時,采用標(biāo)準(zhǔn)k-ω湍流模型,在邊界層邊緣和自由剪切層時,采用k-ε湍流模型,從而提高了計算的準(zhǔn)確性[3]。RNG k-ε由于考慮了平均流動中的旋流,能更好地處理高應(yīng)變率和流線彎曲較大的流動。為對兩種常用湍流模型在水輪機(jī)水動力數(shù)值模擬中的準(zhǔn)確性進(jìn)行比較,本文在相同計算模型、計算網(wǎng)格和參數(shù)設(shè)置下,分別采用k-ωSST湍流模型和RNG k-ε湍流模型進(jìn)行數(shù)值模擬計算,并將計算結(jié)果與水槽試驗結(jié)果進(jìn)行對比分析。
利用Ansys建立設(shè)計0.8 m直徑潮流能葉輪的水動力計算模型,如圖3所示,計算域劃分為旋轉(zhuǎn)域和靜止域,旋轉(zhuǎn)域為包裹葉輪及其周圍的圓柱形區(qū)域,靜止域為距上游3D、下游6D、半徑2.5D的圓柱形區(qū)域,以使葉輪處尾流充分發(fā)展并盡量減小流域壁面效應(yīng)。在Fluent中,采用多參考坐標(biāo)系模型(MRF)進(jìn)行求解。
圖3 葉輪計算域
圖4 葉輪網(wǎng)格示意圖
網(wǎng)格質(zhì)量和邊界條件的選擇對計算結(jié)果具有很大的影響。利用網(wǎng)格劃分工具ICEM CFD進(jìn)行網(wǎng)格劃分,對旋轉(zhuǎn)域采用四面體非結(jié)構(gòu)網(wǎng)格,并在葉輪壁面附近進(jìn)行網(wǎng)格加密以保證葉輪壁面Y+值處于合適的范圍;對靜止域采用六面體結(jié)構(gòu)網(wǎng)格劃分以減小網(wǎng)格數(shù)量,加快收斂速度和計算精度。圖4為網(wǎng)格示意圖,其中旋轉(zhuǎn)域網(wǎng)格單元數(shù)為163.4萬,靜止域網(wǎng)格單元數(shù)為18.8萬,計算總單元數(shù)為182.2萬左右。在邊界條件的設(shè)置上,設(shè)置入口流體速度為恒值速度,采用湍流強(qiáng)度值和水利直徑定義湍流參數(shù),設(shè)置出口邊界條件為自由流出口或壓力出口邊界,葉輪壁面設(shè)置為無滑移壁面邊界條件。
4.1模型水池試驗
利用CFD中的k-ωSST湍流模型和RNG k-ε湍流模型以及葉素動量理論對設(shè)計0.8 m直徑葉輪進(jìn)行數(shù)值模擬計算,并與模型試驗結(jié)果進(jìn)行了對比,從而從計算準(zhǔn)確性和適用性兩方面對CFD數(shù)值模擬方法和葉素動量理論方法進(jìn)行對比,并分析兩種數(shù)值模擬方法的適用范圍。
為驗證葉素動量理論和CFD方法對水輪機(jī)性能分析的準(zhǔn)確性,本文對上述水平軸潮流能水輪機(jī)模型在天津大學(xué)船舶拖曳水池進(jìn)行了葉輪動力特性測試試驗。水池長140 m,寬7 m,水深3 m,試驗裝置主要由葉輪、傳動裝置、轉(zhuǎn)速/扭矩傳感器、發(fā)電機(jī)、控制電路等部分組成,如圖5所示。在試驗過程中將模型固定在拖車上,利用拖車拖動模型在水中運(yùn)動,保持拖車拖曳速度為2.0 m/s,通過改變轉(zhuǎn)軸輸出端的負(fù)載改變?nèi)~輪旋轉(zhuǎn)速度,以將轉(zhuǎn)速控制在230±10 rpm,250±10 rpm,280±10 rpm,300±10 rpm,350±10 rpm,400±10 rpm附近,并通過轉(zhuǎn)速扭矩傳感器測得輸出轉(zhuǎn)速、轉(zhuǎn)矩值,從而獲得不同尖速比下的轉(zhuǎn)速、轉(zhuǎn)矩值。
由于扭矩傳感器安裝在傳動系統(tǒng)的后端,因此,測量的扭矩數(shù)值受傳動系統(tǒng)的摩擦損失的影響,其數(shù)值會略小于葉輪輸出的扭矩值。為準(zhǔn)確地考察葉輪的輸出扭矩,減小傳動系統(tǒng)摩擦損失影響,設(shè)計了傳動系統(tǒng)摩擦損失測量裝置對傳動系統(tǒng)的摩擦阻力矩進(jìn)行了測量。
利用測量的數(shù)據(jù)計算不同尖速比下葉輪的功率系數(shù)。葉輪的功率系數(shù)按下式估算:
式中:Cp為葉輪的功率系數(shù),即能量轉(zhuǎn)換效率;Po為拖曳試驗時扭矩傳感器輸出的機(jī)械功率與傳動系統(tǒng)摩擦力矩產(chǎn)生的損失功率之和;Pi為水對葉輪作用所產(chǎn)生的輸入功率,按式(9)估算:
式中:ρ為水的密度;v為模型相對水的運(yùn)動速度;S為葉輪掃略面積。
圖5 試驗裝置
4.2結(jié)果分析
4.2.1準(zhǔn)確性為比較分析CFD方法和BEM方法在水平軸水輪機(jī)數(shù)值中的準(zhǔn)確性,本文首先對水輪機(jī)的CFD分析中最常用的兩種湍流模型k-ωSST模型和RNG k-ε模型進(jìn)行分析對比以確定適合的湍流模型,再將確定后的湍流模型計算結(jié)果與BEM方法進(jìn)行對比分析,驗證兩種數(shù)值計算方法的計算準(zhǔn)確性。在對水輪機(jī)進(jìn)行數(shù)值模擬計算時,為保證計算結(jié)果的可對比性,數(shù)值模擬計算采用的速度和尖速比值均與水槽試驗值保持一致。
在CFD分析中,為保證結(jié)果的可對比性,除湍流模型及對應(yīng)參數(shù)存在差異外,兩種湍流計算中的計算模型、網(wǎng)格劃分、參數(shù)設(shè)置等均保持一致。計算結(jié)果和試驗結(jié)果對比如圖6所示,兩種湍流模式對水平軸水輪機(jī)性能的預(yù)測均具有一定的準(zhǔn)確性。在最優(yōu)設(shè)計尖速比附近,兩者與實驗值的誤差均在1%以下,且k-ωSST模型計算精度較RNG k-ε模型稍高;在小尖速比下,RNG k-ε湍流模型計算結(jié)果較k-ωSST模型更吻合實驗值,但兩種的計算誤差均隨著尖速比的減小而增大,在TSR=4.77時,兩者誤差分別為4.9%和11.5%;當(dāng)尖速比大于最優(yōu)設(shè)計尖速比,k-ωSST和RNG k-ε的計算結(jié)果變化趨勢相同,且均在最大尖速比8.69下達(dá)到最大誤差,分別為11.3%和10.2%。整體上,在各尖速比下RNG計算值均小于SST計算值;SST湍流模型的整體精度高于RNG模型,但在小尖速比下精度低于RNG湍流模型。在下述與BEM計算方法的對比中,在小尖速比下采用RNG計算值、其他尖速比下采用SST計算值作為CFD計算值進(jìn)行對比分析。
圖6 RNG和SST計算值比較
圖7為BEM數(shù)值模擬方法與CFD數(shù)值模擬方法計算結(jié)果值與實驗值的對比。如圖所示,兩種數(shù)值模擬方法與試驗結(jié)果吻合,在設(shè)計尖速比附近計算精度最高。BEM、CFD、試驗測量的最佳尖速比均在0.65附近,其對應(yīng)的的最佳功率系數(shù)分別為0.415 9,0.396 6,0.395 9,BEM,CFD方法對最佳功率系數(shù)的預(yù)報誤差分別為5.04%和0.15%。在偏離設(shè)計尖速比時,兩種計算方法的誤差均增大:隨著尖速比減小,CFD計算誤差逐漸大于BEM計算誤差;隨著尖速比增大,BEM計算逐漸偏離試驗值且偏差較大。
圖7 試驗值與數(shù)值模擬值的比較
隨著尖速比偏離最佳尖速比,采用CFD計算得到的誤差均增大,其原因可能是:(1)試驗環(huán)境和CFD數(shù)值模擬環(huán)境的差異,如雷諾數(shù)、湍流強(qiáng)度等;(2)CFD計算方法本身的計算誤差:能較好地計算壓差阻力,但對粘性摩擦阻力的計算準(zhǔn)確性偏低;(3)試驗本身的誤差,如壁面效應(yīng)、自由表面的影響等。4.2.2適用性在本文中,葉素動量理論算法的實現(xiàn)基于Matlab編程,單一工況的計算CPU時間在1 s左右,采用葉素動量理論的時間成本主要集中在對葉素動量理論的MATLAB編程上,對所有工況計算的CPU時間在10 s內(nèi)。采用模塊化編程,可方便地對包括設(shè)計尖速比、葉片直徑、葉片數(shù)等葉片外形參數(shù)和計算參數(shù)等進(jìn)行更改計算。且編程計算對計算機(jī)資源要求較低,內(nèi)存占據(jù)少。CFD計算過程涉及三維建模、網(wǎng)格劃分、Fluent計算,時間成本主要集中在后兩步。網(wǎng)格劃分所需時間依賴于網(wǎng)格類型的選擇、網(wǎng)格精度、個人操作熟練度等,需要較多的時間和精力;在Fluent計算過程上,采用8核CPU、內(nèi)存3.43G的計算機(jī)做8核并行計算單一工況所需時間為2 h左右。一旦葉片外形參數(shù)發(fā)生改變,則需重新進(jìn)行三維建模以及后續(xù)過程,若只更改計算工況,則只需在Fluent中更改參數(shù)進(jìn)行計算。
在潮流能水輪機(jī)設(shè)計初期,需要進(jìn)行大量運(yùn)算對不同設(shè)計方案進(jìn)行計算對比,確定最佳葉片數(shù)、設(shè)計尖速比等設(shè)計參數(shù)。若采用CFD方法,則需對每一個方案重復(fù)建模、網(wǎng)格劃分、Fluent分析過程,需要耗費(fèi)大量的時間,采用葉素動量理論,可以在保證一定精度的基礎(chǔ)上節(jié)省時間。初步設(shè)計方案確定后,需要對單一設(shè)計結(jié)果進(jìn)行多工況的計算驗證和詳細(xì)分析,獲取包括流場分布、尾流反作用等詳細(xì)信息,而這些信息無法從葉素動量理論中計算獲得,需采用更復(fù)雜、精度更高的CFD方法。
本文結(jié)合小尺寸水輪機(jī)模型試驗,對比分析了葉素動量理論、CFD方法在水平軸水輪機(jī)性能分析中的準(zhǔn)確性和適用性,結(jié)果說明:
(1)采用葉素動量理論和CFD方法均能較為準(zhǔn)確地分析水輪機(jī)的性能,且CFD的計算精度要高于葉素動量理論的計算精度。
(2)在最佳尖速比附近,兩種數(shù)值方法計算精度最高;在小尖速比下,葉素動量理論的計算精度高于k-ωSST模型的精度,但低于RNG k-ε模型,建議采用RNG k-ε模型的CFD方法進(jìn)行計算;在大尖速比下,BEM計算結(jié)果與試驗偏差較大。
(3)BEM方法計算成本低、精度適中,適合水輪機(jī)設(shè)計初期對性能的分析;CFD計算成本高,計算精度較高,結(jié)果詳細(xì),適合設(shè)計后期對設(shè)計方案的驗證分析。
[1]戴慶忠.潮流能發(fā)電及潮流能發(fā)電裝置[J].東方電機(jī),2010(2):51-66.
[2]ASBahaj,WGJ Batten,GMcCann.Experimental Verification ofNumerical Predictions for the Hydrodynamic Performance ofHorizontal Axis Marine Current Turbine[J].Renewable Energy,2007(32):2479-2490.
[3]王樹杰,盛傳明,袁鵬,等.潮流能水平軸水輪機(jī)湍流模型研究初探[J].中國海洋大學(xué)學(xué)報,2014,44(5):95-100.
[4]Ju Hyun Lee,SunhoPark,DongHwan Kim,et al.Computational methods for performance analysis ofhorizontal axis tidal stream turbines[J].Applied Energy,2012,98:512-523.
[5]盛其虎,孫科,張學(xué)偉,等.潮流能水輪機(jī)水動力分析方法[J].海洋技術(shù)學(xué)報,2014,33(4):105-111.
[6]Prandtl L,Tietjens O K G.Applied Hydroand Aeromechanics,based on Lectures ofI.Prandtl,PhD[M].Mineola:Dovel Publications.com,1957.
[7]Ju Hyun Lee,Dong Hwan Kim,Shin Hyung Rhee,et al.Comptuational and Experimental Analysis for Horizontal Axis Marine Current Turbine Design[C]//Second International Symposiumon Marine Propulsorsamp’ll.Hamburg,Gemany,2011.
[8]盛其虎,趙亞東,張亮.水平軸潮流能水輪機(jī)設(shè)計與數(shù)值模擬[J].哈爾濱工程大學(xué)學(xué)報,2014,35(4):389-394.
[9]MJ Lawson,YLi,D C Sale.Development and Verification ofa Computational Fluid Dynamics Model ofa Horizontal-Axis Tidal Current Turbine[R].Netherlands,2011.
[10]Pape A,Lecanu J.3DNavier-Stokes Computations ofa Stall-regulated Wind Turbine[J].Wind Energy,2004,7(4):309-324.
Comparative Analysis on the Numerical Simulation Methods for Horizontal Axis Tidal Current Turbines
LIAO Wei,WANG Bing-zhen,ZHANG Wei,DUAN Yun-qi
National Ocean Technology Center,Tianjin 300112,China
This paper conducts a comparative analysis on the accuracy and applicability of the blade element momentum(BEM)theory and computational fluid dynamics(CFD)method,the two most widely used numerical simulation methods for hydrodynamic performance analysis of horizontal axis tidal current turbine.To verify the computational results,an open water test is carried out in a towing tank with a self-designed turbine.The results show that both the two methods can provide a satisfactory prediction for the experimental turbine performance,and a higher accuracy is acquired with CFD.With increasing tip speed ratio(TSR),higher deviation is observed with the BEM theory,not applicable for prediction anymore.For low TSR,in contrast,it is suggested that the CFD method of the RNG model be adopted for analysis and calculation.Furthermore,the use of BEM is suited to the comparative analysis of the design scheme under early stage,while CFD is for post-design validation and indepth analysis.
horizontal axis tidal current turbine;numerical simulation;blade element momentum theory;computational fluid dynamic(CFD)method;towing tank test
P743
A
1003-2029(2016)01-0113-05
10.3969/j.issn.1003-2029.2016.01.018
2015-03-14
海洋公益性行業(yè)科研專項經(jīng)費(fèi)資助項目(201205019-3)
廖微(1990-),在讀碩士,主要研究方向為港口、海岸與近海工程。E-mail:lanmengxinfei555@163.com