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

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

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

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

1、仿真分析模型

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

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

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

　　2.1、閥門的流量特性

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

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

圖2 閥門的流量特性

　　2.2、計算實例

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

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

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

圖5 閥前壓力

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

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

3、局部流動模擬分析

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

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

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

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

圖7 計算模型

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

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

圖8 計算網(wǎng)格

　　3.2、邊界條件

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

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

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

圖9 進口速度曲線

　　對此曲線進行二次擬合，則進口速度表達式為

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

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

　　3.3、動態(tài)計算結(jié)果討論

　　閥芯運動過程在不同瞬時的速度和壓力場分布見圖10。

圖10 不同時刻速度和壓力分布圖

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

4、結(jié)論

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

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

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

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

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