，，，

(武漢第二船舶設(shè)計研究所，武漢 430205)

船舶自流注水系統(tǒng)閥門流動噪聲預(yù)測

方超，蔡標(biāo)華，馬士虎，俞健

(武漢第二船舶設(shè)計研究所，武漢 430205)

針對船舶自流注水系統(tǒng)閥門噪聲問題，利用CFD技術(shù)對系統(tǒng)閥門附近流場進(jìn)行3種壓差工況的數(shù)值計算，分別使用閥門噪聲經(jīng)驗公式和聲學(xué)數(shù)值計算法對閥門流動噪聲進(jìn)行分析。計算結(jié)果與測量結(jié)果對比表明，2種方法計算誤差均較小，驗證了閥門噪聲經(jīng)驗公式和聲學(xué)數(shù)值計算法對閥門噪聲預(yù)測的可靠性，在工程實際中可用于閥門噪聲的初步檢測。

閥門；流噪聲；CFD；聲學(xué)計算

自流注水系統(tǒng)利用舷外和艙內(nèi)壓力差對水艙進(jìn)行自流注水，在各種船舶海水系統(tǒng)中有著廣泛的應(yīng)用。自流注水系統(tǒng)閥門噪聲問題突出，系統(tǒng)優(yōu)化設(shè)計的過程中閥門噪聲預(yù)測尤其重要。目前國內(nèi)外學(xué)者針對船舶水動力噪聲預(yù)測做了大量分析研究[1-4],但少見關(guān)于船用閥門噪聲預(yù)測研究的文獻(xiàn)報道。根據(jù)產(chǎn)生機(jī)理不同閥門噪聲分為機(jī)械振動噪聲、氣蝕噪聲和流體噪聲[5-8]。其中機(jī)械振動噪聲一般可通過優(yōu)化結(jié)構(gòu)設(shè)計予以消除。氣蝕噪聲和流體噪聲是由流體介質(zhì)流動所引起，合稱為流動噪聲，二者相互耦合。針對自流注水系統(tǒng)(見圖1)的閥門流動噪聲預(yù)測問題，選擇自流注水系統(tǒng)3種典型的工況，建立模型利用CFD技術(shù)進(jìn)行流場計算，提取流場計算結(jié)果，分別采用VDMA經(jīng)驗公式和數(shù)值計算的方法對系統(tǒng)閥門噪聲進(jìn)行計算，并將計算結(jié)果與試驗結(jié)果進(jìn)行對比，驗證2種閥門噪聲預(yù)測方法的可靠性。

圖1 船舶自流注水系統(tǒng)原理

1 閥門流動噪聲預(yù)測方法介紹

閥門管路的流動噪聲，目前主要有2類預(yù)測方法：傳統(tǒng)的經(jīng)驗公式法；聲學(xué)數(shù)值計算的方法。

1)經(jīng)驗公式法。在已知管路的幾何尺寸和流動參數(shù)情況下，利用相應(yīng)公式估算得到噪聲聲壓級。國際電氣工程標(biāo)準(zhǔn)IEC60534、西德機(jī)器制造業(yè)協(xié)會(VDMA)均針對閥門噪聲預(yù)測給出了估算公式，2種估算公式結(jié)構(gòu)形式和考慮因素大體相同。選擇VDMA公式進(jìn)行自流注水系統(tǒng)閥門噪聲估算。

Lw=40+10lgKv+18lgΔPv-5lgρf+ΔLa，cav(1)

(2)

式中，xF=(p1-p2)/(p1-p2)，Zy為調(diào)節(jié)閥開度為Y時的Z值；球閥Z值一般取0.2。

2)數(shù)值計算法。管內(nèi)流動噪聲源主要來自壁面的偶極子源和雷諾應(yīng)力的四極子源[9]。在高馬赫數(shù)下，四極子聲源項的量級與主要聲源項相同，不可忽略；而在低馬赫數(shù)下，偶極子聲源的貢獻(xiàn)遠(yuǎn)遠(yuǎn)大于四極子聲源，四極子聲源在計算中可忽略不計。本例中流動為低馬赫數(shù)，數(shù)值計算時四極子聲源可以忽略。

