穿流塔板上兩相競(jìng)爭(zhēng)性流動(dòng)特性

張家佳， 王亦飛， 毛霞君， 章岳峰， 于廣鎖（華東理工大學(xué)煤氣化與能源化工教育部重點(diǎn)實(shí)驗(yàn)室，上海市煤氣化工程技術(shù)研究中心，上海 200237）

穿流塔板上兩相競(jìng)爭(zhēng)性流動(dòng)特性

張家佳， 王亦飛， 毛霞君， 章岳峰， 于廣鎖
（華東理工大學(xué)煤氣化與能源化工教育部重點(diǎn)實(shí)驗(yàn)室，上海市煤氣化工程技術(shù)研究中心，上海 200237）

對(duì)一種穿流式固閥塔板進(jìn)行了流體力學(xué)冷態(tài)實(shí)驗(yàn)考察。采用常溫、常壓下的空氣和水作為實(shí)驗(yàn)介質(zhì)，測(cè)試了不同操作條件下塔板的干板壓降、全板壓降和清液層高度，比較了不同的塔板結(jié)構(gòu)參數(shù)對(duì)流體力學(xué)性能參數(shù)的影響。實(shí)驗(yàn)結(jié)果表明，塔板壓降和清液層高度均隨著孔中心距和開孔率的減小而增大，而隨著孔徑的增大而增大。氣液兩相在穿流塔板上為逆流流動(dòng)，且隨著氣液流量的變化呈現(xiàn)出4種不同的流動(dòng)形態(tài)。建立了清液層高度的經(jīng)驗(yàn)式，且計(jì)算值與實(shí)驗(yàn)測(cè)量值吻合良好。推導(dǎo)了描述穿流塔板獨(dú)特的流動(dòng)形態(tài)的參數(shù)即氣體流通分率的數(shù)學(xué)表達(dá)式，計(jì)算結(jié)果表明，該分率隨著孔中心距的減小而減小，并隨著孔徑和開孔率的增大而減小。

穿流塔板；流體力學(xué)；多相流；固閥塔板

在分離設(shè)備中，板式塔由于造價(jià)低、應(yīng)用范圍廣而具有較大的優(yōu)勢(shì)［1］。在溢流式塔板上，出口安定區(qū)兩相的相間接觸減少，對(duì)于氣相夾帶飛灰的洗滌過程，固體易于在此沉積。近年來學(xué)者們將關(guān)注的重點(diǎn)轉(zhuǎn)移到穿流塔板上。穿流塔板是一種沒有降液管的多孔塔板，氣液兩相交替流過塔板開孔區(qū)的獨(dú)特流動(dòng)形式為穿流塔板提供了“自清潔”特性［2］，在含有固相的結(jié)晶和洗滌等過程中，穿流塔板表現(xiàn)出良好的抗堵性能。

在前期的研究中，人們對(duì)穿流塔板的研究主要集中于傳質(zhì)效率，鮮有對(duì)塔板上的流體力學(xué)特性研究。Xu等［3］在穿流塔板上對(duì)甲醇-水和甲醇-異丙醇精餾系統(tǒng)的塔板效率進(jìn)行了實(shí)驗(yàn)研究。Hutchinson等［4］研究了不同負(fù)荷和回流比對(duì)甲醇-水體系的分離效率的影響。Furzer等［5］研究了不同塔板數(shù)的穿流篩孔板的傳質(zhì)效率。Miyahara等［6］研究了穿流塔板上泡沫區(qū)和轉(zhuǎn)變區(qū)的液相傳質(zhì)系數(shù)和氣液相界面積。Garcia等［7］提出了穿流塔板用于精餾過程的設(shè)計(jì)和分析的理論方法。Domingues等［8］提出了新方法嵌套Aspen Plus 12.1模擬軟件，對(duì)溢流浮閥塔板和穿流塔板的全塔效率進(jìn)行了預(yù)測(cè)。Kister等［9］將穿流塔板與溢流式篩孔、浮閥和泡罩塔板進(jìn)行了傳質(zhì)性能的比較。雖然也有一些研究［10-13］對(duì)穿流塔板和其他一系列塔板進(jìn)行比較，但是所研究的穿流塔板的鼓泡元件為篩孔，對(duì)固閥的研究較為少見。

