彭利坤， 屈 鐸， 許文奇， 陳 佳

(海軍工程大學(xué) 動(dòng)力工程學(xué)院，武漢 430033)

大型艦船、潛艇上常設(shè)有縱傾均衡調(diào)水系統(tǒng)以調(diào)節(jié)艦船首尾平衡，其管系中的電液球閥啟閉規(guī)律無(wú)法進(jìn)行精確控制，而輸水過(guò)程中球閥突然關(guān)閉時(shí)會(huì)產(chǎn)生劇烈的水擊現(xiàn)象[1]，不僅導(dǎo)致管系的振動(dòng)噪聲，甚至可能造成管路破裂、閥門(mén)等設(shè)備損壞。

針對(duì)突然關(guān)閥、停泵等引起的流體管路水擊問(wèn)題，王福軍等采用基于特征線法的熱流體仿真平臺(tái)Flowmaster對(duì)某泵站進(jìn)行水力過(guò)渡過(guò)程計(jì)算，并對(duì)泵后閥門(mén)關(guān)閉規(guī)律和空氣閥布置方案進(jìn)行優(yōu)化[2]；由于特征線法難以考慮管道空間效應(yīng)等局限性[3]，CFD(Computational Fluid Dynamics)逐漸被應(yīng)用于水擊壓強(qiáng)計(jì)算[4-5]。Nikpour等[6]采用標(biāo)準(zhǔn)k-ε湍流模型模擬水擊現(xiàn)象，得到了與實(shí)驗(yàn)數(shù)據(jù)基本一致的結(jié)果。張飛等[7]采用短時(shí)傅里葉變換方法，研究了水泵水輪發(fā)電機(jī)組球閥動(dòng)水關(guān)閉前、后及過(guò)程中機(jī)組及球閥的主要頻率成分變化情況。對(duì)于典型的球閥控制流體管路系統(tǒng)，因艦船空間狹小，單相調(diào)壓水箱、蓄能器等傳統(tǒng)的水擊防護(hù)措施不易在艦船上實(shí)施，蔡標(biāo)華[8]采用AMEsim對(duì)艦船首尾移水系統(tǒng)進(jìn)行仿真與實(shí)驗(yàn)分析，認(rèn)為延長(zhǎng)球閥關(guān)閉時(shí)長(zhǎng)能有效抑制水擊；郭蘭蘭等[9-10]采用FLUENT滑移網(wǎng)格技術(shù)對(duì)球閥關(guān)閉水擊壓強(qiáng)進(jìn)行計(jì)算，認(rèn)為相比勻速、勻加速關(guān)閥方式，無(wú)閥腔球閥模型的勻減速關(guān)閥方式對(duì)水擊現(xiàn)象的抑制效果較好。在液壓泵站過(guò)渡過(guò)程中常采用兩階段關(guān)閉液壓閥進(jìn)行水擊防護(hù)，并取得了較好的效果[11]。然而，對(duì)于球閥控制流體管路的水擊防護(hù)，僅有“快開(kāi)(閥)慢關(guān)(閥)”、“先快后慢關(guān)閥”[12]等定性分析，尤其是兩階段線性關(guān)閥具體控制規(guī)律的水擊防護(hù)研究尚未見(jiàn)報(bào)道。因此，以艦船均衡調(diào)水系統(tǒng)為原型，搭建了具有雙向調(diào)水功能的實(shí)驗(yàn)系統(tǒng)，通過(guò)CFD數(shù)值模擬與實(shí)驗(yàn)研究相結(jié)合的方法，重點(diǎn)分析了先快后慢的兩階段線性關(guān)閥的水擊防護(hù)效果，對(duì)球閥關(guān)閉規(guī)律進(jìn)行優(yōu)化，以期為球閥關(guān)閉時(shí)低噪聲控制提供指導(dǎo)。

1 CFD水擊壓強(qiáng)計(jì)算

