蘇文濤，鄭智穎，李小斌，李鳳臣，蘭朝鳳，2

(1.哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院，150001哈爾濱;2.哈爾濱理工大學(xué)電氣與電子工程學(xué)院，150080哈爾濱)

混流式水輪機(jī)在偏工況運(yùn)行時，通常會出現(xiàn)內(nèi)部流動不穩(wěn)定，機(jī)組振動和噪聲增大現(xiàn)象，除了固有的機(jī)械因素外，主要是由于流致振動造成，其中包括30%～60%額定負(fù)荷時尾水管渦帶誘發(fā)大幅壓力脈動，50%額定負(fù)荷時輪葉片進(jìn)口Kármán渦脫流，這樣不僅降低水輪機(jī)運(yùn)行效率，嚴(yán)重時還會造成轉(zhuǎn)輪相鄰葉片中間區(qū)域出現(xiàn)葉道渦，并伴隨空化現(xiàn)象［1］.偏工況運(yùn)行中，在轉(zhuǎn)輪區(qū)主要出現(xiàn)泡狀或片狀的翼型空化，空泡潰滅即為轉(zhuǎn)輪破壞的主要因素.當(dāng)轉(zhuǎn)輪中至少有3個葉片表面出現(xiàn)氣泡時，表明已經(jīng)發(fā)生空化，轉(zhuǎn)輪中空化程度主要由空化數(shù)σ來描述.實(shí)驗(yàn)證明，水輪機(jī)壓力脈動的強(qiáng)度隨空化數(shù)的減小而迅速提高，而運(yùn)行效率則隨之下降［2］.針對混流式水輪機(jī)中內(nèi)部的不穩(wěn)定流動和空化現(xiàn)象，眾多學(xué)者展開了實(shí)驗(yàn)研究和數(shù)值模擬研究.實(shí)驗(yàn)研究主要針對水輪機(jī)流道的各個部分，采用粒子圖像測速儀(PIV)來進(jìn)行流場的測量和分析.Iliescu等［3］采用 PIV技術(shù)對混流式水輪機(jī)尾水管內(nèi)汽、液兩相流動進(jìn)行測量，研究空化尾水管渦帶結(jié)構(gòu).Palafox 等［4］和 Uzol等［5］分別采用PIV實(shí)驗(yàn)研究了軸流式水輪機(jī)葉頂間隙內(nèi)的湍流流動，該軸流式水輪機(jī)比轉(zhuǎn)速較低，且葉頂泄漏較為復(fù)雜，在下游將形成二次流.Katz等［6-7］基于對軸流式水輪機(jī)的PIV測量結(jié)果，引入湍動能組分平衡的觀點(diǎn)，分析了葉頂間隙流動及其尾跡的流動特征.針對偏工況運(yùn)行時的混流式水輪機(jī)尾水管內(nèi)流動，李丹等［8］、王軍等［9］和陳次昌等［10］對其二維流場和渦帶運(yùn)動特征展開了研究，分析渦帶的周期性運(yùn)動特征.實(shí)驗(yàn)研究已經(jīng)揭示了水輪機(jī)內(nèi)復(fù)雜流場的形成原因及其與流場邊界之間的作用結(jié)果;然而，由于三維流動伴隨空化這一現(xiàn)象的復(fù)雜性，目前在眾多發(fā)生流動空化的場合中，實(shí)驗(yàn)研究仍存在很多困難.如導(dǎo)葉之間與轉(zhuǎn)輪葉片之間的流場測試十分復(fù)雜，對于具有復(fù)雜三維流場的水輪機(jī)，采用數(shù)值模擬方法研究空化非定常流動成為觀測空化現(xiàn)象發(fā)展的重要技術(shù).

