郭朋華,王逸遠,張大禹,鄒瀚森

(西安交通大學能源與動力工程學院,710049,西安)

隨著人類對深海探索需求的不斷增加,出現(xiàn)了由深海聲納、無人潛航器等組成的深海觀測網(wǎng)[1-3],傳統(tǒng)的陸基和船基電纜供電方式存在使用成本高、電纜易被切斷的問題,因此利用深海能源實現(xiàn)相關(guān)設(shè)備長周期的原位供電十分必要。深海中蘊含著多種能量形式,諸如溫差能、鹽差能、生物質(zhì)能等[4-5],其中只有海流能具備足夠大的能量密度以滿足相關(guān)設(shè)備的電力需求[6]。

深海環(huán)境下流速很低,遠小于常規(guī)海流能發(fā)電技術(shù)所采用的1.5 m·s-1以上的額定設(shè)計流速[7]。傳統(tǒng)升力型水輪機低流速下的自啟動能力較差,幾乎很難在深海區(qū)域?qū)崿F(xiàn)自啟動[8-9];阻力型垂直軸水輪機雖然具有較好的自啟動能力,但是其能量利用系數(shù)較低[10-14]。阿基米德水輪機是一種極具潛力的阻力型水平軸水輪機,具有良好的自啟動性能且能量利用系數(shù)接近升力型水平軸水輪機[15],可以有效利用深海超低流速的海流能。

在實際應用中,水輪機通常采用陣列形式布置,輪機間的耦合效應可能會提高水輪機的自啟動能力和能量利用系數(shù)[16]。傳統(tǒng)升力型水輪機陣列的機組間干擾較大,排布方式對功率影響顯著。Li等[17]研究了達里厄水輪機之間的水動力耦合規(guī)律,發(fā)現(xiàn)串聯(lián)安裝的水輪機陣列會大大降低下游葉片的功率系數(shù),而并聯(lián)安裝時可以通過減少間距進而增加機群的總功率系數(shù)。袁鵬等[18]通過動量損失的方法對潮流能水輪機陣列進行研究,結(jié)果表明,水輪機并聯(lián)間距小于4倍輪機直徑時,可以更好地利用單位面積的潮流能。相較水平軸水輪機,垂直軸水輪機受湍動能影響較小,輪機之間耦合增益效果更顯著。Antheaume等[19]研究了不同數(shù)量的垂直軸水輪機并聯(lián)排布的影響,發(fā)現(xiàn)合理的排布方式能提高單臺水輪機40.5%的能量利用系數(shù)。陣列中水輪機不同的旋轉(zhuǎn)方向?qū)啓C間的耦合效應也有影響。安琪[20]分析了不同排布方案對水輪機陣列功率及性能的影響規(guī)律,發(fā)現(xiàn)旋轉(zhuǎn)方向相反的水輪機可以提高陣列整體性能。Bai等[21]數(shù)值研究了不同排布方式下的水平軸水輪機轉(zhuǎn)向和間距對性能的影響,結(jié)果表明,直線網(wǎng)格陣列中反向旋轉(zhuǎn)對水輪機的影響大于交錯網(wǎng)格陣列的影響,并且反向旋轉(zhuǎn)的影響隨著水輪機下游距離的增大而迅速減小。部分學者基于對尾跡區(qū)的流場分析,認為交錯排布陣列可以利用被加速的尾流,有可能提高水輪機的能量利用系數(shù)。于書帆等[22]數(shù)值研究了雙水輪機在不同軸間距、相對位置角和不同尖速比情況下的各項性能參數(shù),與孤立輪機相比,雙水輪機交錯排布的能量利用系數(shù)大于孤立輪機的能量利用系數(shù),但由于尾流渦列區(qū)較長,尾流速度恢復所需距離也較長。

綜上所述,水輪機的陣列布置可能帶來耦合增益效應,從而提升整體機群的功率輸出,然而迄今為止沒有文獻對阿基米德水輪機陣列進行研究。因此,本文通過數(shù)值模擬的方式開展阿基米德水輪機的陣列耦合增益研究,為該類型輪機陣列的優(yōu)化布置策略提供參考。

1 數(shù)值模型

1.1 幾何模型建立

本文采用的阿基米德水輪機如圖1所示,本項目組前期對孤立輪機進行了優(yōu)化設(shè)計,其具體幾何參數(shù)如表1所示。

表1 阿基米德水輪機幾何參數(shù)

