沈 娜,韓鳳琴

(華南理工大學(xué)廣州學(xué)院電氣工程學(xué)院,廣東廣州510800)

射流與旋轉(zhuǎn)水斗時空非定常三維流動研究

沈 娜,韓鳳琴

(華南理工大學(xué)廣州學(xué)院電氣工程學(xué)院,廣東廣州510800)

針對噴嘴射流與旋轉(zhuǎn)水斗時空非定常流動關(guān)系進(jìn)行研究,將二維“動畫解析法”完善至三維空間,旨在更準(zhǔn)確地追及射流干涉的影響。分析了射入當(dāng)前水斗入射射流分支的幾何結(jié)構(gòu)以及偏離射流分支的后緣面,三維可視化描述了自由射流和旋轉(zhuǎn)水斗的復(fù)雜非定常流動,得到了自由射流與旋轉(zhuǎn)水斗之間非定常流動的相互作用關(guān)系,進(jìn)而提出了動態(tài)水斗切入射流做瞬態(tài)復(fù)雜運(yùn)動時,射流及水斗背面發(fā)生干涉現(xiàn)象的機(jī)理。分析結(jié)果表明,干涉引起水輪機(jī)性能變化正比于單位轉(zhuǎn)速、射流膨脹率的增加,單位轉(zhuǎn)速變化對于水斗干涉的影響要遠(yuǎn)大于射流直徑膨脹變化帶來的影響。

非定常流動；沖擊式水輪機(jī)；噴嘴射流；旋轉(zhuǎn)水斗；干涉現(xiàn)象

0 引 言

與反擊式水輪機(jī)如混流式、軸流式水輪機(jī)不同,沖擊式水輪機(jī)的流動在時間和空間上都體現(xiàn)出典型的非定常流動特性[1-3]。由于射流和轉(zhuǎn)動水斗的相對位置是時刻變化的,使得沖擊式水輪機(jī)內(nèi)部流動非常復(fù)雜,這也是計算流體力學(xué)在沖擊式水輪機(jī)上的應(yīng)用遠(yuǎn)遠(yuǎn)落后于其它類型水輪機(jī)的原因。由于對沖擊式水輪機(jī)復(fù)雜的內(nèi)部流動機(jī)理缺乏深入了解,準(zhǔn)確地預(yù)測沖擊式水輪機(jī)的性能存在一定的困難,IEC標(biāo)準(zhǔn)中也缺乏這部分內(nèi)容。

過去,沖擊式水輪機(jī)的水力設(shè)計主要是依賴于經(jīng)驗(yàn)及模型試驗(yàn)。隨著現(xiàn)代計算機(jī)技術(shù)以及物理模型的不斷改進(jìn),對沖擊式水輪機(jī)內(nèi)部復(fù)雜流動進(jìn)行數(shù)值模擬也成為可能[4,5]。韓鳳琴、肖業(yè)祥等人在考慮射流膨脹率的前提下二維地對水膜在水斗中的流動進(jìn)行了數(shù)值模擬,并利用相對軌跡二次元地討論了自由射流和旋轉(zhuǎn)水斗背面的干涉問題[6-10]。

然而,沖擊式水輪機(jī)水斗是一個帶有缺口并且曲率急劇變化的復(fù)雜曲面,自由射流想通過水斗缺口較光滑地進(jìn)入旋轉(zhuǎn)水斗。但由于旋轉(zhuǎn)水斗的位置是隨時間時刻變化的,缺口對自由射流的分割也呈現(xiàn)出典型的三維瞬態(tài)現(xiàn)象[11]。文獻(xiàn)12中作者致力數(shù)值模擬前期理論研究與程序開發(fā)。本文將二維“動畫解析法”完善至三維空間,完成了自由射流和旋轉(zhuǎn)水斗間復(fù)雜的非定常流動的數(shù)值解析,并將三維可視化結(jié)果與試驗(yàn)結(jié)果進(jìn)行了比較。文中還討論了單位轉(zhuǎn)速和射流半徑膨脹率對射流與水斗流動相互作用的影響,進(jìn)而提出了動態(tài)水斗切入射流做瞬態(tài)復(fù)雜運(yùn)動時,射流及水斗背面發(fā)生干涉現(xiàn)象的機(jī)理。

