，

(昆明理工大學(xué)土木工程學(xué)院，云南 昆明 650500)

貫流式水輪機三維非定常流DES數(shù)值模擬

歐陽振宇，張立翔

(昆明理工大學(xué)土木工程學(xué)院，云南 昆明 650500)

0 引言

貫流式水輪機具有幾何結(jié)構(gòu)特殊、運行工況多變的特點，其內(nèi)部流場是非常復(fù)雜的黏性、三維非定常湍流流動，而且內(nèi)部流動的工質(zhì)可能處于兩相甚至多相狀態(tài)。因此，為了進(jìn)一步改善貫流式水輪機水力性能,縮短水輪機設(shè)計周期,降低設(shè)計的風(fēng)險，有必要對貫流式水輪機內(nèi)部湍流流動進(jìn)行研究，以便深刻認(rèn)識湍流流動的發(fā)生、發(fā)展規(guī)律，以及對后續(xù)流動的影響。

目前，國內(nèi)外研究者對貫流式水輪機水力性能進(jìn)行研究主要集中在2個方面：一是經(jīng)驗設(shè)計與模型試驗相結(jié)合，通過反復(fù)試驗，最終得到較優(yōu)的設(shè)計，但是試驗具有周期長、費用高，且易受外界干擾等不足；二是數(shù)值模擬研究，隨著計算流體力學(xué)理論與計算機軟、硬件的快速發(fā)展，利用計算流體力學(xué)技術(shù)研究葉輪機械內(nèi)部流場和優(yōu)化設(shè)計已成為趨勢[1-5]。數(shù)值模擬研究方法突出的優(yōu)點在于周期短、成本低。

針對某水電站貫流式模型水輪機，以標(biāo)準(zhǔn)的k-ε模型計算的三維定常模擬的收斂結(jié)果作為初始場，采用分離渦模型和滑移網(wǎng)格方法進(jìn)行內(nèi)部湍流場非定常的模擬，獲得了水輪機湍流場的瞬態(tài)流動特征。結(jié)果表明,DES湍流模型能夠準(zhǔn)確模擬貫流式水輪機不同瞬時內(nèi)部流場，研究成果為水力誘發(fā)的水力發(fā)電機組振動提供理論依據(jù)。

1 數(shù)值模擬方法

1.1 計算對象及網(wǎng)格設(shè)計

以某總裝機容量為3×40MW的水電站貫流式水輪機機組的單機建立模型，電站為徑流式開發(fā)，以發(fā)電為主，混凝土閘壩最大壩高 49 m，水庫正常蓄水位 325 m，相應(yīng)庫容 1.217×108 m3。水輪機燈泡體部分如圖1所示，其上有固定導(dǎo)葉、16個活動導(dǎo)葉以及轉(zhuǎn)輪5個葉片。系統(tǒng)選取的計算區(qū)域為在正常工作開度下從進(jìn)水口到尾水管出口的整個流道，模型順河向為Z軸方向，指向下游為正；橫河向為Y軸方向，指向右岸為正；X軸為垂直方向，向上為正，模型及坐標(biāo)如圖2所示。全流道由3部分連接而成，分別為進(jìn)水口部分、轉(zhuǎn)輪室部分和尾水管部分，計算網(wǎng)格采用適應(yīng)性較強的非結(jié)構(gòu)網(wǎng)格。由于水輪機內(nèi)部湍流的發(fā)生主要從活動導(dǎo)葉出口至轉(zhuǎn)輪室出口和尾水管前部，故在該區(qū)間對網(wǎng)格適當(dāng)加密，總計網(wǎng)格節(jié)點數(shù)共約165萬個，網(wǎng)格單元數(shù)共約618萬個。

圖1 燈泡體

圖2 流道網(wǎng)格

1.2 湍流模型和離散方法

(1)

d為(S-A)模型中計算單元離壁面的距離;CDES為常數(shù)，取 0.65;Δ為所有計算單元中心與最臨近計算單元中心距離的最大值，即

Δ=max(Δx,Δy,Δz)

(2)