對大型混流式水輪機(jī)機(jī)組內(nèi)部空化流動數(shù)值研究的主要目的是從宏觀和微觀兩方面揭示流動穩(wěn)定性問題.目前，對于葉片式流體機(jī)械內(nèi)部流場的數(shù)值計算，仍以雷諾平均數(shù)值模擬(RANS)方法為主，而該方法中雷諾應(yīng)力基于線性模型計算，使某些包含細(xì)節(jié)信息的高頻分量被濾除，數(shù)值模擬不能體現(xiàn)真實(shí)的多尺度湍流流動.大渦模擬(LES)的思想是將小尺度的渦從流動中分離出來，通過小尺度渦和大尺度渦之間的相互影響來描述不同尺度渦結(jié)構(gòu)的相互作用.數(shù)學(xué)表達(dá)上，可在大渦流場的運(yùn)動控制方程中添加相應(yīng)的附加應(yīng)力項(xiàng)來體現(xiàn)小渦的影響，該附加應(yīng)力項(xiàng)即為亞格子(SGS)應(yīng)力項(xiàng).一般認(rèn)為湍流中小尺度渦的結(jié)構(gòu)趨于各向同性，基于不同渦特征即可構(gòu)造不同的SGS模型.在水輪機(jī)數(shù)值模擬方面，Zhang等［11］對某一混流式水輪機(jī)進(jìn)行了三維LES，并將計算結(jié)果與壓感片實(shí)驗(yàn)獲得的結(jié)果作了比較，得到了令人滿意的結(jié)果，但其研究僅限于葉片部位;Liu等［12］采用LES研究了混流式水輪機(jī)尾水管渦帶的形態(tài)，通過對比不同負(fù)荷下的渦帶形態(tài)，發(fā)現(xiàn)越偏離設(shè)計工況，漩渦強(qiáng)度越大;黃劍峰等［13］使用LES動態(tài)亞格子湍流模型和先進(jìn)的滑移網(wǎng)格技術(shù)對混流式水輪機(jī)全流道動靜干涉的三維非定常湍流場進(jìn)行了動態(tài)LES，計算結(jié)果表明LES對復(fù)雜流場的仿真效果更加準(zhǔn)確.針對上游尾跡對入口流動的擾動，Wang等［14］對某一混流式水輪機(jī)的單流道模型進(jìn)行了LES模擬，發(fā)現(xiàn)強(qiáng)尾跡將造成葉片前緣處的大面積流動分離，并沿展向形成主渦.Ma等［15］采用混流式水輪機(jī)單流道模型，驗(yàn)證了LES中SGS模型對計算結(jié)果的影響，并和實(shí)驗(yàn)進(jìn)行了對比，發(fā)現(xiàn)適應(yīng)壁面的渦黏模型和實(shí)驗(yàn)結(jié)果符合最好.而最近，RANS-LES混合湍流模型也被逐漸應(yīng)用在水輪機(jī)模擬中［16-17］，較好的模擬結(jié)果體現(xiàn)了LES自身特性帶來的優(yōu)勢.

由于水輪機(jī)全流道的幾何復(fù)雜性，上、下游部件之間流動干擾，以及旋轉(zhuǎn)部件與靜止部件之間的流動干涉，使全流道計算及分析相對困難;另外，由于空化的產(chǎn)生是由于液體在低壓區(qū)發(fā)生相變，如何準(zhǔn)確地捕捉到空化區(qū)域也成為了必要問題.目前，對于水輪機(jī)的數(shù)值研究，很少進(jìn)行全特性驗(yàn)證.本文用LES方法，并結(jié)合空化模型，對某一大型混流式水輪機(jī)在不同活動導(dǎo)葉開度下的特性及在不同工況點(diǎn)的空化流場進(jìn)行全通道三維數(shù)值計算，對比分析單相流與空化兩相流計算結(jié)果，并對外特性和內(nèi)流場進(jìn)行實(shí)驗(yàn)結(jié)果對比，驗(yàn)證針對偏工況運(yùn)行的混流式水輪機(jī)的基于空化模型的大渦模擬方法.

1 計算模型與數(shù)值模擬方法

1.1 水輪機(jī)模型與幾何參數(shù)

對某大型混流式水輪機(jī)的模型進(jìn)行數(shù)值模擬，幾何模型如圖1所示.從流動入口到流動出口依次分為蝸殼、固定導(dǎo)葉、活動導(dǎo)葉、轉(zhuǎn)輪及尾水管5個通道部件，其中入口段前面增加了直管段部分，用來穩(wěn)定蝸殼的入口流動狀態(tài).該模型的具體幾何參數(shù)參見表1.

表1 模型水輪機(jī)基本幾何參數(shù)

圖1 混流式水輪機(jī)幾何模型

1.2 控制方程

