周云龍，張　祎

(東北電力大學(xué) 能源與動(dòng)力工程學(xué)院，吉林 吉林 132012)

?

CFD技術(shù)在固液兩相流漿體水擊壓強(qiáng)計(jì)算中的應(yīng)用

周云龍，張祎

(東北電力大學(xué) 能源與動(dòng)力工程學(xué)院，吉林 吉林 132012)

摘要：采用計(jì)算流體力學(xué)(CFD)動(dòng)網(wǎng)格方法的彈性光順?lè)ㄅc局部網(wǎng)格重劃法對(duì)閥門關(guān)閉過(guò)程進(jìn)行了動(dòng)態(tài)模擬，對(duì)閥門突然關(guān)閉引起的固液兩相流漿體水擊壓強(qiáng)進(jìn)行了數(shù)值計(jì)算。通過(guò)改變閥門關(guān)閉速度分析了閥門關(guān)閉時(shí)間對(duì)水擊壓強(qiáng)產(chǎn)生的影響進(jìn)行了分析，并比較了不同初始流速對(duì)水擊壓強(qiáng)的影響。結(jié)果表明：延長(zhǎng)閥門關(guān)閉時(shí)間將減小漿體水擊壓強(qiáng)的最大值，關(guān)閉時(shí)間超過(guò)2s后減小效果可忽略不計(jì)；閥門關(guān)閉時(shí)間一定時(shí)，初始流速越大漿體水擊壓強(qiáng)最大值也越大。

關(guān)鍵詞：固液兩相流；漿體水擊；動(dòng)網(wǎng)格；數(shù)值模擬

壓管道中閥門的快速關(guān)閉(開啟)以及泵等設(shè)備的突然啟動(dòng)(停止)會(huì)導(dǎo)致管道中流場(chǎng)的急劇變化，從而引起管道中某個(gè)截面上介質(zhì)的流速、流向的激烈變化，這種發(fā)生在管道中的水力瞬變現(xiàn)象稱為水擊現(xiàn)象，是管道中最常見的水力瞬變現(xiàn)象之一。伴隨水擊現(xiàn)象的產(chǎn)生，管路中的壓強(qiáng)可能升高到正常工作壓強(qiáng)的幾倍至幾十倍。對(duì)水擊現(xiàn)象的處理不當(dāng)輕則會(huì)引起噪音、管路振動(dòng)，重則會(huì)損壞管路元件，甚至引起管道破裂、跑漿跑水等嚴(yán)重水力事故，因此在管路系統(tǒng)的設(shè)計(jì)與安裝過(guò)程對(duì)水擊壓強(qiáng)的計(jì)算必不可少[1-4]。

隨著近年來(lái)計(jì)算機(jī)技術(shù)的發(fā)展，水擊現(xiàn)象的研究方法逐漸由實(shí)驗(yàn)法、圖解法向數(shù)值模擬方法發(fā)展[5-12]。目前，對(duì)單相液體水擊現(xiàn)象的研究已經(jīng)比較成熟，但對(duì)固液兩相流漿體水擊的研究計(jì)算較為少見，并主要以理論分析為主。韓文亮[13]等推導(dǎo)出了偽均質(zhì)固液兩相流水擊現(xiàn)象的基本方程，建立偽均質(zhì)固液兩相流水擊的數(shù)值模型，并用特征線法對(duì)其進(jìn)行了分析計(jì)算；周云龍[14-15]等對(duì)長(zhǎng)距離漿體輸送管道停泵水錘進(jìn)行了分析研究，同時(shí)推導(dǎo)出了考慮含氣量的漿體水擊的基本方程，并運(yùn)用到了工程實(shí)際中；王偉英[16]等人推導(dǎo)出了含粗顆粒物料的固液兩相流漿體在垂直管道中的水擊波速公式，并與實(shí)驗(yàn)結(jié)果進(jìn)行了比較。

采用計(jì)算流體力學(xué)(CFD)中的動(dòng)網(wǎng)格方法定義管道閥門系統(tǒng)，可以對(duì)閥門運(yùn)動(dòng)的瞬態(tài)過(guò)程進(jìn)行模擬，同時(shí)能夠?qū)崿F(xiàn)流場(chǎng)的可視化，便于進(jìn)一步分析研究，因此越來(lái)越多的被運(yùn)用于水擊現(xiàn)象的模擬分析之中[17-20]。

本文運(yùn)用計(jì)算流體力學(xué)方法(CFD)對(duì)管道中閥門突然關(guān)閉引起的固液兩相流漿體水擊現(xiàn)象進(jìn)行了數(shù)值計(jì)算，并分析了閥門關(guān)閉時(shí)間以及初始流速對(duì)漿體水擊壓強(qiáng)的影響。