1.1 實(shí)驗(yàn)系統(tǒng)簡(jiǎn)介

實(shí)驗(yàn)系統(tǒng)原理如圖1所示，通過(guò)控制空氣瓶?jī)?nèi)壓縮空氣的壓力來(lái)控制調(diào)水的速度；通過(guò)控制空氣瓶?jī)?nèi)空氣的輸送方向?qū)崿F(xiàn)雙向調(diào)水；其中管路1、管路2分別安裝數(shù)字步進(jìn)球閥Va(最快啟閉時(shí)間3 s，即啟閉最大速度30 °/s)、電液球閥Vb(啟閉時(shí)間約0.5 s)。根據(jù)相似性原理和實(shí)驗(yàn)室具體情況搭建的實(shí)驗(yàn)系統(tǒng)臺(tái)架如圖2所示。

圖1 實(shí)驗(yàn)系統(tǒng)原理示意圖Fig.1 Schematic of experimental system

圖2 實(shí)驗(yàn)系統(tǒng)整體實(shí)物圖Fig.2 Physical map of experimental system

實(shí)驗(yàn)系統(tǒng)中的數(shù)字球閥和電液球閥僅驅(qū)動(dòng)裝置不同，與流場(chǎng)接觸的內(nèi)部結(jié)構(gòu)基本相同，如圖3所示，閥芯外壁與閥體內(nèi)壁構(gòu)成了閥腔。

1.2 數(shù)值模擬方法

由球閥閥芯球體、閥腔、閥芯通道的幾何關(guān)系，構(gòu)建球閥一定開(kāi)度下時(shí)的剖面與閥口投影示意圖如圖4所示，推導(dǎo)得到球閥流通面積A與開(kāi)度θ的關(guān)系

(1)

1-左閥體；2-閥芯；3-右閥體；4-閥腔；5-密封圈；6-進(jìn)口管道；7-閥芯通道；8-出口管道

(a)

(b)圖4 球閥剖面與閥口投影示意圖Fig.4 Cutaway view and projection drawing of valve

根據(jù)實(shí)驗(yàn)系統(tǒng)實(shí)際管長(zhǎng)，取閥前管長(zhǎng)l1、閥后管長(zhǎng)l2建立球閥全開(kāi)時(shí)的流場(chǎng)幾何模型，同時(shí)針對(duì)球閥結(jié)構(gòu)與運(yùn)動(dòng)特點(diǎn)，將整個(gè)流域分為進(jìn)口管路、閥芯通道、閥腔和出口管路四部分建模。采用結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行網(wǎng)格劃分。

首先對(duì)網(wǎng)格模型進(jìn)行了網(wǎng)格無(wú)關(guān)性模擬。針對(duì)0.5 s時(shí)長(zhǎng)單線性勻速關(guān)閥工況，在20萬(wàn)網(wǎng)格的基礎(chǔ)上，每次增加2萬(wàn)網(wǎng)格對(duì)其進(jìn)行仿真，結(jié)果發(fā)現(xiàn)，當(dāng)網(wǎng)格增加到約28萬(wàn)時(shí)，其最大水擊壓強(qiáng)的計(jì)算結(jié)果與網(wǎng)格約為26.2萬(wàn)的計(jì)算結(jié)果相差僅為0.3%，所以采用約26.2萬(wàn)的網(wǎng)格模型。該模型網(wǎng)格質(zhì)量在0.3以上，并對(duì)運(yùn)動(dòng)邊界進(jìn)行了適當(dāng)加密，其網(wǎng)格如圖5所示。

