傾斜塔板漏液的三維數(shù)值模擬及實(shí)驗(yàn)

張平 王學(xué)平 李雅俠

摘 ? ? ?要：塔體受到橫向載荷作用，板式塔的塔板傾斜，塔板上液層出現(xiàn)非等高。在冷模試驗(yàn)塔裝置中以空氣-水為實(shí)驗(yàn)介質(zhì)，對(duì)傾斜塔板的漏液性能進(jìn)行了實(shí)驗(yàn)研究。實(shí)驗(yàn)表明，塔板傾斜方向與液相流動(dòng)方向相同時(shí)，傾斜塔板漏液小于水平塔板，傾角越大漏液越少，塔板氣相操作下限減小;塔板傾斜方向與液相流動(dòng)方向相反時(shí)，傾斜塔板漏液大于水平塔板，傾角越大漏液越多，塔板氣相操作下限增大。利用計(jì)算流體力學(xué)的雙歐拉模型對(duì)傾斜塔板的漏液狀態(tài)進(jìn)行了三維數(shù)值模擬，模擬結(jié)果同實(shí)驗(yàn)結(jié)果吻合較好。模擬結(jié)果顯示，傾斜塔板在液層非等高剖面上發(fā)生不均勻漏液，傾角大的塔板在液層薄的區(qū)域發(fā)生氣相溝流現(xiàn)象。

關(guān) ?鍵 ?詞：計(jì)算流體力學(xué);漏液;數(shù)值模擬;氣含率;傾斜塔板

中圖分類號(hào)：TQ053.5 ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼： A ? ? ? 文章編號(hào)： 1671-0460（2019）12-2875-04

Abstract： The tray of the plate tower was inclined and the liquid layer on the tray was non-equal when the tower body was subjected to lateral load. Weeping of inclined trays was studied experimentally using air-water as experimental medium in cold model test tower. The results showed that the weeping of the inclined tray was smaller than that of the horizontal tray， the weeping of tray reduced with the increase of inclination angle and the lower limit of gas phase operation of tray decreased when the inclined direction of the tray was the same as the liquid flow direction. The weeping of the inclined tray was larger than that of the horizontal tray; the larger the inclination angle， the more the weeping， and the lower limit of gas phase operation of tray was increased when the inclination direction of the tray was opposite to the direction of liquid flow. The two-Euler model of computational fluid dynamics （CFD） was used to simulate the weeping of inclined trays. The simulation results were in good agreement with the experimental results. The simulation results indicated that uneven weeping occurred on non-contour profiles in the inclined tray. Gas phase channeling appeared in the region of thin liquid layer with large inclination tray.

Key words： Computational fluid dynamics; Weeping; Numerical simulation; Gas holdup; Inclined tray

板式塔在操作狀態(tài)時(shí)，塔板上存留一定高度的液層，塔體受到橫向載荷發(fā)生撓曲時(shí)，塔板上液層在液體流動(dòng)方向上出現(xiàn)不等高[1]。塔頂撓曲最大，液層非等高從塔底開始向上逐層加劇。張平[2]研究了這種液層不等高對(duì)塔板流場(chǎng)的影響，發(fā)現(xiàn)塔板回流區(qū)位置發(fā)生了改變。均勻氣相通過非等高的液層時(shí)，液層厚的區(qū)域比液層薄的區(qū)域更易發(fā)生漏液，造成了塔板沿液體流動(dòng)方向發(fā)生不均勻漏液。塔板漏液是一種導(dǎo)致塔板效率降低的縱向返混現(xiàn)象，實(shí)驗(yàn)證明，塔板漏液不均勻分布對(duì)塔板效率有很大影響[3]趙培等[4]用實(shí)驗(yàn)方法研究了十字旋閥塔板的漏液性能，獲得了十字旋閥塔板漏液點(diǎn)孔速的關(guān)聯(lián)式。劉靜[5]以泡罩立體篩板為研究對(duì)象得到漏液點(diǎn)氣速與堰高和液體流量的關(guān)系。在0.5m×0.5m的矩形篩孔塔板的漏液操作狀態(tài)的CFD模擬中，孫玉春[6]進(jìn)行了初步的嘗試，模擬結(jié)果與L&B方法[7]的計(jì)算結(jié)果基本吻合。Ali[8，9]利用歐拉-歐拉雙流體模型數(shù)值模擬了矩形塔板和圓形塔板的干板壓降和塔板漏液，模擬結(jié)果與實(shí)驗(yàn)結(jié)果基本一致。