1 水斗表面與射流截面的離散

沖擊式水輪機(jī)水斗曲面復(fù)雜,其邊界由以下幾部分組成：分水刃、內(nèi)邊緣、側(cè)邊緣、外邊緣、直壁和缺口。本文采用了三維邊界貼體網(wǎng)格對水斗表面進(jìn)行離散化,如圖1所示。由于水斗左右兩邊關(guān)于分水刃對稱,為了簡化研究,選取分水刃右半邊水斗作為研究對象。沿著邊界曲線,水斗的斜率在缺口和直壁的連接處突然發(fā)生變化,這一連接處的結(jié)點(diǎn)稱為缺口底部。缺口的存在希望使非定常射流能夠較光滑地流入旋轉(zhuǎn)水斗。當(dāng)通過打開噴針開度來增加流量時,射流半徑隨之增加,水流外邊緣將非常接近缺口底部,射流與缺口之間可能發(fā)生干涉,水斗缺口對先行水斗背面也可能產(chǎn)生干涉影響。

圖1 沖擊式水斗的復(fù)雜邊界

為了探究旋轉(zhuǎn)水斗和自由射流間復(fù)雜的非定常相互作用,將射流的截面在考慮射流膨脹率的情況下將其離散成有641個節(jié)點(diǎn)的邊界貼體網(wǎng)格,如圖2所示。每一個射流節(jié)點(diǎn)都是由jR和kθ兩個參數(shù)確定的。每個射流單元是由四個相鄰點(diǎn)來表征。

圖2 邊界貼體網(wǎng)格對射流的離散

圖2中的水斗缺口將射流分成入射射流分支和偏離射流分支。由于缺口的幾何形狀在空間上是三維的,那么缺口和射流的交點(diǎn)就應(yīng)該在jR-和kθ-兩方向上尋找。通過連接這些交點(diǎn),就可以得到被缺口切割的三維邊界。沿著缺口的射流節(jié)點(diǎn)組成入射射流分支的出口移動邊界,同時它也組成了偏離射流分支的入口移動邊界。缺口的幾何形狀對入射射流分支和偏離射流分支有很大的影響。

2 射流與水斗相互作用數(shù)值計算

2.1 非定常射流與旋轉(zhuǎn)水斗的時空離散

為了模擬直線前行的射流和旋轉(zhuǎn)水斗之間的非定常相互作用,利用動畫解析法來進(jìn)行流動的數(shù)值計算。將以角速度Ω旋轉(zhuǎn)的具有NB枚水斗的水輪機(jī)中相鄰兩枚水斗間NdivF等分割,各個微小角度旋轉(zhuǎn)之間的狀態(tài)稱為“一個時空步長”,如圖3所示。引入時空離散參數(shù)iF表示每個時空步長。當(dāng)旋轉(zhuǎn)水斗的分水刃尖端剛開始接觸自由射流最內(nèi)緣代表流線的瞬間定義為時空同期離散的起點(diǎn),即iF=0。這樣,非定常射流與旋轉(zhuǎn)水斗在時間與空間上進(jìn)行了離散。

圖3 射流與水斗的時空同期離散

用動畫解析法離散的空間步長ΔθF為

ΔθF=2π/(NB·NdivF)

(1)

離散的時間步長ΔtF為

ΔtF=2π/(Ω·NB·NdivF)

(2)

在用動畫解析法模擬射流與水斗的關(guān)系時,定義了3個水斗,即：研究對象為當(dāng)前水斗、與當(dāng)前水斗相鄰并先接收到射流的為先行水斗、與當(dāng)前水斗相鄰較當(dāng)前水斗后接收到射流的為后續(xù)水斗,如圖4所示。對于研究射流和水斗的關(guān)系時,從當(dāng)前水斗接觸到射流開始。當(dāng)前水斗剛接觸第一滴射流時,后續(xù)水斗還處在射流區(qū)之外,如圖4a所示。隨著射流逐漸的進(jìn)入當(dāng)前水斗,當(dāng)前水斗移動某一位置時,后續(xù)水斗剛好接觸到第一點(diǎn)射流,如圖4b所示。

