基于固液雙介質(zhì)的導(dǎo)水機構(gòu)雙葉柵流場研究

劉福春 張婷 劉超 李翔宇

摘?要：基于Pro/E模型，Gambit網(wǎng)格，采用固液兩相流介質(zhì)，用計算流體力學(xué)（Computational Fluid Dynamics）方法分析了某電站導(dǎo)水機構(gòu)雙列葉柵的位置關(guān)系是否會導(dǎo)致水輪機水力損失的變化;對雙葉柵5種位置關(guān)系的壓力、速度分布及水力損失進行水力計算，證明雙葉柵的位置關(guān)系會對導(dǎo)水機構(gòu)流場產(chǎn)生較大作用。CFD計算結(jié)果表明：在最優(yōu)工況和額定工況下，雙葉柵相對位置為較大時，水力損失較小。

關(guān)鍵詞：固液兩相流;導(dǎo)水機構(gòu);流場分析;設(shè)計優(yōu)化

1 緒論

對現(xiàn)有文獻進行總結(jié)發(fā)現(xiàn)，水輪機導(dǎo)水機構(gòu)固定導(dǎo)葉和活動導(dǎo)葉的位置分布、導(dǎo)葉數(shù)量、導(dǎo)葉轉(zhuǎn)角、導(dǎo)葉的正負(fù)曲率都對水輪機性能有影響[1-5]。但在現(xiàn)有的水輪機導(dǎo)水機構(gòu)的仿真模擬研究中，對于活動導(dǎo)葉和固定導(dǎo)葉相對位置對導(dǎo)水機構(gòu)流暢的研究較少，也不能普遍適用于所有現(xiàn)存的水輪機，所以對于某新建電站的導(dǎo)水機構(gòu)流場仿真模擬仍具有實際意義[6-18]。本文對某電廠混流式機組水輪機導(dǎo)水機構(gòu)進行CFD計算，希望能對機組運行中的水輪機水力振動及損失有所參考。

2 計算方法和流動模型

2.1 水輪機的基本參數(shù)

轉(zhuǎn)輪直徑1.7m，最大水頭238m，最小水頭220m，額定水頭231m;導(dǎo)葉20個。泥沙顆粒直徑為0.5mm，體積濃度為5%。

2.2 計算模型及網(wǎng)格

使用Pro/E軟件對水輪機由蝸殼進口至尾水管出口進行三維建模。模型如圖1、2所示。用前處理軟件Gambit進行混合網(wǎng)格劃分。

2.3 水流控制方程

2.3.1 流體相的基本方程

2.3.2 固粒相的基本方程

固粒相粒徑為Dpk的粒子所滿足的連續(xù)方程和動量方程為：

對連續(xù)方程和動量方程在貼體坐標(biāo)下采用有限體積法在模型中進行分割。控制方程分割時，壓力相采用二階中心差分格式，其他項采用二階迎風(fēng)格式，采用壓力速度耦合算法[19]。

2.4 邊界條件給定

進口邊界為蝸殼進口截面流量給定，出口邊界為尾水管末端截面。模型外壁為組合封閉面，內(nèi)部為流場，模型壁面區(qū)域使用標(biāo)準(zhǔn)函數(shù)法處理[19]。

3 計算結(jié)果與分析

3.1 導(dǎo)葉的分布關(guān)系及導(dǎo)葉相對水力損失公式

將水輪機模型建立坐標(biāo)系，將導(dǎo)水機構(gòu)重心作為坐標(biāo)原點，將模型的右、前方向分別作為作為X、Y軸，與導(dǎo)水機構(gòu)垂直方向作為Z軸，根據(jù)導(dǎo)葉間隔18°，固定導(dǎo)葉不動，將活動導(dǎo)葉沿Z軸依次旋轉(zhuǎn)3.6°，得到以下方案。方案1雙葉柵位置為0°;方案2雙葉柵位置為3.6°;方案3雙葉柵位置為72°;方案4雙葉柵位置為10.8°;方案5雙葉柵位置為14.4°。

根據(jù)導(dǎo)水機構(gòu)的水力學(xué)特性，水力損失如下式所示：

3.2 結(jié)果與分析

本文分別對額定工況和最優(yōu)工況工況進行了計算。從導(dǎo)水機構(gòu)進口到出口，壓力、速度分布比較均勻，流線清晰，壓力和速度呈負(fù)相關(guān)，壓力在逐漸降低，而速度在升高。

CFD計算后，對五種方案下壓力和流速場分析表明，雙葉柵擺放位置不同對流場影響較大。在雙葉柵位置為14.4°附近時水力損失較小，壓力場和速度場均勻分布，不容易產(chǎn)生壓力突變區(qū)，一定程度上減弱了空化空蝕。

4 結(jié)論