采用LES方法(Smagorinsky-Lilly模型)對混流式水輪機(jī)內(nèi)部流動進(jìn)行數(shù)值模擬計算.對不可壓縮牛頓流體流動，用過濾算子對Navier-Stokes方程進(jìn)行變量過濾，可得基于LES的動量守恒方程為

式中:ui為i方向流動速度，xi為笛卡爾坐標(biāo)系下x、y、z3 個分量，p為壓力，ρ為流體的密度，υ為流體的運(yùn)動黏度，τij為SGS應(yīng)力項(xiàng).

采用Smagorinsky-Lilly模型有

其中CS為Smagorinsky系數(shù)，Δ為LES過濾尺度，Sij為可解尺度的變形率張量，δij為克羅內(nèi)克符號.Smagorinsky-Lilly模型中，SGS應(yīng)力可表述為

式中，渦黏性μt為亞格子尺度的湍動黏度.

Smagorinsky常數(shù)CS可根據(jù)Kolmogorov常數(shù)CK獲得:

通常為減小SGS應(yīng)力的擴(kuò)散，提高計算過程收斂性，特別是在近壁區(qū)，CS取值可適當(dāng)減小.考慮到全通道中壁面對流動的強(qiáng)烈作用，在近壁區(qū)可以按照 Van Driest模型來修正CS值［18］，即

其中y+為網(wǎng)格節(jié)點(diǎn)到壁面的量綱一的距離，參數(shù)A+=25，CS0=0.1.

為更準(zhǔn)確模擬空化流動的發(fā)生情況，計算中除使用單相流動模型之外，還引入了描述空泡運(yùn)動的空化模型.商業(yè)CFD軟件ANSYS中目前已經(jīng)集成了基于Rayleigh-Plesset方程的幾種空化模型，如 Singhal模型、Schnerr-Sauer模型和 Zwart-Gerber-Belamri(ZGB)模型［19］，這幾種模型均描述了空化過程中的相變，考慮了流場中空泡的密度變化與其運(yùn)動特征的聯(lián)系.薛瑞等［20］和劉厚林等［21］針對不同的湍流場，對3種空化模型進(jìn)行了比較研究，發(fā)現(xiàn)ZGB模型的計算收斂性較佳，計算穩(wěn)定.故文中計算選用ZGB空化模型.

1.3 網(wǎng)格劃分與邊界條件

根據(jù)水輪機(jī)的幾何模型，依次按照蝸殼、固定導(dǎo)葉、活動導(dǎo)葉、轉(zhuǎn)輪、尾水管5個部分分別劃分網(wǎng)格，然后組合為模型水輪機(jī)的全流道網(wǎng)格，其中在相鄰部件之間，均設(shè)置了網(wǎng)格交界面(interface).網(wǎng)格劃分的最終結(jié)果如圖2所示，針對具有翼型結(jié)構(gòu)的導(dǎo)葉和轉(zhuǎn)輪部件，在其近壁面處均進(jìn)行了邊界層網(wǎng)格加密，以滿足壁面流動計算的要求.

圖2 網(wǎng)格劃分細(xì)節(jié)圖

對計算結(jié)果進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證，選定設(shè)計工況對不同節(jié)點(diǎn)數(shù)的網(wǎng)格進(jìn)行了模擬驗(yàn)證，通常無關(guān)性驗(yàn)證以模型實(shí)驗(yàn)效率為標(biāo)準(zhǔn)，以核對數(shù)值計算得到的水力效率.該計算工況下，活動導(dǎo)葉開度a=18 mm、單位轉(zhuǎn)速n11=60.6 r/min(即最優(yōu)工況)，模型實(shí)驗(yàn)效率為η=94.8%.計算發(fā)現(xiàn)，在網(wǎng)格數(shù)超過一定數(shù)量后，隨網(wǎng)格節(jié)點(diǎn)數(shù)目的增加，計算得到的水輪機(jī)水力效率值趨于穩(wěn)定，并和模型實(shí)驗(yàn)的效率相吻合，故當(dāng)計算域網(wǎng)格數(shù)量達(dá)到829萬時即可滿足要求.該工況下葉片表面的y+值見圖2，可以看出，y+值符合大渦模擬的計算要求.水輪機(jī)各個部件的網(wǎng)格參數(shù)如表2所示.