根據(jù)調(diào)水系統(tǒng)的工況，調(diào)水時(shí)兩端水箱的壓差為4 bar，其中一端與空氣瓶連通，而另一端連通大氣，因此，設(shè)置參考?jí)毫? bar，進(jìn)口設(shè)置為總壓4 bar的壓力進(jìn)口；出口設(shè)置為0 bar的壓力出口。采用FLUENT滑移網(wǎng)格技術(shù)和標(biāo)準(zhǔn)k-ε湍流模型模擬球閥關(guān)閉瞬態(tài)過(guò)程，同時(shí)定義交界面以便各流域進(jìn)行數(shù)據(jù)交換，并通過(guò)UDF(User Defined Function)實(shí)現(xiàn)對(duì)閥芯運(yùn)動(dòng)規(guī)律的精確控制。

(a)閥芯通道及進(jìn)出口管路

(b)閥腔圖5 幾何模型網(wǎng)格Fig.5 Mesh of geometrical model

球閥由全開(kāi)到完全閉合是一個(gè)動(dòng)態(tài)過(guò)程，管道中流體為非穩(wěn)態(tài)湍流流動(dòng)，其連續(xù)性方程和動(dòng)量方程為

(2)

(3)

CFD計(jì)算時(shí)采用標(biāo)準(zhǔn)k-ε湍流模型來(lái)使上述方程封閉。

1.3 數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

根據(jù)實(shí)驗(yàn)系統(tǒng)特性對(duì)模型各參數(shù)選取如表1所示。

表1 模型參數(shù)Tab.1 Parameters of model

對(duì)球閥關(guān)閉過(guò)程進(jìn)行瞬態(tài)數(shù)值模擬，監(jiān)測(cè)閥前截面的平均壓力，也即水擊壓強(qiáng)。研究水擊問(wèn)題需考慮水的可壓縮性，根據(jù)流體體積模量表達(dá)式(4)[13]、流體密度ρ的表達(dá)式(5)及水擊波波速c的表達(dá)式(6)，在Fluent中采用UDF定義三者之間的聯(lián)系以考慮水的可壓縮性。

(4)

(5)

(6)

式中:K為水的體積模量;V為水的體積；p為壓強(qiáng)。

以0.5 s時(shí)長(zhǎng)勻速關(guān)閉電液球閥工況的水擊壓強(qiáng)數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比如圖6所示。

圖6 水擊壓強(qiáng)計(jì)算值與電液球閥實(shí)驗(yàn)值對(duì)比Fig.6 Water hammer pressure of experiment and simulation

從圖6中可知，實(shí)驗(yàn)得到的最大水擊壓強(qiáng)為7.31 bar，計(jì)算所得的最大水擊壓強(qiáng)為7.722 bar，就最大水擊壓強(qiáng)而言，計(jì)算值與實(shí)驗(yàn)值誤差為5.64%。仿真采用單線性勻速關(guān)閥；實(shí)驗(yàn)中電液球閥動(dòng)作是電磁閥打開(kāi)、液壓推動(dòng)閥芯旋轉(zhuǎn)關(guān)閉的，加速階段很短，可近似看作勻速關(guān)閥。實(shí)驗(yàn)和數(shù)值模擬中關(guān)閥規(guī)律上存在的這點(diǎn)不同是造成誤差的主要原因之一。數(shù)值模擬的水擊壓強(qiáng)在閥門(mén)完全關(guān)閉后先上下波動(dòng)，最后穩(wěn)定在4 bar，這與實(shí)驗(yàn)結(jié)果吻合較好，說(shuō)明考慮水的可壓縮性進(jìn)行關(guān)閥水擊現(xiàn)象的數(shù)值模擬是可行的。

1.4 CFD計(jì)算結(jié)果分析

