大口徑中心型蝶閥在實(shí)際應(yīng)用中由于需承受較大水壓作用，其閥板上往往利用筋板增加強(qiáng)度，從而減少閥板變形。利用流體動力學(xué)分析軟件，以RANS 方程為控制方程，采用標(biāo)準(zhǔn)的紊流模型(k-ε 模型)和壓力耦合方程組的半隱式方法(SIMPLE 算法)，研究在100%開度下閥板表面三種典型筋板分布形式(“﹟”形、“◇“形和“//”形)對流經(jīng)閥門流體流動特性的影響規(guī)律。數(shù)值計(jì)算結(jié)果顯示，在滿足閥板強(qiáng)度要求前提下，采用對稱性“//“形筋板分布有利于減少閥門流體阻力和動力水頭損失，增加流量。

　　蝶閥是工業(yè)和農(nóng)業(yè)領(lǐng)域常見的閥門，主要起切斷和節(jié)流作用。蝶閥啟閉件是一個(gè)圓盤形的閥板，在閥體內(nèi)繞自身軸線旋轉(zhuǎn)，從而達(dá)到啟閉或調(diào)節(jié)的目的。近年來，隨著計(jì)算機(jī)技術(shù)及計(jì)算流體動力學(xué)軟件的發(fā)展，越來越多的人利用數(shù)值計(jì)算的方法分析蝶閥內(nèi)部三維流動特性，如壓力流速分布、分離流動區(qū)域等。袁新明等數(shù)值模擬研究了閥門的阻力特性;諸葛偉林等采用基于非結(jié)構(gòu)、非交錯(cuò)網(wǎng)格的有限體積法求解用兩方程模型封閉的雷諾平均N-S方程組，對蝶閥的三維分離流動進(jìn)行了數(shù)值模擬，得出蝶閥的流動阻力系數(shù)隨著蝶閥關(guān)閉角度的增大呈指數(shù)性增長的結(jié)論;沈新榮等針對自行開發(fā)的一類安裝配流板的新型電動蝶閥模型，進(jìn)行了三維數(shù)值模擬和實(shí)驗(yàn)研究，對不同開度下電動蝶閥的三維湍流流動進(jìn)行了數(shù)值計(jì)算;劉健等對大口徑蝶閥運(yùn)用商用流體計(jì)算軟件FLUENT，對其不同開度情況下的流場形式進(jìn)行了三維數(shù)值模擬分析;劉華坪等利用動網(wǎng)格技術(shù)對管路中常見的閥門進(jìn)行了動態(tài)數(shù)值模擬;HuangandKim使用三維數(shù)值模擬技術(shù)分析了在蝶閥內(nèi)部的不可壓縮的流體的流動狀況，并給出了分析后的速度場圖和應(yīng)力分布圖;LinandSchohl分析了CFD在蝶閥領(lǐng)域的應(yīng)用;黃國權(quán)等利用CFD軟件FLUENT對中心型蝶閥流場進(jìn)行數(shù)值模擬，定性給出了閥門在不同開度下重要部位的受力情況，比較直觀地給出了閥門在不同工況下流道內(nèi)部的速度分布，得到了渦流的形成過程，及速度對渦流形成及擴(kuò)展的影響;類似的研究工作還有許多。這些研究對于指導(dǎo)蝶閥的設(shè)計(jì)、改善其流動狀況、減小流動阻力具有非常重要的意義。

　　對于大口徑蝶閥，由于其在應(yīng)用時(shí)需承受較大的流體壓力，其閥板容易產(chǎn)生變形而降低工作可靠性。因此，往往需要設(shè)計(jì)加強(qiáng)筋增加閥板強(qiáng)度。目前，加強(qiáng)筋的分布形式主要有“#”形、“◇”形和“//”形三種(圖1)。

圖1 蝶閥閥板加強(qiáng)筋的三種典型結(jié)構(gòu)

　　工程實(shí)際應(yīng)用過程中發(fā)現(xiàn)，不同形狀加強(qiáng)筋對流經(jīng)閥門流體流動特性的影響不一樣。為揭示其影響的內(nèi)在規(guī)律，利用流體動力學(xué)分析軟件，以RANS方程為控制方程，采用標(biāo)準(zhǔn)的紊流模型(k-ε模型)和壓力耦合方程組的半隱式方法(SIMPLE算法)，給出了在100%開度下閥板表面三種典型的筋板分布形式對閥門流動特性的影響規(guī)律。

1、數(shù)值計(jì)算條件設(shè)置

　　1.1、控制方程

　　選擇單相流體的不可壓縮三維粘性流動模型，采用不可壓縮流動的雷諾方程組與k-ε湍流模型構(gòu)成封閉的方程組來描述，如式(1)所示。

(1)

　　式中：ρ為密度; 為壓力;k為湍流動能;ε為湍流動能耗散率; 為速度矢量;xi和xj為坐標(biāo)參數(shù);v為湍流粘度;vt為湍流粘性系數(shù);Gk和Gε為由平均速度梯度引起的湍流動能及湍流動能耗散率的產(chǎn)生項(xiàng);cμ、σε、σk、c1、c2為常數(shù)，在標(biāo)準(zhǔn)k-ε模型中，cμ=0.09，σε=1.3，σk=1.0。

　　1.2、計(jì)算條件

　　采用不可壓縮流動的雷諾時(shí)均方程組，湍流模型采用標(biāo)準(zhǔn)k-ε模型;所有方程中的對流項(xiàng)均用二階迎風(fēng)格式離散，離散方程的求解采用壓力耦合方程組的半隱式方法(SIMPLE算法)。連續(xù)性方程和動量方程收斂殘差標(biāo)準(zhǔn)均為10-3。進(jìn)口邊界條件為速度進(jìn)口，進(jìn)口速度設(shè)1.0m/s。出口邊界條件定為自由流動，其他邊界條件均為固壁;計(jì)算中忽略重力對流場的影響;全流場計(jì)算定常流動，得到蝶閥內(nèi)流場的詳細(xì)分布情況。選取蝶閥及其前后一段管道作為計(jì)算域，流動方向?yàn)閤軸正向。采用非結(jié)構(gòu)網(wǎng)格，將蝶閥閥板及尾跡區(qū)域的網(wǎng)格局部加密，以確保網(wǎng)格質(zhì)量和求解精度。

