呂 超，劉艷龍，陳緒鑫，殷宏鑫，孫銘赫

（東北大學(xué)秦皇島分?？刂乒こ虒W(xué)院，河北 秦皇島 066004）

氫氣作為能源，具有能量密度大、無(wú)毒無(wú)味無(wú)腐蝕、便于存儲(chǔ)和運(yùn)輸及零碳排放的優(yōu)點(diǎn)，因此具有很廣闊的應(yīng)用前景［1］。2019年3月，氫能首次被寫(xiě)進(jìn)我國(guó)《政府工作報(bào)告》之中［2］，這表明我國(guó)對(duì)于氫能非常重視。氫能的有效利用，氫氣的運(yùn)輸是重中之重，其中管道運(yùn)輸具有運(yùn)載量大、效率高、經(jīng)濟(jì)效益高等優(yōu)勢(shì)，因此具有良好的市場(chǎng)前景［3］。

在氫氣的管道運(yùn)輸中，如果管道內(nèi)部有污染物堆積或管道存在較大的加工誤差，則會(huì)造成管道內(nèi)壁形成凸起。氫氣管道內(nèi)壁的粗糙程度對(duì)氫氣的流場(chǎng)具有較大的影響，尤其是在直角彎管處，如果管道較為粗糙，就會(huì)加劇氫氣輸送的不穩(wěn)定性，造成安全隱患，而管道內(nèi)壁形成的凸起是造成管道內(nèi)部粗糙程度高的重要原因之一，因此研究氫氣管道內(nèi)部凸起對(duì)流場(chǎng)的影響非常必要。在管道受損分析方面，單會(huì)娜［4］對(duì)管道產(chǎn)生缺陷的原因進(jìn)行了分析與檢測(cè)，提出不銹鋼管道的焊縫表面易受到特定介質(zhì)的腐蝕。Pluvinage［5］提出了利用缺口修正失效評(píng)估圖對(duì)氫氣管道缺陷進(jìn)行評(píng)定的方法。趙博鑫等［6］對(duì)氫氣管道的鋼管和材料進(jìn)行了適應(yīng)性分析，指出管道焊縫內(nèi)部有夾雜物等缺陷極易造成氫腐蝕而破壞管道。李強(qiáng)林等［7］對(duì)氣體管道內(nèi)部凹坑缺陷進(jìn)行了分析與評(píng)價(jià)，分析了有關(guān)氣體管道出現(xiàn)凹坑缺陷的原因?？梢?jiàn)在氣體輸送管道方面，國(guó)內(nèi)外均著重研究管道的缺陷問(wèn)題，而對(duì)于管道內(nèi)部凸起則研究較少，尤其是微小凸起對(duì)管道的影響的流場(chǎng)分析。從文獻(xiàn)［7］可知，氣體管道不僅會(huì)有凹坑缺陷，在表面也存在凸起，雖然這些凸起形狀各異，但根據(jù)平滑程度大致近似于錐形、圓柱形以及圓柱與球冠的組合體等形狀。本文運(yùn)用Fluent 模擬氫氣流動(dòng)時(shí)管道內(nèi)壁不同形狀的凸起對(duì)管內(nèi)流場(chǎng)的影響。

1 模型建立

1.1 整體物理模型

以歐洲大規(guī)模使用的氫氣管道為研究對(duì)象，結(jié)合我國(guó)實(shí)際情況，設(shè)計(jì)氫氣管道的尺寸為：管道直徑600 mm，直角彎管兩邊接口的直線管道長(zhǎng)度均為6 000 mm，呈90°分布，其中心截面均在z=0平面上，3種凸起均位于直角彎管凸側(cè)內(nèi)壁。運(yùn)用SCDM軟件進(jìn)行建模，z=0的直角彎管二維截面如圖1所示。

圖1 z=0的直角管道二維截面Fig.1 Two-dimensional section diagram of right angle pipe with z=0

1.2 凸起處的微小模型

3 種典型凸起均為對(duì)稱形狀，其對(duì)稱軸均通過(guò)直角彎管圓心。凸起的二維截面如圖2所示。凸起的底部寬度均為18 mm，最大高度為11.95 mm，且在管壁的同一位置，均對(duì)稱于同一對(duì)稱軸。數(shù)字1、2、3、4、5、6、7代表的直線是本次研究使用Fluent模擬設(shè)定的監(jiān)測(cè)線，將它們分別命名為監(jiān)測(cè)線1～監(jiān)測(cè)線7。在數(shù)字8 代表的曲線上按照等距原則設(shè)定20 個(gè)監(jiān)測(cè)點(diǎn)。利用這些監(jiān)測(cè)線、監(jiān)測(cè)點(diǎn)獲取凸起上各點(diǎn)的壓力值和速度大小。

