楊元龍，鄭 文

(1.中國艦船研究設計中心，武漢430064;2.武漢船舶設計研究所，武漢430064)

船用大壓降給水管道多級孔板設計及節(jié)流特性研究

楊元龍1，鄭 文2

(1.中國艦船研究設計中心，武漢430064;2.武漢船舶設計研究所，武漢430064)

為減少大壓降給水管道振動和噪聲,確保船舶汽水系統(tǒng)的穩(wěn)定循環(huán)特性,以某船用大壓降給水管道為原型,采用CFD數(shù)值模擬的方法分析六級同心孔板和六級偏心孔板的壓降特性、流場結(jié)構(gòu)及湍流脈動規(guī)律,計算結(jié)果表明,六級偏心節(jié)流孔板適用于大壓降短管距的船舶給水管路系統(tǒng)。數(shù)值計算的結(jié)果與理論設計吻合較好。

多級節(jié)流孔板;給水管道;汽蝕

節(jié)流孔板廣泛應用于船舶汽水系統(tǒng),作為船舶系統(tǒng)管網(wǎng)內(nèi)流體壓力和流量匹配的關鍵裝置,其主要運行原理是由于孔板限流產(chǎn)生局部阻力,促使孔板上下游產(chǎn)生壓降,進而達到節(jié)流降壓的設計目的。對于大壓降系統(tǒng)管路,由于節(jié)流孔板前后阻塞壓差的存在,布置單級節(jié)流孔板極易誘發(fā)流體汽蝕而導致管路發(fā)生劇烈振動和噪音［1-2］。對于具有給水循環(huán)流量大、溫度高、壓降大及管道距離短等運行和結(jié)構(gòu)特點的船用給水管路系統(tǒng),必須采用多級節(jié)流孔板來降低管路系統(tǒng)壓力［3］。

目前國內(nèi)外學者針對節(jié)流孔板結(jié)構(gòu)的理論設計和數(shù)值優(yōu)化做了大量分析研究［4-7］。但少見關于船舶大壓降給水管道節(jié)流孔板設計與分析的文獻報道。文中基于某大型陸上試驗測量數(shù)據(jù)作為計算邊界條件,根據(jù)多級節(jié)流孔板的設計方法和原則,對船用大壓降給水管道多級節(jié)流孔板進行結(jié)構(gòu)設計,采用CFD數(shù)值模擬的手段對六級同心孔板和六級偏心孔板的節(jié)流特征開展研究,并詳細分析了給水管路孔板局部流域的壓力、流速和湍動能等關鍵水力參數(shù)的變化特性及影響因素,討論適用于船舶大壓降給水管道的多級節(jié)流孔板結(jié)構(gòu)形式和布置方式。

1 數(shù)值方法敏感性驗證

以某船用給水管道單級節(jié)流孔板為原型,采用ICEM CFD軟件進行物理建模,基于O網(wǎng)格技術(shù)實現(xiàn)非結(jié)構(gòu)化網(wǎng)格劃分,利用CFX 12軟件對網(wǎng)格模型進行前處理、計算求解和后處理。同時,為減少離散解的產(chǎn)生,降低方程計算的偽擴散性,采用全隱式耦合算法迭代計算求解流速和壓力的耦合過程。

計算邊界參數(shù)設置為:進口壓力為4.5 MPa,給水溫度為104.5℃,出口壓力為0.13 MPa,管壁為絕熱條件和無滑移邊界條件,給水物性參數(shù)的計算標準均來自IAPWS IF97數(shù)據(jù)庫。

結(jié)合來自于試驗數(shù)據(jù)的計算邊界條件,開展單級節(jié)流孔板壓降特性分析,以驗證數(shù)值方法的敏感性和準確性。給水管道單級節(jié)流孔板前后流域壓力變化規(guī)律見圖1、2,在節(jié)流孔板上游區(qū)域壓力為4.5 MPa,給水流經(jīng)孔板后,在節(jié)流孔板的局部阻力作用下,促使流體壓力快速降低,導致節(jié)流孔板下游區(qū)域給水壓力約為0.1MPa,同時數(shù)值計算的壓力曲線與試驗測量數(shù)據(jù)基本吻合,驗證了CFD數(shù)值計算方法的可靠性和計算結(jié)果的準確性。

