為了研究泵系統(tǒng)調(diào)閥過程的瞬態(tài)特性和內(nèi)流機(jī)理，在一維分析軟件Flowmaster中建立了包含管路、閥門和泵在內(nèi)的仿真模型，并以三維簡(jiǎn)化閘閥為模型，采用Fluent6.2進(jìn)行計(jì)算，對(duì)開啟過程的非定常內(nèi)部流動(dòng)進(jìn)行數(shù)值模擬研究.采用動(dòng)網(wǎng)格的方法分析了閥門開啟過程中閥芯運(yùn)動(dòng)引起的流場(chǎng)變形。結(jié)果表明：直線特性和對(duì)數(shù)特性的調(diào)節(jié)閥都具有快開特性，即流量變化對(duì)閥門的相對(duì)開度相當(dāng)敏感，當(dāng)閥門開度為10%～20%時(shí)，水擊壓力迅速下降;而通過內(nèi)部流態(tài)分析可知，在閥門開度較小的工況下，閥后流場(chǎng)紊亂，造成較大的水力損失，使阻力系數(shù)值增加，當(dāng)閥門開度小于50%時(shí)，穩(wěn)態(tài)和瞬態(tài)工況下閥門的阻力系數(shù)值有較大的區(qū)別由分析可知，研究閥門開啟過程的瞬態(tài)特性，以及建立內(nèi)部流態(tài)模型，都不能完全按照通常的穩(wěn)態(tài)理論進(jìn)行，尤其對(duì)閥門開度較小的工況，應(yīng)對(duì)其進(jìn)行一定程度的修正，以保證計(jì)算結(jié)果的正確性。

　　閥門作為管道系統(tǒng)中一種阻力可變的節(jié)流元件，通過改變其開度，可以改變管道系統(tǒng)的工作特性，從而實(shí)現(xiàn)調(diào)節(jié)流量和改變壓力的目的。它既是一種調(diào)節(jié)元件，同時(shí)也是一種控制元件，是實(shí)現(xiàn)管道系統(tǒng)安全經(jīng)濟(jì)輸送的重要設(shè)備.使用閥門調(diào)節(jié)流量和控制瞬變壓力的關(guān)鍵因素是閥門的工作特性.由于閥門的工作特性受管道系統(tǒng)和工作狀態(tài)的影響很大，因而必須針對(duì)具體系統(tǒng)對(duì)閥門的特性進(jìn)行具體分析.

　　通常，對(duì)于在瞬態(tài)操作條件下工作的閥門，其設(shè)計(jì)也是利用了穩(wěn)態(tài)的結(jié)果，采取通常的設(shè)計(jì)方法，且在對(duì)泵系統(tǒng)進(jìn)行計(jì)算時(shí)，把閥門的動(dòng)態(tài)模型處理為靜態(tài)模型進(jìn)行計(jì)算，而瞬態(tài)性能與穩(wěn)態(tài)性能有很大的區(qū)別，因此很有必要研究閥門調(diào)節(jié)過程中的工作特性。但是，目前對(duì)閥門瞬態(tài)工作特性的研究工作，大多針對(duì)泵系統(tǒng)進(jìn)行水錘分析，閥門只是作為其中的一個(gè)元件;對(duì)于閥門內(nèi)部流態(tài)和外部工作特性的研究工作，由于非定常計(jì)算的難度和工作量較大，僅在試驗(yàn)方面有一定的進(jìn)展，而在數(shù)值模擬方面至今未見如何解決閥門調(diào)節(jié)過程中精確定義其邊界問題的方法。

　　文中應(yīng)用外特性仿真分析以及內(nèi)部流態(tài)分析軟件，針對(duì)閥門瞬態(tài)工作過程的調(diào)節(jié)特性進(jìn)行研究，并對(duì)閥門的瞬態(tài)和穩(wěn)態(tài)阻力系數(shù)進(jìn)行對(duì)比分析，為分析閥門的動(dòng)態(tài)特性提供一定的參考依據(jù)。