表2 計算域網(wǎng)格結(jié)點(diǎn)數(shù)及單元數(shù)

邊界條件設(shè)置中，流動入口設(shè)置為計算域進(jìn)口(蝸殼前直管段);采用總壓進(jìn)口條件，該值通過模型實(shí)驗(yàn)中所采用的水頭和進(jìn)口水流速度來給定;流動出口設(shè)置為尾水管的出口，采用壓力出口;體壁面均采用無滑移邊界條件.

1.4 求解方法

選取ANSYS中的CFX求解器進(jìn)行水輪機(jī)內(nèi)部全流道三維模擬計算.針對旋轉(zhuǎn)部件(轉(zhuǎn)輪)與靜止部件(活動導(dǎo)葉出口、尾水管進(jìn)口)的交界面，坐標(biāo)變換模型均采用瞬態(tài)動靜干涉(Transient Rotor-Stator)模式.對流項(xiàng)離散格式選用 High resolution格式，同時計算的每個時間步長中保證轉(zhuǎn)輪旋轉(zhuǎn)1°，計算迭代100次，收斂的殘差標(biāo)準(zhǔn)為1×10-5.考慮到流場中近壁區(qū)的局部低雷諾數(shù)性質(zhì)，在LES求解中，近壁區(qū)的流動采用標(biāo)準(zhǔn)的壁面函數(shù)來處理.

2 計算結(jié)果及分析

數(shù)值模擬中，選取與Su等［22］的實(shí)驗(yàn)測量相同的流動工況，即在偏離最優(yōu)導(dǎo)葉開度a=18 mm的兩條活動導(dǎo)葉開度線a=10 mm和a=14 mm上，分別選取5個工況點(diǎn)進(jìn)行LES(單相流和空化流動)計算，具體計算工況分布如表3所示.

表3 計算工況點(diǎn)參數(shù)

2.1 內(nèi)流場和外特性驗(yàn)證

為對比LES與RANS對混流式水輪機(jī)內(nèi)部流場和空化性能預(yù)測的精度，Su等［23］對某一偏工況下的流場進(jìn)行了三維全通道單相流LES，結(jié)果表明，LES比RANS能很好地再現(xiàn)流場中的空化渦.為了描述水輪機(jī)內(nèi)部流場的特征，下面給出對活動導(dǎo)葉開度14 mm時單位轉(zhuǎn)速為65.4 r/min、單位流量為0.346 7 m3/s的計算結(jié)果，計算啟用空化模型，即為多相流LES計算，并將結(jié)果和未啟用的空化模型的單相LES進(jìn)行比較.

圖3為兩種計算方法對無葉區(qū)附近流場的預(yù)測情況，選取沿導(dǎo)葉高度方向的中間平面，并以靠近固定導(dǎo)葉側(cè)的測量線端點(diǎn)為x軸起始點(diǎn)，圖中給出了基于平均速度場的流線及該測量線上量綱一的速度U/U0分布，其中U為當(dāng)?shù)亟^對速度，U0為根據(jù)蝸殼進(jìn)口流量計算得到的平均速度.

圖3 無葉區(qū)流場數(shù)值模擬與實(shí)驗(yàn)結(jié)果對比(a=14 mm)

定量分析表明，活動導(dǎo)葉之后的尾跡流動在導(dǎo)葉壓力面一側(cè)(轉(zhuǎn)輪側(cè))流速較高，而在吸力面一側(cè)(固定導(dǎo)葉側(cè))流速較低，這是翼型擾流的顯著特征.通過對比文獻(xiàn)［22］中PIV實(shí)驗(yàn)的流場分布，發(fā)現(xiàn)兩種方法獲得的活動導(dǎo)葉尾跡平均流場流線走勢與PIV實(shí)驗(yàn)結(jié)果一致，LES空化計算獲得的測量線上速度值的大小與PIV實(shí)驗(yàn)獲得的結(jié)果符合得更好.