圖4 先行水斗、當(dāng)前水斗和后續(xù)水斗

2.2 作用在當(dāng)前水斗上的射流節(jié)點(diǎn)

選擇有18個水斗,幾何比轉(zhuǎn)速為B/Dref=0.35(B為水斗內(nèi)表面寬度,Dref為轉(zhuǎn)輪參考直徑)的沖擊式水輪機(jī)轉(zhuǎn)輪作為非定常流動分析的數(shù)值模型。在最優(yōu)單位速度nDH=40 r/min的情況下,分水刃尖端與射流速度比Ust/C0=0.628,無因次噴針行程Sn/Rt=1.4(Sn為噴針行程,Rt為噴嘴半徑)。設(shè)定兩相鄰水斗間離散的動畫單元數(shù)為NdivF=40,射流三維沿射流半徑j(luò)R方向的離散數(shù)NdivR0=40,沿圓周方向kθ的離散數(shù)NdivR0=20。射流半徑的膨脹率設(shè)定為kRj=0.1%。iF為每一個時空步長。

通過曲面微分幾何及在微曲面上流體節(jié)點(diǎn)的投影可以得到射流節(jié)點(diǎn)在旋轉(zhuǎn)水斗上的真實(shí)三維位置,射流與水斗交點(diǎn)的計算規(guī)則如圖5所示[13]。

圖5 射流與水斗交點(diǎn)的計算規(guī)則

通過改進(jìn)的動畫解析法,對不同的動畫單元,作用在水斗表面的三維射流節(jié)點(diǎn)如圖6所示。

圖6 不同動畫單元射流在當(dāng)前水斗上的三維落點(diǎn)

圖6中列出了每隔20個動畫單元入射射流分支在當(dāng)前水斗的入射情況,水斗上的藍(lán)色網(wǎng)格即表示自由射流在當(dāng)前水斗上的入射面。隨著動畫單元數(shù)的增大,進(jìn)入水斗中的射流逐漸增多,直至整個射流進(jìn)入當(dāng)前水斗。由圖中可以看到,在iF到達(dá)40左右時,整個射流已經(jīng)進(jìn)入水斗。從iF=40到80,都是整個射流射入當(dāng)前水斗。需要注意的是,當(dāng)iF推進(jìn)到40以后,后續(xù)水斗將開始接觸并逐漸攔截射流,由iF=90左右,不再是整個射流進(jìn)入當(dāng)前水斗。隨著動畫單元數(shù)的繼續(xù)增大,當(dāng)前水斗中接到的射流越來越少。到iF=127以后,當(dāng)前水斗上已沒有射流射入其中。

圖7 不同動畫單元射流后緣面的三維示意

2.3 自由射流后緣面的數(shù)值計算

旋轉(zhuǎn)水斗缺口將自由射流非定常地分為入流和偏流兩股分支。偏流分支的后緣面由缺口在每個動畫單元分割自由射流的內(nèi)移動邊界的疊加得來,所以水斗缺口形狀對射流后緣面起著決定性的作用。

當(dāng)iF比較大的時候,后緣面的精確形狀對射流進(jìn)入水斗有很大的影響。為此我們利用動畫解析法得到了在不同動畫單元射流偏流分支的三維后緣面形狀,如圖7所示。

2.4 射流和水斗的相互作用

射流和水斗之間存在復(fù)雜時空非定常相互作用。文獻(xiàn)8利用“動畫單元解析法”二維地表示了自由射流被當(dāng)前水斗分割為入射射流分支和偏離射流分支的情況,如圖8所示。但是二維視圖還不足以準(zhǔn)確分析射流與水斗間的復(fù)雜非定常相互作用程度。

圖8 水斗對射流的二維非定常分割