圖1　簡(jiǎn)化后的管道幾何模型

圖2　網(wǎng)格劃分及閥門運(yùn)動(dòng)過(guò)程

2研究對(duì)象與基本參數(shù)

本文研究的對(duì)象是上游為恒定水頭，下游為閘閥的管道系統(tǒng)，圖1為簡(jiǎn)化后的管道幾何模型。管道直徑為0.1 m，閥芯前管道長(zhǎng)度為30 m，閥芯后管道長(zhǎng)度為2 m。管路上游恒定壓強(qiáng)為70 560 Pa，下游出口壓力為標(biāo)準(zhǔn)大氣壓101.3 kPa。管路中的流體介質(zhì)為水與灰分混合的偽均質(zhì)流漿體，漿體中灰分的體積濃度Cv為9.09%，其中水的密度ρw為1×103kg/m3，灰分密度ρs為2×103kg/m3，漿體的密度由公式(1)計(jì)算[14]：

ρm=ρsCv+ρw(1-Cv) .

(1)

將圖1中的流體區(qū)域作為計(jì)算域，利用Gambit軟件對(duì)該計(jì)算域進(jìn)行非結(jié)構(gòu)網(wǎng)格劃分，并對(duì)閥體以及閥門附近部分管道進(jìn)行了局部細(xì)化。采用動(dòng)網(wǎng)格中的彈性光順?lè)ㄅc局部網(wǎng)格重劃法實(shí)現(xiàn)閥門關(guān)閉的動(dòng)態(tài)過(guò)程，閥體與上下游管道接觸面設(shè)置為交界面(interface)，接觸部分允許流體流過(guò)，不接觸部分設(shè)置為固體壁面(wall)，不允許流體流過(guò)，閥體以及閥門附近網(wǎng)格劃分以及閥體運(yùn)動(dòng)過(guò)程如圖2所示。

3數(shù)值模擬計(jì)算與分析

3.1漿體水擊控制方程

描述漿體水擊現(xiàn)象的控制方程為由守恒定律推導(dǎo)出來(lái)的連續(xù)性方程與運(yùn)動(dòng)方程[13，20]：

(1)

(2)

式中:am為漿體水擊波速，ρm為漿體密度。

3.2邊界條件與解析條件

本文使用Fluent軟件進(jìn)行數(shù)值模擬計(jì)算，進(jìn)出口邊界條件分別定義為壓力入口及壓力出口，采用壓力基非穩(wěn)態(tài)求解器求解。閥門部分為運(yùn)動(dòng)邊界，運(yùn)動(dòng)規(guī)律為勻速關(guān)閉。

3.3模擬計(jì)算結(jié)果與分析

圖3-圖5分別為入口初始流速為4.5 m/s，閥門1 s關(guān)閉時(shí)水擊現(xiàn)象發(fā)生后不同閥門開度時(shí)管道中軸向速度變化云圖?？梢钥闯鲩y門開始關(guān)閉后，由于閥體的作用，閥門處的流體速度與流向會(huì)發(fā)生劇烈變化，緊貼閥門處流體流速由于閥體的作用首先開始滯止，這將引發(fā)水擊現(xiàn)象。

圖3 開度為0.2圖4 開度為0.5

圖5　開度為0.8

圖6　不同閥門關(guān)閉時(shí)間水擊壓強(qiáng)變化

圖7　不同閥門關(guān)閉時(shí)間水擊壓強(qiáng)最大值

圖8　不同初始流速水擊壓強(qiáng)變化

3.3.1閥門關(guān)閉時(shí)間對(duì)漿體水擊現(xiàn)象的影響

圖6為同一模型不同閥門關(guān)閉時(shí)間下水擊現(xiàn)象發(fā)生后管道閥門處水擊壓強(qiáng)的變化情況。

從圖中可以看出，改變閥門關(guān)閉時(shí)間將對(duì)水擊壓強(qiáng)的最大值產(chǎn)生巨大影響。這是因?yàn)殚y門的關(guān)閉過(guò)程并不是瞬時(shí)完成的，而是需要一個(gè)關(guān)閉過(guò)程，而水擊現(xiàn)象發(fā)生后水擊壓強(qiáng)將以水擊波的形式向管道上游傳播，當(dāng)水擊波傳播到閥門上游管路元件或水池時(shí)，將會(huì)發(fā)生反射，形成一個(gè)減壓波向反方向傳播。