本文測(cè)量了穿流固閥塔板的塔板壓降、清液層高度等塔板流體力學(xué)參數(shù)，并對(duì)清液層高度提出關(guān)聯(lián)式。描述了穿流塔板上兩相逆流接觸的4種不同流動(dòng)形態(tài)。定義了描述穿流塔板獨(dú)特流動(dòng)形態(tài)的特征參數(shù)即塔板開孔區(qū)氣體流通分率，并推導(dǎo)了該氣體流通分率的計(jì)算模型。研究了氣液兩相流量、塔板孔徑、孔中心距和開孔率等結(jié)構(gòu)參數(shù)對(duì)塔板的流體力學(xué)性能的影響。

1　研究方法

本文在內(nèi)徑500 mm的有機(jī)玻璃塔內(nèi)對(duì)穿流固閥塔板進(jìn)行了流體力學(xué)特性研究。塔板上兩相流動(dòng)形式如圖1所示，實(shí)驗(yàn)裝置圖如圖2所示，塔板結(jié)構(gòu)參數(shù)見表1，其中dh為孔徑，th為孔中心距，n為孔數(shù)，AO為開孔率。采用常溫常壓的空氣和水為測(cè)試流體?？諝庥晒娘L(fēng)機(jī)經(jīng)氣體均布器進(jìn)入塔底，向上經(jīng)過塔板后從塔頂氣體出口排出。水由儲(chǔ)罐經(jīng)過液體均布器進(jìn)入塔頂，向下經(jīng)過塔板后從塔底返回儲(chǔ)罐以供循環(huán)。

調(diào)節(jié)空氣和水的流量以保證塔板處于正常的操作狀態(tài)，塔板壓降由壓差傳感器采集，清液層高度由U型液位計(jì)測(cè)量。塔板上圓形固閥呈正三角形排列，實(shí)驗(yàn)中分別對(duì)4塊塔板進(jìn)行測(cè)試，以考察孔徑、孔中心距離和開孔率等結(jié)構(gòu)參數(shù)對(duì)塔板流體力學(xué)性能的影響。

圖1　穿流固閥塔板的流動(dòng)形態(tài)Fig.1 Flow pattern of dual-flow fixed valve tray

圖2　實(shí)驗(yàn)流程圖Fig.2 Experiment flow chart

表1　穿流塔板結(jié)構(gòu)參數(shù)Table 1 Dual-flow trays specifications

2　結(jié)果與討論

2.1兩相流型

氣液兩相在穿流塔板上逆流流動(dòng)，液體穿過閥孔流入下一級(jí)塔板，與下一塔板的泡沫區(qū)和噴射區(qū)接觸，為兩相提供了更大的相接觸面積，其有效傳遞空間更大。穿流塔板上的兩相競(jìng)爭(zhēng)性地流過同一開孔區(qū)域的特征，表現(xiàn)為塔板上兩相混合的動(dòng)態(tài)變化過程，即在某一閥孔中的不同時(shí)間內(nèi)氣液兩相交替流過，在某一瞬間，一部分閥孔為液體占據(jù)，另一部分為氣體流動(dòng)通道。