2、結(jié)果分析

　　2.1、閥板無量綱流體阻力比較

　　在相同計(jì)算條件下，對閥板三種典型的加強(qiáng)筋形狀產(chǎn)生的無量綱流體阻力進(jìn)行比較，如圖2所示，數(shù)值為正說明產(chǎn)生的是阻力，數(shù)值越大則閥板對流體的阻力越大。從計(jì)算結(jié)果可以看出：三種加強(qiáng)筋導(dǎo)致閥板受到的流體阻力大小不同，阻力從大到小依次為“#”形、“◇”形和“//”形。進(jìn)一步從閥板表面壓力分布情況以及流體流經(jīng)閥板后流場的結(jié)構(gòu)和形態(tài)改變情況進(jìn)行深入分析，揭示導(dǎo)致差異的本質(zhì)原因。

圖2 相同條件下不同加強(qiáng)筋形狀產(chǎn)生的無量綱流體阻力

　　2.2、壓力分布情況比較

　　圖3為三種典型加強(qiáng)筋結(jié)構(gòu)形成的閥板表面三維壓力分布云圖和截面壓力分布云圖。

圖3 壓力分布情況

　　對于閥板而言，三種情況形成的高壓區(qū)基本上集中于閥板前部迎水面處，以中間逐漸向兩側(cè)減小，直至閥板背面形成低壓區(qū)，并且壓力達(dá)到最小。而對于閥板上的加強(qiáng)筋而言，壓力分布形態(tài)則有所不同：“//”形加強(qiáng)筋的高壓區(qū)集中于前端面，低壓區(qū)集中于后端面;“◇”形加強(qiáng)筋的高壓區(qū)集中于菱形的前角區(qū)域附近，而壓力最小值位于菱形上下兩側(cè)尖角后側(cè)區(qū)域附近;“#”形加強(qiáng)筋的高壓區(qū)較為分散，而低壓區(qū)主要分布于豎直筋條的后緣區(qū)域。加強(qiáng)筋上壓力分布形態(tài)的差異是造成閥門流體阻力不同的根本原因之一。

　　2.3、流體速度矢量及流線分布情況比較

　　從流體流經(jīng)三種加強(qiáng)筋形成的速度矢量及流線分布圖(圖4)可以看出：“//”形加強(qiáng)筋對流體的流動影響最小，流線進(jìn)過加強(qiáng)筋后流動方向基本上未發(fā)生明顯的改變，流動分離現(xiàn)象并不明顯，流阻較小;“#”形加強(qiáng)筋對流體的流動影響最大，流線經(jīng)過加強(qiáng)筋后被截?cái)�，并在加�?qiáng)筋內(nèi)部形成卡門漩渦，對流體阻礙作用最大。而流線經(jīng)過“◇”形加強(qiáng)筋后會在上下兩側(cè)尖角之后形成一定的分離，但是對流體的阻礙作用沒有“#”形加強(qiáng)筋的強(qiáng)烈。

圖4 流體速度矢量及流線分布情況

3、結(jié)論

　　本文通過數(shù)值模擬，研究了閥板在100%開度下其表面三種典型的筋板分布形式對流經(jīng)閥門流體的流動特性的影響規(guī)律。通過對無量綱流體阻力進(jìn)行比較以及表面壓力分布情況和流體速度矢量及流線分布情況等方面的分析。得出在滿足閥板強(qiáng)度要求的前提下，采用對稱性“//”形筋板分布有利于減少閥門的流體阻力和動力水頭損失并且增加流量的結(jié)論。