圖4給出了不同活動導(dǎo)葉開度下LES計算結(jié)果與實(shí)驗(yàn)結(jié)果的對比，包含了啟用和未啟用空化模型時的情況.可見，單位轉(zhuǎn)速與單位流量的關(guān)系與模型實(shí)驗(yàn)符合很好，表明采用空化計算的LES能夠很好地預(yù)測水輪機(jī)模型的外特性，圖4(b)中活動導(dǎo)葉開度a=18 mm的計算結(jié)果作為對比.從圖4(a)和(b)可以看出，采用空化兩相LES計算的結(jié)果與實(shí)驗(yàn)更加吻合，而單相LES的計算結(jié)果與實(shí)驗(yàn)值有些許偏離.在上述工況下，同一開度下的水輪機(jī)單位流速與單位流量曲線接近線性變化關(guān)系，且3個開度下的特性曲線基本平行，即該范圍內(nèi)，水輪機(jī)運(yùn)行時的實(shí)際轉(zhuǎn)速與全流道內(nèi)流量也應(yīng)成線性關(guān)系.

圖4 不同活動導(dǎo)葉開度下的外特性對比

圖5給出了活動導(dǎo)葉開度a=10 mm時通過空化LES空化計算水輪機(jī)中真實(shí)流量-轉(zhuǎn)速(Q～N)與效率 -轉(zhuǎn)速(η-N)的關(guān)系曲線，其中a=18 mm的計算結(jié)果作為對比.由于在水頭不變的情況下，流動中n11與N成線性變化關(guān)系，圖5的結(jié)果也同時驗(yàn)證了圖4的結(jié)論，即在上述工況變化范圍中，n11～Q11與Q～N均呈現(xiàn)線性關(guān)系.額定工況下，計算得混流式水輪機(jī)的水力效率接近95%，且水力效率隨著轉(zhuǎn)輪轉(zhuǎn)速的降低而減小.當(dāng)轉(zhuǎn)速N＜800 r/min時，η迅速下降.在小開度a=10 mm工況，水力效率下降較小，但運(yùn)行效率總體較低，約87%左右.

為了驗(yàn)證空化LES對運(yùn)行效率計算的準(zhǔn)確程度，圖6給出了效率曲線與模型水輪機(jī)實(shí)驗(yàn)獲得綜合特性曲線的對比.圖中數(shù)據(jù)點(diǎn)來自于圖5，且將實(shí)際轉(zhuǎn)輪轉(zhuǎn)速轉(zhuǎn)換成了單位轉(zhuǎn)速.

圖6中，計算得到的效率點(diǎn)(η，n11)由圓圈標(biāo)示，圓圈中心標(biāo)示了效率點(diǎn)的位置.可以看出，效率點(diǎn)均非常接近于模型實(shí)驗(yàn)的開度線a=10 mm和a=18 mm，即空化計算的效率曲線與綜合特性曲線相符合.另外，對比模型綜合特性曲線可知，設(shè)計工況處效率曲線陡峭，而偏工況時，效率曲線比較平緩.可見，啟用空化模型的LES可以很好地再現(xiàn)并預(yù)測水輪機(jī)的外特性.

圖5 不同活動導(dǎo)葉開度下空化兩相流計算的轉(zhuǎn)速－流量與轉(zhuǎn)速－效率曲線

圖6 使用空化LES計算得到不同開度下效率曲線與綜合特性曲線對比圖

2.2 葉道渦預(yù)測及其演化

葉道渦是轉(zhuǎn)輪內(nèi)部發(fā)生空化之后形成的渦結(jié)構(gòu)，其潰滅時會對葉片表面產(chǎn)生嚴(yán)重的侵蝕，所以準(zhǔn)確預(yù)測葉道渦的生成與發(fā)展對水輪機(jī)正常運(yùn)轉(zhuǎn)具有重要意義.為了對葉道渦形態(tài)進(jìn)行可視化，ANSYS軟件提供了許多渦結(jié)構(gòu)特征提取方法，如螺旋度(Helicity)方法、渦流強(qiáng)度(Swirling strength)方法、渦量(Vorticity)方法、Q方法和λ2方法等.

用渦量準(zhǔn)則顯示轉(zhuǎn)輪區(qū)的渦量分布.通過給定渦量的閾值，渦量準(zhǔn)則可以正確判定出漩渦存在的區(qū)域，下面分別對啟用和不啟用空化模型時的工況進(jìn)行分析.針對活動導(dǎo)葉開度a=10 mm、n11=59.9 r/min工況，先分別計算出空化LES和單相LES模擬結(jié)果的平均速度場，再采用0.5倍平均渦量 ＜|Ω|＞作為閾值進(jìn)行葉道渦顯示，如圖7所示.可以看出，空化LES可以捕捉到更多的葉道渦，且更加符合模型實(shí)驗(yàn)獲得的葉道渦形態(tài)［22］;但單相LES則不能很好地捕捉到沿流向的葉道渦運(yùn)動.