1、仿真分析模型

　　針對(duì)典型的泵裝置系統(tǒng)，建立仿真模型如圖1所示，供水裝置系統(tǒng)分為很多個(gè)計(jì)算模塊，其中進(jìn)出口壓力大小設(shè)為大氣壓力值，閥門的運(yùn)動(dòng)特性由閥門控制器來確定，圖中黑色圓點(diǎn)代表節(jié)點(diǎn)。

圖1 泵裝置系統(tǒng)仿真模型

2、閥門的調(diào)節(jié)特性

　　2.1、閥門的流量特性

　　閥門的流量特性通常以相對(duì)流量NQ(某一開度時(shí)的流量與全開時(shí)流量之比)與相對(duì)開度φ(某一開度時(shí)閥桿行程或轉(zhuǎn)角與全開閥桿行程或轉(zhuǎn)角之比)的關(guān)系來表示，并以流量系數(shù)KQ(相對(duì)流量與閥門的最大通過能力的乘積)來衡量。

　　通常，閥門說明書上提供的流量系數(shù)是以清水為介質(zhì)，閥門前后的壓差為0.1MPa，流體的密度為1000kg/m3。這種規(guī)定條件下的流量特性稱為閥門的固有流量特性。常見的調(diào)節(jié)閥固有流量特性有快開、直線和等百分比特性等3種.如圖2所示，曲線1，2，3分別是理想的快開特性、直線特性和等百分比特性閥門的特性曲線。但是在對(duì)閥門進(jìn)行調(diào)節(jié)的過程中，由于閥門開度的改變，會(huì)對(duì)水力產(chǎn)生一定的擾動(dòng)作用，并對(duì)水擊有一定的反射作用，這稱為閥門的動(dòng)態(tài)特性，它與管路系統(tǒng)和調(diào)閥規(guī)律都有關(guān)系。下面針對(duì)閥門調(diào)節(jié)所引起的管路負(fù)載快速變化的過程，對(duì)泵裝置系統(tǒng)的水力變化進(jìn)行分析，以研究閥門調(diào)節(jié)過程中，泵裝置系統(tǒng)特性以及閥門的動(dòng)態(tài)特性。

圖2 閥門的流量特性

　　2.2、計(jì)算實(shí)例

　　水平輸水管道計(jì)算模型見圖1。兩種特性(直線特性和等百分比特性)的閥門安裝在距離心泵600m處。設(shè)閥門調(diào)節(jié)時(shí)間為10s。閥門在水力改變過程中的流量特性和閥前壓力計(jì)算結(jié)果見圖3-5。

圖3 直線閥固有特性與動(dòng)態(tài)特性

圖4 等百分比閥固有特性與動(dòng)態(tài)特性

圖5 閥前壓力

　　計(jì)算過程為閥門開啟的過程，由計(jì)算結(jié)果可知，即便安裝在短距離管線上，直線特性和等百分比特性的調(diào)節(jié)閥都具有快開特性，亦即閥在其開啟的初始行程和關(guān)閉的末段行程中，其流量變化對(duì)閥門的相對(duì)開度相當(dāng)敏感，因此閥門在此時(shí)的調(diào)節(jié)作用很明顯。分析其原因：當(dāng)管路系統(tǒng)中的管道長(zhǎng)度足夠長(zhǎng)時(shí)，管路的摩阻遠(yuǎn)遠(yuǎn)大于閥門的阻力，可忽略不計(jì)，于是各種不同特性閥門的動(dòng)態(tài)特性都向快開特性靠近，且管道越長(zhǎng)快開特性越明顯，而當(dāng)管道長(zhǎng)度較小時(shí)，閥門的動(dòng)態(tài)特性和靜態(tài)特性較接近。如圖5所示，在閥開啟行程的前10%～20%，閥門的水擊壓力下降特別迅速，減壓幅值也很大.由于壓力的瞬變會(huì)帶來一系列的問題，因此在泵裝置系統(tǒng)中，應(yīng)注意閥門的正確操作。