圖1 阿基米德水輪機幾何參數(shù)示意圖Fig.1 Geometrical parameters of the Archimedes turbine

1.2 計算域建立

利用ANSYS SpaceClaim建立計算域,坐標原點位于螺旋1葉片底部。正對來流方向的流域橫截面尺寸為12d×6d,螺旋葉片前沿距入口2d,后端距出口10d,具體尺寸如圖2所示。

(a)并聯(lián)陣列

(b)串聯(lián)陣列

1.3 計算設(shè)置及邊界條件

采用四面體網(wǎng)格劃分,并對阿基米德螺旋葉片周圍進行局部網(wǎng)格加密,邊界層為20層,調(diào)整第一層邊界層的厚度,以得到接近1的y+值。在計算設(shè)置上,采用二階迎風算法,湍流模型為SSTk-ω模型。選取水作為流體介質(zhì),密度為998.2 kg·m-3,動力黏度為10-3Pa·s,入口速度為0.5 m·s-1,出口設(shè)置為壓力出口,相對壓力為0。速度分量和湍流性質(zhì)均采用10-5的收斂標準。

水輪機的能量利用系數(shù)和啟動力矩系數(shù)是水輪機重要的性能特征,尤其是部署在深海的水輪機,良好的啟動力矩系數(shù)可以在一定程度上簡化控制系統(tǒng)和額外的機械裝置。能量利用系數(shù)Cp、啟動力矩系數(shù)Ct,s和葉尖速比λ計算式為

(1)

(2)

(3)

式中:T為力矩,N·m;ω為轉(zhuǎn)速,rad·s-1;ρ為流體密度,kg·m-3;R為水輪機半徑,m;V為來流速度,m·s-1。

1.4 數(shù)值模型無關(guān)性驗證

為了保證數(shù)值模擬結(jié)果的準確性,本文進行了網(wǎng)格無關(guān)性驗證、時間步長無關(guān)性驗證和計算區(qū)域無關(guān)性驗證。對單個阿基米德螺旋輪機進行不同網(wǎng)格尺寸劃分,以性能指標Cp、Ct,s作為判斷指標,網(wǎng)格數(shù)量及計算結(jié)果如表2所示。當網(wǎng)格數(shù)量為133萬時,Cp、Ct,s的波動幅度低于1.5%,繼續(xù)增大網(wǎng)格數(shù)量對計算結(jié)果影響較小。時間步長的影響如表3所示,當時間步長為0.01 s時,Cp、Ct,s的波動幅度低于1%,縮小時間步長對計算結(jié)果影響不大。對不同尺寸的計算區(qū)域進行無關(guān)性驗證結(jié)果如表4所示,選擇圖2所示計算區(qū)域尺寸為參考值。當計算區(qū)域縮放比例為1.1時,Cp、Ct,s的波動幅度低于1%,擴大計算區(qū)域面積對計算結(jié)果影響不大。故在保證計算精度的前提下本文依據(jù)圖2的計算區(qū)域,選擇133萬級別的網(wǎng)格和0.01 s的時間步長進行數(shù)值模擬。

表2 網(wǎng)格無關(guān)性驗證

表3 時間步長無關(guān)性驗證

表4 計算區(qū)域無關(guān)性驗證

1.5 實驗驗證

為了進一步驗證數(shù)值模型的準確性,對阿基米德水輪機進行了循環(huán)水槽實驗。實驗裝置示意圖如圖3所示,實驗段來流面積為0.5 m×0.8 m,長為1 m,實驗流速為0.5 m·s-1。阿基米德水輪機的螺桿水平置于實驗段的中心,阻塞比為17.7%,采用鈍體水槽修正方法[23-24]補償阻塞比帶來的影響,結(jié)果對比如圖4所示,數(shù)值模擬結(jié)果與修正后的實驗數(shù)據(jù)吻合較好,驗證了本文數(shù)值模擬方法的可靠性。

圖3 實驗裝置示意圖Fig.3 Diagram of experimental device

圖4 數(shù)值模擬結(jié)果與實驗數(shù)據(jù)對比Fig.4 Validation of Cp between experiment and simulation

1.6 速度控制策略