對(duì)于實(shí)驗(yàn)系統(tǒng)中的球閥，由式(1)可知，當(dāng)閥芯由全開(kāi)旋轉(zhuǎn)80°時(shí)，球閥恰好完全關(guān)閉。對(duì)實(shí)驗(yàn)系統(tǒng)球閥管路模型在有效關(guān)閥時(shí)長(zhǎng)內(nèi)采用先快后慢的兩階段線性關(guān)閥進(jìn)行CFD數(shù)值仿真，且本節(jié)中數(shù)值模擬的所有工況中整個(gè)關(guān)閥(球閥閥芯從全開(kāi)旋轉(zhuǎn)90°)時(shí)長(zhǎng)均采用0.5 s。其中，“以270 °/s速度快關(guān)43.75%”工況的球閥模型縱截面各時(shí)刻的壓力云圖和速度云圖分別如圖7、圖8所示。從圖7和圖8中可以看出，球閥閥芯在由全開(kāi)到關(guān)閉的過(guò)程中，由于閥芯的節(jié)流作用，流速發(fā)生了劇烈的變化，造成了嚴(yán)重的壓力損失，而且進(jìn)口管道和出口管道的壓差會(huì)隨著球閥的關(guān)閉而增大。

(a)0.4 s

(b)0.5 s

(c)0.6 s

(d)0.7 s圖7 各時(shí)刻壓力云圖Fig.7 Pressure contour at different time

(a)0.4 s

(b)0.5 s

(c)0.6 s

(d)0.7 s圖8 各時(shí)刻速度云圖Fig.8 Velocity contour at different time

圖9 快關(guān)階段角度不同的關(guān)閥規(guī)律Fig.9 Closing law with different angle at quick-closing stage

圖10 快關(guān)階段角度不同時(shí)水擊壓強(qiáng)曲線Fig.10 Water hammer pressure of different closing law with different angle at quick-closing stage

表2 快關(guān)階段角度不同時(shí)的最大水擊壓強(qiáng)Tab.2 Maximum pressure of different closing law with different angle at quick-closing stage

快關(guān)階段速度均為270 °/s、快關(guān)階段角度不同的關(guān)閥規(guī)律如圖9所示，其對(duì)應(yīng)的水擊壓強(qiáng)曲線如圖10所示，其最大水擊壓強(qiáng)如表2所示。從圖10和表2可知：

(1)采用先快后慢的兩階段線性關(guān)閥時(shí)，閥門(mén)完全關(guān)閉后的壓力振蕩較勻速關(guān)閥時(shí)有了明顯改善;

(2)較線性關(guān)閥而言，采用先快后慢的兩階段線性規(guī)律進(jìn)行關(guān)閥時(shí)，水擊壓強(qiáng)峰值將提前出現(xiàn)，當(dāng)快關(guān)階段角度太大時(shí)，反而導(dǎo)致水擊壓強(qiáng)變大，其中“以270 °/s速度快關(guān)全角度56.25%”的最大水擊壓強(qiáng)達(dá)到9.593 bar；當(dāng)快關(guān)階段角度較小時(shí)，則對(duì)水擊的抑制效果不佳，其中“以270 °/s速度快關(guān)31.25%”的最大水擊壓強(qiáng)為7.083 bar;

(3)就本模型而言，關(guān)閥時(shí)長(zhǎng)均為0.5 s時(shí)采用“以270 °/s速度快關(guān)43.75%”的關(guān)閥方式的水擊防護(hù)效果較好，其最大水擊壓強(qiáng)僅為6.884 bar。

在上述結(jié)果(3)的基礎(chǔ)上研究快關(guān)角度均為43.75%時(shí)快關(guān)階段速度不同工況下的水擊壓強(qiáng)，具體關(guān)閥規(guī)律如圖11所示，其對(duì)應(yīng)的關(guān)閥水擊壓強(qiáng)如圖12所示，不同關(guān)閥規(guī)律下的最大水擊壓強(qiáng)如表3所示。

表3 快關(guān)階段速度不同時(shí)的最大水擊壓強(qiáng)Tab.3 Maximum pressure of different closing law with different velocity at quick-closing stage

從圖11、圖12和表3中可知：快關(guān)階段速度越快，水擊壓強(qiáng)最大值出現(xiàn)時(shí)刻越早；當(dāng)快關(guān)角度均為全行程的43.75%時(shí)，最大水擊壓強(qiáng)隨快關(guān)階段速度增加而減小，且其衰減程度變緩。實(shí)際應(yīng)用中關(guān)閥速度不可能無(wú)限增大。

