呂強(qiáng) 張文斌 郭勇 陳如華 王曉暉 匡開林

摘要：為解決用于水表井中的遠(yuǎn)傳設(shè)備供電不足，對水表管道結(jié)構(gòu)進(jìn)行適當(dāng)改造，引入垂直型螺旋式水輪機(jī)組對鋰電池進(jìn)行充電，使得設(shè)備處于最佳運(yùn)行狀態(tài)。本文將針對90°和120°兩組不同包角的葉片在不同工況下進(jìn)行全流道數(shù)值模擬，分析內(nèi)流狀態(tài)并進(jìn)行選優(yōu)。結(jié)果表明，當(dāng)轉(zhuǎn)速由100r/min到1000r/min過程中，90°包角渦輪出現(xiàn)反轉(zhuǎn)，但是回收效率逐漸增加，120°包角渦輪正轉(zhuǎn)，且功率呈拋物線趨勢。在轉(zhuǎn)速約為500r/min時，120°包角渦輪達(dá)到最優(yōu)工況，在此工況下120°包角渦輪內(nèi)部流動湍動能較小，能量耗散弱，運(yùn)行更加穩(wěn)定。

關(guān)鍵詞：垂直型螺旋式水輪機(jī);包角;數(shù)值模擬;湍動能

中圖分類號：TK73文獻(xiàn)標(biāo)志碼：A

0 引言

目前自動抄表系統(tǒng)在供電、天然氣等領(lǐng)域已經(jīng)廣泛應(yīng)用，自來水自動抄表的應(yīng)用也在普及，影響其推廣的主要原因就是設(shè)備供電問題。大口徑水表安裝在水表井中，環(huán)境潮濕，交流電源無法引入。在水表井中用電設(shè)備有數(shù)據(jù)采集模塊和遠(yuǎn)傳發(fā)射終端兩部分，由于安裝環(huán)境的制約，數(shù)據(jù)采集模塊和遠(yuǎn)傳終端供電只能采用內(nèi)置干電池供電。

遠(yuǎn)傳終端是將采集模塊采集到的數(shù)據(jù)運(yùn)用GPRS技術(shù)發(fā)送到互聯(lián)網(wǎng)，通過對應(yīng)的網(wǎng)頁和軟件顯示出來，做為收費(fèi)的依據(jù)。遠(yuǎn)傳終端設(shè)備內(nèi)部裝有一組電池，在設(shè)計時一方面需要盡量加大電池的容量，另一方面需要減小設(shè)備的功耗。

為解決以上問題帶來的不便，本文針對自來水表這一流量監(jiān)測設(shè)備對管道進(jìn)行簡易改進(jìn)，引入垂直型螺旋式水輪發(fā)電機(jī)對蓄電池進(jìn)行不間斷供電，確保各種遠(yuǎn)傳監(jiān)測設(shè)備處于最佳運(yùn)行狀態(tài)。本研究針對兩種不同包角的渦輪采用ANYSIS CFX數(shù)值模擬軟件進(jìn)行全流道數(shù)值模擬，分析不同包角的葉片在不同工況下回收效率，為管道式超微型水輪機(jī)組選擇最優(yōu)渦輪，對未來管道供電系統(tǒng)的研發(fā)工作有一定指導(dǎo)作用。

所設(shè)計發(fā)電裝置裝配圖如圖1所示，該裝置由線圈槽、磁鐵槽、渦輪和螺紋主體四大部分組成，管道來流沖擊渦輪驅(qū)動渦輪旋轉(zhuǎn)，通過連接軸帶動磁鐵槽中磁鐵來切割磁感線產(chǎn)生感應(yīng)電動勢，從而產(chǎn)生電流來為儲蓄電池供電，而螺紋主體則是一種密封裝置，內(nèi)置傳動軸傳輸渦輪產(chǎn)生的扭矩，該裝置結(jié)構(gòu)簡單，便于安裝和維護(hù)，滿足在不同環(huán)境下自來水管道用電設(shè)備的需求。

1葉輪幾何造型

本文研究對象為一直徑為D=50mm的自來水管道，長度H=150mm，管道內(nèi)壓力P=0.5MPa，流量Q=25m/h。其葉輪設(shè)計參數(shù)為：葉輪直徑d=16mm，葉輪寬度b=15mm，葉片數(shù)Z=3，葉片包角分為90°和120°，所設(shè)計渦輪物理模型如圖2及圖3所示。

2 數(shù)值計算

本文采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，分別對兩種不同包角的渦輪進(jìn)行了網(wǎng)格無關(guān)性研究，兩種渦輪在網(wǎng)格總數(shù)均大于300萬時其效率較為穩(wěn)定，波動較小，因此將兩種模型的總網(wǎng)格數(shù)控制在300萬左右。本文90°包角模型總網(wǎng)格數(shù)量為326萬4489，其節(jié)點(diǎn)數(shù)量為57萬6836，120°包角模型總網(wǎng)格數(shù)量為327萬2654，其節(jié)點(diǎn)數(shù)量為57萬8275。

利用ANYSIS-CFX軟件對兩種不同包角的渦輪模型進(jìn)行全流道數(shù)值模擬，整個計算域如圖4所示。分析類型為穩(wěn)態(tài)，湍流模型選用k-ε模型。采用速度進(jìn)口，壓力出口，輸送介質(zhì)為25℃清水，計算收斂標(biāo)準(zhǔn)為10，動靜交界面采用Frozen rotor，固壁采用無滑移邊界條件。