當(dāng)減壓波到達(dá)閥門位置時(shí)，閥門并沒(méi)有完全關(guān)閉(間接水擊)，減壓波所附加的壓力降將使水擊現(xiàn)象引起的壓強(qiáng)躍升值變小。閥門關(guān)閉時(shí)間越長(zhǎng)，單位時(shí)間內(nèi)閥門的開度變化越小，反射回來(lái)的減壓波抵消的壓強(qiáng)份額越大，實(shí)際水擊壓強(qiáng)越小。當(dāng)閥門在反射的減壓波到達(dá)閥門位置之前已完全關(guān)閉(直接水擊)，水擊現(xiàn)象將不受到減壓波的影響，這樣的水擊壓強(qiáng)最大值較大。

圖7為閥門處水擊壓強(qiáng)最大值隨閥門關(guān)閉時(shí)間變化的情況?？梢钥闯?，隨著閥門關(guān)閉時(shí)間的延長(zhǎng)，水擊壓強(qiáng)的最大值將急劇下降，但當(dāng)閥門關(guān)閉時(shí)間超過(guò)2 s后，閥門關(guān)閉時(shí)間對(duì)水擊壓強(qiáng)最大值的影響可忽略不計(jì)。在工程實(shí)際中可以合理設(shè)置閥門關(guān)閉時(shí)間，避免直接水擊現(xiàn)象的發(fā)生，減小水擊現(xiàn)象的危害。

3.3.2初始流速對(duì)漿體水擊壓強(qiáng)的影響

本文對(duì)閥門關(guān)閉時(shí)間為1s的同一模型，不同初始流速的水擊現(xiàn)象進(jìn)行了數(shù)值模擬計(jì)算。圖8為不同初始流速閥門處水擊壓強(qiáng)的變化情況，可以看出隨著初始流速的增大，閥門處的水擊壓強(qiáng)也隨之增大。

這是因?yàn)樗畵衄F(xiàn)象發(fā)生后，由于閥體的作用閥門處的流體速度將發(fā)生劇烈變化，這種流體流速、流向的改變將在管道中附加水擊壓強(qiáng)，即水擊壓強(qiáng)的實(shí)質(zhì)是管道中流體動(dòng)能部分或全部變化為勢(shì)能時(shí)引起的壓強(qiáng)變化，水擊壓力即是流體的慣性力，因此改變管道初始流速會(huì)對(duì)水擊壓強(qiáng)產(chǎn)生影響。

4結(jié)論

本文采用計(jì)算流體力學(xué)中的動(dòng)網(wǎng)格方法對(duì)閥門[21-22]關(guān)閉引起的有壓管道漿體水擊現(xiàn)象進(jìn)行了動(dòng)態(tài)模擬，并分析了閥門關(guān)閉時(shí)間以及初始流速對(duì)漿體水擊壓強(qiáng)的影響，得出了以下結(jié)論：

(1)閥門關(guān)閉時(shí)間會(huì)對(duì)水擊壓強(qiáng)的最大值產(chǎn)生影響，水擊壓強(qiáng)的最大將隨著閥門關(guān)閉時(shí)間的延長(zhǎng)而減小，但關(guān)閉時(shí)間超過(guò)2 s后，閥門關(guān)閉時(shí)間對(duì)漿體水擊壓強(qiáng)的影響可以忽略不計(jì)。

(2)在閥門關(guān)閉時(shí)間一定時(shí)，由于水擊現(xiàn)象附加壓力的實(shí)質(zhì)為管道中流體的慣性力，因此初始流速會(huì)對(duì)漿體水擊壓強(qiáng)產(chǎn)生影響，初始流速越快漿體水擊壓強(qiáng)的最大值越大。

參考文獻(xiàn)

[1]王學(xué)芳，葉宏開.工業(yè)管道中的水錘[M].北京： 科學(xué)出版社，1991.

[2]劉光臨，蔣勁，沈正.陽(yáng)邏電廠冷卻水系統(tǒng)水錘事故分析[J].華中電力，1999，12(2)：24-26.

[3]賈佑?xùn)|，安增琴，李萍，等.一起停泵水錘事故分析及其防止[J].萊鋼科技，2007(3)：74-75.

[4]秋元德三.水撃作用と圧力脈動(dòng)[M].日本：日本工業(yè)新聞社，1988.

[5]楊玲霞，李樹慧，侯詠梅，等.水擊基本方程的改進(jìn)[J].水利學(xué)報(bào)，2007，38(8)：948-952.

[6]Wylie E B，Streeter V L.Fluid Transients[M].New York:McGraw-Hill，1978.