以上研究者對(duì)塔板漏液的研究，都認(rèn)為塔板是水平的。由于塔體受到風(fēng)載荷等橫向載荷作用，傾斜塔板上產(chǎn)生非等高的液層，如圖1所示。

本文采用能改變塔板傾斜角度的實(shí)驗(yàn)塔裝置，利用水-空氣為實(shí)驗(yàn)操作介質(zhì)，測(cè)定了不同溢流強(qiáng)度、不同塔板傾角下，漏液隨空塔氣速不同的變化量。利用歐拉-歐拉雙流體方法中標(biāo)準(zhǔn)k-ε湍流模型，對(duì)傾斜塔板在漏液狀態(tài)下的狀況進(jìn)行了數(shù)值模擬，探索了傾斜塔板沿液相流動(dòng)方向不均勻漏液的分布情況。

1 ?實(shí)驗(yàn)塔裝置及測(cè)定方法

1.1 ?實(shí)驗(yàn)裝置

測(cè)量?jī)A斜塔板漏液的冷模實(shí)驗(yàn)塔裝置的結(jié)構(gòu)示意圖如圖2所示。

實(shí)驗(yàn)塔與底座利用螺栓連接，實(shí)驗(yàn)塔兩側(cè)各有一個(gè)調(diào)節(jié)螺釘，可以改變塔體傾斜角度。試驗(yàn)中，空氣首先由氣泵輸送，由轉(zhuǎn)子流量計(jì)測(cè)速后由塔底進(jìn)入塔內(nèi)，自下而上通過塔中的各層塔板;同時(shí)水箱中的水由水泵加壓輸送，經(jīng)轉(zhuǎn)子流量計(jì)測(cè)速后從塔頂注入塔內(nèi)，氣液兩相在塔板上接觸，空氣繼續(xù)上升從塔頂排空。實(shí)驗(yàn)裝置中用于測(cè)定漏液的傾斜塔板的結(jié)構(gòu)尺寸如表1所示。

1.2 ?實(shí)驗(yàn)步驟

實(shí)驗(yàn)前，將塔體左右兩側(cè)的調(diào)節(jié)螺釘調(diào)到同一高度，然后調(diào)高左側(cè)的螺釘，使塔體向右傾斜，使測(cè)定塔板的傾斜方向與液相流動(dòng)方向相反。記錄螺釘?shù)恼{(diào)整高度，塔底收集漏液的閥門處于打開狀態(tài)。試驗(yàn)塔正常操作開始后，保持液相流量和塔板傾角不變，降低氣泵流量，使塔板產(chǎn)生漏液，關(guān)閉塔底收集漏液的閥門，2 min后打開收集漏液閥門，計(jì)量漏液的體積。保持塔板漏液狀態(tài)及塔板傾角不變的情況下，改變液相流量，測(cè)量漏液。在氣相和液相流量都不變的情況下，調(diào)節(jié)螺釘?shù)母叨?，改變塔板的傾角，測(cè)量漏液。

2 ?CFD數(shù)值模擬

2.1 ?數(shù)學(xué)模型

當(dāng)氣相以一定速度通過塔板上一定厚度的液層時(shí)，氣液兩相是互相貫穿的連續(xù)介質(zhì)，故選用歐拉-歐拉方法，采用標(biāo)準(zhǔn)k-ε湍流模型封閉連續(xù)方程和動(dòng)量方程。

2.2 ?網(wǎng)格劃分及邊界條件

冷模試驗(yàn)塔的塔板直徑為380 mm，可以忽略液面落差的影響，故傾斜塔板數(shù)值模擬沒有設(shè)定液體入口和液體出口，模擬時(shí)塔板上方“PATCH”一定高度的水相。傾斜塔板的計(jì)算模型及邊界條件如圖3所示，X軸方向代表冷模試驗(yàn)塔液相流動(dòng)方向。

