張忠孝

(東港市水利事務(wù)服務(wù)中心，遼寧 東港 118300)

0 引 言

泵站是水利工程不可或缺的一部分，可以實(shí)現(xiàn)農(nóng)業(yè)灌溉水、工業(yè)用水、生活用水的長距離運(yùn)輸，保障社會經(jīng)濟(jì)的可持續(xù)發(fā)展[1-6]。但是在泵站運(yùn)行過程中，進(jìn)水口附近區(qū)域時(shí)常會出現(xiàn)漩渦，由于漩渦會產(chǎn)生吸氣現(xiàn)象，從而導(dǎo)致空化、機(jī)組振動、機(jī)組運(yùn)行效率降低以及卷入漂浮物等問題，對于泵站的長期安全與穩(wěn)定運(yùn)行造成負(fù)面影響[7-8]。

影響泵站進(jìn)水口水力特性的因素較多，如進(jìn)水口方向、角度、進(jìn)水流量、進(jìn)水口布置形式、淹沒深度等[9-13]。其中，進(jìn)水口淹沒深度是較重要的因素之一，淹沒深度會導(dǎo)致泵站存在死水區(qū)域，當(dāng)來水經(jīng)過死水區(qū)域時(shí)由于流速分布變得極不均勻，就容易形成漩渦。淹沒深度不僅關(guān)系著泵站的長期安全運(yùn)行，而且還與工程成本直接相關(guān)，淹沒深度過大或者過小均不能達(dá)到良好的工程效益。因此，有必要針對復(fù)雜工況下的泵站進(jìn)水口水力特性進(jìn)行研究。

本文利用數(shù)值模擬方式，對0.8、1和1.2 m3/s共3種流量以及-0.1、0、0.1、0.3、0.5、0.7、0.9和1.1 m共8種淹沒深度下的泵站進(jìn)水口漩渦水力特性進(jìn)行模擬分析，以期能為泵站設(shè)計(jì)和運(yùn)行提供參考。

1 基本理論

在流體力學(xué)當(dāng)中，主要存在連續(xù)方程、動量守恒方程以及能量守恒方程等三大控制方程，由于自然界水體屬于單向流動，不涉及化學(xué)反應(yīng)，因而本文僅考慮流體的連續(xù)性和動量守恒。選定控制方程后，需要進(jìn)行計(jì)算區(qū)域的離散化，常用的控制方程離散化方法有有限差分法、有限單元法、有限體積法等，其中有限體積法具有運(yùn)算速度快的優(yōu)點(diǎn)，故本文采用有限體積法。前人研究表明，Realizablek-ε模型在進(jìn)行泵站前池?cái)?shù)值模擬時(shí)，漩渦范圍與實(shí)際情況更加符合，因而本文選用Realizablek-ε模型作為數(shù)值模擬時(shí)的湍流模型[14]。

2 數(shù)值模型

2.1 幾何模型及邊界條件

采用UG10軟件構(gòu)建泵站進(jìn)水口幾何模型。該進(jìn)水口包括5臺泵站機(jī)組(流道頂板厚度h1=1.2 m)，將建立好的泵站幾何模型導(dǎo)入Fluent數(shù)值模擬軟件中進(jìn)行網(wǎng)格劃分。在Fluent軟件中，泵站出口處設(shè)置為自由出流邊界條件，泵站進(jìn)口處選用流量進(jìn)口，自由表面選用剛蓋假定，泵站前池、流道邊壁、引渠等結(jié)構(gòu)均采用無滑移的邊界條件[15]。

2.2 網(wǎng)格劃分

采用Mesh軟件對泵站進(jìn)行網(wǎng)格劃分，將泵站劃分為四面體網(wǎng)格結(jié)構(gòu)。為了計(jì)算精確，在邊墩等位置處采取網(wǎng)格加密處理，正常網(wǎng)格的大小尺寸為0.2 m的四面體，加密網(wǎng)格的大小尺寸為0.02 m的四面體，整個(gè)模型網(wǎng)格數(shù)量達(dá)到202×104個(gè)，經(jīng)網(wǎng)格劃分后的模型見圖1。

圖1 泵站進(jìn)水口數(shù)值模型

網(wǎng)格劃分太稀會增加模型計(jì)算的離散誤差，導(dǎo)致仿真計(jì)算結(jié)果與實(shí)際結(jié)果相差較大；網(wǎng)格劃分太密，會導(dǎo)致模擬計(jì)算時(shí)間大大延長，影響計(jì)算分析效率，因而必須進(jìn)行網(wǎng)格無關(guān)性分析。以前池水力損失作為泵站模型網(wǎng)格劃分無關(guān)性的特征參數(shù)，分別計(jì)算得到前池水力損失隨網(wǎng)格數(shù)量的變化曲線，見圖2。從圖2中可以看到，隨著網(wǎng)格數(shù)量的增加，前池水力損失呈先增大后逐漸穩(wěn)定的變化特征，當(dāng)網(wǎng)格數(shù)量>200×104個(gè)時(shí)，前池水力損失為0.021 5 m，之后便不再隨網(wǎng)格數(shù)量的變化而變化。因此，本文劃分的網(wǎng)格數(shù)量為202×104個(gè)具有無關(guān)性。

圖2 前池水力損失隨網(wǎng)格數(shù)量變化

2.3 工況設(shè)計(jì)

為了更全面地了解各種工況下泵站進(jìn)水口的漩渦特性，設(shè)計(jì)0.8、1和1.2 m3/s這3種不同流量，每種流量下進(jìn)水口淹沒深度分別為-0.1、0、0.1、0.3、0.5、0.7、0.9和1.1 m。以1#機(jī)組進(jìn)水口為研究對象，在每種工況下，從上往下選取表面層、進(jìn)水口頂板、進(jìn)水口頂板以下0.2 m共3個(gè)水平斷面進(jìn)行漩渦分析。