圖2 典型凸起的二維截面Fig.2 Two-dimensional cross-section diagram of typical bulges

2 網(wǎng)格劃分及模擬參數(shù)設(shè)置

2.1 網(wǎng)格劃分

為減少計(jì)算量，提高計(jì)算精度，對(duì)原來(lái)的幾何模型進(jìn)行簡(jiǎn)化，取原始管道模型的中軸面建立二維管道模型。在劃分網(wǎng)格時(shí)采用四邊形的結(jié)構(gòu)化網(wǎng)格劃分方法，劃分完成過(guò)后將其轉(zhuǎn)化為非結(jié)構(gòu)網(wǎng)格導(dǎo)入Fluent求解器進(jìn)行計(jì)算。

利用ICEM 所劃分出的網(wǎng)格（取凸起附近的一部分管道示意）如圖3所示。圖3（a）的網(wǎng)格數(shù)量約為65萬(wàn)，圖3（b）和圖3（c）網(wǎng)格數(shù)量約為70萬(wàn)，經(jīng)檢驗(yàn)網(wǎng)格質(zhì)量均在0～1 之間且平均網(wǎng)格質(zhì)量在0.9以上。

圖3 幾何模型的結(jié)構(gòu)化網(wǎng)格劃分Fig.3 Structured mesh division of geometric model

2.2 模擬參數(shù)設(shè)置

2.2.1 假設(shè)條件

忽略氫氣流動(dòng)對(duì)管道幾何尺寸造成的微變、氫氣流過(guò)凸起表面及其附近對(duì)管道內(nèi)的溫度變化、氫氣管道內(nèi)的密度變化、氫氣管道內(nèi)的壁厚。

2.2.2 模擬參數(shù)設(shè)置

本次模擬的氫氣流速為8 m/s，在圓形管內(nèi)流動(dòng)，特征長(zhǎng)度L取圓管的直徑d。推導(dǎo)得圓管內(nèi)流動(dòng)的雷諾數(shù)計(jì)算公式為［8］

式中：V為截面氫氣的平均速度；v為流體的運(yùn)動(dòng)黏度。

氫氣的運(yùn)動(dòng)黏度為241×10-6m2/s，圓管直徑為0.6 m，代入式（1）得出氫氣流動(dòng)的雷諾數(shù)Re=19 917＞2 320，管道內(nèi)的風(fēng)流狀態(tài)為湍流。因此選擇標(biāo)準(zhǔn)k-ε模型計(jì)算流體，常數(shù)取值見(jiàn)表1。壓力-速度耦合算法選擇SIMPLE算法。

表1 標(biāo)準(zhǔn)k-ε特定常數(shù)取值Tab.1 Standard k-ε specific constant values

氫氣管道中的溫度為300 K，F(xiàn)luent 初始?xì)鈮簽?.013 25 kPa，其他參數(shù)設(shè)置見(jiàn)表2。

表2 Fluent模型主要參數(shù)設(shè)置Tab.2 Main parameter settings of Fluent model

3 結(jié)果與討論

3.1 壓力分布

直角管道凸起處的壓力大小與分布影響管道的強(qiáng)度和使用壽命，某部分壓力過(guò)大或者分布不均可能會(huì)造成細(xì)微裂縫，嚴(yán)重影響輸送安全。因此對(duì)凸起表面壓力分布進(jìn)行了模擬仿真，氫氣管道直角彎管凸起處的壓力分布云圖如圖4所示。

根據(jù)模擬過(guò)程中凸起不同位置（Y坐標(biāo)值）的壓力的變化，繪制變化曲線，如圖5所示。

分析圖4和圖5可以得出結(jié)論：凸起附近較大的壓力場(chǎng)均在右側(cè)靠近監(jiān)測(cè)線的區(qū)域，且錐形凸起附近的壓力峰值最大，圓柱形凸起和組合體凸起的壓力峰值幾乎相等。隨著縱坐標(biāo)的增大，凸起左側(cè)輪廓線壓力逐漸增大，但壓力普遍較小，且均為負(fù)壓，凸起右側(cè)輪廓線上壓力會(huì)先上升然后驟降，使得此曲線的上頂點(diǎn)處附近壓力值較小，監(jiān)測(cè)線3 的平均壓力最小，圓柱形凸起頂部壓力先下降而后突然上升，然后趨于穩(wěn)定，曲線8 上壓力值大小呈現(xiàn)出拋物線形式的變化，開(kāi)口向上，壓力先降后升，且曲線末端（縱坐標(biāo)最大位置）壓力值小于始端。