圖7 平均流場的葉道渦可視化

針對葉道渦的演化過程，下面將進(jìn)一步對渦結(jié)構(gòu)的變化進(jìn)行說明.如圖8所示，計算工況為活動導(dǎo)葉開度14 mm，單位轉(zhuǎn)速61.0 r/min，進(jìn)行空化LES計算.

圖8 空化LES計算獲得的葉道渦演化過程

圖8中給出了水輪機(jī)轉(zhuǎn)輪內(nèi)部葉道渦生成、發(fā)展和潰滅的一個典型演化過程，空化渦使用Q方法進(jìn)行提取，并和相同工況下的模型實(shí)驗(yàn)進(jìn)行了比較.可見，葉道渦位于葉片通道之間，并初生于葉片進(jìn)水邊附近，隨時間變化，轉(zhuǎn)輪順時針旋轉(zhuǎn)，空化渦體積逐漸增長(1～3)，達(dá)到最大之后，又逐漸減小直到消失(4～5).需要注意的是，從1～5，在葉片的出水邊(模型實(shí)驗(yàn)未監(jiān)測葉片出水邊)附近，有大量空化渦沿著展向分布，這是由于葉道渦在不斷生成積累過程中，隨主流流動而被帶到了出水邊.部分空化渦將被進(jìn)一步帶到流動通道的下游，在尾水管中形成碎渦結(jié)構(gòu).

文獻(xiàn)［22］已通過實(shí)驗(yàn)發(fā)現(xiàn)，轉(zhuǎn)輪在偏工況下運(yùn)行時，其中葉道渦的變化呈現(xiàn)高頻特征，且其演化頻率為轉(zhuǎn)頻的整數(shù)倍，該倍數(shù)一般為轉(zhuǎn)輪葉片數(shù).在上述計算中，設(shè)置時間步長為轉(zhuǎn)輪每轉(zhuǎn)動1°的時間，可以得到葉道渦變化的一個周期為轉(zhuǎn)輪轉(zhuǎn)動 24°，對應(yīng)時間為 0.005 2 s，約 192.31 Hz.而該流動條件下，單位轉(zhuǎn)速為61.0 r/min，轉(zhuǎn)輪實(shí)際轉(zhuǎn)速 757.8 r/min，故轉(zhuǎn)頻為 12.63 Hz，可見葉道渦演化的周期正好約為轉(zhuǎn)頻的15倍，和葉片數(shù)相同.于是，上述數(shù)值模擬的結(jié)果和模型實(shí)驗(yàn)的結(jié)果相符合，即使用空化LES計算能夠很好地預(yù)測水輪機(jī)偏工況運(yùn)行時的葉道渦運(yùn)動.

2.3 尾水管渦帶演化

尾水管渦帶是水輪機(jī)中另一種空化現(xiàn)象，從轉(zhuǎn)輪出口的中心開始，向肘管發(fā)展.由于轉(zhuǎn)輪出口水流的周期性運(yùn)動，該渦帶也在不斷旋轉(zhuǎn)運(yùn)動，隨不同工況的變化渦帶具有不同的形態(tài)，其運(yùn)動會引發(fā)尾水管壁面的周期性壓力脈動.

針對活動葉開度10 mm、單位轉(zhuǎn)速64.9 r/min工況分別使用空化LES和單相LES計算，提取尾水管中空化渦帶結(jié)構(gòu)，如圖9所示.將25℃下水的汽化壓力作為等值面提取渦形態(tài)，可見，空化LES計算得到的渦帶形態(tài)與模型實(shí)驗(yàn)觀測到的結(jié)果接近，該工況下尾水管渦帶呈現(xiàn)螺旋狀.另外在該工況下單相LES出現(xiàn)了分岔的渦帶形態(tài)，說明單相流計算不能很好地處理氣液密度變化和壓力變化的關(guān)系，而空化兩相流則考慮了空泡密度與流場運(yùn)動的關(guān)系，故能獲得較為滿意的空化形態(tài).