　　由圖3，4還可知，閥門在t=10s時(shí)即調(diào)節(jié)結(jié)束之后，流量的變化會(huì)滯后一定的時(shí)間，兩種調(diào)節(jié)方案均在t=15s左右達(dá)到穩(wěn)定，其原因是由于泵的機(jī)組慣量所帶來的滯后性，且隨著機(jī)組慣量的增大，其滯后性會(huì)迅速增大，這一點(diǎn)在泵裝置系統(tǒng)中具有普遍性。

3、局部流動(dòng)模擬分析

　　3.1、閥門模型建立及網(wǎng)格劃分

　　雖然動(dòng)網(wǎng)格方法可以實(shí)現(xiàn)閥門的啟閉過程，但是由于閥門運(yùn)動(dòng)速度較快時(shí)，網(wǎng)格的更新速度加快，對(duì)網(wǎng)格的質(zhì)量要求較高，從而增加了計(jì)算時(shí)間。圖6給出了實(shí)際計(jì)算過程中閥芯網(wǎng)格更新的過程。由圖可知，在閥門開啟過程中，網(wǎng)格不斷地被拉伸更新，網(wǎng)格密度變得逐漸稀疏，這對(duì)于高精度的數(shù)值模擬造成了很大的障礙。因此，如何保證網(wǎng)格質(zhì)量且減少計(jì)算量，成為一個(gè)需要解決的關(guān)鍵問題。

圖6 網(wǎng)格變形圖

　　為解決這個(gè)問題，可參考區(qū)域動(dòng)態(tài)滑移法，將閥門分為幾個(gè)計(jì)算區(qū)域，網(wǎng)格拉伸區(qū)和網(wǎng)格靜止區(qū)域。圖7為簡(jiǎn)化的閘閥三維模型，動(dòng)靜區(qū)域間采用滑移面連接的方法。在程序計(jì)算的過程中，單獨(dú)對(duì)網(wǎng)格更新區(qū)域采用含有移動(dòng)邊界的N-S方程離散，而靜止區(qū)域采用N-S方程離散。

圖7 計(jì)算模型

　　將計(jì)算區(qū)域分區(qū)后，閥芯運(yùn)動(dòng)區(qū)域單獨(dú)進(jìn)行網(wǎng)格更新，這在計(jì)算中可以保證靜止區(qū)域原有的網(wǎng)格分布質(zhì)量，故大大提高了計(jì)算效率和精度。

　　對(duì)簡(jiǎn)化模型進(jìn)行網(wǎng)格劃分，閥芯在計(jì)算過程中，網(wǎng)格不斷地被拉伸和壓縮，因此，對(duì)閥芯部分加密，劃分結(jié)果如圖8所示。

圖8 計(jì)算網(wǎng)格

　　3.2、邊界條件

　　閥芯的運(yùn)動(dòng)采用動(dòng)網(wǎng)格技術(shù)實(shí)現(xiàn)，邊界新的位置由Fluent自動(dòng)執(zhí)行更新。在使用動(dòng)網(wǎng)格時(shí)，只需給定初始網(wǎng)格和使用用戶自定義函數(shù)定義運(yùn)動(dòng)區(qū)域的運(yùn)動(dòng)即可。

　　其邊界定義：采用速度進(jìn)口和自由出流邊界條件，在壁面處采用無滑移邊界條件，近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù).閥芯的運(yùn)動(dòng)方式采用動(dòng)邊界文件(udf)來控制。閥芯采用直線的啟閉規(guī)律，其運(yùn)動(dòng)速度為0.025m/s，運(yùn)動(dòng)時(shí)間為2s;計(jì)算時(shí)間步長(zhǎng)為0.001s。