圖4 直角彎管凸起附近壓力分布Fig.4 Pressure distribution near the bulges of right-angle bend pipe

圖5 3種凸起不同監(jiān)測(cè)位置的壓力變化Fig.5 Pressure changes at different monitoring positions of three bulges

3.2 速度分布

氣體對(duì)凸起表面的沖擊也是必須要考慮的因素，這涉及氫氣流動(dòng)時(shí)的穩(wěn)定性，能直接反映沖擊情況的便是速度大小分布與方向，管道直角彎管處內(nèi)壁凸起附近氣體流速云圖如圖6所示。氣體速度分布情況如圖7所示。

分析圖6和圖7可知：凸起表面的氫氣流動(dòng)速度均非常低，且在凸起左右兩側(cè)靠近管壁處形成低速區(qū)域，背風(fēng)側(cè)形成速度死區(qū)，組合體形凸起速度死區(qū)最??；氫氣流動(dòng)方向一開(kāi)始平行于管壁，后逐漸變化，到錐形凸起頂點(diǎn)或圓柱體凸起和組合體凸起右側(cè)頂點(diǎn)時(shí)與管壁大約呈60°夾角，隨著繼續(xù)流動(dòng)又逐漸平行于管壁。

圖6 直角彎管凸起附近速度分布Fig.6 Velocity distribution near the bulges of right-angle bend pipe

圖7 凸起附近速度矢量分布Fig.7 Distribution of velocity vector near bulges

3.3 剪切力分布

剪切力的大小體現(xiàn)了流體對(duì)管道內(nèi)部的沖擊情況，本文根據(jù)模擬過(guò)程中凸起不同位置對(duì)剪切力變化的影響，繪制了變化曲線，如圖8所示。

圖8 3種凸起不同監(jiān)測(cè)位置的剪切力變化Fig.8 Shear force changes at different monitoring positions of the three bulges

分析圖8可知，凸起迎風(fēng)側(cè)剪切力變化規(guī)律是：自下向上剪切力逐漸增大，在中上段驟增，使迎風(fēng)側(cè)頂點(diǎn)處剪切力較大；圓柱形凸起頂部隨著Y坐標(biāo)的增大，剪切力最初驟降，后趨于平緩；而較為平緩的組合體凸起頂部曲面上，剪切力隨著縱坐標(biāo)的增大呈波浪式變化，總體是下降趨勢(shì)，變化的劇烈程度較小。

4 結(jié)論

（1）錐形凸起相較于其他2種凸起壓力峰值最大，圓柱形凸起壓力峰值最低；凸起迎風(fēng)側(cè)附近區(qū)域形成較大的壓力場(chǎng)，迎風(fēng)側(cè)頂部壓力較?。煌蛊鹎骓敳繅毫o(wú)突變現(xiàn)象，斜平面頂部壓力先突變后趨于平緩；凸起背風(fēng)側(cè)壓力隨著位置改變變化不明顯。

（2）3 種凸起表面相近區(qū)域氫氣流動(dòng)速度差幾乎相同，凸起表面氫氣流動(dòng)速度較低并在凸起兩側(cè)形成低速區(qū)域，背風(fēng)側(cè)區(qū)域面積較大，形成速度死區(qū)，組合體形凸起速度死區(qū)最?。粴錃饬鬟^(guò)時(shí)速度方向與管壁之間的夾角變化是從0°到約60°再到0°。

（3）剪切力分析結(jié)果表明：氫氣對(duì)迎風(fēng)側(cè)頂部沖擊較大，其他段無(wú)明顯變化，組合體形凸起相對(duì)于其他2 種凸起剪切力峰值較小，流體對(duì)其沖擊力度較小；頂部形狀變化越平緩，流體對(duì)其沖擊越平穩(wěn)；背風(fēng)側(cè)幾乎無(wú)沖擊。

（4）圓柱形凸起壓力峰值小，而迎風(fēng)側(cè)頂部受到的沖擊較大，使得其對(duì)氫氣流動(dòng)阻礙作用較小而且容易被流體沖蝕，使凸起形狀發(fā)生改變，更有利于流體流動(dòng)，產(chǎn)生良性循環(huán)。