（1）本文對某電站水輪機在兩種典型工況下進行了CFD仿真計算，得到了不同狀態(tài)下的導(dǎo)水機構(gòu)壓力、速度場及水力損失。

（2）雙葉柵相對位置較小時容易使壓力場局部增大，容易導(dǎo)致導(dǎo)葉局部磨損。所以，建議在一定程度上擴大分布圓的半徑，降低導(dǎo)葉高度，盡量增大導(dǎo)葉的開度。

參考文獻：

[1]沈春和，桂林.低比轉(zhuǎn)速水輪機活動導(dǎo)葉固液兩相數(shù)值模擬研究[J].貴州電力技術(shù)，2015，18（05）：26-30+33.

[2]李玲娟.混流式水輪機內(nèi)部固液兩相流數(shù)值模擬及磨損預(yù)估[D].西華大學(xué)，2015.

[3]黃劍峰，張立翔，姚激，龍立焱.水輪機泥沙磨損兩相湍流場數(shù)值模擬[J].排灌機械工程學(xué)報，2016，34（02）：145-150.

[4]華紅.含沙水中錦屏二級水電站水輪機內(nèi)部流動及性能預(yù)測研究[D].西華大學(xué)，2013.

[5]劉梅，宋文武，王興林，徐志堅.二重葉柵不同環(huán)量配比對水輪機組性能影響研究[J].中國農(nóng)村水利水電，2015（08）：184-188.

[6]柯強，劉小兵，周倩倩.混流式水輪機導(dǎo)葉相對位置對其水力性能的影響[J].大電機技術(shù)，2016（04）：29-35.

[7]王文全，李偉忠，閆妍.水輪機活動導(dǎo)葉勻速開關(guān)過程的三維湍流數(shù)值模擬[J].排灌機械工程學(xué)報，2017，35（02）：133-137.

[8]陳維勤，鈴木敏曉，陳梁年.混流式水輪機固定導(dǎo)葉的設(shè)計研究[J].水電站機電技術(shù)，2017，40（01）：8-11.

[9]范江峰，尤莉莎.基于CFD的不同曲率導(dǎo)葉對水輪機內(nèi)流場影響分析[J].水利科技與經(jīng)濟，2017，23（10）：22-27.

[10]肖玉紅.混流式水輪機蝸殼·導(dǎo)葉及轉(zhuǎn)輪內(nèi)部水流的二維與三維數(shù)值模擬比較分析[J].安徽農(nóng)業(yè)科學(xué)，2011，39（28）：17613-17615+17627.

[11]段銘鈺，吳承亮，黃逸哲.軸流泵全流道三維定常流場特性分析研究[J].裝備制造技術(shù)，2016（06）：28-31.

[12]文哲男，巴德純，鄧文娟，岳向吉.ISIGHT在混流式水輪機活動導(dǎo)葉優(yōu)化設(shè)計中的應(yīng)用[J].機械設(shè)計與制造，2017（05）：5-8.

[13]米紫昊，李太江，李勇.混流式水輪機導(dǎo)水機構(gòu)損傷三維數(shù)值模擬分析[J].人民黃河，2014，36（06）：129-131.

[14]楊光葉.混流式水輪機活動導(dǎo)葉CFD計算分析[J].現(xiàn)代制造技術(shù)與裝備，2014（05）：13-14.

[15]劉萬江，魏顯著，韓秀麗，許彬.高水頭混流式水輪機導(dǎo)葉位置關(guān)系對水力性能的影響[J].大電機技術(shù)，2013（06）：46-48+51.

[16]陳元林，覃大清.導(dǎo)葉雙列葉柵CFD數(shù)值計算及結(jié)果分析方法研究[J].大電機技術(shù)，2013（04）：38-41.

[17]劉勝柱，梁武科，郭鵬程，羅興.導(dǎo)葉雙列葉柵位置對水力性能的影響[J].水利水電技術(shù)，2004（01）：63-64+67-95.

[18]齊學(xué)義，周慧利，高志遠.含沙水流水輪機兩列導(dǎo)葉相對位置對活動導(dǎo)葉磨損的影響[J].蘭州理工大學(xué)學(xué)報，2013，39（01）：37-41.

[19]廖翠林，王福軍.X型葉片水輪機轉(zhuǎn)輪流場的數(shù)值模擬[J].水力發(fā)電學(xué)報，2008（03）：141-144.

作者簡介：劉福春（1990—?），男，漢族，山東德州人，碩士研究生，中級工程師，流體機械及其自動化;張婷（1990—?），女，漢族，山東聊城人，碩士研究生，中級工程師，電氣工程及其自動化;劉超（1989—?），男，漢族，河北保定人，碩士研究生，中級工程師，能源與動力工程;李翔宇（1992—?），女，漢族，河北張家口人，碩士研究生，助理工程師，電氣工程。