　　由于在閥門調(diào)節(jié)過程中，進(jìn)口速度處于變化的過程，因此，首先通過對(duì)泵系統(tǒng)裝置在Flowmaster中仿真分析，對(duì)計(jì)算得出的數(shù)據(jù)進(jìn)行擬合，得出的方程作為計(jì)算邊界，閥門進(jìn)口速度曲線見圖9。

圖9 進(jìn)口速度曲線

　　對(duì)此曲線進(jìn)行二次擬合，則進(jìn)口速度表達(dá)式為

υ=-1.1133t2+3.9283t+0.0423， (1)

　　其相關(guān)系數(shù)為R2=0.998。

　　3.3、動(dòng)態(tài)計(jì)算結(jié)果討論

　　閥芯運(yùn)動(dòng)過程在不同瞬時(shí)的速度和壓力場(chǎng)分布見圖10。

圖10 不同時(shí)刻速度和壓力分布圖

　　由于閥門本身的結(jié)構(gòu)特點(diǎn)，進(jìn)口段的通流面積較大，而經(jīng)過閥芯時(shí)通流面積突然縮小.這種結(jié)構(gòu)特點(diǎn)對(duì)流動(dòng)的影響，從靜壓分布上表現(xiàn)為進(jìn)口段靜壓力相對(duì)較大，表明速度較小;出口端靜壓力較小，表明速度較大。而在閥芯位置處，截面的急劇縮小使流速迅速增大，由圖也可知，該位置的靜壓力值在整個(gè)流道內(nèi)是最低的。在出口段靠近閥芯處，有一個(gè)較大范圍的旋渦區(qū)，這是由于流體流向的突然改變產(chǎn)生的，會(huì)引起較大的能量損失。

　　隨著閥門開度的不斷增大，由于其通流面積不斷地增大，閥芯附近的壓力變化范圍逐漸增大，而閥后的負(fù)壓區(qū)域逐漸減小。從速度矢量圖上也能反映這個(gè)變化，t=1.40s時(shí)流動(dòng)的旋渦區(qū)相對(duì)于t=0.05s時(shí)的減小了，但旋渦伴在閥門開啟的整個(gè)過程中都是存在的。

　　3.4、瞬態(tài)和穩(wěn)態(tài)阻力系數(shù)的比較

　　推得通過閥門的流量和損失系數(shù)的關(guān)系為

　　式中：Δh為壓力水頭，m;g為重力加速度，m/s²;A為管道截面積，m2;ξ為閥門的損失系數(shù).

　　由式(2)可知，當(dāng)閥門開度變化時(shí)，損失系數(shù)的變化導(dǎo)致了通過閥門流量的變化;這個(gè)流量的變化又對(duì)管路產(chǎn)生壓力變化。所以，確定一個(gè)正確的閥門損失系數(shù)與開度變化的關(guān)系，對(duì)于系統(tǒng)動(dòng)態(tài)模型的建立是很重要的。

　　在閥芯開啟過程以及固定的情況下，利用公式(3)計(jì)算出閥門的損失系數(shù)，計(jì)算數(shù)據(jù)見表1，表中φ為閥門相對(duì)開度。

表1 穩(wěn)態(tài)和瞬態(tài)阻力系數(shù)表

　　式中：Δp為被測(cè)閥門的壓力損失，即閥門的前后壓差;υav為流體在管道內(nèi)的平均流速;ρ為流體密度。

　　為了定量分析閥芯運(yùn)動(dòng)過程和穩(wěn)態(tài)過程的區(qū)別，設(shè)其阻力系數(shù)絕對(duì)誤差為

　　式中：ξs為瞬態(tài)計(jì)算阻力系數(shù);ξw為穩(wěn)態(tài)計(jì)算阻力系數(shù)。

　　將計(jì)算得到的數(shù)據(jù)以圖表的形式表示，見圖11。