為了與實(shí)驗(yàn)對(duì)比，模擬塔板的結(jié)構(gòu)尺寸與表1相同，塔板上高度為150 mm，塔板下高度為50 mm。模型利用FLUENT軟件進(jìn)行數(shù)值模擬，采用非結(jié)構(gòu)六面體貼體網(wǎng)格。經(jīng)過網(wǎng)格無關(guān)性檢驗(yàn)，最后確定X軸方向、Y軸方向的網(wǎng)格大小為2 mm，Z軸方向的網(wǎng)格大小為5 mm。篩孔處的網(wǎng)格進(jìn)行了細(xì)化，網(wǎng)格尺寸最大為1.5 mm，整個(gè)計(jì)算域的網(wǎng)格總數(shù)大約為140萬左右。采用空氣和水為模擬工作介質(zhì)，氣體入口采用速度入口邊界條件，壓力出口作為氣體出口的邊界條件，出口靜壓為0，無滑移壁面邊界作為壁面條件。

2.3 ?數(shù)值模擬過程

SIMPLEC算法應(yīng)用到壓力-速度耦合方程中，動(dòng)量、湍動(dòng)能和湍動(dòng)耗散率方程的離散都采用二階迎風(fēng)格式，計(jì)算精度為10-3。模擬開始時(shí)，計(jì)算域的空氣體積分率為1，即全部為空氣，然后在塔板上方賦予一定高度的水，然后從氣體入口開始計(jì)算。監(jiān)測(cè)氣體注入的時(shí)間，提取塔板下部的液相體積分率，同時(shí)提取篩孔處液相的質(zhì)量流率。

3 ?結(jié)果和討論

3.1 ?傾斜塔板漏液的實(shí)驗(yàn)結(jié)果

在塔板溢流強(qiáng)度Lw=1.62×10-3m3·s-1·m-1時(shí)，塔板漏液量隨空塔氣速動(dòng)能因子變化曲線如圖4所示，圖4（a）表示塔板傾斜方向和板上液體流動(dòng)方向相同，圖4（b）表示方向相反。方向相同時(shí)，傾斜塔板漏液量小于水平塔板，傾角越大漏液量越少，傾角2?和4?塔板最大減少分別8.07%和15.07%。方向相反時(shí)，傾斜塔板漏液量大于水平塔板，傾角越大漏液量越多，傾角2?和4?塔板最大增加分別36.85%和46.81%。在FS=0.9 m·s-1（kg·m-3）0.5時(shí)，漏液量的減少變緩，這時(shí)塔板有氣相溝流現(xiàn)象。

為了保證板式塔的正常操作，規(guī)定塔板漏液量不能大于液相流量的10%，當(dāng)溢流強(qiáng)度Lw=1.62×10-3 m3·s-1·m-1時(shí)，冷模塔的操作下限對(duì)應(yīng)的漏液為4.16×10-2 kg·s-1，水平塔板的氣相操作下限為0.25 m·s-1。方向相同時(shí)，2°傾角塔板和4°傾角塔板氣相操作下限分別0.20和0.17 m·s-1;方向相反時(shí)，2°傾角塔板和4°傾角塔板氣相操作下限分別0.31 m·s-1和0.37 m·s-1。

FS=0.44 m·s-1（kg·m-3）0.5時(shí)，塔板漏液量隨塔板溢流強(qiáng)度的變化曲線如圖5所示，圖5（a）表示塔板傾斜方向和板上液體流動(dòng)方向相同，圖5（b）表示方向相反?？梢钥闯觯迓┮毫侩S著溢流強(qiáng)度的增加而增大。方向相同時(shí)，隨著溢流強(qiáng)度的增加，水平塔板和傾斜塔板的漏液之差逐漸減小，方向相反時(shí)，漏液之差逐漸加大。

3.2 ?CFD模擬結(jié)果

當(dāng)Fs=0.44 m?s-1（kg?m-3）0.5時(shí)，氣相注入時(shí)間t=0.5 s時(shí)，在液層非等高剖面內(nèi)的漏液情況如圖6所示。圖中可以觀察到，沿塔板X軸方向，水平塔板液層等高并產(chǎn)生均勻漏液，傾斜塔板液層非等高發(fā)生不均勻漏液。傾角4?的塔板在液層薄的區(qū)域有氣相溝流現(xiàn)象。

FS=0.44 m·s-1（kg·m-3）0.5時(shí)，水平塔板漏液的CFD模擬結(jié)果同實(shí)驗(yàn)和H&M[11]公式的對(duì)比如圖7所示?？梢钥闯?，CFD的模擬結(jié)果同實(shí)驗(yàn)結(jié)果以及公式計(jì)算結(jié)果是基本吻合的。CFD模擬和實(shí)驗(yàn)測(cè)量的結(jié)果比公式計(jì)算結(jié)果小，因?yàn)闅馑僭龃螅幸徊糠致┮罕粴庀鄮Щ厮濉?/p>