數(shù)值計算流體噪聲目前有2種途徑：基于FW-H積分形式方程的解法；先結(jié)合CFD計算與Lighthill理論求解壁面偶極子，再使用邊界元法求解聲學(xué)Helmholtz方程。積分求解FW-H方程的優(yōu)點是計算量和計算格式要求相對較低，但也存在只能計算遠(yuǎn)場輻射，內(nèi)存計算誤差大和不能考慮結(jié)構(gòu)和聲學(xué)裝置的影響等缺點；邊界元法的優(yōu)點是能考慮聲傳播的問題，可以進(jìn)行復(fù)雜聲場的計算，不足是只考慮了壁面的偶極子聲源。

管內(nèi)流場噪聲計算需要考慮管壁對聲波的反射吸收等聲傳播問題，采用積分求解FW-H方程法誤差較大，因此宜采用邊界元法。邊界元法以Lighthill方程為控制方程，對控制方程進(jìn)行傅里葉變換得到Helmholtz方程。

2pa6k2pa=-

(3)

通過CFD流場瞬態(tài)計算，提取管道閥門壁面壓力脈動信號，導(dǎo)入聲學(xué)計算軟件轉(zhuǎn)換成等效的流體聲源，并映射到聲學(xué)邊界元網(wǎng)格上，求解Helmholtz方程可得到管道內(nèi)聲場特性。

2 閥門瞬態(tài)流場計算

經(jīng)驗公式法和數(shù)值計算法預(yù)測閥門噪聲都是基于CFD計算結(jié)果，流場計算結(jié)果對于閥門噪聲預(yù)測的精確度有重要的影響。采用經(jīng)驗公式結(jié)合流場仿真計算結(jié)果預(yù)測閥門噪聲基于瞬態(tài)計算結(jié)果，本文分別在3種壓差工況下對球閥進(jìn)行瞬態(tài)流場數(shù)值計算。

2.1 幾何模型和計算工況

球閥三維數(shù)值模型進(jìn)出口管徑為100 mm，閥芯通徑為80 mm，3種典型工況見表1。對計算幾何模型進(jìn)行簡化處理，計算流場時不考慮壁面厚度。利用ICEM軟件對幾何模型進(jìn)行網(wǎng)格劃分，網(wǎng)格總數(shù)為263萬，滿足計算的精度要求。

表1 閥門仿真計算工況

流場邊界條件均為壓力進(jìn)口和壓力出口，管內(nèi)壁表面和閥芯內(nèi)表面為固體壁面邊界，采用無滑移條件。流體介質(zhì)為水，不可壓縮，流動過程無熱量交換，溫度為20 ℃ ，計算過程中忽略重力影響。

2.2 流場計算結(jié)果

首先利用RANS模型對流場進(jìn)行定場模擬，流場穩(wěn)定后，再用LES模型進(jìn)行非定場計算。計算時間步長設(shè)為5×10-5s,計算時間步數(shù)設(shè)為200步，殘差收斂標(biāo)準(zhǔn)10-3。模擬中得到最高頻率為10 kHz,滿足工程實際中關(guān)注的頻率范圍[10]。

3種工況流場壓力和速度分布見圖2，流動噪聲的產(chǎn)生與流體壓力和速度密切相關(guān)。

圖2 3種工況下壓力和速度分布

對于工況一，系統(tǒng)流量為496.8 m3/h，由于閥門的節(jié)流作用，最高流速達(dá)到了61.62 m/s，閥后壓力迅速降低，存在較大的負(fù)壓區(qū)，負(fù)壓低于 20 ℃時海水飽和蒸汽壓力，產(chǎn)生空化現(xiàn)象，誘發(fā)嚴(yán)重的空化噪聲。閥門處流場的劇烈變化會對閥門結(jié)構(gòu)產(chǎn)生強(qiáng)烈的沖擊，導(dǎo)致閥門及管壁振動并向外輻射噪聲。工況二、三因注水艙背壓的原因流量減小，分別為367.9 m3/h和197.7 m3/h，同時由于水艙背壓的存在，空化現(xiàn)象消失。