[7]劉曉平，陳燕國(guó)，孫東坡，等.管道水擊數(shù)值模擬方法分析[J].黃河水利職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2007，19(3)：7-11.

[8]孫新嶺，雍歧衛(wèi)，蔣明，等.有限元法在管道水擊計(jì)算中的應(yīng)用探索[J].后勤工程學(xué)院學(xué)報(bào)，2006(1)：67-69.

[9]A.Imre，I.Bama，G，et al.Theoretical study of flashing and water hammer in a supercritical water cycle during pressure drop[J].Nuclear Engineering and Design，2010，240(6)：1569-1574.

[10] 唐曉寅，葉江濤.長(zhǎng)管水擊最大水擊壓強(qiáng)的解析[J].后勤工程學(xué)院學(xué)報(bào)，2004(2)：32-36.

[11] M.S.Ghidaoui，M.Zhao，D.A.McInnis，etal.A review of water hammer theory and practice [J].Applied Mechanics Reviews，2005，58(1/6):49.

[12] 侯詠梅.水擊理論與計(jì)算研究[D].鄭州：鄭州大學(xué)，2003.

[13] 韓文亮，柴宏恩，韓軍.偽均質(zhì)固液兩相流水擊的數(shù)值模擬[J].有色金屬，2000，52(1):37-41.

[14] 周云龍，李書芳，辛凱.長(zhǎng)輸漿體管路停泵水錘的研究[J].流體機(jī)械，2013，41(8)：26-30.

[15] 周云龍，孫斌，段曉寧等.復(fù)合管道內(nèi)液-固兩相流漿體水擊壓強(qiáng)計(jì)算[J].2004，25(2)：251-254.

[16] 王英偉，曹華德，曹斌等.含粗顆粒固體物料在垂直管流中水擊規(guī)律[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2014，22(5)：909-915.

[17] 華曄，廖偉麗.CFD技術(shù)在管道閥門水擊計(jì)算中的應(yīng)用[J].電網(wǎng)與清潔能源，2009，25(3):72-75.

[18] 劉華坪，陳浮，馬波.基于動(dòng)網(wǎng)格與UDF技術(shù)的閥門流場(chǎng)數(shù)值模擬[J].汽輪機(jī)技術(shù)，2009，50(2):106-108.

[19] 郭蘭蘭，耿介，石碩，等.基于UDF方法的閥門變速關(guān)閉過(guò)程中的水擊壓強(qiáng)計(jì)算研究[J].山東大學(xué)學(xué)報(bào)：理學(xué)版，2014，49(3)：27-30.

[20] 秦宏波，白曉寧，胡壽根.基于CFD的管道固液兩相流輸送的數(shù)值計(jì)算及實(shí)驗(yàn)對(duì)比[J].水力采煤與管道運(yùn)輸，2001，(2)：3-5.

[21] 韓文亮，柴宏恩，韓軍.偽均質(zhì)固液兩相流水擊的數(shù)值模擬：Ⅰ—理論[J].有色金屬，2000，52(1)：37-41

[22] 孫斌，劉彤.納米流體在內(nèi)置扭帶管的傳熱數(shù)值模擬[J].東北電力大學(xué)學(xué)報(bào)，2015，35(2)：10-17.

Application of CFD Method in Computing of Slurry Water Hammer on Liquid-Solid Tow-Phase Flow

ZHOU Yun-long，ZHANG Wei

(Energy Resource and Power Engineering College，Northeast Dianli University，Jilin Jilin 132012)

Abstract：By using spring-base smoothing technology and local remeshing technology of dynamic mesh in CFD method to calculate the pipe valve closing，and the slurry water hammer pressure in pipeline with slurry flows carrying solid particles caused by the suddenly close of the valve was calculated.By changing closing speed of the valve，the influence of valve’s closure time on slurry water hammer pressure were analyzed，and the influence of different initial velocity of the slurry on slurry water hammer were compared.Results shows that the increase of valve’s closure time will reduce the slurry water hammer’s max pressure on liquid-solid two-phase flow systems，the Influence is neglected when valve’s closure time is mors than 2 seconds；when valve closing time is fixed，the increase of velocity of the slurry will increase the slurry water hammer’s max pressure.

Key words：Liquid-solid tow-phase flow；Slurry water hammer；Dynamic mesh；Numerical simulation

中圖分類號(hào)：TB126

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1005-2992(2016)01-0056-04

作者簡(jiǎn)介：周云龍(1960-)，男，吉林省扶余市人，東北電力大學(xué)能源與動(dòng)力工程學(xué)院教授，博士，博士生導(dǎo)師，主要研究方向：液-固兩相流.

收稿日期：2015-12-10