由本文2.2節(jié)數(shù)值計算得到的射入面及2.3節(jié)的后緣面,得到三維空間不同動畫單元下射入面(圖中藍(lán)色面)和后緣面(圖中紅色面)之間的間隙,如圖9所示。當(dāng)iF=10時,上游射流的小部分進(jìn)入當(dāng)前水斗(作用在當(dāng)前水斗的內(nèi)表面),而下游射流(作用在當(dāng)前水斗的背面)的后緣面則流出當(dāng)前水斗射向先行水斗。在射流入射面和后緣面之間存在一點(diǎn)很小的空隙。當(dāng)iF=20時,更多的射流進(jìn)入當(dāng)前水斗,而后緣面也以更快的速度流出。當(dāng)iF=30時,幾乎所有的射流進(jìn)入當(dāng)前水斗,此時射入面和后緣面之間的空間增大。當(dāng)iF=40時,所有的射流進(jìn)入當(dāng)前水斗,而被當(dāng)前水斗缺口分割的偏離射流分支的完全后緣面射入先行水斗,此時射入面和后緣面之間有較寬的間隙。入射面和后緣面之間的間隙能直觀表示射流與水斗背面的干涉情況,如圖中iF=10及iF=20所示,間隙越小,射流與水斗后緣面越容易發(fā)生干涉。本沖擊式水輪機(jī)的設(shè)計良好,干涉情況不嚴(yán)重。如果設(shè)計不佳時,射流后緣面和入射面的位置更加接近甚至可能互換,表示干涉現(xiàn)象趨于嚴(yán)重。且三維空隙形狀比二維真實(shí),可以更加準(zhǔn)確地預(yù)測特性。

圖9 在不同動畫單元自由射流的三維非定常分割

2.5 內(nèi)部流態(tài)觀察試驗(yàn)與數(shù)值模擬比較

本研究在開發(fā)沖擊式水輪機(jī)水斗表面非定常自由水膜流動數(shù)值解析方法的同時,還與合作單位委托華中科技大學(xué)進(jìn)行了沖擊式水輪機(jī)的模型試驗(yàn),試驗(yàn)觀測到的水斗內(nèi)部水膜流態(tài)分布如圖10所示。圖10a表示了iF=40時射流全部射入當(dāng)前水斗,與圖10b的數(shù)值仿真結(jié)果有較好的吻合度。比起試驗(yàn)需要經(jīng)費(fèi)、技術(shù)等條件,數(shù)值仿真顯現(xiàn)令人滿意優(yōu)勢。如：對于當(dāng)前水斗來說,由于有先行水斗和后續(xù)水斗的干擾,真實(shí)的部分水膜會被擋住而看不見,數(shù)值仿真可視化可以幫助解決這個問題,一般看不見的部分可以由數(shù)值仿真可視化結(jié)果來表現(xiàn),這對于干涉現(xiàn)象描述起到非常好的作用。該自編軟件可快捷在PC機(jī)上運(yùn)行。

圖10 iF=40水斗內(nèi)的水膜流態(tài)

3 射流后緣面和水斗背面的干涉

如果入射射流分支的入射面和偏離射流分支的后緣面之間的空隙太狹窄,當(dāng)前水斗背面有可能滲入偏離射流而發(fā)生干涉,將導(dǎo)致水輪機(jī)的水力效率下降。比轉(zhuǎn)速nDH以及射流半徑膨脹率kRj都會影響射流后緣面與水斗背面的干涉情況,本節(jié)討論nDH及kRj對水力性能的影響程度。

3.1 自由射流的膨脹

在射流流動的過程中,射流外表面與大氣的摩擦產(chǎn)生切向應(yīng)力,從而導(dǎo)致射流外表面流速減小并發(fā)生動能損失。射流流速和動能持續(xù)減小導(dǎo)致射流直徑逐漸擴(kuò)大產(chǎn)生射流膨脹。

為了便于研究自由射流膨脹率對水斗內(nèi)部流動能量轉(zhuǎn)化的影響,引入了自由射流半徑膨脹率kRj,其定義如下

kRj=(Rj-R0)/(Xj-X0)×100%

(3)

式中,R0、X0是射流在縮流位置的射流半徑及縮流位置在靜止坐標(biāo)中沿X軸的位置；Rj、Xj是射流在膨脹以后的任意斷面射流半徑及該斷面在靜止坐標(biāo)系下的位置。

圖11 自由射流的膨脹

3.2 單位轉(zhuǎn)速和射流膨脹率對干涉的影響