2 多級節(jié)流孔板設計

根據(jù)參考文獻［1］中的設計方法,各級節(jié)流孔板的壓降按照幾何級數(shù)遞減設計原則,即第一級節(jié)流孔板的壓降為Δp1,第二級節(jié)流孔板的壓降為Δp1/2,第三級節(jié)流孔板壓降為Δp1/22,…,第n級節(jié)流孔板的壓降為Δp1/2n-1,最后一級孔板的節(jié)流輸出壓力為設計目標壓力。

圖1 孔板壓力變化

圖2 孔板壓力分布

針對大壓降的船用給水管路,為防止節(jié)流孔板發(fā)生汽蝕現(xiàn)象而導致管道劇烈振動和噪聲,在節(jié)流孔板設計過程中,采用多級降壓的方式,確保各級孔板節(jié)流后壓力高于流體飽和壓力。根據(jù)參考文獻［1］和［8］的設計原則和方法,結(jié)合試驗測量數(shù)據(jù):給水流量為50 t/h,溫度為104.5℃,密度為955 kg/m3,進口壓力為4.5 MPa,背壓為0.13 MPa,對船用給水系統(tǒng)管道多級節(jié)流孔板進行結(jié)構(gòu)設計,計算參數(shù)詳見表1。

表1 節(jié)流孔板計算 mm

3 六級孔板節(jié)流特性分析

3.1 六級同心孔板節(jié)流特性分析

基于上述六級孔板結(jié)構(gòu)設計參數(shù),將1～6號孔板的間距設置為55,55,55,40,40 mm,進行六級孔板幾何建模,見圖3。

如圖4所示,給水管路及節(jié)流孔板采用結(jié)構(gòu)化網(wǎng)格方式進行網(wǎng)格劃分,管壁利用O網(wǎng)格方法添加邊界層網(wǎng)格,并進行局部網(wǎng)格細化?；诰W(wǎng)格無關解的計算分析,計算域共計230 000個網(wǎng)格單元。

圖3 六級節(jié)流孔板結(jié)構(gòu)

圖4 網(wǎng)格模型

六級同心節(jié)流孔板前后流線分布規(guī)律見圖5、6。由圖5、6可見,第一級孔板上游給水以較低流速流過孔板,孔板的節(jié)流作用導致第一級節(jié)流孔板下游流速急劇升高,形成流體射流現(xiàn)象。沿著流體流動方向,節(jié)流孔板開孔直徑依次增大,且流體流動未達到完全發(fā)展,導致流體射流程度減弱。

圖5 同心孔板流動分布規(guī)律

圖6 同心孔板流速分布

六級同心節(jié)流孔板前后流速矢量分布見圖7。由圖7可知,流體快速穿過六級節(jié)流孔板過程中誘發(fā)的射流效應,對各級孔板之間的流域產(chǎn)生卷吸作用,導致各級孔板之間流體形成了局部渦流區(qū)域,極易導致管道和孔板的低頻振動。

六級同心節(jié)流孔板湍動能分布規(guī)律見圖8。由圖8可知,管路上游的節(jié)流孔板區(qū)域流體流動湍動能較高,最大值達到380 m2/s2,增強湍流流體微團之間的動量傳遞,導致管路上游孔板區(qū)域強烈脈動。沿著流體流動的方向,流體射流效應減弱,流體流速不斷下降,流速梯度減小,促使流體湍流脈動強度逐漸降低。

圖7 同心孔板流速矢量分布規(guī)律

圖8 同心孔板湍動能分布

六級節(jié)流孔板壓力和流速的變化見圖9。

圖9 同心孔板壓力和流速分布曲線

