何興建，李 翔，李 軍

[1.中國(guó)特種設(shè)備檢測(cè)研究院，北京 100013；2.中國(guó)石油大學(xué)(北京)，北京 102249]

異徑管是石油化工裝置管道中常用的連接配件，可起到改變管徑調(diào)節(jié)介質(zhì)流速、改變管道走向而管內(nèi)介質(zhì)流線都能順暢,減少其壓降和溫降[1]，且可增加管道柔性減緩管道應(yīng)力等多方面重要作用[2]，常用異徑管的地方是泵的進(jìn)出口、調(diào)節(jié)閥的進(jìn)出口、溫度計(jì)擴(kuò)大管左右和再沸器進(jìn)料分配管后[3]?，F(xiàn)階段異徑管的設(shè)計(jì)與選用主要根據(jù)國(guó)內(nèi)外標(biāo)準(zhǔn)進(jìn)行[4-5]。在實(shí)際使用中由于異徑管沖蝕穿孔、應(yīng)力腐蝕、管件成形中在內(nèi)壁面產(chǎn)生的溝痕等原始缺陷，以及邊緣連接環(huán)縫區(qū)的開(kāi)裂而導(dǎo)致的異徑管失效事故時(shí)有發(fā)生[6-7]。沖蝕穿孔是異徑管局部失效的主要原因之一。

近年來(lái)，國(guó)內(nèi)外學(xué)者運(yùn)用實(shí)驗(yàn)、理論、數(shù)值計(jì)算等方法對(duì)異徑管沖蝕腐蝕磨損、流場(chǎng)以及應(yīng)力等進(jìn)行了部分研究，如偶國(guó)富等通過(guò)二維建模對(duì)異徑管氣液兩相沖蝕磨損進(jìn)行了研究，得出流體流向?yàn)椤靶∵M(jìn)大出”時(shí)對(duì)腐蝕產(chǎn)物保護(hù)膜的影響要比“大進(jìn)小出”小的多[8]；陳孫藝等對(duì)異徑彎管的無(wú)力矩環(huán)向應(yīng)力進(jìn)行了計(jì)算分析，得出內(nèi)壓在異徑彎管中引起的環(huán)向應(yīng)力大小及分布與異徑管大端圓面平均半徑、小端圓面平均半徑、壁厚、錐底角及經(jīng)向彎角,管截面的圓周角有關(guān)[9]。但對(duì)尺寸與沖蝕情況之間關(guān)系的研究較少，為此作者通過(guò)對(duì)不同尺寸異徑管進(jìn)行沖蝕模擬研究，得出沖蝕情況與異徑管尺寸之間的關(guān)系，以期為異徑管的選取、制造及安全檢測(cè)提供參考意見(jiàn)。

1 數(shù)學(xué)模型

1.1 數(shù)學(xué)模型建立

1.1.1 氣相湍流模型

流場(chǎng)在異徑管的變徑區(qū)域變化比較大，為能較好的模擬變徑區(qū)域的真實(shí)情況，采用適合流動(dòng)類型比較廣泛的RNGk-ε模型預(yù)測(cè)異徑管的流場(chǎng)。并在商用軟件Fluent流場(chǎng)模擬平臺(tái)上進(jìn)行數(shù)值模擬研究。

其方程如下。

(1)

(2)

式中：Gk是由于平均速度梯度引起的湍動(dòng)能所產(chǎn)生；ρ是流體密度；Gb是由浮力引起的湍動(dòng)能所產(chǎn)生；C1ε和C3ε是經(jīng)驗(yàn)參數(shù)(分別取1.44和0.09)；YM指可壓湍動(dòng)能的膨脹對(duì)總耗散率的影響；σk是湍動(dòng)能對(duì)應(yīng)的普朗特?cái)?shù)(默認(rèn)值為1.0)；σε是湍動(dòng)能耗散率所對(duì)應(yīng)的普朗特?cái)?shù)(默認(rèn)值為1.3)。

1.1.2 氣固兩相流模型

