朱多彪，白曉飛，沈 云

1.國(guó)家知識(shí)產(chǎn)權(quán)局專(zhuān)利局專(zhuān)利審查協(xié)作湖北中心，湖北武漢 430200 2.湖北宏源電力設(shè)計(jì)咨詢(xún)有限公司，湖北武漢 430000 3.中水淮河規(guī)劃設(shè)計(jì)研究有限公司，安徽合肥 230000

基于大渦模擬的貫流式水輪機(jī)尾水渦帶研究

朱多彪1，白曉飛2，沈 云3

1.國(guó)家知識(shí)產(chǎn)權(quán)局專(zhuān)利局專(zhuān)利審查協(xié)作湖北中心，湖北武漢 430200 2.湖北宏源電力設(shè)計(jì)咨詢(xún)有限公司，湖北武漢 430000 3.中水淮河規(guī)劃設(shè)計(jì)研究有限公司，安徽合肥 230000

基于大渦模擬方法并采用Smagorinsky-Lilly動(dòng)態(tài)亞格子模型，分別計(jì)算了燈泡貫流式水輪機(jī)的6個(gè)不同工況的非定常流場(chǎng)，得到了每個(gè)工況尾水管段尾水管渦帶信息。研究發(fā)現(xiàn)：尾水管中心區(qū)域的脈動(dòng)強(qiáng)度大于外圍。水頭對(duì)尾水渦帶形態(tài)影響顯著，隨著水頭的降低，渦帶形狀由雙螺旋變成單螺旋，直至平直狀。

大渦模擬；貫流式水輪機(jī)；壓力脈動(dòng)；渦帶

1963年，美國(guó)氣象學(xué)者 Smagorinsky[1]為研究天氣預(yù)報(bào)問(wèn)題首次提出了大渦模擬（Large Eddy Simulation）方法[2]。該方法認(rèn)為湍流運(yùn)動(dòng)是由不同尺度的渦組成，通過(guò)濾波計(jì)算，將湍流的瞬時(shí)運(yùn)動(dòng)分解成小尺度渦和大尺度渦兩部分。其中，大尺度渦擁有主導(dǎo)的湍流動(dòng)能，通過(guò)控制方程直接求解；小尺度（亞格子尺度）渦主要起耗散作用，在高雷諾數(shù)下趨于各向同性，受邊界條件影響較小，可用通用模型進(jìn)行模擬。因此，亞格子模型是大渦數(shù)值模擬控制方程封閉的前提，是大渦模擬方法的核心。Smagorinsky提出了第一個(gè)亞格子模型，Lilly[3]改進(jìn)并發(fā)展了Smagrinsky-Lilly動(dòng)態(tài)亞格子模型，其中Smagorinsky模型系數(shù)SC根據(jù)解析尺度提供的信息動(dòng)態(tài)計(jì)算，由當(dāng)時(shí)當(dāng)?shù)氐乃俣葓?chǎng)動(dòng)態(tài)處理，故而能動(dòng)態(tài)地反映湍流信息，且與實(shí)驗(yàn)結(jié)果符合良好。

目前，LES在湍流的數(shù)值計(jì)算中獲得了大量的關(guān)注，我國(guó)的研究者也取得了豐碩的成果。茅媛婷[4]將LES方法應(yīng)用于混流式水輪機(jī)的研究，取得了與實(shí)驗(yàn)吻合較好的結(jié)果，夏晨宇[5]將大渦模擬方法應(yīng)用于貫流式水輪機(jī)壓力脈動(dòng)的研究，然而大渦模擬方法用于貫流式水輪機(jī)尾水渦帶的研究并不多見(jiàn)。本文以燈泡貫流式水輪機(jī)為研究對(duì)象，重點(diǎn)針對(duì)尾水渦帶進(jìn)行了研究，為大渦模擬的應(yīng)用和貫流式水輪機(jī)尾水渦帶的研究做了積極探索。

1 幾何模型與計(jì)算方法

1.1 計(jì)算區(qū)域與計(jì)算工況