單位轉(zhuǎn)速nDH(=nDref/H1/2)表現(xiàn)了水斗圓周速度Uref(～nDref)和射流速度Cj(～H1/2)的比率。所以兩射流表面間空隙的形狀直接受單位轉(zhuǎn)速nDH的影響。在保持轉(zhuǎn)速n不變的情況下,降低水頭H,射流速度Cj隨之變小,此時nDH增大,間隙隨之變窄,射流與水斗背面發(fā)生干涉的可能性加劇,如表1所示。除單位轉(zhuǎn)速對干涉有影響,射流膨脹率也會產(chǎn)生影響。

表1中清晰地展示了兩射流分支面上的數(shù)值可視化效果。當(dāng)nDHopt=40 rpm,kRj=0.1%時,即使當(dāng)iF=20時空隙也足以避免這種干涉；但當(dāng)nDH=45 rpm,kRj=0.3%時,即使當(dāng)iF=40時空隙也不足以避免這種干涉。

3.3 水輪機(jī)效率損失

本文3.2節(jié)中數(shù)值預(yù)測了水斗背面與射流間的非定常干涉現(xiàn)象,同時也計算得到了在不同工況和射流半徑膨脹率條件下射流干涉引起的水輪機(jī)效率損失,如表2所示。

表1 空隙的三維表現(xiàn)

表2 不同射流半徑膨脹率和比轉(zhuǎn)速時的效率損失

由表2中的結(jié)果可見nDH的增大直接導(dǎo)致了效率的下降。如在相同的射流半徑膨脹率kRj=0.3%下,nDH/nDHopt=1.25時的效率下降大約為7.9%,而nDH/nDHopt=1.5時的效率下降大約為39.1%,單位轉(zhuǎn)速增大1.2倍時效率下降了5倍。由此可見,單位轉(zhuǎn)速及射流半徑膨脹率的增加都會引起沖擊式水輪機(jī)水力性能的下降,單位轉(zhuǎn)速變化對于水斗的影響要遠(yuǎn)大于射流直徑膨脹變化帶來的影響。

4 結(jié) 論

本文數(shù)值模擬了旋轉(zhuǎn)水斗與自由射流之間的復(fù)雜非定常相互作用,得到如下結(jié)論：

(1)通過使用邊界貼體網(wǎng)格來離散水斗復(fù)雜的曲面和自由射流的過流截面,并對非定常射流和旋轉(zhuǎn)水斗在時間和空間上進(jìn)行離散,數(shù)值模擬了沖擊式水輪機(jī)的時空非定常三維流動。

(2)偏離射流分支后緣面和入射射流分支射入面有間隙形成。兩股射流分支之間的間隙呈現(xiàn)出復(fù)雜三維形狀,并且與時間相關(guān)。入射面和后緣面之間的間隙可以直觀表示射流與水斗背面的干涉情況,間隙越小越容易發(fā)生干涉。

(3)三維表現(xiàn)了射入當(dāng)前水斗的射流節(jié)點(diǎn)、入射射流的射入面和偏離射流后緣面的動態(tài)特性,得到了時空瞬態(tài)變化時射流及旋轉(zhuǎn)水斗非定常流動的相互作用。由干涉引起的水輪機(jī)性能變化正比于單位轉(zhuǎn)速、射流膨脹率的增加,單位轉(zhuǎn)速變化對于水斗干涉的影響要遠(yuǎn)大于射流直徑膨脹變化帶來的影響。

[1]肖業(yè)祥, 鄭愛玲, 韓鳳琴, 等. CFD法研究多噴嘴沖擊式水輪機(jī)的射流干涉[J].華南理工大學(xué)學(xué)報: 自然科學(xué)版, 2007, 35(3): 66- 70．

[2]韓鳳琴, 肖業(yè)祥, 久保田 喬. 小型水輪機(jī)渦殼三維粘性流動解析[J]. 農(nóng)業(yè)機(jī)械學(xué)報, 2004, 35(2)： 40- 43.

[3]肖業(yè)祥, 韓鳳琴, 久保田 喬. 非定常水斗表面自由水膜的流動機(jī)理[J]. 華南理工大學(xué)學(xué)報: 自然科學(xué)版, 2006, 34(4): 75- 79, 90．