為了獲取尾水渦帶的演化過程，對活動導(dǎo)葉開度14 mm、單位轉(zhuǎn)速61.0 r/min的工況進(jìn)行空化LES計算，如圖10所示.圖中顯示為從尾水管開始的錐管段渦帶部分，事實(shí)上渦帶發(fā)展是從轉(zhuǎn)輪下方的泄水錐開始，為了顯示清楚渦帶的旋轉(zhuǎn)情況，特意截取尾水管中的一段空化渦帶.圖10表征了尾水管渦帶旋轉(zhuǎn)一周的演化過程，其所用時間約為轉(zhuǎn)輪旋轉(zhuǎn)2.5周的時間，可見，渦帶的演化是一個低頻過程，即約為轉(zhuǎn)頻的0.4倍，和文獻(xiàn)［22］的結(jié)論相符合.

圖9 尾水管渦帶形態(tài)對比

圖10 空化LES計算獲得的空化渦帶演化過程

由形態(tài)結(jié)構(gòu)可看出，從泄水錐生成的空化渦被攜帶到下游時，會隨著轉(zhuǎn)輪的旋轉(zhuǎn)而變成螺旋狀結(jié)構(gòu)，并不斷增長，獲得的空化渦結(jié)構(gòu)和模型實(shí)驗(yàn)比較符合.而中心渦帶在增長的同時，在接近尾水管壁面的部分有破碎渦生成，可以推測，部分破碎渦來源于轉(zhuǎn)輪出水邊集結(jié)的空化渦.這些破碎渦在渦帶旋轉(zhuǎn)一周的時間中，將從下游的肘管處瀉出，之后，渦帶形態(tài)將開始下一個周期變化.從機(jī)理上看，由于空化LES計算特別考慮了空化區(qū)域的壓力、密度變化和流場運(yùn)動的關(guān)系，故能更好地處理流場中壓力梯度變化劇烈的區(qū)域.可見通過LES可以預(yù)測尾水管入口處壓力梯度較高的區(qū)域，和實(shí)際流動情況相符合.

綜上所述，采用空化模型的LES不僅可以準(zhǔn)確地獲得混流式水輪機(jī)內(nèi)部三維流場，更能很好地捕捉到空化區(qū)域的形態(tài)特征，尤其針對葉道渦和尾水管渦帶的形態(tài)演化，對水輪機(jī)的穩(wěn)定運(yùn)行具有重要的意義.

3 結(jié)論

1)LES可以預(yù)測尾水管入口存在壓力梯度較高的區(qū)域，且能捕捉到尾水管錐管段復(fù)雜的渦帶旋轉(zhuǎn)流動，而RANS則不能準(zhǔn)確預(yù)測.

2)采用空化兩相流LES計算的外特性結(jié)果與實(shí)驗(yàn)對比吻合度高，而單相流計算結(jié)果與實(shí)驗(yàn)值偏差較大.

3)采用空化兩相流的LES可以很好地捕捉到轉(zhuǎn)輪葉片之間沿流向發(fā)展的葉道渦，與實(shí)驗(yàn)結(jié)果相吻合，驗(yàn)證了基于空化模型的LES方法在預(yù)測偏工況下混流式水輪機(jī)性能的可行性與準(zhǔn)確性.

4)空化LES獲得的葉道渦演化過程為高頻運(yùn)動，而尾水管渦帶為低頻運(yùn)動，和模型實(shí)驗(yàn)結(jié)果相符合.

［1］D?RFLER P，SICK M，COUTU A.Flow-induced pulsation and vibration in hydroelectric machinery［M］.London:Springer，2013.

［2］WEI X Z，SU W T，LI X B，et al.Effect of blade perforation on Francis hydro-turbine cavitation characteristics ［J］.Journal of Hydraulic Research，2014，52(3):412-420.

［3］ILIESCU M S，CIOCAN G D，AVELLAN F.Analysis of the cavitating draft tube vortex in a francis turbine using particle image velocimetry measurements in two-phase flow［J］.ASME Journal of Fluid Engineering，2008，130:21105.

［4］PALAFOX P，OIDFIELD M L G，LAGRAFF J E，et al.PIV maps of tip leakage and secondary flow fields on a low-speed turbine blade cascade with moving end wall［J］.Physics of Fluids，2008，130:11001.