3種工況下部分流動參數(shù)計算結(jié)果見表2。

表2 3種工況下流場計算結(jié)果

3 閥門流噪聲計算結(jié)果與分析

3.1 經(jīng)驗公式計算結(jié)果

根據(jù)流場計算結(jié)果，結(jié)合式(2)計算得到3種工況下自流注水系統(tǒng)閥門流噪聲總聲壓級。VDMA閥門水動力噪聲聲壓級估算公式計算的閥門噪聲為空氣噪聲，測試條件與標(biāo)準(zhǔn)化試驗方案相同，測點設(shè)在調(diào)節(jié)閥出口法蘭的平面內(nèi)，距管道外壁1 m處，計算結(jié)果為流體動力噪聲總聲壓級。

經(jīng)驗公式計算和試驗測量結(jié)果對比見表3。

表3 VDMA經(jīng)驗公式噪聲計算結(jié)果

工況一公式估算結(jié)果比試驗結(jié)果大2.62 dB，工況二和工況三公式估算結(jié)果比試驗結(jié)果分別小13.00 dB和10.97 dB。VDMA閥門水動力噪聲聲壓級估算公式只考慮了在封閉管道系統(tǒng)中的流體動力學(xué)噪聲，而并不包括調(diào)節(jié)閥內(nèi)部可動部件產(chǎn)生的響聲，以及在固體材料中傳播的聲音，或由于某種原因反射和共振造成的聲音放大效應(yīng)。目前的空氣噪聲試驗測量技術(shù)手段無法將多種噪聲源混合的噪聲分開測量，試驗測量空氣噪聲包含閥門及管壁振動輻射噪聲，因此理論上試驗所測空氣噪聲相對式(1)估算噪聲較大，考慮流場數(shù)值計算誤差的因素，可認(rèn)為估算結(jié)果與真實情況基本吻合。

3.2 聲學(xué)數(shù)值計算結(jié)果

將管道壁面壓力脈動導(dǎo)入聲學(xué)計算軟件Virtual Lab中計算流動噪聲，計算得到不同工況下監(jiān)測面聲壓分布見圖3。監(jiān)測面設(shè)在調(diào)節(jié)閥進(jìn)口法蘭處，監(jiān)測面中心處設(shè)置監(jiān)測點，監(jiān)測點聲壓頻域分布見圖4。

由計算結(jié)果可知，工況一、二、三總聲壓級分別為207.4，191.4，191.3 dB，閥門流動噪聲隨閥門前后壓差的減小而減小。當(dāng)壓差減小到一定程度后，流體介質(zhì)流速降低，閥門處流體不再發(fā)生空化，空化噪聲消失，繼續(xù)減小閥門壓差對流動噪聲的抑制作用減弱，與式(1)預(yù)測趨勢相吻合。

圖3 3種工況下聲壓分布云圖(100 Hz)

圖4 3種工況下監(jiān)測點聲壓頻域分布

數(shù)值計算結(jié)果與試驗測量結(jié)果對比見圖5，3種工況下運(yùn)用水聽器測得監(jiān)測點處流動噪聲總聲壓級分別為198.3，186.5，180.3 dB，相對計算值誤差分別為9.1，4.9，11.0 dB，數(shù)值計算結(jié)果與試驗測量結(jié)果隨閥門前后壓差減小變化趨勢基本一致，考慮流場和聲場的計算誤差等因素，數(shù)值計算結(jié)果和試驗測量結(jié)果基本吻合。

圖5 不同壓差下監(jiān)測點總聲壓級計算值和測量值比較

4 結(jié)論

1)2種方法均可以較好地預(yù)測閥門流動噪聲。將計算結(jié)果與試驗測量結(jié)果對比，結(jié)果表明2種方法計算誤差均較小，2種方法預(yù)測閥門噪聲具有可靠性，可用于工程中閥門噪聲的初步檢測，具有較大的工程實用價值。

