船舶典型充汽管路水動(dòng)力特性數(shù)值預(yù)測(cè)

楊元龍，李堅(jiān)波，孫寶芝，楊龍濱，王興剛中國(guó)艦船研究設(shè)計(jì)中心，湖北武漢40064海軍駐上海江南造船（集團(tuán)）有限責(zé)任公司軍事代表室，上海09哈爾濱工程大學(xué)動(dòng)力與能源工程學(xué)院，黑龍江哈爾濱5000

船舶典型充汽管路水動(dòng)力特性數(shù)值預(yù)測(cè)

楊元龍1，李堅(jiān)波2，孫寶芝3，楊龍濱3，王興剛1
1中國(guó)艦船研究設(shè)計(jì)中心，湖北武漢430064
2海軍駐上海江南造船（集團(tuán)）有限責(zé)任公司軍事代表室，上海201913
3哈爾濱工程大學(xué)動(dòng)力與能源工程學(xué)院，黑龍江哈爾濱150001

建立船舶典型充汽管路三維計(jì)算模型，運(yùn)用計(jì)算流體力學(xué)方法進(jìn)行充汽管路流體水動(dòng)力特性數(shù)值研究，捕捉流體壓力、流速、密度、湍流動(dòng)能等關(guān)鍵參數(shù)分布規(guī)律，基于流致振動(dòng)預(yù)測(cè)機(jī)理，揭示與流致振動(dòng)密切相關(guān)的流體流動(dòng)能量。計(jì)算結(jié)果顯示：充汽管路彎頭部位產(chǎn)生由外壁向內(nèi)壁發(fā)展的二次流動(dòng)現(xiàn)象，導(dǎo)致彎頭區(qū)域呈現(xiàn)內(nèi)壁壓力小、流速高、密度小、湍流動(dòng)能及耗散率高等情況；基于流體流動(dòng)能量分布規(guī)律，充汽管路彎頭區(qū)域能量出現(xiàn)最大值，且彎頭內(nèi)壁流動(dòng)能量高于外壁，可以預(yù)測(cè)充汽管路彎頭內(nèi)壁側(cè)承受流致振動(dòng)破壞較嚴(yán)重，實(shí)際充汽管路檢測(cè)數(shù)據(jù)驗(yàn)證了數(shù)值預(yù)測(cè)結(jié)果。

充汽管路；水動(dòng)力特性；流致振動(dòng)

0　引言

船舶充汽系統(tǒng)管路是實(shí)現(xiàn)從增壓鍋爐向蒸汽蓄熱器供汽的關(guān)鍵管路，具有充汽負(fù)荷高、熱沖擊大、周期性頻繁波動(dòng)的特點(diǎn)［1-2］。在周期性大負(fù)荷擾動(dòng)的充汽過(guò)程中，充汽系統(tǒng)管路的壓力場(chǎng)和速度場(chǎng)產(chǎn)生劇烈的高頻率脈動(dòng)，不僅在管路近壁區(qū)域激發(fā)很高的壁面剪切應(yīng)力，更增強(qiáng)流體湍流脈動(dòng)程度，引起顯著的流體動(dòng)力學(xué)效應(yīng)，從而誘發(fā)充汽管路產(chǎn)生劇烈的流致振動(dòng)現(xiàn)象［3-5］。因此研究周期性大負(fù)荷擾動(dòng)下充汽管路水動(dòng)力特性對(duì)船用蒸汽動(dòng)力系統(tǒng)控制及安全穩(wěn)定運(yùn)行具有重要意義。

目前國(guó)內(nèi)外學(xué)者針對(duì)蒸汽系統(tǒng)管路流場(chǎng)特征做了大量研究［6-9］，但在公開(kāi)文獻(xiàn)中關(guān)于船舶充汽系統(tǒng)管路水動(dòng)力及流致振動(dòng)特性的研究較少。本文將以實(shí)際船舶典型充汽管路為原型，采用CFD計(jì)算流體軟件，模擬船舶典型充汽管路的流動(dòng)，捕捉壓力、密度、流速等關(guān)鍵參數(shù)的分布規(guī)律，結(jié)合充汽管路流致振動(dòng)的誘發(fā)機(jī)理，揭示與流致振動(dòng)密切相關(guān)的流體流動(dòng)能量，明晰充汽管路的“薄弱”部位，最終摸索出大負(fù)荷擾動(dòng)條件下船舶充汽系統(tǒng)管路流體水動(dòng)力及流致振動(dòng)變化特性的預(yù)測(cè)方法。