基于圓柱坐標(biāo)系建立三維葉片，進(jìn)而建立燈泡貫流式水輪機(jī)全流道模型，轉(zhuǎn)輪直徑D1=7.4m，輪轂直徑d=3.1m，葉片數(shù)Z=5，導(dǎo)葉數(shù)16。計(jì)算控制域即水輪機(jī)的全流道流場(chǎng)空間，在進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證的基礎(chǔ)上，采用非結(jié)構(gòu)三維網(wǎng)格技術(shù)劃分網(wǎng)格。水輪機(jī)發(fā)電機(jī)額定轉(zhuǎn)速n=83.33r/min，旋轉(zhuǎn)方向?yàn)橄蛳掠雾槙r(shí)針。本文按照水頭高低，選擇了運(yùn)行范圍內(nèi)從最高到最低水頭的6個(gè)不同工況（詳見(jiàn)表1），進(jìn)行了三維紊流非定常流場(chǎng)計(jì)算。

1.2 計(jì)算求解方法

1.2.1 濾波方法

大渦模擬第一步就是濾波，把流動(dòng)分成小尺度量和大尺度量。濾波運(yùn)算通過(guò)消除湍流中小尺度脈動(dòng)，濾掉比過(guò)濾網(wǎng)格小的漩渦，從而得到可解尺度運(yùn)動(dòng)的控制方程。過(guò)濾方法如下：

表1 計(jì)算工況

1.2.2 計(jì)算方法

為了節(jié)省計(jì)算時(shí)間，首先用RANS方法定常計(jì)算出一個(gè)初始流場(chǎng)，再改用LES方法計(jì)算非定常的流動(dòng)狀態(tài)，LES計(jì)算直到湍流變得統(tǒng)計(jì)地穩(wěn)定。

非定常階段采用大渦模擬計(jì)算方法，亞格子模型采用Smagorinsky-Lilly動(dòng)態(tài)模型，用以模擬小尺度渦的運(yùn)動(dòng)。由于水輪機(jī)的進(jìn)口壓力比較穩(wěn)定且易于獲得，因此采用壓力進(jìn)出口邊界條件，壓差為水頭，考慮吸出高度。固壁采用無(wú)滑移邊界，動(dòng)靜干涉選擇滑移網(wǎng)格模型?？紤]到轉(zhuǎn)輪頻率的關(guān)系，采用 0.002s 為時(shí)間步長(zhǎng)，即每個(gè)時(shí)間步轉(zhuǎn)輪轉(zhuǎn)過(guò)1°。離散化方法采用有界中心差分格式，壓力差值選擇PRESTO!格式，壓力速度耦合選擇PISO算法。

2 尾水管渦帶特性分析

通過(guò)設(shè)置等壓面，選擇合適的壓強(qiáng)，顯示尾水渦帶形態(tài)。

圖1 （a）顯示工況1尾水渦帶以泄水錐為起點(diǎn)，由互相對(duì)稱(chēng)的兩條螺旋線(xiàn)交錯(cuò)形成“麻花”狀。該工況下轉(zhuǎn)輪出流帶有環(huán)量，水流旋轉(zhuǎn)著進(jìn)入尾水管，為渦帶的形成創(chuàng)造了條件，但渦帶處于尾水管中心區(qū)域，并沒(méi)有發(fā)生偏心。此時(shí)旋轉(zhuǎn)的水流在尾水管內(nèi)受離心力的作用，而附著在尾水管壁面，導(dǎo)致流速放緩、湍流粘滯性上升，形成“貼壁效應(yīng)”[6]，故而尾水管內(nèi)渦帶外圍區(qū)域壓力脈動(dòng)不強(qiáng)，尾水管壁面也不會(huì)受到強(qiáng)烈的水力振動(dòng)。工況2也是雙螺旋交織形成的渦帶，與工況1不同的是，兩條螺旋的對(duì)稱(chēng)性不明顯，大小相差較大形成主次之分。大螺旋與工況1的相似，較粗壯，隨著流場(chǎng)向下游發(fā)展而衰減變細(xì)直至消失。小螺旋與大螺旋交織在一起，相對(duì)瘦弱且螺旋狀不明顯，并且比大螺旋較早消失。圖1（c）顯示工況3的尾水渦帶形態(tài)與工況2相似，也是由兩條小渦帶交織形成。不同的是，其中由一條呈螺旋狀而另一條較為平直，并且兩條渦帶大小相當(dāng)，交織更加緊密并有相融為一體的趨勢(shì)。由圖1（d）可以看出，工況4的尾水渦帶僅由一條螺旋帶組成。渦帶起源于泄水錐，較為粗壯而隨流場(chǎng)向下游衰減，但衰減緩慢拖延綿長(zhǎng)。工況5的水頭較低，因此全流道整體壓差較小，尾水管并沒(méi)有出現(xiàn)螺旋狀的渦帶，而僅因?yàn)樾顾F的導(dǎo)流作用而形成了一條平直的渦帶。工況6水頭也較低，因此與工況5相似，尾水渦帶仍為一條直筒狀的平直渦帶，且較為短小，衰減較快，很快便消失。