以往文獻中關(guān)于水輪機陣列研究,通常將陣列中的輪機設(shè)置為相同的轉(zhuǎn)速以便計算,然而實際中由于不同位置處的水輪機來流速度不同,因此陣列中各輪機的最優(yōu)工作轉(zhuǎn)速也應該是不同的。本文通過計算得到不同流速下最佳功率輸出所對應的轉(zhuǎn)速,采用多項式擬合的方法獲得孤立輪機的最佳功率曲線。利用Fluent中用戶自定義函數(shù)UDF對陣列中的孤立輪機進行轉(zhuǎn)速控制,通過迭代計算使轉(zhuǎn)速匹配到最佳功率,計算流程如圖5所示,其中系數(shù)隨時間變化:當t≤10時,r=0.5;當1050時,r=0.1。

圖5 速度控制框圖Fig.5 Flow chart of the speed control method

2 并聯(lián)與串聯(lián)陣列耦合增益研究

2.1 并聯(lián)陣列

2.1.1 并聯(lián)陣列能量利用系數(shù)

對并聯(lián)的阿基米德水輪機陣列,分別進行同向和反向旋轉(zhuǎn)布置,計算結(jié)果如圖6所示。

阿基米德水輪機并聯(lián)陣列中單個輪機的Cp均有不同程度的提高,同向和反向旋轉(zhuǎn)布置的輪機陣列的Cp在距離Dx=1.3d時均達到峰值,分別較孤立輪機的Cp提升了5%、7%。隨著距離Dx的增加,在1.5d之后兩種旋轉(zhuǎn)方向的水輪機的Cp均逐漸減小,直至接近孤立輪機的Cp。

(a)同向旋轉(zhuǎn)

(b)反向旋轉(zhuǎn)

并聯(lián)陣列間距Dx=1.3d時的流場如圖7所示,圖中上方輪機為螺旋2,下方輪機為螺旋1,兩臺并聯(lián)輪機的周圍形成高速區(qū),中間流域的加速效果尤為顯著。由于阿基米德水輪機具有一定的軸向長度,所以受阻塞效應影響,阿基米德水輪機會在入流面和側(cè)面持續(xù)受到位于中間流域加速區(qū)的影響,從而獲得高于孤立輪機的性能參數(shù)。此外,根據(jù)角動量

(a)同向旋轉(zhuǎn)

(b)反向旋轉(zhuǎn)

守恒,旋轉(zhuǎn)輪機的周圍會產(chǎn)生與之轉(zhuǎn)向相反的尾流,同樣由于阿基米德水輪機的軸向長度,相反轉(zhuǎn)向的尾流會出現(xiàn)在水輪機的側(cè)面和尾流面。反向旋轉(zhuǎn)陣列中的輪機尾流會與相鄰輪機產(chǎn)生嚙合效應,所以在輪機間距較小時,反向旋轉(zhuǎn)陣列的能量利用系數(shù)大于同向旋轉(zhuǎn)陣列的能量利用系數(shù)。

隨著Dx的不斷增加,這種增益效果逐漸減小,當Dx=4d時,Cp較孤立輪機僅提高2%;當Dx小于1.3d時,輪機之間的干涉影響逐漸顯著。

2.1.2 并聯(lián)陣列啟動力矩系數(shù)

并聯(lián)陣列的Ct,s如圖8所示,相比反向安裝,同向安裝陣列的自啟動性能更為優(yōu)越,其Ct,s大于孤立輪機,同樣是在間距為Dx=1.3d附近達到峰值,螺旋2的Ct,s較孤立輪機提升了19%。隨著距離Dx的不斷增大,Ct,s逐漸接近孤立輪機。綜合考慮Ct,s、Cp,并聯(lián)陣列的間距應該選擇在Dx=1.3d左右。

(a)同向安裝

(b)反向安裝

2.2 串聯(lián)陣列

2.2.1 串聯(lián)陣列能量利用系數(shù)

對串聯(lián)的阿基米德水輪機陣列,同樣進行同向和反向旋轉(zhuǎn)布置,結(jié)果如圖9所示。

(a)同向旋轉(zhuǎn)

(b)反向旋轉(zhuǎn)

圖10 串聯(lián)陣列同向旋轉(zhuǎn)速度云圖Fig.10 Velocity contour of the serial array with the same rotation direction

阿基米德水輪機串聯(lián)放置時,隨著間距的不斷增加,上游輪機在Dz=3d時恢復到孤立輪機的Cp,此后基本不受下游輪機的影響。下游輪機的Cp隨距離的增加先快速上升,當間距Dz=5d時恢復到孤立輪機Cp的70%,隨后上升趨勢明顯變緩,當間距Dz增加到10d時,Cp恢復到孤立輪機的80%,對應的流場尾跡速度云圖如圖10所示。雖然在Dz=5d處的流場仍存在較強的尾跡渦,但是下游阿基米德輪機已經(jīng)能夠恢復至較高的Cp,這是由于阿基米德水輪機運行尖速比較低,受湍動能的影響較小。而在較高尖速比運行下的升力型水平軸水輪機中,串聯(lián)陣列要獲得類似的恢復效果通常間距要在8d以上[25],以上結(jié)果表明阿基米德水輪機可以組成更為密集的串聯(lián)陣列,從而最大化陣列的能量密度。