[4]宋文武, 符杰, 周敏, 等. 多噴嘴沖擊式水輪機(jī)內(nèi)部流動及部件振形的研究[J]. 大電機(jī)技術(shù), 2008, 6: 41- 45．

[5]賀麗萍, 王茜蕓, 劉萬江, 等. 沖擊式水輪機(jī)內(nèi)部流態(tài)的數(shù)值模擬[C]∥第十九次中國水電設(shè)備學(xué)術(shù)討論會論文集, 大連: 黑龍江科學(xué)技術(shù)出版社： 2013, 11: 119- 124．

[6]沈娜, 韓鳳琴, 久保田喬. 基于微分幾何水斗內(nèi)移動網(wǎng)格的時空非定常流動[J]. 華南理工大學(xué)學(xué)報: 自然科學(xué)版, 2012, 40(2): 39- 44.

[7]韓鳳琴, 久保田 喬,劉潔. 數(shù)值模擬沖擊式水輪機(jī)內(nèi)部非定常流動 [J]. 華中理工大學(xué)學(xué)報, 2000, 28(11): 14- 16．

[8]周晶林, 韓鳳琴, 吳迪, 等. 兩相流射流的最小半徑數(shù)值研究[J]. 水電能源科學(xué), 2007, 25(5)： 110- 114．

[9]肖業(yè)祥, 沈娜, 韓鳳琴, 等. 水斗曲面三維非正交貼體網(wǎng)格數(shù)值生成[J]. 工程熱物理學(xué)報, 2008, 29(4): 595- 598.

[10]HAN Fengqin, SHEN Na, LI Longxiang, et al. Unsteady separation of Jet branch by rotating cutout of Pelton bucket[J] . Science China: Technological Sciences, 2011, 54(2): 302- 310.

[11]沈娜, 韓鳳琴. 沖擊式水輪機(jī)數(shù)值流動解析研究[J]. 水力發(fā)電, 2015, 41(12): 68- 71,75．

[12]沈娜. 靜止系射流與旋轉(zhuǎn)系移動水膜流間時空相互干涉[D]. 廣州： 華南理工大學(xué), 2011.

[13]韓鳳琴, 劉平安, 肖業(yè)祥, 等. 復(fù)雜急曲率變化水斗貼體坐標(biāo)的曲面擬合[J]. 水電能源科學(xué), 2005, 23(4): 68- 69．

[14]KUBOTA T. Observation of Jet Interference in 6-nozzle Pelton Turbine[J]. Journal of Hydraulic Research, 1989, 27(6): 753- 755．

(責(zé)任編輯 高 瑜)

Research of 3D Space-Time Unsteady Flow between Jet Flow and Rotating Buckets

SHEN Na, HAN Fengqin

(School of Electrical Engineering, Guangzhou College of South China University of Technology,Guangzhou 510800, Guangdong, China)

The unsteady flow relationship between injection jet and rotating buckets is studied in order to more accurately investigate the jet interference effects through improving the cartoon analysis method from 2D to 3D. The geometry of the impinging branch landed on relevant bucket surface and the deflected branch toward the preceding bucket are numerically analyzed, and the complex unsteady flow between free jet and rotating buckets are described with visualized 3D image. The interaction relationship between free jet and rotating buckets is obtained, and then the mechanism of interference phenomenon between free jet and rotating bucket with the transient complex movement when the dynamic bucket cutting into the free jet is proposed. The analysis results show that the performance change of turbine caused by the interference is proportional to the increase of unit speed and jet expansion, the interference caused by unit speed is much greater than that by free jet expansion．

unsteady flow; Pleton turbine; injection jet; rotating bucket; interference phenomenon

2016- 06- 02

國家自然科學(xué)基金資助項(xiàng)目(50879026)；浙江省重大科技專項(xiàng)(2008C11057)；華南理工大學(xué)廣州學(xué)院優(yōu)秀青年骨干教師項(xiàng)目(B1130007)

沈娜(1984—),女,湖北武漢人,講師,博士,主要從事新能源發(fā)電技術(shù)研究．

TK735

A

0559- 9342(2017)03- 0088- 06