設(shè)定異徑管的最大入口含塵顆粒濃度為ρi(顆粒)=20 g/m3，對(duì)應(yīng)的顆粒相的體積分率為ε=ci/ρp?10%，模擬空間屬稀相范疇，可以忽略顆粒之間的碰撞[10]，可以使用DPM模型。顆粒相的計(jì)算采用穩(wěn)態(tài)藕合求解；用顆粒隨機(jī)軌道模型計(jì)算顆粒軌道時(shí)，跟蹤計(jì)算顆粒沿軌道的質(zhì)量和動(dòng)量的變化，將這些物理量引入到隨后的氣相流場(chǎng)計(jì)算中，既考慮氣體對(duì)顆粒相的作用，又考慮顆粒對(duì)氣體的作用，交替求解顆粒相與氣相的控制方程，直到計(jì)算結(jié)果都達(dá)到收斂標(biāo)準(zhǔn)。

1.1.3 沖蝕模型

在模擬中所采用的DPM模型中沖蝕速率公式如下。

(3)

其中C(dp)為顆粒直徑的函數(shù)，設(shè)定為一常數(shù)1.8×e-9，α為顆粒對(duì)壁面的沖擊角，f(α)為沖擊角的函數(shù)，v是顆粒相對(duì)壁面的速度，b(v)是此相對(duì)速度的函數(shù)，在模擬過(guò)程中速度指數(shù)函數(shù)為常數(shù)2.6。

1.2 邊界條件與網(wǎng)格劃分

1.2.1 幾何結(jié)構(gòu)

變徑管的結(jié)構(gòu)尺寸為入口管道直徑取為D=100 mm，長(zhǎng)度為1 000 mm，出口管直徑取D=70 mm，長(zhǎng)度取為800 mm，變徑角度為θ=10°、20°、30°、45°、60°、70°、90°，坐標(biāo)原點(diǎn)取在變徑開(kāi)始中心，具體尺寸見(jiàn)圖1。

圖1 異徑管結(jié)構(gòu)尺寸與橫截面網(wǎng)格劃分

1.2.2 入口邊界條件

氣相入口速度取v=10 m/s，入口ρ(顆粒)=1、5、10、15、20 g/m3。顆粒相入口邊界條件為將顆粒入口處的射流源設(shè)為面源，顆粒均勻地分布在整個(gè)入口截面的網(wǎng)格上，由每一個(gè)網(wǎng)格中心射入，設(shè)定顆粒的入口速度與氣相入口速度相同。計(jì)算過(guò)程中顆粒無(wú)質(zhì)量變化，無(wú)熱量傳遞，所以只設(shè)定顆粒密度和粒徑。模擬粉料選用單一粒徑煤粉，顆粒堆積密度為973.43 kg/m3，顆粒粒徑為80 μm。

1.2.3 出口與壁面邊界條件

氣相出口邊界條件按充分發(fā)展的管流條件處理，所有變量在出口截面處軸向梯度?D/?z為零。顆粒相出口邊界條件為完全逃逸。氣相流場(chǎng)在壁面采用無(wú)滑移邊界條件，近壁網(wǎng)格采用標(biāo)準(zhǔn)壁面函數(shù)近似處理。顆粒對(duì)壁面的碰撞恢復(fù)系數(shù)取為完全反彈。

1.2.4 網(wǎng)格劃分

氣相流場(chǎng)選用SIMPLE算法求解。計(jì)算網(wǎng)格采用結(jié)構(gòu)化網(wǎng)格，對(duì)應(yīng)網(wǎng)格節(jié)點(diǎn)數(shù)分別為294 000～350 000個(gè)，橫截面網(wǎng)格具體情況見(jiàn)圖1。

2 結(jié)果與討論

2.1 尺寸對(duì)沖蝕情況的影響

2.1.1 對(duì)沖蝕部位的影響

6種尺寸異徑管在入口ρ(顆粒)=15 g/m3時(shí)的沖蝕嚴(yán)重位置圖見(jiàn)圖2。由圖2可以看出隨尺寸變化的規(guī)律為：在小角度時(shí)，沖蝕嚴(yán)重區(qū)域比較均勻，其形狀為斑點(diǎn)狀，分布在0.1～0.4 L之間壁面區(qū)域；隨角度的上升，沖蝕區(qū)域(20°)向0.1～0.5 L之間擴(kuò)展，嚴(yán)重區(qū)連成片狀；隨角度繼續(xù)上升，嚴(yán)重沖蝕區(qū)域片狀形狀更加明顯，嚴(yán)重區(qū)域集中在出口管與變徑區(qū)域相連處，如在60°時(shí)，只在0.1～0.2 L之間且?guī)缀踹B成環(huán)狀,在90°時(shí)，嚴(yán)重沖蝕區(qū)域完全連成環(huán)狀，其寬度為0.1～0.5 L。(注：嚴(yán)重沖蝕區(qū)域定義為最大沖蝕率的50%～100%)