在本實(shí)驗(yàn)中，對(duì)于不同的氣液流量，塔板上有4種不同的流動(dòng)狀態(tài)，如圖3和圖4所示。在低氣速和清液層高度較小時(shí)，氣體通過一部分閥孔無規(guī)則地鼓泡，而且鼓泡是間斷的，即塔板上液體間斷地阻塞氣體流動(dòng)。在低氣速和清液層高度略大時(shí)，出現(xiàn)蜂窩狀泡沫，這是塔板上充滿一層很薄的、湍動(dòng)程度很低的蜂窩狀泡沫層，氣液兩相無規(guī)則地通過塔板。在氣速和清液層高度適中時(shí)，塔板上出現(xiàn)均勻泡沫層，此時(shí)泡沫層開始湍動(dòng)，在其中形成無規(guī)則的波并使液體從某些孔中滴漏下去。在高氣速和清液層高度較大情況下形成循環(huán)泡沫區(qū)，這時(shí)泡沫激烈湍動(dòng)，在泡沫層中形成循環(huán)（中心處向上，壁面附近向下），液體主要在塔壁附近滴漏。

圖3　穿流塔板上氣液接觸狀態(tài)Fig.3 Gas-liquid flow pattern on the dual-flow tray

圖4　泡沫層周期振動(dòng)的類型Fig.4 Type of periodic vibration on the froth layer