1　物理結(jié)構(gòu)及網(wǎng)格模型

圖1為船舶典型充汽系統(tǒng)管路結(jié)構(gòu)圖。充汽管外徑為0.15m，管路總長(zhǎng)為11.2m，其包含6個(gè)直角彎頭，彎曲半徑為0.45m。

如圖2所示，充汽系統(tǒng)管路采用六面體結(jié)構(gòu)化網(wǎng)格進(jìn)行網(wǎng)格劃分，筒壁處基于O型網(wǎng)格技術(shù)手段添加邊界層網(wǎng)格，網(wǎng)格近壁尺寸y+為30～40，保證近壁網(wǎng)格尺度滿足壁面率要求?；诰W(wǎng)格敏感性分析，計(jì)算域共計(jì)300 000個(gè)網(wǎng)格單元。

圖1　充汽管路結(jié)構(gòu)Fig.1 Structure of charging pipe

圖2　充汽管路網(wǎng)格模型Fig.2 Gridmodelof charging pipe

2　計(jì)算邊界條件

為準(zhǔn)確預(yù)測(cè)船舶充汽管路的水動(dòng)力變化特性，基于實(shí)際船舶典型充汽系統(tǒng)管路的流動(dòng)和傳熱參數(shù)，設(shè)置邊界計(jì)算條件：過(guò)熱蒸汽進(jìn)口壓力為6.18 MPa，進(jìn)口溫度為430℃，出口壓力為5.6 MPa，管壁滿足絕熱（管路包覆保溫層，且外表面溫度低）和無(wú)滑移邊界條件。過(guò)熱蒸汽介質(zhì)參數(shù)的計(jì)算標(biāo)準(zhǔn)均來(lái)自IAPWSIF97數(shù)據(jù)庫(kù)［10］。

3　充汽管路水動(dòng)力特性預(yù)測(cè)

3.1水動(dòng)力特性分析

圖3為充汽管路壓力變化云圖。由圖可知，沿著充汽管路中蒸汽介質(zhì)流動(dòng)方向，蒸汽壓力逐漸降低，由于在充汽過(guò)程中，管路內(nèi)表面存在一定粗糙度，導(dǎo)致過(guò)熱蒸汽與管路產(chǎn)生摩擦阻力。同時(shí)，充汽管路上有6個(gè)直角彎頭，對(duì)蒸汽介質(zhì)產(chǎn)生局部阻力。因此，在摩擦阻力和局部阻力的綜合作用下，充汽管路的壓力不斷降低。如圖3所示，直角彎頭的外壁壓力明顯高于內(nèi)壁壓力，主要是由于在離心力的作用下，推擠彎頭中心區(qū)域的蒸汽介質(zhì)至外壁，導(dǎo)致彎頭外壁區(qū)域壓力明顯高于內(nèi)壁。

圖3　充汽管路壓力變化云圖Fig.3 Pressure contoursof charging pipe

圖4所示為充汽管路蒸汽介質(zhì)流速矢量分布規(guī)律?？梢钥闯觯刂淦苈方橘|(zhì)流動(dòng)方向，蒸汽流速逐漸升高。從充汽管路6號(hào)彎頭處橫剖面流速矢量圖可以看出，由于離心力的作用誘發(fā)彎頭橫剖面產(chǎn)生二次流動(dòng)，致使流速較高的蒸汽快速流向彎頭內(nèi)壁面，且彎頭中心區(qū)域流體向外壁擴(kuò)散，彎頭內(nèi)壁附近流體橫掠中心區(qū)域向外壁流動(dòng)。計(jì)算結(jié)果與文獻(xiàn)［6］中的數(shù)值模擬結(jié)果吻合較好。

圖4　充汽管路流速矢量分布規(guī)律Fig.4 Velocity vector distributionsof charging pipe