圖2 沖蝕嚴(yán)重區(qū)域圖[ρi(顆粒)=15 g/m3]

2.1.2 對(duì)沖蝕最大值的影響

對(duì)7種尺寸異徑管的不同入口ρ(顆粒)進(jìn)行對(duì)比研究，結(jié)果見(jiàn)圖3。從圖3中可知：在低濃度時(shí)(1 g/m3),異徑管沖蝕率隨角度上升而緩慢上升，不同尺寸異徑管的沖蝕率相差不大；隨著入口ρ(顆粒)上升，曲線有顯著變化，入口ρ(顆粒)為 10 g/m3時(shí)，隨著尺寸的變化，沖蝕率先快速上升，后在30°時(shí)上升幅度變緩慢，在70°之后沖蝕率開(kāi)始下降，并在90°時(shí)沖蝕率下降為8×10-8kg/m2·s，小于70°時(shí)的1.13×10-7kg/m2·s，也小于45°時(shí)的8.75×10-8kg/m2·s,且隨入口ρ(顆粒)的上升，此種趨勢(shì)變得更加明顯。這是因?yàn)轭w粒撞擊壁面角度隨變徑區(qū)域角度上升而增大，而壁面沖蝕嚴(yán)重程度隨入射角度的上升的趨勢(shì)為先上升后下降，并存在最嚴(yán)重的撞擊角度[11]。同時(shí)入口ρ(顆粒)越大，差異累積越嚴(yán)重，故在入口ρ(顆粒)為20 g/m3時(shí)，曲線先上升后下降的趨勢(shì)最明顯。此外流場(chǎng)在變徑區(qū)域轉(zhuǎn)向，進(jìn)而帶動(dòng)顆粒轉(zhuǎn)向，導(dǎo)致顆粒與壁面的實(shí)際撞擊角度并不是異徑管的變徑角度。

角度/(°)圖3 沖蝕率隨角度的變化(v=10 m/s)

2.2 入口ρ(顆粒)對(duì)不同尺寸異徑管的沖蝕影響

為分析入口ρ(顆粒)與異徑管沖蝕嚴(yán)重程度之間的關(guān)系，選定入口速度為10 m/s,ρ(顆粒)分別為1、5、10、15、20 g/m3，并對(duì)異徑管沖蝕嚴(yán)重的部位進(jìn)行取點(diǎn)研究，結(jié)果見(jiàn)圖4。

ρ(顆粒)/(g·m-3)圖4 沖蝕率隨入口濃度的變化(v=10 m/s)

由圖4可以看出：沖蝕率隨著入口ρ(顆粒)的上升而增加，在7種不同尺寸異徑管中，沖蝕率的上升都呈現(xiàn)線性上升的趨勢(shì)。這是因?yàn)樵诠潭ǖ娜肟谒俣认骂w粒的能量以及與壁面碰撞的次數(shù)已經(jīng)確定，ρ(顆粒)的上升只會(huì)增加入射顆粒的數(shù)量，顆粒數(shù)量只與沖蝕率成線性的關(guān)系，故出現(xiàn)沖蝕率隨ρ(顆粒)增加而線性增長(zhǎng)的現(xiàn)象。

此外在低濃度時(shí)，沖擊嚴(yán)重程度主要由變徑角度與顆粒數(shù)量來(lái)決定，可知存在一最嚴(yán)重沖蝕的變徑角度。但隨著ρ(顆粒)上升，加之在大角度時(shí)，顆粒在變徑區(qū)域的ρ(顆粒)將進(jìn)一步加大，ρ(顆粒)的增大導(dǎo)致流體湍動(dòng)能增大，進(jìn)而增大顆粒對(duì)壁面的撞擊次數(shù)，出現(xiàn)70°尺寸異徑管在高濃度20 g/m3時(shí)沖蝕率大于60°尺寸異徑管的現(xiàn)象。

3 異徑管使用、壽命及檢測(cè)意見(jiàn)

實(shí)際工業(yè)中，在滿足生產(chǎn)技術(shù)能力與工藝要求的前提下，應(yīng)盡量使用小角度的異徑管，以減小其嚴(yán)重沖蝕。在特殊情況下，需采用大角度的異徑管時(shí)，應(yīng)盡量避開(kāi)最嚴(yán)重沖蝕變徑角度。