DES的核心思想是在流動平緩區(qū)域采用RANS求解，在具有大尺度回流和分離等渦脫落區(qū)域采用LES求解。在每個計算單元中，首先應(yīng)用RANS求解湍流流動，當(dāng)RANS計算得到的湍流尺度大于當(dāng)?shù)鼐W(wǎng)格尺度時，進(jìn)一步轉(zhuǎn)換為LES的亞格子模型求解。

系統(tǒng)在DES Spalart-Allmaras模型求解中，RANS湍流應(yīng)力的封閉采用單方程Spalart-Allmaras湍流模型，空間參數(shù)離散采用二階迎風(fēng)格式。LES求解采用Smagorinsky亞格子模型，空間參數(shù)離散采用具有二階精度的中心差分格式。采用有限體積法和非交錯網(wǎng)格對瞬態(tài)控制方程進(jìn)行離散，時間項采用二階全隱式格式，源項和擴散項采用二階中心格式，對流項采用二階迎風(fēng)格式。壓力和速度的耦合求解采用SIMPLEC算法。整個計算過程由CFD軟件ANSYS FLUENT完成。

1.3 初始條件和邊界條件

與單部件的流動模擬相比，全流道所有部件整體模擬耦合計算的邊界條件更容易給定，在動靜部件間不會產(chǎn)生不準(zhǔn)確的邊界條件，計算時只需指定進(jìn)口和出口的邊界條件即可，計算結(jié)果與實際情況更接近。為了研究葉片旋轉(zhuǎn)對水輪機內(nèi)部流場結(jié)構(gòu)的影響，在模型水輪機正常工作轉(zhuǎn)速下進(jìn)行計算。給定速度進(jìn)口和尾水管自由出流邊界條件，在壁面采用無滑移邊界條件，近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)。初始計算參數(shù)設(shè)置的情況如表1所示。

表1 初始計算參數(shù)

采用標(biāo)準(zhǔn)k-ε模型的定常計算結(jié)果作為非定常計算的初始場。分離渦模擬計算中時間步長為0.0001s，計算收斂后，時間步向前推進(jìn)，同時轉(zhuǎn)輪網(wǎng)格相應(yīng)轉(zhuǎn)動到新的位置，開始進(jìn)行新時間步上的計算，共計算0.016s，轉(zhuǎn)輪轉(zhuǎn)動12.96°。

1.4 動靜區(qū)域計算方法

定常湍流計算采用多參考系(MRF)模型，非定常湍流計算采用滑移網(wǎng)格模型。當(dāng)轉(zhuǎn)動部件和非轉(zhuǎn)動部件之間的交互作用很強時應(yīng)使用滑移網(wǎng)格模型，滑移網(wǎng)格模型是模擬動靜干擾非定常流場最精確的方法。MRF模型計算出流場的結(jié)果可作為瞬態(tài)滑移網(wǎng)格模型計算的初始條件。計算中，將水輪機動靜干擾較為強烈的轉(zhuǎn)輪區(qū)域設(shè)置為滑移網(wǎng)格區(qū)域。

2 結(jié)果及分析

2.1 水輪機內(nèi)部壓力場及湍動粘度分布

在水輪機正常工作的工況下，對模型水輪機內(nèi)的非定常流場進(jìn)行模擬，記錄下壓力場和湍動粘度隨時間的變化。圖3為壓力場在不同時刻的壓強分布。從圖3可以看出，隨著轉(zhuǎn)輪的轉(zhuǎn)動，水輪機內(nèi)部壓力場在該時間段內(nèi)不斷增大，在燈泡體內(nèi)壓強增大較為明顯，隨后壓力的增大沿轉(zhuǎn)輪至尾水管發(fā)展。