圖5～圖6給出了充汽管路湍流動(dòng)能和湍流耗散率變化云圖。由圖可知，沿著流動(dòng)方向，充汽管路湍流脈動(dòng)動(dòng)能和耗散率逐漸增大。由于湍流動(dòng)能與蒸汽流速成正比關(guān)系，沿著充汽管路中蒸汽介質(zhì)流動(dòng)方向，蒸汽介質(zhì)壓力不斷降低，流速升高，促使湍流動(dòng)能和耗散率逐漸增大，且在二次流動(dòng)作用下，彎頭內(nèi)壁區(qū)域的湍流動(dòng)能和耗散率明顯高于外壁，增強(qiáng)湍流流體微團(tuán)之間的動(dòng)量傳遞，導(dǎo)致彎管區(qū)域湍流劇烈脈動(dòng)。

圖5　充汽管路湍流動(dòng)能分布規(guī)律Fig.5 Turbulence kinetic energy distributions of charging pipe

圖6　充汽管路湍流耗散率分布規(guī)律Fig.6 Turbulence eddy dissipation distributionsof charging pipe

充汽管路壁面剪切應(yīng)力變化規(guī)律如圖7所示。在充汽管路直管區(qū)域，壁面剪切應(yīng)力不是很大，但在彎頭部位壁面剪切應(yīng)力突然增大，且沿著充汽介質(zhì)流動(dòng)方向，流速不斷升高，剪切應(yīng)力峰值也逐漸增大。同時(shí)，由于彎頭內(nèi)壁湍流脈動(dòng)更劇烈，彎頭內(nèi)壁側(cè)剪切應(yīng)力最大值約是外壁的3倍。

圖7　充汽管路壁面剪切應(yīng)力分布規(guī)律Fig.7 Wall shear distributionsof charging pipe

3.2流致振動(dòng)預(yù)測(cè)分析

圖8～圖9示出了充汽管路流速分布曲線和云圖。充汽管路內(nèi)蒸汽流速呈先快速升高再急劇降低的變化規(guī)律，究其原因主要是沿著流動(dòng)方向壓力不斷減小，促使流速快速升高，并由于彎管部位二次流動(dòng)的作用，導(dǎo)致流速急劇升高。因此在兩相鄰彎管間，呈流速由高到低的周期性變化規(guī)律，且由于彎管沿程連續(xù)作用，導(dǎo)致蒸汽流速快速升高。

圖８　充汽管路流速變化曲線Fig.8 Velocity curve of charging pipe

圖9　流速局部分布規(guī)律Fig.9 Local velocity distributions

圖10～圖11給出了充汽管路內(nèi)蒸汽密度變化曲線和云圖。由圖可知，沿著充汽管路蒸汽流動(dòng)方向，蒸汽密度呈先快速降低再逐漸升高的分布規(guī)律。主要是由于沿著流動(dòng)方向，壓力降低，蒸汽流速升高，在質(zhì)量流量守恒條件下，可壓縮過(guò)熱蒸汽的密度不斷降低。但由于彎管區(qū)域流速急劇升高，進(jìn)而導(dǎo)致彎管部位密度快速減小。因此在兩相鄰彎管間，出現(xiàn)明顯的蒸汽密度由小到大的周期性變化過(guò)程，且由于彎管沿程連續(xù)作用，導(dǎo)致蒸汽密度不斷減小。

圖10　充汽管路內(nèi)蒸汽密度變化曲線Fig.10 Steam density curve of charging pipe

圖11　充汽管路內(nèi)蒸汽密度變化云圖Fig.11 Steam density contoursof charging pipe

流致振動(dòng)的誘發(fā)機(jī)理主要包括湍流激振、流體彈性不穩(wěn)定、漩渦脫落和聲共振。文獻(xiàn)［5］指出，湍流激振、流體彈性不穩(wěn)定是流體誘導(dǎo)振動(dòng)的主要機(jī)理。對(duì)于流動(dòng)機(jī)制復(fù)雜的船舶充汽系統(tǒng)管路，湍流誘導(dǎo)充汽管路的壓力和速度產(chǎn)生隨機(jī)性的脈動(dòng)，而脈動(dòng)的壓力和速度場(chǎng)產(chǎn)生的湍流激振力不斷供給管子能量，管子吸收能量產(chǎn)生振動(dòng)且引起顯著的流體動(dòng)力效應(yīng)。因此，基于湍流激振、流體彈性不穩(wěn)定2種機(jī)理，誘發(fā)充汽管路流致振動(dòng)的激振力與流體流動(dòng)的能量成正比關(guān)系，其中流體流動(dòng)能量等于流體平均密度ρ乘以流速U的平方值［3-5］，即ρU2。