3計算結(jié)果分析

3.1外特性分析

采用ANYSIS CFX對兩種不同包角的模型從轉(zhuǎn)速100r/min到1000r/min進(jìn)行數(shù)值模擬，因?yàn)樽詠硭艿缐毫νǔ］^為恒定，所以我們給定進(jìn)口速度為3.5m/s，出口壓力為0.5MPa。對兩組不同包角渦輪的扭矩展開分析討論，其結(jié)果如圖5和圖6所示，圖5為兩種渦輪的輸出扭矩特性，從圖中可以看出，90°包角渦輪出現(xiàn)反轉(zhuǎn)，扭矩隨轉(zhuǎn)速增加而增大，而120°包角渦輪隨轉(zhuǎn)速增加扭矩呈下降趨勢，最后停止轉(zhuǎn)動且出現(xiàn)反轉(zhuǎn)，說明兩種模型由于設(shè)計結(jié)構(gòu)或者安裝位置的原因出現(xiàn)葉片背面受力大于工作面受力，在運(yùn)行過程中背面承受壓力對渦輪造成較大能量損失。

圖6是兩種不同包角渦輪的輸出功率特性曲線，從圖中可以看出90°包角渦輪雖然出現(xiàn)反轉(zhuǎn)，但是回收功率卻隨轉(zhuǎn)速逐漸增加;120°包角渦輪回收功率隨轉(zhuǎn)速呈拋物線分布，轉(zhuǎn)速大約為500r/min時，回收功率達(dá)到最高。

3.2 內(nèi)流場分析

90°包角渦輪與120°包角渦輪在轉(zhuǎn)速為500r/min時得扭矩和回收功率相近，因此我們對兩組渦輪在轉(zhuǎn)速為500r/min時展開內(nèi)流分析及討論：

3.2.1壓力分布

圖7和圖8為90度和120度包角渦輪的壓力分布云圖，從渦輪局部壓力分布來看，90度包角渦輪壓力集中在葉片的背面，造成葉片背面所受壓力始終大于工作面所受壓力，所以渦輪工作開始就出現(xiàn)了反轉(zhuǎn)現(xiàn)象，且隨著轉(zhuǎn)速的增加，葉片背面所受壓力逐漸大于工作面所受壓力，渦輪功率逐漸增加。120度包角渦輪壓力集中在渦輪的工作面，在轉(zhuǎn)速為900r/min之前，葉片工作面所受壓力始終大于葉片背面所受壓力，渦輪處于正轉(zhuǎn)狀態(tài)，當(dāng)轉(zhuǎn)速大于900r/min時壓力大小分布相反，此時也出現(xiàn)反轉(zhuǎn)現(xiàn)象。

3.2.2速度矢量分析

圖9與圖10為兩種模型的速度矢量圖，來流由上而下沖擊渦輪，從圖中可以看到來流沖擊到葉片之后的流動狀態(tài)，兩種渦輪都在來流的沖擊下在工作面產(chǎn)生了較強(qiáng)的漩渦。90度包角渦輪由于是反轉(zhuǎn)狀態(tài)，故產(chǎn)生的漩渦強(qiáng)度大，能量耗散強(qiáng)度大，雖然功率隨轉(zhuǎn)速增加而增大，但是造成了管道能量的損失。120度包角渦輪流動則顯得更為紊亂，這是由于增大了包角之后，葉片對流體的束縛能力增強(qiáng)，但是又造成了工作面內(nèi)流體形成二次流，轉(zhuǎn)速增加，葉片工作面與背面的受力則逐漸趨于相等直至出現(xiàn)反轉(zhuǎn)的狀態(tài)，所以該包角渦輪的功率呈拋物線形分布，在轉(zhuǎn)速約為500r/min之后效率持續(xù)下降直至反轉(zhuǎn)。

4 總結(jié)

通過對兩種不同包角渦輪在不同轉(zhuǎn)速下進(jìn)行全流道數(shù)值模型，分析了其外特性及內(nèi)部流動狀態(tài)，從而確定了垂直型螺旋式渦輪包角對性能的影響，主要得出以下結(jié)論：

（1）在來流一定的情況下，90度包角渦輪工作面受力小于葉片背面受力，造成渦輪反轉(zhuǎn)，但效率隨轉(zhuǎn)速逐漸增加。120度包角渦輪則隨轉(zhuǎn)速增加葉片工作面受力與背面逐漸趨于平衡，最終出現(xiàn)反轉(zhuǎn)，效率在500r/min時最優(yōu)。

（2）120°包角渦輪相對于90°包角葉片對流體則表現(xiàn)出更強(qiáng)的約束力，但由于結(jié)構(gòu)為螺旋結(jié)構(gòu)，增大包角造成流體在工作面內(nèi)形成二次流，產(chǎn)生較強(qiáng)的漩渦，耗散了能量，而90度包角渦輪由于反轉(zhuǎn)造成了更大的能量損失。

（3）兩種渦輪在選型時應(yīng)根據(jù)實(shí)際情況，當(dāng)轉(zhuǎn)速小于500r/min時則選取120度包角渦輪，此時能量損失較弱，效率較高;當(dāng)轉(zhuǎn)速大于500r/min時，則選取90包角渦輪，對應(yīng)輸出功率大。