圖4為工況下不同時刻水輪機在Y軸截面湍動粘度分布。圖4中，隨著時間的推移，水輪機整體湍動能有先增大，后減小的趨勢。在t=0s時刻，水輪機進(jìn)水口湍動粘度較大，在其他區(qū)域的湍流強度不大，表明在水輪機開啟不久，進(jìn)口處湍流運動劇烈；在t=0.008 s時刻，進(jìn)口湍動粘度減小，沿水流方向湍動粘度增大，特別在葉輪區(qū)域及尾水管出口處湍動粘度急劇增大，表明進(jìn)口的湍流運動向前推進(jìn)，葉輪區(qū)域及尾水管出口湍流運動劇烈且復(fù)雜；在t=0.016s時刻，水輪機內(nèi)部流場整體湍動粘度減小，進(jìn)口和出口湍動粘度減小明顯，葉輪及葉輪出口湍動粘度有減小的趨勢，表明水輪機運行趨于穩(wěn)定。

圖3 水輪機Y軸截面不同瞬時壓強分布

圖4 水輪機Y軸截面不同瞬時湍動粘度分布

2.2 水輪機葉輪區(qū)域的渦量分布變化及分析

由圖5和圖6可知，隨著時間的推移，水輪機轉(zhuǎn)輪室進(jìn)出口處的渦量都在增加，在進(jìn)口處渦量增加幅度相對較少，渦量由靠近內(nèi)環(huán)壁面處開始向整個平面均勻擴散，表明水輪機在啟動過程中葉輪室進(jìn)口處渦街會經(jīng)歷一個發(fā)展壯大、失穩(wěn)破碎和均勻分布的演化過程；在轉(zhuǎn)輪室出口處，渦量增加明顯，渦量不斷向壁面方向擴散、聚集和增長，渦量分布變化明顯，表明轉(zhuǎn)輪的轉(zhuǎn)動對尾水管渦量分布及數(shù)量的影響明顯，旋轉(zhuǎn)能夠引起湍流流動結(jié)構(gòu)的劇烈變化。另一方面，比較同一時刻轉(zhuǎn)輪室進(jìn)出口的渦量變化可知，進(jìn)口處的渦量穩(wěn)定趨勢較出口處明顯，表明流體由進(jìn)口向出口的流經(jīng)過程，流向渦對在不同瞬時經(jīng)過相互碰撞、破碎、合并和脫落等過程充分發(fā)展之后，對旋轉(zhuǎn)的轉(zhuǎn)輪產(chǎn)生強動靜干擾，產(chǎn)生的渦帶進(jìn)一步延伸到轉(zhuǎn)輪葉道中可能誘發(fā)葉片的共振，造成破壞，影響機組穩(wěn)定運行。這種復(fù)雜的非定常特性，對理解機組運行過程中水輪機導(dǎo)水機構(gòu)內(nèi)的實際流動狀態(tài)具有重要意義。

圖5 水輪機轉(zhuǎn)輪室進(jìn)口不同瞬時渦量分布

3 結(jié)束語

通過建立包含貫流式水輪機的全流道模型，應(yīng)用FLUENT軟件，采用有限體積法和DES湍流模型對其內(nèi)部流場進(jìn)行三維非定常數(shù)值模擬，并分析了內(nèi)部流場壓力和湍動能粘度，轉(zhuǎn)輪室進(jìn)出口的渦量變化及影響。渦結(jié)構(gòu)在水輪機的激流振動中起著決定性作用，利用分離渦模型可以準(zhǔn)確地再現(xiàn)三維動態(tài)渦結(jié)構(gòu)，為研究水輪機的渦激振動及其對機組出力擺動的影響提供了理論依據(jù)，對水電站的流固耦合研究有參考意義。數(shù)值模擬結(jié)果表明,DES結(jié)合了RANS方法及LES方法，能高效和準(zhǔn)確地模擬高旋運動的復(fù)雜渦旋流場，表現(xiàn)出LES方法準(zhǔn)確模擬大尺度漩渦分離的特點，顯示出了優(yōu)越性，具有廣闊的應(yīng)用前景。由于DES方法結(jié)合了RANS與LES各自方法的優(yōu)點,而且利用對方的長處有效地彌補了自身的不足，在當(dāng)前有限的計算條件下，成為可以準(zhǔn)確而高效地模擬復(fù)雜流道三維非定常湍流流動的一種可行的辦法。

[1] 宋文武，胡 杰.高水頭貫流式水輪機的全流場解析[J].西華大學(xué)學(xué)報(自然科學(xué)版),2010,29(2):169-173.