塔板向左傾斜，F(xiàn)S=0.22 m·s-1（kg·m-3）0.5，Lw=1.08×10-3 m3·s-1·m-1時(shí)，沿著塔板X軸方向篩孔的漏液分布如圖8所示?？梢钥闯?，水平塔板漏液沿著X軸方向變化不大，相鄰篩孔漏液平均相差為0.51%，可以認(rèn)為是均勻漏液。傾斜塔板沿著X軸方向，漏液量逐漸減小，相鄰篩孔漏液相差最小為2.56%，最大為21.56%，發(fā)生不均勻漏液。

4 ?結(jié)論

（1）塔板傾斜方向與板上液相流動(dòng)方向相同時(shí)，傾斜塔板漏液量小于水平塔板，傾角越大漏液量越少，塔板氣相操作下限減小。塔板傾斜方向與板上液相流動(dòng)方向相反時(shí)，傾斜塔板漏液量大于水平塔板，塔板氣相操作下限增大，傾角越大漏液量越多。

（2）傾斜塔板的漏液量隨塔板溢流強(qiáng)度增加而增大。隨著溢流強(qiáng)度的增加，出現(xiàn)液面落差，塔板傾斜方向與板上液相流動(dòng)方向相同時(shí)，水平塔板和傾斜塔板的漏液之差逐漸減小;塔板傾斜方向與板上液相流動(dòng)方向相反時(shí)，漏液之差逐漸加大。

（3）沿著液相流動(dòng)方向，水平塔板基本呈均勻漏液狀態(tài)，傾斜塔板發(fā)生不均勻漏液，尤其在液層最薄處，漏液量和其它篩孔相差很大。三維數(shù)值模擬可以定量分析每個(gè)篩孔的漏液分布情況，對(duì)研究?jī)A斜塔板的不均勻漏液有較好的實(shí)際意義。

參考文獻(xiàn)：

[1] 張平， 王翠華， 吳劍華. 板式塔工作狀態(tài)下塔頂撓度的計(jì)算[J]. 化學(xué)工程， 2011， 39（10）： 37-39.

[2]張平， 李雅俠， 秦然，等. 傾斜塔板液相三維流場(chǎng)的數(shù)值模擬[J]. 高?；瘜W(xué)工程學(xué)報(bào)， 2014， 28（6）： 1217-1222.

[3] 曾愛武， 劉福善， 許松林，等. 篩板不均勻漏液的規(guī)律及其影響[J]. 高校化學(xué)工程學(xué)報(bào)， 1996， 10（1）： 80-83.

[4]趙培， 施坤明， 熊丹柳，等. 十字旋閥塔板漏液性能[J]. 化工進(jìn)展， 2010， 29（2）： 223-227.

[5] 劉靜， 褚雅志， 王領(lǐng)， 等. 泡罩立體篩板的流體力學(xué)性能[J]. 化學(xué)工程， 2017， 45（7）： 18-21.

[6]孫玉春，涂鄭禹， 夏君. 篩孔塔板漏液操作狀態(tài)下氣液兩相流場(chǎng)的CFD模擬研究[J]. 化學(xué)工程師， 2011， 187（4）， 23-25.

[7]LOCKETT M J， BAMIK S. Weeping from sieve trays[J]. IEC Proc. Des. Dev.， 1986， 25（2）： 561-569.

[8] ALI Z， SEYYED H H， RAHHAR R， CFD study of weeping rate in the rectangular sieve trays[J]. Journal of the Taiwan Institute of Chemical Engineers， 2013， 44（1）： 27-33.

[9]ALI Z， SEYYED H H， RAHBAR R. CFD and experimental studies of liquid weeping in the circular sieve tray columns[J]. Chemical Engineering Research and Design， 2013， 91（12）： 2333-2345.

[10]VAN BATEN J M， KRISHNA R. Modelling sieve tray hydraulics using computational fluid dynamics[J].Chemical Engineering Science， 2000， 77： 143-151.

[11]KRISHNA R. CFD Simulations of Sieve Tray Hydrodynamics[J]. Trans.ICHEME， 1999， 77（10）： 639-646.