3 結(jié)果分析

3.1 渦量值與水位關(guān)系

不同工況下渦量隨水位的變化關(guān)系見圖3。

圖3 不同工況下渦量隨水位變化特征曲線

從圖3中可知，當(dāng)淹沒深度為-0.1和0 m時(shí)，隨著水位的上升，渦量值逐漸增大；當(dāng)淹沒深度≥0.1 m后，隨著水位上升，渦量呈先增大后減小的變化特征，最大渦量值出現(xiàn)在水位1.2 m處。這說明在水位較淺的表層區(qū)域，渦流強(qiáng)度隨著水深增大而逐漸增大，渦流強(qiáng)度在頂板附近達(dá)到最大值；此后隨著水深增加，渦流強(qiáng)度逐漸減小，這主要是因?yàn)楫?dāng)泵站存在淹沒深度時(shí)，在進(jìn)水口處就會形成死水區(qū)域，當(dāng)外面流水進(jìn)入這些死水區(qū)域后，會造成這個(gè)區(qū)域水流流速紊亂，表面層作為水體跟空氣的交界面，在交界面處存在明顯的剪切作用，并跟隨水體深度生不斷的能量交換，當(dāng)來流運(yùn)動到隔墩附近時(shí)就會出現(xiàn)脫流現(xiàn)象，從而導(dǎo)致在交界面出現(xiàn)角動量，為渦流的產(chǎn)生創(chuàng)造了條件；當(dāng)水位達(dá)到一定深度后，此時(shí)渦量耗散能量逐漸大于積累的能量，因而渦量值會逐漸減低。

3.2 渦量值與流量關(guān)系

不同流量下渦核區(qū)域渦量值見表1。從表1中可以看到，當(dāng)流量一定時(shí)，統(tǒng)一淹沒深度下，進(jìn)水口頂板處的渦量值最大，其次為進(jìn)水口頂板以下0.2 m，最小的為表面層處，這與上文分析結(jié)果一致；在同一位置處，隨著淹沒深度的增加，渦量值呈先增大后減小的變化特征，最大渦量值出現(xiàn)在淹沒深度為0.5 m時(shí)，這說明在淹沒深度為0.5 m時(shí)，泵站進(jìn)水口死水區(qū)域更容易形成渦流并有利于渦流的持續(xù)發(fā)展。對比相同位置和淹沒深度下的不同流量的渦量值可以發(fā)現(xiàn)，隨著流量的增大，渦量最大值逐漸增大，這是因?yàn)榱髁吭酱?，水的流速越大，動能越大，漩渦強(qiáng)度越大，因而渦量值也越大。

表1 不同工況渦量值與流量的關(guān)系

3.3 最大渦量值與淹沒深度的關(guān)系

不同工況下最大渦量值隨淹沒深度的變化關(guān)系見圖4。從圖4中可以看到，相同流量下，當(dāng)淹沒深度小于0.5 m時(shí)，隨著淹沒深度的增加，泵站進(jìn)水口處的最大渦量值逐漸增大；當(dāng)淹沒深度大于0.5 m時(shí)，最大渦量值隨著淹沒深度的增加逐漸減小，最大渦量值出現(xiàn)在淹沒深度0.5 m處，隨著流量的增大，最大渦量值也逐漸增大，當(dāng)泵站進(jìn)水口寬高比一定時(shí)，存在臨界淹沒深度，當(dāng)淹沒深度達(dá)到這一數(shù)值時(shí)，進(jìn)水口將發(fā)生較為嚴(yán)重的漩渦現(xiàn)象。

圖4 不同工況最大渦量值與淹沒深度的關(guān)系

3.4 討 論

通過以上分析可知，泵站進(jìn)水口漩渦與流量和淹沒深度息息相關(guān)，特別是淹沒深度對渦量值的影響最為明顯；在不同的流量下，最大渦量值均出現(xiàn)在流道頂板區(qū)域附近，且當(dāng)淹沒深度為0.5 m時(shí)，渦量值最大。從理論上講，為了減小漩渦對泵站進(jìn)水口水力特性的影響，應(yīng)將淹沒深度設(shè)置為小于0.5 m或者大于0.5 m，但是當(dāng)淹沒深度小于0.5 m時(shí)，泵站進(jìn)水口處的流速分布又會變得異常紊亂、極不均勻；當(dāng)淹沒深度為0.9 m左右時(shí)，流速分布最為均勻，因此從改善流速分布和減小渦量綜合考慮，將淹沒深度設(shè)置為0.9 m最為合適。

4 結(jié) 語

采用數(shù)值模擬方式，對不同流量和淹沒深度下的泵站進(jìn)水口漩渦水力特性進(jìn)行分析，結(jié)論如下：

1) 淹沒深度一定時(shí)，不同流量工況下，最大渦量值均出現(xiàn)在流道頂板附近區(qū)域；隨著流量的增加，漩渦值逐漸增大。

2) 流量一定時(shí)，隨著淹沒深度的增加，最大渦量值呈先增大后減小的變化特征；當(dāng)淹沒深度為0.5 m時(shí)，渦量值最大。

3) 從流速分布均勻性和控制渦量大小綜合考慮，認(rèn)為淹沒深度設(shè)置為0.9 m時(shí)，泵站進(jìn)水口水力特性最佳。

4) 受限于研究數(shù)據(jù)和資料的不足，本文僅對流量和淹沒深度兩個(gè)因素進(jìn)行了模擬分析，關(guān)于其他因素的綜合影響還需要在今后做進(jìn)一步研究。