在清液層高度較大且氣速較高的情況下，充氣液體在整個(gè)塔板上呈現(xiàn)擺動(dòng)狀態(tài)，此時(shí)可以分為兩種振動(dòng)方式：第1種方式在清液層中等高度時(shí)出現(xiàn)，它有一個(gè)波節(jié)直徑，這時(shí)液體在塔板上從一側(cè)向另一側(cè)來回?cái)[動(dòng)，同時(shí)還伴隨著圍繞塔板中心的旋轉(zhuǎn)運(yùn)動(dòng)，液體在波峰處周期性地滴漏；第2種方式在清液層高度較小時(shí)出現(xiàn)，它有一個(gè)波節(jié)圓，這種波對(duì)塔體中心線是對(duì)稱的，而且有連續(xù)的波峰和波谷。這兩種振動(dòng)方式的振幅均隨氣速和持液量的增加而增加。

2.2壓降

塔板壓降是影響氣體分布性能的主要因素，在低壓條件下尤為顯著［14］。本實(shí)驗(yàn)中塔板的干板壓降（pd）如圖5所示。結(jié)果表明，干板壓降隨著孔中心距、孔徑和開孔率的減少而增大，且隨著氣體表觀流速（G）的增大而增大。

圖5　塔板的干板壓降Fig.5 Dry-plate pressure drop of trays

穿流塔板上的壓力平衡方程可寫為

其中：ht為塔板的總壓降，以液層高度表示；hdG為氣體通過x分率的閥孔的壓降；hL為氣體穿過泡沫層的壓降；hLd為液體通過1－x分率閥孔的壓降；h′L為推動(dòng)液體穿過1－x分率閥孔的液體靜壓力。

hdG一般形式為［15］：

其中：ρL為水密度；ρG為蒸汽密度；g為重力加速度；u2x為穿過x分率的閥孔的氣體線性速度；Cv為孔流系數(shù)。

孔流系數(shù)Cv的計(jì)算模型存在很多的不同形式。對(duì)于穿流塔板上的穿孔流動(dòng)，Leibson［16］最早以關(guān)聯(lián)圖的形式提出了孔流系數(shù)的函數(shù)表達(dá)，歸納為關(guān)聯(lián)式如下：

其中：tt為塔板的板厚；dh為閥孔孔徑；Ah為開孔面積。

穿流塔板上的流體動(dòng)力學(xué)特征表現(xiàn)為液體的動(dòng)量轉(zhuǎn)化為克服靜壓梯度，以使液體向下流動(dòng)的靜壓頭。塔板上兩相的振蕩是由液體重力與氣體壓降和表面張力總和的動(dòng)態(tài)競(jìng)爭(zhēng)所導(dǎo)致。相比篩孔塔板，固閥塔板一般采用較大的孔徑，故閥孔處的表面張力可忽略不計(jì)，固閥上的力的競(jìng)爭(zhēng)為液體重力和氣體慣性。對(duì)于此復(fù)雜過程，本文關(guān)注氣體流動(dòng)部分：

圖6～圖8給出了在相同的實(shí)驗(yàn)條件下，不同塔板的壓降隨著氣體和液體流量（L）的變化規(guī)律。由圖可知，塔板壓降隨著液體流量的增加而增加，其隨氣體流量的變化亦有相同的趨勢(shì)。由塔板的比較可得，塔板壓降隨著孔中心距和開孔率的減小而增大，隨著孔徑的增大而增大。特別地，如圖7所示，孔徑對(duì)于塔板壓降的影響并不顯著。

圖6　不同孔中心距的塔板（Tray 1#，2#）下，氣液兩相流量對(duì)塔板壓降的影響Fig.6 Effect of gas and liquid flow on total pressure drop of tray 1#and tray 2#

圖7　不同孔徑的塔板（Tray 2#，3#）下，氣液兩相流量對(duì)塔板壓降的影響Fig.7 Effect of gas and liquid flow on total pressure drop of tray 2#and tray 3#

2.3清液層高度

Kotschering等［17］提出了穿流塔板上的當(dāng)量清液層高度的關(guān)聯(lián)式，Xu等［3］對(duì)該關(guān)聯(lián)式做了修正，該修正方程如下所示：

然而，該關(guān)聯(lián)式并沒有包含孔中心距的影響，根據(jù)實(shí)驗(yàn)結(jié)果，對(duì)該關(guān)聯(lián)式修正如下：

其中：L為液體流量；Us為表觀氣速；ML為液相摩爾質(zhì)量；指數(shù)n由式（7）計(jì)算。

圖8　不同開孔率塔板（Tray 1#，4#）下，氣液兩相流量對(duì)塔板壓降的影響Fig.8 Effect of gas and liquid flow on total pressure drop of tray 1#and tray 4#

由式（6）所得的計(jì)算值與實(shí)驗(yàn)測(cè)量值的對(duì)照如圖9所示，該關(guān)聯(lián)式的標(biāo)準(zhǔn)偏差為14.3%。

圖9　清液層高度的實(shí)驗(yàn)值與計(jì)算值的對(duì)照Fig.9 Parity plot for the experimental and calculated clear liquid height

清液層高度隨著氣體和液體流量的變化趨勢(shì)如圖10～圖12所示。清液層高度隨著氣體流量的增大而增大，并對(duì)液體流量的增大呈現(xiàn)出近似線性增大的規(guī)律。由不同塔板間的比較可得，清液層高度隨著孔中心距和開孔率的減小而增大，而隨著孔徑的增大而增大。在低負(fù)荷條件下，泡沫層高度顯著減少，只有較少的孔內(nèi)存在兩相的接觸，氣液兩相接觸時(shí)間縮短，此時(shí)不同塔板的清液層高度幾乎沒有區(qū)別。

圖10　不同孔中心距塔板（Tray 1#，2#）下，氣液兩相流量對(duì)清液層高度的影響Fig.10 Effect of gas and liquid flow on the clear liquid height of tray 1#and tray 2#

圖11　不同孔中心距塔板（Tray 2#，3#）下，氣液兩相流量對(duì)清液層高度的影響Fig.11 Effect of gas and liquid flow on the clear liquid height of tray 2#and tray 3#

圖12　不同孔中心距塔板（Tray 1#，4#）下，氣液兩相流量對(duì)清液層高度的影響Fig.12 Effect of gas and liquid flow on the clear liquid height of tray 1#and tray 4#

2.4開孔區(qū)的氣體流通分率

穿流塔板上的“漏液”是產(chǎn)生于塔板上的湍流波和液體表面的振蕩，液體以連續(xù)水柱而非簡(jiǎn)單的水滴形式穿過閥孔。在某一時(shí)刻，氣液兩相可能通過同一孔流動(dòng)或者分別穿過不同的孔流動(dòng)。

氣體和液體流動(dòng)的自由通道面積總和等于塔板的開孔面積，即

其中：AFG為氣體流通面積；AFL為液體流通面積。

自由通道面積對(duì)于每一塔板是固定值，其值等于一段時(shí)間（ts）內(nèi)氣體和液體流通面積的平均值，其中ts為一段充足的時(shí)間以保證以下積分收斂：

氣液兩相穿過自由面積主要有兩相同時(shí)通過某一孔或者分別在不同的孔內(nèi)流動(dòng)，所有孔內(nèi)振蕩性的交替流過氣體和液體，或兩者的組合形式。

由此，塔板上開孔區(qū)的氣體流通分率定義如下：

由于沒有直接的測(cè)量方法用于該參數(shù)的測(cè)量，故其必須由塔板壓降等其他參數(shù)推導(dǎo)得出。由實(shí)驗(yàn)測(cè)量得到塔板壓降和清液層高度后，式（4）可寫為：

可得出

氣體流通分率受各參數(shù)影響的結(jié)果如圖13～15所示。由圖可知，該氣體流通分率隨著氣體流量的增大而增大，隨著液體流量的增大而減小。由不同塔板間的比較可知，塔板上開孔區(qū)的氣體流通分率隨著孔中心距的增大而增大，而隨著孔徑和開孔率的減小而增大。但孔徑對(duì)于該分率的影響并不顯著（圖14）。如圖15所示，當(dāng)塔板的開孔率已足夠大時(shí)，氣體流通分率將不再隨液體流量的變化而變化。

3　結(jié) 論

本文通過對(duì)孔徑、孔中心距和開孔率等結(jié)構(gòu)參數(shù)對(duì)穿流式固閥塔板的流體力學(xué)性能的影響進(jìn)行研究，得出塔板壓降和清液層高度對(duì)于塔板結(jié)構(gòu)參數(shù)的變化規(guī)律，即隨著孔中心距和開孔率的減小而增大，而隨著孔徑的增大而增大。對(duì)清液層高度與塔板結(jié)構(gòu)參數(shù)和操作變量進(jìn)行了關(guān)聯(lián)，擬合值與實(shí)驗(yàn)測(cè)量值吻合良好。對(duì)穿流塔板的特征參數(shù)即氣體流通分率進(jìn)行了定義并推導(dǎo)計(jì)算式，計(jì)算結(jié)果表明，塔板上氣體流通分率隨著孔中心距的減小而減小，并隨著孔徑和開孔率的增大而減小。且在各個(gè)塔板結(jié)構(gòu)變量中，孔徑對(duì)氣體流通分率的影響最不顯著。

圖13　孔中心距對(duì)氣體流通分率的影響Fig.13 Effect of center distance on the fraction of holes passing gas

圖14　孔徑對(duì)氣體流通分率的影響Fig.14 Effect of hole’s diameter on the fraction of holes passing gas

圖15　開孔率對(duì)氣體流通分率的影響Fig.15 Effect of open area on the fraction of holes passing gas

［1］ LUO Na，QIAN Feng，YE Zhencheng，et al.Estimation of

mass-transfer efficiency for industrial distillation columns［J］. Industrial&Engineering Chemistry Research，2012，51（7）：3023-3031.

［2］ GAICIA J A，F(xiàn)AIR J R.A fundamental model for the prediction of distillation sieve tray efficiency.1.Database development［J］.Industrial&Engineering Chemistry Research，2000，39（6）：1809-1817.

［3］ XU Z P，AFACAN A，CHUANG K T.Efficiency of dualflow trays in distillation［J］.The Canadian Journal of Chemical Engineering，1994，72（4）：607-613.

［4］ HUTCHINSON M，BADDOUR R.Ripple trays：A new tool for vapor-liquid contacting［J］.Chemical Engineering Progress，1956，52：503-508.

［5］ FURZER I A，DUFFY G J.Mass transfer on sieve plates without downcomers［J］.Chemical Engineering Journal，1977，14（3）：217-224.

［6］ MIYAHARA T，KURIHARA M，ASODA M，et al.Gasliquid interfacial area and liquid-phase mass transfer coefficient in sieve plate columns without doencomer operating at high gas velocities［J］.Journal of Chemical Engineering of Japan，1990，23（3）：280-285.

［7］ GAICIA J A，F(xiàn)AIR J R.Distillation sieve trays without downcomers：Prediction of performance characteristics［J］. Industrial&Engineering Chemistry Research，2002，41（6）：1632-1640.

［8］ DOMINGUES T L，SECCHI A R，MENDES T F.Overall efficiency evaluation of commercial distillation columns with valve and dualflow trays［J］.AIChE Journal，2010，56（9）：2323-2330.

［9］ KISTER H Z.Distillation Design［M］.New York：McGraw-Hill，1992.

［10］ BILLET R.Separation tray without downcomers［J］. Chemical Engineering&Technology，2001，24（11）：1103-1112.

［11］ SHOUKRY E，CERMAK J，KOLAR V.On the hydrodynamics of sieve plates without downcomers：I. Dynamics nature of the gas liquid dispersions［J］.Chemical Engineering Journal，1974，8（1）：27-40.

［12］ SHOUKRY E，KOLAR V.On the hydrodynamics of sieve plates without downcomers：Ⅱ.A model of plate operation and the phenomenon of multiplicity of steady state［J］. Chemical Engineering Journal，1974，8（1）：41-51.

［13］ FURZER I A.Froth heights on dualflow trays with a ternary azeotropic system of ethyl acetate-ethanol-water［J］. Industrial&Engineering Chemistry Research，2000，39（5）：1430-1436.

［14］ WIJIN E F.On the lower operating range of sieve and valve trays［J］.Chemical Engineering Journal，1998，70（2）：143-155.

［15］ 賴文衡.關(guān)于新型垂直篩板塔板壓降問題［J］.化學(xué)工程，1983，11（5）：38-41.

［16］ LEIBSON I，KELLEY R E，BULLINGTON L A.How to design perforated trays［J］.Petrol Refine，1957，36：127-128.

［17］ KOTSCHERING N A，OLESKI W M，DILMAN W W. Investigation of dualflow tray behavior under distillation conditions［J］.Khimicheskaya Promyshlennost，1960，7：591-593.

Competitive Characteristics of Two Phase Flow on the Dual-Flow Tray

ZH ANG Jia-jia， WANG Yi-fei， MAO Xia-jun， ZHANG Yue-feng， YU Guang-suo
（Key Laboratory of Coal Gasification and Energy Chemical Engineering of Ministry of Education，
Shanghai Engineering Research Center of Coal Gasification，East China University of Science and Technology，Shanghai 200237，China）

This paper addresses an experimental investigation in the hydrodynamic behavior of a new type of dual-flow fixed valve tray.Water/air system was used at atmospheric pressure and ambient temperature.The dry pressure drop，total pressure drop and clear liquid height were measured and the comparisons between different tray geometries were made.The results show that the total pressure drop and clear liquid height have the same trend，which is increased with the decreased center distance of holes，increased hole’s diameter and decreased open area.Correlation for clear liquid height is proposed，and the agreement of the experimental and calculated data is demonstrated.The results show that gas and liquid flow counter-currently through the tray holes and four main hydrodynamic regimes are distinguished at different gas/liquid loads.A characteristic parameter named the fraction of holes passing gas is defined and deduced.The calculated results show that the fraction of holes passing gas is reduced by the decreased center distance of holes，while the increased hole’s diameter and open area leads to the less fraction.

dual-flow tray；hydrodynamics；multiphase flow；fixed valve tray

TQ546

A

1006-3080（2016）01-0009-06 DOI：10.14135/j.cnki.1006-3080.2016.01.002

2015-06-28

國(guó)家高技術(shù)研究發(fā)展計(jì)劃（863計(jì)劃）（2012AA053101）

張家佳（1986-），女，安徽蚌埠人，主要研究氣液兩相流傳熱強(qiáng)化。

王亦飛，E-mail：wangyf@ecust.edu.cn