圖12～13給出了充汽管路流體流動(dòng)能量分布曲線和云圖。從圖中可以發(fā)現(xiàn)，蒸汽流速和密度的綜合作用，導(dǎo)致在兩相鄰彎管間出現(xiàn)明顯的流體流動(dòng)能量由大到小的周期性變化過(guò)程，且由于彎管沿程連續(xù)作用，導(dǎo)致流體流動(dòng)能量不斷增大?；谏鲜龅牧髦抡駝?dòng)預(yù)測(cè)分析機(jī)理，彎管部位流體流動(dòng)能量更大，預(yù)示充汽管路的彎管區(qū)域承受流致振動(dòng)破損更嚴(yán)重。

圖12　充汽管路流體流動(dòng)能量變化曲線Fig.12 Fluid energy curve of charging pipe

圖13　流體流動(dòng)能量分布規(guī)律Fig.13 Fluid flow energy distributions

圖14顯示了充汽管路彎管部位流致振動(dòng)預(yù)測(cè)數(shù)據(jù)與實(shí)際彎頭檢測(cè)結(jié)果的對(duì)比。由于彎管區(qū)域二次流動(dòng)的影響，促使彎管內(nèi)側(cè)流速高，流體流動(dòng)能量大，從而導(dǎo)致充汽管路內(nèi)部區(qū)域承受流致振動(dòng)破損更加嚴(yán)重，另外彎頭部位流致振動(dòng)數(shù)值預(yù)測(cè)與實(shí)際船舶充汽管路檢測(cè)結(jié)果基本吻合，驗(yàn)證了數(shù)值方法預(yù)測(cè)充汽管路流致振動(dòng)分布規(guī)律的可行性和準(zhǔn)確性。

圖14　流體流動(dòng)能量局部分布規(guī)律Fig.14 Local fluid energy distributions

4　結(jié)論

本文通過(guò)建立船舶典型充汽管路三維計(jì)算模型，運(yùn)用CFD方法開(kāi)展充汽管路系統(tǒng)水動(dòng)力特性的數(shù)值研究，明晰了壓力、流速、密度、湍流動(dòng)能等關(guān)鍵參數(shù)的分布規(guī)律，并基于流致振動(dòng)預(yù)測(cè)機(jī)理揭示了與流致振動(dòng)密切相關(guān)的流體流動(dòng)能量。主要結(jié)論如下：

1）在離心力作用下，充汽管路彎頭部位產(chǎn)生由外壁向內(nèi)壁發(fā)展的二次流動(dòng)現(xiàn)象；

2）二次流動(dòng)導(dǎo)致彎頭區(qū)域內(nèi)壁壓力小、流速高、密度小、湍流動(dòng)能及耗散率高，外壁壓力大、流速低、密度大、湍流動(dòng)能及耗散率低；

3）基于流體流動(dòng)能量分布規(guī)律，充汽管路彎頭部位流動(dòng)能量出現(xiàn)最大值，且彎頭內(nèi)壁流動(dòng)能量高于外壁，可以預(yù)測(cè)充汽管路彎頭內(nèi)壁側(cè)承受流致振動(dòng)破壞較嚴(yán)重，實(shí)際充汽管路破損數(shù)據(jù)驗(yàn)證了本文CFD預(yù)測(cè)結(jié)果。

［1］胡繼敏，金家善，嚴(yán)志騰.汽源參數(shù)變化對(duì)儲(chǔ)汽筒充汽系統(tǒng)變工況性能的影響［J］.海軍工程大學(xué)學(xué)報(bào)，2012，24（5）：102-107. HU Jimin，JIN Jiashan，YAN Zhiteng.Influence of steam supply parameters variation on performances of steam accumulator system［J］.Journal of Naval University of Engineering，2012，24（5）：102-107.

［2］孫寶芝，郭家敏，雷雨，等.船用蒸汽蓄熱器非平衡熱力過(guò)程［J］.化工學(xué)報(bào)，2013，64（增刊1）：59-65. SUN Baozhi，GUO Jiamin，LEIYu，et al.Non-equilibrium thermodynamics process of marine steam accumulator［J］.CIESC Journal，2013，64（Supp 1）：59-65.