2.2.2 串聯(lián)陣列啟動力矩系數(shù)

不同放置間距下的Ct,s如圖11所示。無論是同向還是反向安裝,上游輪機在Dz=3d時Ct,s均出現(xiàn)最小值,隨后逐步恢復至孤立輪機的Ct,s。下游輪機的Ct,s隨間距的變化情況與Cp類似,在間距Dz=5d之前快速上升,隨后增長趨勢變緩。特別需要注意的是,在距離Dz小于2d時下游輪機的Ct,s極低,對于同向安裝間距為1.5d、2d時甚至出現(xiàn)了Ct,s為負數(shù)的工況,此時上游輪機的尾流效應會嚴重影響下游輪機的自啟動性能。

(a)同向安裝

(b)反向安裝

3 交錯排布水輪機陣列性能研究

通過對并聯(lián)陣列和串聯(lián)陣列的研究,可以發(fā)現(xiàn)上游輪機的尾流效應和相鄰輪機間的阻塞效應是影響水輪機Cp、Ct,s的主要因素。因此,可以采用交錯的陣列排布方式,一方面可以減少尾流低速區(qū)的影響,另一方面可以有效利用阻塞的耦合增益效應。

本文以典型的三角陣列排布方式為研究對象,以垂直于來流方向間距Dx和平行于來流方向間距Dz為研究變量,分析了兩種不同旋轉(zhuǎn)方案下三角形陣列的耦合增益效應,兩種旋轉(zhuǎn)方案如圖12所示。

(a)同向旋轉(zhuǎn)

(4)

同向、反向旋轉(zhuǎn)陣列中各輪機Cp分別如圖13、圖14所示。對于同向旋轉(zhuǎn)布置,當列間距Dx小于d時,過小的間距使得下游輪機處于尾流效應明顯的區(qū)域,此時上游輪機尾跡的不利影響超過了輪機間的耦合增益效應,使得下游輪機的Cp均低于孤立輪機;當Dx=1.2d時,只有在行間距Dz小于3d的情況下才會高于孤立輪機。隨著列間距Dx的繼續(xù)增大,下游輪機的Cp均高于孤立輪機,當Dx=2d時,Cp出現(xiàn)峰值,較孤立輪機提升6%,此時尾流效應影響較小,而下游輪機之間以及下游與上游輪機之間均存在較強的耦合增益效應。反向旋轉(zhuǎn)布置時下游輪機Cp的變化趨勢與同向旋轉(zhuǎn)相似,但峰值出現(xiàn)在Dx=1.2d處,較孤立輪機提升8%。

(a)螺旋2能量利用系數(shù)

(b)螺旋3能量利用系數(shù)

(a)螺旋2能量利用系數(shù)

(b)螺旋3能量利用系數(shù)

(a)同向旋轉(zhuǎn)

(b)反向旋轉(zhuǎn)

Fig.15 Diagram of coupling gain boundary of triangular arrays

4 結(jié) 論

(1)受益于耦合增益效應,并聯(lián)陣列條件下輪機的功率輸出和自啟動性能均有提升。反向旋轉(zhuǎn)時耦合增益最大,在間距為1.3d時Cp提升至孤立輪機的107%;同向旋轉(zhuǎn)時輪機具有更高的Ct,s,在間距為1.3d時可以提升至孤立輪機的119%。

(2)受尾跡效應影響,串聯(lián)陣列中上游輪機的Cp、Ct,s略有降低,在間距大于5d時恢復至接近孤立輪機的水平;下游輪機的Cp、Ct,s均出現(xiàn)明顯下降,在間距大于5d之后恢復趨勢變緩,Cp恢復至孤立輪機Cp的70%。

(3)受耦合增益和尾跡效應的綜合作用,三角形陣列在保持一定間距的條件下,可以使Cp得到8%的提升。上游輪機和下游輪機進行反向旋轉(zhuǎn)布置既可以獲得較大的功率提升,又能夠?qū)崿F(xiàn)緊湊布置。