由各工況的渦帶形態(tài)圖可以看出，水頭與尾水渦帶形態(tài)密切相關(guān)，在高水頭的工況下會(huì)形成雙螺旋“麻花狀”渦帶，隨著水頭的減小，逐漸演變成為單螺旋渦帶，當(dāng)水頭進(jìn)一步減小時(shí)渦帶則呈平直狀。然而在工況變化過(guò)程中葉片開(kāi)度、水頭、流量、出力等因素都有變化，僅有水頭和出力的變化較有規(guī)律，水頭與出力的變化如何影響渦帶形態(tài)有待進(jìn)一步深入研究。

綜合來(lái)看，由于燈泡貫流式機(jī)組的進(jìn)水流道和尾水管為直線(xiàn)型，流道布置采用水力軸心全對(duì)稱(chēng)型式，所以各工況下的渦帶形態(tài)都表現(xiàn)出圓周對(duì)稱(chēng)的形態(tài)，且處于尾水管軸線(xiàn)中心區(qū)域，降低了渦帶偏心并撞擊尾水管壁面的可能性。與采用缺乏對(duì)稱(chēng)性的蝸殼和彎肘型尾水管布置相比，燈泡機(jī)組提高了運(yùn)行的安全性。當(dāng)高水頭時(shí)，轉(zhuǎn)輪為環(huán)量出口，水流旋轉(zhuǎn)著進(jìn)入尾水管，旋轉(zhuǎn)的環(huán)量產(chǎn)生了阻止水流脫離邊壁的離心力，有利于減少擴(kuò)散管中的能量損失，但會(huì)產(chǎn)生螺旋狀渦帶，渦帶中心區(qū)的脈動(dòng)強(qiáng)度大于低水頭工況的平直渦帶。因此，水輪機(jī)運(yùn)行中水頭高低的選擇值得慎重權(quán)衡。

3 結(jié)論

高水頭時(shí)，尾水管內(nèi)有兩條螺旋線(xiàn)交錯(cuò)形成“麻花狀”渦帶，隨著水頭降低，渦帶逐漸演變成為單螺旋線(xiàn)；低水頭時(shí)工況，渦帶則成平直狀。水頭決定了尾水渦帶的形狀，水頭越高螺旋越明顯。尾水管中心區(qū)域的脈動(dòng)強(qiáng)度沿軸線(xiàn)向下游衰減，至尾水管中段流場(chǎng)開(kāi)始趨于平穩(wěn)。尾水管渦帶一直處于沿軸線(xiàn)小半徑區(qū)域，并沒(méi)有發(fā)生偏心，也不會(huì)撞擊尾水管壁面，因此，不會(huì)對(duì)廠房造成較大振動(dòng)。

[1]Smagorinsky J. General circulation experiments with the Primitive equations[J]. Mon. Weath. Rev.，1963，91（3）：99–164

[2]張兆順，崔桂香，許春曉.湍流大渦數(shù)值模擬的理論和應(yīng)用[M].北京：清華大學(xué)出版，2008.

[3]D. K. Lilly. A Proposed Modification of the Germano Subgrid-Scale Closure Model. Physics of Fluids， 4：633–635， 1992.

[4]茅媛婷，鄭源，周大慶，等.混流式水輪機(jī)尾水管壓力脈動(dòng)特性分析[C].第四屆全國(guó)水力機(jī)械及其系統(tǒng)學(xué)術(shù)會(huì)議論文集，2011：373-376.

[5]夏晨宇，李龍，朱多彪，等.基于大渦模擬的貫流式水輪機(jī)壓力脈動(dòng)分析[J].中國(guó)農(nóng)村水利水電，2015，（2）：143-147.

[6]朱多彪，李龍，沈云.基于升力法的貫流式水輪機(jī)葉片設(shè)計(jì)及可行性分析[J].水電能源科學(xué)，2013，07：158-161.

TK7

A

1674-6708（2015）150-0147-02