從圖9可以看出,第一級節(jié)流孔板上游流體的流速為6 m/s,流經(jīng)第一級孔板的開孔時,流速快速升高至約120m/s,流體速度相應增大20倍。在第一級孔板下游流域,流體流速逐漸下降,直至給水管路出口。從圖9還可以發(fā)現(xiàn),由于第一級孔板區(qū)域的流速急劇增大,流體的靜壓能轉(zhuǎn)換為動壓能,導致在第一級節(jié)流孔板局部流域的流體壓力急劇降低。在第一級孔板下游流域,隨著各級孔板的開孔直徑增大,流體速度逐漸減小,促使流體動壓能轉(zhuǎn)換為靜壓能,流體壓力不斷升高,壓力變化見圖10。

但流體壓力不會完全上升至管路孔板上游流體壓力,主要是由于各級孔板下游形成大尺寸的流體漩渦(見圖7),渦流的能量耗散轉(zhuǎn)變?yōu)榱黧w的熱力學能,因此六級孔板節(jié)流后壓力不能上升至管路孔板的上游流體壓力。

圖10 同心孔板壓力分布云圖

從圖9可見,在第一級節(jié)流孔板與第六級孔板之間的流域形成了長度約為200 mm的局部負壓區(qū)域,其流體壓力低于給水溫度為104.5℃對應的飽和壓力0.118 8 MPa,導致流體發(fā)生汽蝕現(xiàn)象,進而引發(fā)流體的阻塞流動,從而易導致給水管路孔板的高頻振動。

3.2 六級偏心孔板節(jié)流特性分析

為優(yōu)化大壓降小管距的給水管路多級孔板的節(jié)流特性,結(jié)合實際船舶管路孔板設計經(jīng)驗,將六級同心孔板設置為六級偏心孔板,開孔的偏心距設計為8 mm,并開展六級偏心孔板節(jié)流特性的計算分析。

六級偏心節(jié)流孔板前后流線分布見圖11。由圖11可知,流體以較高流速穿過第一級節(jié)流孔板的下游區(qū)域時,由于節(jié)流孔板的偏心結(jié)構(gòu),改善流體的流動水力結(jié)構(gòu),流體并未產(chǎn)生嚴重的射流現(xiàn)象。

圖11 偏心孔板流動分布

六級偏心孔板前后區(qū)域流速矢量分布見圖12。沿著給水管路中流體的流動方向,部分相鄰的兩級孔板之間形成小尺寸的漩渦。究其原因主要是第一級至第六級偏心孔板結(jié)構(gòu)的逐級導流作用,規(guī)避節(jié)流孔板區(qū)域流體射流和射流現(xiàn)象的發(fā)生,降低各級孔板之間渦流的誘發(fā)概率,減小流體局部渦流的尺度,因此減緩了管道和孔板的低頻振動。

六級偏心節(jié)流孔板湍動能分布見圖13、14。

圖12 偏心孔板流速矢量分布

圖13 偏心孔板湍動能橫截面分布規(guī)律

圖14 偏心孔板湍動能縱截面分布規(guī)律

由圖13、14可知,在偏心孔板的導流作用下,降低流體微團的湍流脈動,且各級偏心節(jié)流孔板之間的流體湍動能分布比較均勻,其中湍動能的最大值約為40 m2/s2,遠小于六級同心孔板的湍動能最高值。

六級偏心孔板流速變化見圖15～17。

圖15 偏心孔板流速分布

由圖可知,由于各級偏心孔板的導流作用,在兩相鄰孔板間出現(xiàn)明顯的流體速度由小到大的周期性變化過程,且由于各級偏心節(jié)流孔板沿程連續(xù)作用,導致流體的速度逐漸降低。流體的最高速度約為53 m/s,降低了流體的沖擊動能,弱化流體對管路和孔板的沖擊,減緩了給水管路的振動。

圖16 偏心孔板流速橫截面分布

圖17 偏心孔板流速縱截面分布

六級偏心孔板的壓力變化見圖18～20。

圖18 偏心孔板壓力分布

圖19 偏心孔板壓力縱截面分布

圖20 偏心孔板壓力橫截面分布