［5］UZOL O，CHOW Y C，KATZ J，et al.Average passage flow field and deterministic stresses in the tip and hub regions of a multistage turbomachine［J］.Experiments in Fluids，2003，125(4):714-725.

［6］ UZOL O，KATZ J.Flow measurement techniques in turbomachinery［M］//Handbook of Experimental Fluids Mechanics.Heidelberg:Springer，2007:917-957.

［7］KATZ J，CHOW Y C，SORANNA F，et al.Experimental Characterization of Flow Structure and Turbulence in Turbomachines［C］//5th International Symposium on Fluid Machinery and Fluids Engineering.Jeju:［s.n.］，2012.

［8］李丹，陳次昌，季全凱，等.PIV測試水輪機(jī)尾水管錐管流場［J］.西華大學(xué)學(xué)報:自然科學(xué)版，2004(S1):228-230.

［9］王軍，付之躍，溫國珍，等.HL220改型裝置不同工況下尾水管進(jìn)口流動狀態(tài)PIV試驗(yàn)研究［J］.工程熱物理學(xué)報，2005，26(S1):93-96.

［10］陳次昌，李丹，季全凱，等.混流式水輪機(jī)尾水管內(nèi)部流場的 PIV測試［J］.機(jī)械工程學(xué)報，2006，42(12):83-88.

［11］ZHANG L X，GUO Y K，WANG W Q.Large eddy simulation of turbulent flow in atrue 3D francis hydroturbine passage with dynamicalfluid-structure interaction［J］.International Journal for Numerical Methods in Fluids，2007，54(5):517-541.

［12］LIU X B，ZENG Y Z，CAO S Y.Numerical prediction of vortex flow in hydraulic turbine draft tube for LES［J］.Journal of Hydrodynamics，Series B，2005，17(4):448-454.

［13］黃劍峰，張立翔，王文全.混流式水輪機(jī)全流道三維非定常湍流場的動態(tài)大渦模擬［J］.中國農(nóng)村水利水電，2009(1):101-108.

［14］WANG W Q，ZHANG L X，YAN Y.Large-eddy simulation of turbulent flow considering inflow wakes in a Francis turbine blade passage［J］.Journal of Hydrodynamics，Series B，2007，19(2):201-209.

［15］MA J M，WANG F J，TANG X L.Comparison of several subgrid-scale models for large-eddy simulation of turbulentflows in water turbine［C］//The 4thInternational Symposium on Fluid Machinery and Fluid Engineering.Beijing:［s.n.］，2008.

［16］NETO A D A，JESTER-ZUERKER J，JUNG A，et al.Evaluation of a Francis turbine draft tube flow at part load using hybrid RANS-LES turbulence modelling［J］.IOP Conf Series:Earth and Environmental Science，2012(15):062010.

［17］SENTYABOV A V，GAVRILOV A A，DEKTEREV A A，et al.Numerical investigation of the vortex core precession in a model hydro turbine with the aid of hybrid methods for computation of turbulent flows［J］.Thermophysics and Aeromechanics，2014， 21(6):707-718.

［18］王福軍.計算流體動力學(xué)分析:CFD軟件原理與應(yīng)用［M］.北京:清華大學(xué)出版社，2004.

［19］ZWART P J，GERBER A G，BELAMRI T.A twophase flow model for predicting cavitation dynamics［C］//Fifth International Conference on Multiphase Flow.Yokohama:ICMF COP，2004:No.152.

［20］薛瑞，張淼，許戰(zhàn)軍，等.對不同空化模型的比較研究［J］.西北水電，2014(2):85-89.

［21］劉厚林，劉東喜，王勇，等.三種空化模型在離心泵空化流計算中的應(yīng)用評價［J］.農(nóng)業(yè)工程學(xué)報，2012，28(16):54-59.

［22］SU W T，LI X B，LI F C，et al.Experimental investigation on the characteristics of hydrodynamic stabilities in Francis hydro-turbine models［J］.Advances in Mechanical Engineering，2014，Article ID 486821.

［23］SU W T，LI F C，LI X B，et al.Assessment of LES performance in simulating complex 3D flows in turbomachines ［J］. Engineering Applications of ComputationalFluid Mechanics， 2012， 6(3):355-364.