2)大壓差工況下可采取增大注水艙背壓減小閥門前后壓差的措施抑制閥門噪聲。3種不同壓差工況下，閥門前后壓差越小，閥門流動噪聲越低，但壓差減小到一定程度后，空化現(xiàn)象消失，減小壓差對閥門噪聲的抑制作用減弱。

3)工程實際中，采用VDMA經(jīng)驗公式法計算閥門噪聲一般是基于臺架試驗測量得到的流動參數(shù)，本文通過流場數(shù)值計算得到相應(yīng)流動參數(shù)，閥門噪聲計算結(jié)果誤差較小。采用VDMA經(jīng)驗公式與流場數(shù)值計算相結(jié)合的方法預(yù)測閥門噪聲未見相關(guān)文獻(xiàn)報道，較傳統(tǒng)方法有操作簡單、成本低的優(yōu)勢。

經(jīng)驗公式法和聲學(xué)數(shù)值計算法預(yù)測閥門噪聲都是基于閥門瞬態(tài)流場計算結(jié)果，流場計算的精確度對閥門噪聲預(yù)測的準(zhǔn)確性具有較大影響，后續(xù)研究中可分別采用不同的湍流模型計算閥門流場，進(jìn)行比較后選擇合適的模型以進(jìn)一步提高2種預(yù)測方法的準(zhǔn)確性。

[1] 劉明靜，馬運(yùn)義.潛艇艏部聲吶流噪聲計算方法研究[J].船海工程，2009,38(5)：46-50.

[2] 孟堃宇.基于大渦模擬的潛艇脈動壓力與流噪聲性能數(shù)值計算[D].上海：上海交通大學(xué)，2011.

[3] 董運(yùn)義，吳崇健，張京偉.基于Lilley寬帶聲源模型的舵葉流噪聲數(shù)值模擬研究[J].船海工程，2009,38(5)：24-31.

[4] 張楠.三維孔腔流激噪聲的大渦模擬與聲學(xué)類比預(yù)報與驗證研究[J].船舶力學(xué)，2010,14(1-2)：181-190.

[5] 蔡標(biāo)華,鄭海斌.海水系統(tǒng)噪聲分析及其減振降噪措施綜述[J].艦船科學(xué)與技術(shù)，2009,31(12)：114-117.

[6] 王典佳，唐浩.基于CFD方法的調(diào)節(jié)閥流動噪聲預(yù)測[J].噪聲與振動控制，2014(5)：106-109.

[7] 封海波.海水管路系統(tǒng)中閥門動態(tài)特性與噪聲控制研究[D].哈爾濱，哈爾濱工程大學(xué)，2003.

[8] 楊元龍，鄭文.船用大壓降給水管道多級孔板設(shè)計及節(jié)流特性研究[J].船舶工程，2015,44(3)：164-168.

[9] 劉曉宙，繆國平，唐春權(quán)，等.粘性流體通過二維振蕩圓柱的聲輻射研究[J].船舶力學(xué)，2002,6(4)：75-88.

[10] 魏志.閥體后90°圓形彎管內(nèi)流場和噪聲的數(shù)值模擬[D].上海：上海交通大學(xué)，2013.

Prediction of Flow Noise in Valve of Self-flowing Water Injection System

FANGChao,CAIBiao-hua,MAShi-hu,YUJian

(Wuhan Second Ship Design and Research Institute, Wuhan 430205, China)

In order to predict the noise of the ship’s water injection system, the flow field of the system was calculated by using the CFD technique under the conditions of three kinds of pressure difference. Based on the numerical results of CFD, the valve noise empirical formula and the acoustics numerical method were used to analyze the valve flow noise. The calculated results were compared with the measured results, showing that the two methods had little error, which verified the reliability of the two methods. The two kinds of valve noise prediction methods can be used in engineering practice for the initial detection of valve noise.

valve; flow noise; CFD; acoustic computation

U661.44

A

1671-7953(2017)06-0108-04

10.3963/j.issn.1671-7953.2017.06.025

2016-12-28

2017-03-10

方超(1991—)，男，碩士生

研究方向：艦船振動與噪聲控制技術(shù)