圖11 瞬態(tài)和穩(wěn)態(tài)的阻力系數(shù)誤差圖

　　由表1可知，隨著閥門開度φ的增大，瞬態(tài)阻力系數(shù)迅速降低，當(dāng)閥門開度到50%時(shí)，穩(wěn)態(tài)和瞬態(tài)阻力系數(shù)值達(dá)到較小的值4左右。由閥門阻力系數(shù)的計(jì)算公式可知，此時(shí)閥芯所受到的阻力較小，即此時(shí)閥門的通流能力較強(qiáng);而在開度為4%時(shí)，穩(wěn)態(tài)和瞬態(tài)阻力系數(shù)值均較大，即閥芯所受到的阻力值較大。同時(shí)比較分析圖11的瞬態(tài)和穩(wěn)態(tài)阻力系數(shù)絕對(duì)誤差值可知，隨著閥門開度的增大，絕對(duì)誤差值逐漸降低，當(dāng)閥門開度達(dá)到50%時(shí)，穩(wěn)態(tài)和瞬態(tài)阻力系數(shù)值基本達(dá)到一致。

　　由以上數(shù)據(jù)分析可知，在分析開啟或關(guān)閉過程中閥門的阻力特性時(shí)，不能完全按照通常的穩(wěn)態(tài)理論進(jìn)行，尤其是對(duì)閥門開度較小的工況，如果閥門開度小于50%時(shí)，應(yīng)對(duì)其進(jìn)行一定程度的修正，以保證計(jì)算結(jié)果的正確性。

4、結(jié)論

　　該研究驗(yàn)證了動(dòng)網(wǎng)格技術(shù)可用于閥門開啟過程中分析，閥門開啟過程中的瞬態(tài)特性，對(duì)于研究閥門動(dòng)態(tài)特性、優(yōu)化閥門結(jié)構(gòu)與強(qiáng)度設(shè)計(jì)，可提供一定的參考依據(jù)。

　　1)在閥門快速調(diào)節(jié)的過程中，閥門動(dòng)態(tài)特性和靜態(tài)特性表現(xiàn)出較大的區(qū)別，且直線特性和等百分?jǐn)?shù)特性的調(diào)節(jié)閥都具有快開特性，即流量變化對(duì)閥門的相對(duì)開度相當(dāng)敏感，在閥門開度10%～20%范圍時(shí)，水擊壓力迅速下降。

　　2)比較分析靜態(tài)和動(dòng)態(tài)的閥門阻力系數(shù)，閥門開度大于50%時(shí)，靜態(tài)和動(dòng)態(tài)阻力系數(shù)值的差別較小;而在開度值小于50%時(shí)，動(dòng)、靜阻力系數(shù)值有較大的區(qū)別，且隨著閥門開度的降低而增大。因此，在建立閥門的動(dòng)態(tài)模型時(shí)應(yīng)當(dāng)加以考慮。

　　3)動(dòng)網(wǎng)格技術(shù)模擬閥門的開啟過程中，水流的瞬態(tài)變化情況，相對(duì)于靜態(tài)仿真，動(dòng)態(tài)仿真更能準(zhǔn)確地模擬出這個(gè)過程，且在閥門的開啟過程中，流場(chǎng)呈現(xiàn)復(fù)雜的非定常特性，閥芯后面在開度較小時(shí)出現(xiàn)很強(qiáng)的旋渦。

　　4)結(jié)合外特性分析軟件Flowmaster和局部流動(dòng)分析軟件Fluent對(duì)閥門的動(dòng)態(tài)特性進(jìn)行了研究，解決了閥門調(diào)節(jié)過程中復(fù)雜的邊界條件的定義問題。對(duì)于在瞬態(tài)工況下運(yùn)行的閥門，其模型建立以及水力性能預(yù)測(cè)，都不能完全按照通常的穩(wěn)態(tài)理論進(jìn)行。