沿著給水的流動方向,流體壓力逐漸降低。主要是由于偏心孔板的導流作用,強化偏心孔板附近區(qū)域的流體能量耗散,導致大部分動壓能轉(zhuǎn)換為熱能,少部分動壓轉(zhuǎn)換為靜壓,因此流體的壓力不斷降低。各級偏心孔板節(jié)流下游流域沒有負壓區(qū)產(chǎn)生,且孔板節(jié)流后的流體壓力均高于給水溫度為104.5℃對應的飽和壓力0.118 8MPa,因此避免流體發(fā)生汽蝕現(xiàn)象,減少船舶給水管路系統(tǒng)的高頻振動。數(shù)值計算壓力曲線與理論設計數(shù)值吻合較好。

4 結(jié)論

1)節(jié)流孔板下游產(chǎn)生渦流,誘發(fā)流體流動能量耗散,促使流體部分動壓能轉(zhuǎn)換為靜壓能,導致節(jié)流孔板下游壓力不能恢復至孔板上游壓力。

2)六級同心節(jié)流孔板會導致孔板下游產(chǎn)生大尺度渦流和汽蝕現(xiàn)象,極易誘發(fā)管路振動,不適用于大壓降小管距的船舶給水管路系統(tǒng)。

3)六級偏心節(jié)流孔板能減小渦流尺度,抑制孔板汽蝕現(xiàn)象,發(fā)揮孔板限流降壓功效,減緩流體沖擊管壁的動量,規(guī)避大壓降小管距的船舶給水管路劇烈振動和噪聲。

［1］張寶峰.多級節(jié)流孔板的設計計算［J］.西北電力技術(shù),2005(5):27-30.

［2］聶 方,張 賢,董 君.凝結(jié)水泵再循環(huán)管道上調(diào)節(jié)閥和節(jié)流孔板的選型計算分析［J］.節(jié)能技術(shù), 2011,29(1):470-473.

［3］李 妍,陸道綱,曾小康.適用于大壓降小間距管道的節(jié)流件設計及分析［J］.核動力工程,2013,34(4): 126-129.

［4］王秋穎.船用蒸汽管道系統(tǒng)振動及抗沖擊特性有限元仿真分析［J］.船海工程,2007,36(5):69-72.

［5］YAN Y,THORPE R B.Flow regime transition due to cavitation in the flow through an orifice［J］.International Journal of Multiphase Flow,1990,16(6):1023-1045.

［6］趙繼剛,劉德榮.化工行業(yè)水泵再循環(huán)管道節(jié)流孔板的選型計算［J］.當代化工,2012,41(11):1232-1234.

［7］靈 敏,李洪林.大壓降短距離管道節(jié)流孔板數(shù)值模擬研究 ［J］.長江大學學報:自然科學版,2011,8 (2):68-70.

［8］電力工業(yè)部.DL/T 5054-1996火力發(fā)電廠汽水管道設計技術(shù)規(guī)定［S］.1996.

Design of the Multistage Orifice Plate and Throttle Characteristics for Marine Feed Pipe with Large Pressure Drop

YANG Yuan-long1,ZHENG Wen2
(1 China Ship Development and Design Center,Wuhan 430064,China; 2Wuhan Ship Development Design Institute,Wuhan 430064,China)

In order to reduce vibration and noise from feed pipe with large pressure drop and ensure stable ship steam and water circle characteristics,multistage throttle orifice need to be designed reasonably.Taking the feed pipe with large pressure drop as the prototype,the CFDmethod is utilized to analyze numerically the pressure drop,flow field structure and turbulent distributions in detail between concentric and eccentric orifice plate.The calculated results show that the six-stage eccentric orifice plate is suitable to the feed pipeline system with large pressure drop and small space.The numerical simulation results agreed with those of theatrical design well.

multistage throttle orifice;feed pipe;cavitation

10.3963/j.issn.1671-7953.2015.03.039

U664.5

A

1671-7953(2015)03-0164-05

2015-02-28

修回日期:2015-03-05

國家自然科學基金資助項目(51309063)

楊元龍(1986-),男,碩士,助理工程師

研究方向:艦船蒸汽動力系統(tǒng)性能仿真及設計

E-mail:long31609@163.com