［3］SUN B Z，YANG Y L.Numerically investigating the influence of tube support plates on thermal-hydraulic characteristics in a steam generator［J］.Applied Thermal Engineering，2013，51（1/2）：611-622.

［4］PITTIGREW M J，TAYLOR C E，SUBASH N. Flow-induced vibration specifications for steam generators and liquid heat exchangers：AECL-11401［R］. Chalk River，Ontario：Fluid Sealing and Dynamics Branch Chalk River Laboratories，1995.

［5］AXISA F，ANDUNESJ，VILLARD B.Overview of numerical methods for predicting flow-induced vibration ［J］.Journalof Pressure Vessel Technology，1988，110 （1）：6-14.

［6］江山，張京偉，吳崇健，等.基于FLUENT的90°圓形彎管內(nèi)部流場(chǎng)分析［J］.中國(guó)艦船研究，2008，3（1）：37-41. JIANG Shan，ZHANG Jingwei，WU Chongjian，et al. Numerical simulation of inner flow in 90°bending duct of circular-section based on FLUENT［J］.Chinese Journal of Ship Research，2008，3（1）：37-41.

［7］查友其，張成，丁建軍，等.船用海水冷卻管路沖擊響應(yīng)計(jì)算［J］.中國(guó)艦船研究，2013，8（4）：68-73. ZHA Youqi，ZHANG Cheng，DING Jianjun，et al. Shock response calculation of marine seawater cooling piping systems［J］.Chinese Journal of Ship Research，2013，8（4）：68-73.

［8］FERNG Y M，TSENG Y S，PEIB S，etal.A two-phase methodology to predict FAC wear sites in the piping system of a BWR［J］.Nuclear Engineering and Design，2008，238（9）：2189-2196.

［9］龐天照，郭偉，盛元平，等.基于流固耦合的U型管流致振動(dòng)數(shù)值分析［J］.中國(guó)艦船研究，2011，6（4）：29-33. PANG Tianzhao，GUOWei，SHENG Yuanping，et al. Numerical simulation of fluid-induced vibration in U-tube based on FSI［J］.Chinese Journal of Ship Research，2011，6（4）：29-33.

［10］WAGNERW，COOPER JR，DITTMANN A，et al. The IAPWS industrial formulation 1997 for the thermodynamic properties of water and steam［J］.Journal of Engineering for Gas Turbines and Power Transactionsof the ASME，2000，122（1）：150-182.

［責(zé)任編輯：田甜］

Num erical prediction on the hydrodynamic characteristics for typical ship charging pipes

YANGYuanlong1，LIJianbo2，SUN Baozhi3，YANG Longbin3，WANGXinggang1
1 China Ship Development and Design Center，Wuhan 430064，China
2 Naval Military Representative Office in Jiangnan Shipyard（Group）Co.，Ltd.，Shanghai201913，China
3 College of Power and Energy Engineering，Harbin Engineering University，Harbin 150001，China

A three dimensional calculation model of typical ship charging pipes is established in this paper，where its hydrodynamic characteristics are investigated using the CFD method.The distribution of key parameters such as pressure，velocity，density，and turbulence kinetic energy are obtained.Based on the prediction mechanism of flow induced vibration，the flow energy，which is closely related to flow induced vibration，is revealed.The calculated results show that the secondary flow phenomenon，occurring from the outer wall to the inner wall，is most prominent at the elbow of charging pipes.This secondary flow could lead to low pressure and density，higher velocity and turbulence kinetic energy with the region. Based on fluid flow-energy distributions，the flow-energy are estimated to reach its maximum value at the elbow；meanwhile，the flow-energy for the inner wall will be higher than that for the outer wall.These hydrodynamic results suggest that the flow induced vibration damage for the inner wall around the elbow could be rather serious，which is then verified with actual test data.

charging pipe；hydrodynamic characteristics；flow induced vibration

U664.5

A

10.3969/j.issn.1673-3185.2015.04.020

2014-09-01網(wǎng)絡(luò)出版時(shí)間：2015-7-29 9:23:57

楊元龍（通信作者），男，1986年生，碩士，助理工程師。研究方向：艦船蒸汽動(dòng)力系統(tǒng)性能仿真及設(shè)計(jì)。E-mail：long31609@163.com