此外在進(jìn)行異徑管的壽命評(píng)估時(shí),主要根據(jù)異徑管的沖蝕影響下凹坑的深度隨時(shí)間的變化，并根據(jù)深度情況進(jìn)行評(píng)估，故可以采用以下公式進(jìn)行：

L深=k×E/ρ×n

其中：L深為沖蝕深度，m；E為沖蝕速率，kg/m2·s；ρ為異徑管材料密度，kg/m3；n為異徑管使用時(shí)間，s；k為由材料確定的待定系數(shù)。

對(duì)工業(yè)中某特定角度的異徑管進(jìn)行檢測(cè)時(shí)，應(yīng)重點(diǎn)對(duì)靠近出口側(cè)的變徑區(qū)域壁面進(jìn)行檢測(cè)；同時(shí)隨ρ(顆粒)的上升，重點(diǎn)測(cè)量的范圍應(yīng)向入口側(cè)延伸；必要時(shí)可先進(jìn)行沖蝕模擬，預(yù)測(cè)沖蝕嚴(yán)重位置，為壁厚測(cè)量提供前期準(zhǔn)備，再進(jìn)行實(shí)際測(cè)量以期提高異徑管檢測(cè)的準(zhǔn)確性。

4 結(jié) 論

(1) 異徑管的沖蝕區(qū)域主要位于變徑區(qū)域壁面，嚴(yán)重的沖蝕部位位于靠近出口的變徑壁面。

(2) 異徑管的嚴(yán)重沖蝕率隨入口濃度的上升線性上升，隨角度上升而呈現(xiàn)先增大后減小的趨勢(shì)。且嚴(yán)重沖蝕區(qū)域的形狀從小角度的斑點(diǎn)狀向大角度的環(huán)狀變化。

(3) 工業(yè)使用異徑管時(shí)應(yīng)避開(kāi)最嚴(yán)重沖蝕變徑角，在進(jìn)行異徑管的檢測(cè)時(shí)，應(yīng)對(duì)靠近出口側(cè)的變徑壁面加密測(cè)點(diǎn)。同時(shí)可減小入口ρ(顆粒)和加厚變徑區(qū)域壁厚的方式減弱異徑管變徑區(qū)域的沖蝕破壞作用。此外根據(jù)作者研究得出的異徑管沖蝕速率可以方便、快捷的預(yù)測(cè)出管道的壽命。

[ 參 考 文 獻(xiàn) ]

[1] 梁志正，等.重整加氫聯(lián)合裝置氛腐蝕及其防護(hù)措施[J].石油化工腐蝕與防護(hù),1995,16(1):24-26.

[2] 陳孫藝，等.異徑管三維實(shí)體模型的計(jì)算機(jī)建模方法[J].工程圖學(xué)學(xué)報(bào),2005(2)：137-141.

[3] 鐘文學(xué).談?wù)勁涔茉O(shè)計(jì)的誤區(qū)[J].化工設(shè)備與管道,2001(4)：33-35.

[4] 陳孫藝，等.異徑管的標(biāo)準(zhǔn)及其失效分析[J].化工設(shè)備與管道,2005,6(45)： 25-28.

[5] 劉一民譯.工廠制造的鍛制鋼對(duì)接焊管件[M].蘭州：蘭州石油機(jī)械研究所,1994：1.

[6] 栗雪勇.加氫裂化裝置煉高含硫油后的腐蝕與防腐[C].//第五屆壓力容器使用管理學(xué)術(shù)會(huì)議論文集.南京：加工高含硫油設(shè)備和管道防腐蝕技術(shù)研討會(huì)，2002：42.

[7] 沈中華.乙二醇裝置六效蒸發(fā)系統(tǒng)腐蝕分析[J].揚(yáng)子石油化工,1999,14(3):10-13.

[8] 偶國(guó)富，等.異徑管沖蝕失效的流固耦合數(shù)值模擬[J].力學(xué)學(xué)報(bào)，2010，42(2)：197-204.

[9] 陳孫藝，等.異徑彎管的無(wú)力矩環(huán)向應(yīng)力解析解[J].壓力容器,2007,24(2)：35-39.

[10] 朱紅鈞，等.Fluent12流體分析及工程仿真[M].北京:清華大學(xué)出版社,2011:41-42.

[11] Finnie I.The mechainism of erosion of ductile metals[M].Newyork:Proc,3Td US Nation,conf,on Appli,mesh,ASME,1958:527-532.