2 實(shí)驗(yàn)結(jié)果分析

軟件設(shè)置了多種兩階段關(guān)閥控制規(guī)律，其中關(guān)閥時(shí)長(zhǎng)均為4 s的兩階段關(guān)閥規(guī)律如圖13所示，CFD仿真計(jì)算的關(guān)閥水擊壓強(qiáng)曲線如圖14所示，實(shí)驗(yàn)值曲線如圖15所示。

圖11 快關(guān)階段速度不同的關(guān)閥規(guī)律Fig.11 Closing law with different velocity at quick-closing stage

圖12 快關(guān)速度不同時(shí)水擊壓強(qiáng)曲線Fig.12 Water hammer pressure of different closing law with different velocity at quick-closing stage

圖13 兩階段關(guān)閥規(guī)律Fig.13 Valve-closing law of two stage

圖14 不同關(guān)閥規(guī)律水擊壓強(qiáng)仿真值Fig.14 Water hammer pressure of different valve-closing law of simulation

圖15 不同兩階段關(guān)閥規(guī)律的水擊壓強(qiáng)實(shí)驗(yàn)值Fig.15 Water hammer pressure of different valve-closing law of experiment

從圖14、圖15中可看出，采用先快后慢的兩階段線性關(guān)閥規(guī)律時(shí)，水擊壓強(qiáng)的峰值出現(xiàn)的時(shí)間提前，且閥門(mén)完全關(guān)閉后的壓力振蕩明顯減小；但如果快關(guān)階段角度太大反而會(huì)導(dǎo)致最大水擊壓強(qiáng)增加，而快關(guān)階段角度太小則造成水擊壓強(qiáng)抑制效果不佳，以上四種關(guān)閥規(guī)律中以快慢速比1.33且快關(guān)50%時(shí)的水擊抑制效果最好，這與CFD仿真結(jié)論基本一致。

同時(shí)將不同關(guān)閥規(guī)律最大水擊壓強(qiáng)的實(shí)測(cè)數(shù)據(jù)與仿真值對(duì)比如表4所示。從中可看出，就最大水擊壓強(qiáng)而言，仿真值與實(shí)驗(yàn)值誤差均在6%以內(nèi)，考慮到調(diào)水壓力難以精確控制等因素，誤差在允許范圍內(nèi)。

表4 最大水擊壓強(qiáng)仿真與實(shí)驗(yàn)值對(duì)比Tab.4 Maximum pressure of experiment and simulation

3 結(jié) 論

(1)本文針對(duì)艦船縱傾均衡系統(tǒng)中球閥控制流體管路在突然關(guān)閥時(shí)產(chǎn)生劇烈水擊問(wèn)題，設(shè)計(jì)搭建了包含開(kāi)關(guān)規(guī)律可精確控制數(shù)字球閥的實(shí)驗(yàn)系統(tǒng)，通過(guò)對(duì)比實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了CFD應(yīng)用于水擊計(jì)算的可行性。

(2)通過(guò)分析CFD計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)可知：當(dāng)關(guān)閥時(shí)長(zhǎng)一定時(shí)，采用先快后慢的兩階段線性關(guān)閥能有效抑制球閥關(guān)閉后的壓力振蕩情況，且水擊壓強(qiáng)峰值出現(xiàn)時(shí)刻提前，但是快關(guān)階段角度太大時(shí)，反而會(huì)增大關(guān)閥過(guò)程水擊壓強(qiáng)；而當(dāng)快關(guān)階段角度較小時(shí)，則造成水擊壓強(qiáng)抑制效果不佳，就本研究而言，以快慢速比1.6左右且快關(guān)角度在45%附近時(shí)水擊防護(hù)效果較好。