[2] 楊栗晶.雙渦輪燈泡貫流式水輪機三維定常流動解析[J].水利水電技術(shù),2012,43(7):123-126.

[3] 姚 激，張立翔.基于滑移網(wǎng)格的小型垂直軸風(fēng)力機氣動性能的數(shù)值模擬[J].機械與電子，2013，(5)：12-15.

[4] 張 偉,陳紅勛.離心葉輪設(shè)計及非設(shè)計工況湍流 RANS模擬的驗證和確認(rèn)[J].水動力學(xué)研究與進(jìn)展(A輯)，2011，26(1):73-84.

[5] 周大慶,吳玉林,劉樹紅.軸流式水輪機模型飛逸過程三維湍流數(shù)值模擬[J].水利學(xué)報,2010,41(2):233-238.

[6] Lin Peifeng,Wu Dichong,Zhu Zefei.Nanoparticle transportation and brownian diffusion in planar jet flow via large eddy simulation[J].Journal of Fluid Dynamics,2013,(4):354-358.

[7] 康小偉,韓 維,王允良.充氣式翼型繞流大渦模擬分析[J].海軍航空工程學(xué)院學(xué)報,2013,(2):137-140.

[8] 周 磊,解茂昭,羅開紅.大渦模擬在內(nèi)燃機中應(yīng)用的研究進(jìn)展[J].力學(xué)學(xué)報,2013,(4):467-482.

[9] Fauconnier D,Langhe C D,Dick E.A family of dynamic finite difference schemes for large-eddy simulation[J].Journal of Computational Physics，2009，228(6)：1830-1861.

[10] 黃劍峰,張立翔,王文全,等.混流式水輪機三維非定常流分離渦模型的精細(xì)模擬[J].中國電機工程學(xué)報，2011,31(26):83-89.

Simulation of Three-dimensional Unsteady Flows in a Tubular Hydro-Turbine on Detached Eddy Simulation

OUYANGZhen-yu,ZHANGLi-xiang

(School of Civil Engineering,Kunming University of Science and Technology,Kunming 650500,China)

以某水電站為研究對象，建立包含貫流式水輪機的全流道數(shù)值模型，基于不可壓縮瞬態(tài)N-S方程，進(jìn)行了三維非定常流研究。將全流道定常計算的結(jié)果作為初始流場，應(yīng)用分離渦模型(DES)和滑移網(wǎng)格技術(shù)，進(jìn)行了水輪機動靜干擾的非定常湍流數(shù)值模擬，得到了全流場的瞬時壓力分布和湍動粘度分布，轉(zhuǎn)輪室進(jìn)出口處的瞬時渦量分布，捕捉到水輪機轉(zhuǎn)輪室進(jìn)出口處復(fù)雜的三維動態(tài)渦結(jié)構(gòu)。研究工作對于探討水輪機渦激振動形成的機理具有一定的參考意義。

貫流式水輪機；計算流體動力學(xué)；三維非定常湍流；分離渦模擬

Based on the transient incompressible NS equations,taking a hydropower station as the research object,the establishment of a numerical model that contains the full passage tubular turbine,conducted a three-dimensional unsteady flow research.The whole calculation runner steady flow field as a result of the initial application of the DES and slip grid technology,were often unsteady numerical simulation of turbulent movement of turbine interference,has been instantaneous pressure distribution and turbulent flow over the kinematic viscosity distribution,import and export of instantaneous vorticity distribution at the runner chamber,turbine runner to capture the complex three-dimensional dynamic vortex chamber inlet and outlet structures.Studies to explore the mechanism for the formation of vortex-induced vibration of turbine has a certain reference value.

tubular hydro-turbine;computational fluid dynamics(CFD);3D unsteady turbulences;DES

2014-01-06

國家自然科學(xué)基金資助項目(51279071)

TV734.1

A

1001-2257(2014)06-0048-04

歐陽振宇(1990-)，男，湖南常德人，碩士研究生，研究方向為水輪機流固耦合；張立翔(1959-)，男，云南石屏人，教授，博士研究生導(dǎo)師，研究方向為流固耦合。