氣相F因子對(duì)規(guī)整填料內(nèi)液膜流動(dòng)的影響

李相鵬，楊立軍，李　鋒

(浙江工業(yè)大學(xué) 機(jī)械工程學(xué)院，浙江 杭州310014)

?

氣相F因子對(duì)規(guī)整填料內(nèi)液膜流動(dòng)的影響

李相鵬，楊立軍，李鋒

(浙江工業(yè)大學(xué) 機(jī)械工程學(xué)院，浙江 杭州310014)

摘要：采用CFD數(shù)值模擬方法研究規(guī)整填料內(nèi)液相參數(shù)改變時(shí)氣相F因子對(duì)氣液兩相流體動(dòng)力學(xué)的影響，建立基于流體力學(xué)軟件FLUENT并與Navier-Stokes方程和VOF方法相結(jié)合的二維模型，湍流結(jié)構(gòu)模擬采用RNG K-ε模型.對(duì)Mellapak 250Y型金屬板波紋填料在不同液相流量和表面張力下氣相F因子對(duì)液膜流體動(dòng)力學(xué)的影響進(jìn)行模擬分析，并與文獻(xiàn)實(shí)驗(yàn)結(jié)果對(duì)比驗(yàn)證.結(jié)果表明：不同液相參數(shù)下，受氣相F因子的影響，液膜逐漸由穩(wěn)定流動(dòng)變?yōu)橐旱?；持液量隨流量的增大而增加，隨表面張力的增加先增大后減??；當(dāng)流量超過臨界值時(shí)，有效傳質(zhì)面積突然增大，而有效傳質(zhì)面積隨表面張力的增加而持續(xù)增大；受氣相F因子和液相參數(shù)的影響，速度分布Vy存在臨界值，當(dāng)液相參數(shù)大于臨界值時(shí)Vy波動(dòng)較大，液面變化較大.

關(guān)鍵詞：氣相F因子；氣液兩相流；規(guī)整填料；液相流量

規(guī)整填料塔具有效率高、通量大、壓降低和放大效應(yīng)小等優(yōu)點(diǎn)，已在諸多工業(yè)分離過程中得到廣泛應(yīng)用[1].目前，對(duì)于規(guī)整填料塔的研究主要集中在傳質(zhì)過程和流體力學(xué)的實(shí)驗(yàn)方面[2]，但對(duì)于復(fù)雜的流場(chǎng)、溫度場(chǎng)和濃度場(chǎng)等的實(shí)驗(yàn)研究，并未形成一個(gè)普遍適用的經(jīng)驗(yàn)?zāi)Ｐ?，而僅僅適用于特定的實(shí)驗(yàn)范圍，因此所得出的預(yù)測(cè)結(jié)果具有一定的局限性[3].隨著計(jì)算機(jī)技術(shù)的廣泛應(yīng)用，計(jì)算流體力學(xué)(CFD)已經(jīng)成為研究流體動(dòng)力學(xué)的重要分析方法[4].Van Gulijk[5]利用CFD軟件對(duì)填料內(nèi)液相流動(dòng)及徑向返混行為進(jìn)行了研究，提出了“Toblerone”三維簡(jiǎn)化模型.在此基礎(chǔ)上Van Baten等[6-8]建立針對(duì)KATAPAK-S型反應(yīng)精餾規(guī)整填料的8孔道模型，模擬液相軸向和徑向的返混情況，得出的返混系數(shù)與實(shí)驗(yàn)值近似.Szulczewaka B等[9]建立了Mellapak 250Y型金屬填料的二維兩相逆流模型，對(duì)多種形式波紋板的氣液兩相逆流流動(dòng)的過程進(jìn)行了模擬，并將氣速、液速和液體的物理性質(zhì)對(duì)氣液相界面面積的影響進(jìn)行模擬研究.Raynal等[10]為研究規(guī)整填料液膜厚度、持液量等，建立了“Z”狀流動(dòng)路線的二維兩相流模型，并進(jìn)一步針對(duì)界面?zhèn)髻|(zhì)進(jìn)行了研究.浙江工業(yè)大學(xué)李相鵬[11]通過模擬研究液相物性對(duì)規(guī)整填料內(nèi)氣液兩相流的影響，重點(diǎn)研究了密度、粘度改變時(shí)對(duì)液膜流動(dòng)的影響；李相鵬等[12]模擬研究不同氣相F因子對(duì)液膜流體動(dòng)力學(xué)的影響，特別是各工況下液膜動(dòng)態(tài)運(yùn)動(dòng)分布、持液量以及有效傳質(zhì)面積的變化情況.

筆者結(jié)合前人的研究結(jié)果，建立填料的二維模型，并結(jié)合Navier-Stokes方程和VOF模型，研究氣液兩相流在填料內(nèi)的流動(dòng)情況，并通過改變液相參數(shù)，研究氣相F因子對(duì)液膜分布狀況、兩相流場(chǎng)的結(jié)構(gòu)變化、速度分布以及有效傳質(zhì)面積的影響情況.

1計(jì)算模型

1.1建模

Mellapak 250Y型金屬板波紋填料參數(shù)如表1所示.根據(jù)結(jié)構(gòu)尺寸，選其中的5個(gè)波紋建立二維氣液兩相流模型，其中模型總長(zhǎng)180 mm，氣相入口9 mm，氣相出口63 mm，液相入口5 mm，液相出口3 mm，波紋板間距為24 mm如圖1(a)所示.氣相由氣相進(jìn)口進(jìn)入，在壓力的作用下穿過曲折的流道，并最終由氣相出口流出，在此過程中，氣相與由液相進(jìn)口進(jìn)入的液相在曲折的流道中不斷地接觸，達(dá)到傳質(zhì)、混合的目的，液相最終由液相出口流出；圖1(b)為局部的網(wǎng)格、尺寸圖，網(wǎng)格為結(jié)構(gòu)化網(wǎng)格，利于跟蹤復(fù)雜的界面變化，同時(shí)在界面區(qū)域?qū)W(wǎng)格進(jìn)行加密.

表1　Mellapak 250Y型填料參數(shù)

圖1　模型結(jié)構(gòu)示意圖Fig.1　Model structure diagram

1.2邊界條件和計(jì)算方法

設(shè)置邊界條件，氣相入口和液相入口為velocity-in，氣相出口和液相出口為pressure-out，壁面無(wú)滑移，計(jì)算方法為VOF多相流模型和RNGK-ε氣液兩相湍流結(jié)構(gòu)相結(jié)合，這種方法能夠更為準(zhǔn)確的跟蹤氣液兩相界面的變化情況，以及液膜的形成和發(fā)展過程.

2模型的驗(yàn)證

為了驗(yàn)證CFD模型和計(jì)算方法的準(zhǔn)確性，以文獻(xiàn)[14]中的復(fù)雜表面液膜流動(dòng)實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，采用相同的試驗(yàn)工況建立模型，并將實(shí)驗(yàn)結(jié)果與模擬結(jié)果相互對(duì)比.根據(jù)文獻(xiàn)[14]中規(guī)整填料的結(jié)構(gòu)形式，選用其中的C型表面結(jié)構(gòu)作為實(shí)驗(yàn)對(duì)象，其表面結(jié)構(gòu)的圓弧半徑為1.587 5 mm，圓弧之間間距為0.2 mm，單節(jié)波紋總長(zhǎng)3.375 mm；以硅油為對(duì)比模擬介質(zhì)，其物理參數(shù)如下：密度969 kg/m3，粘度0.088 5 Pa·s，表面張力0.021 4 N/m，單位寬度液相流量為1.979×10-6m3/(s·m).如圖2所示，將流動(dòng)穩(wěn)定時(shí)的液膜狀態(tài)與CFD模擬結(jié)果對(duì)比發(fā)現(xiàn)：兩者契合度較好，誤差較小，絕對(duì)值偏差均小于10%，平均絕對(duì)偏差為3.06%，所建立的模型符合要求.

圖2　模擬結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比Fig.2　Comparison of simulation result and experimental result

3液相流量改變時(shí)氣相F因子對(duì)液膜流動(dòng)的影響

液膜流動(dòng)狀態(tài)主要是由氣相因子決定的，氣相作用在液膜表面使得液膜發(fā)生變化，但氣相F因子的作用會(huì)受到液相流量的影響，因此在此以空氣-甲醇為模擬對(duì)象，研究單位寬度液相流量改變時(shí)逆向氣載對(duì)氣液兩相流體動(dòng)力學(xué)特性的影響，單位寬度液相流量Q分別為1.5×10-4，3.0×10-4，4.5×10-4，6.0×10-4m3/(s·m)，氣相進(jìn)口速度取2 m/s(氣相F因子1.66 Pa0.5).

3.1對(duì)液膜流速度場(chǎng)分布的影響

由圖3可得規(guī)整填料內(nèi)的氣液流動(dòng)狀況，液膜速度沿波紋板表面呈周期性變化，靠近液膜的氣體在曳力的作用下一同下降，并與其他氣體產(chǎn)生剪切作用，從而在液膜附近形成回流.

圖3　氣液兩相流分布云圖Fig.3　Contours of gas-liquid

圖4為不同單位寬度液相流量時(shí)，氣相F因子作用下液膜自由界面沿填料表面分布圖.從圖4可得：?jiǎn)挝粚挾纫合嗔髁縌從1.5×10-4～4.5×10-4m3/(s·m)時(shí)，在氣相F因子作用下液膜表面分布基本一致，液膜貼合在填料表面，厚度在6.0×10-4m左右，只有微小波動(dòng)，且隨著時(shí)間的變化不大，只有在波紋板波峰和波谷處液流發(fā)生積聚，尤其在波紋板波谷處液流積聚比較明顯，4.5×10-4m3/(s·m)即為單位液相流量的峰值，當(dāng)小于峰值時(shí)，液膜流動(dòng)穩(wěn)定，而當(dāng)液體流量超過峰值時(shí)，液膜則會(huì)發(fā)生劇烈波動(dòng)，以Q=6.0×10-4m3/(s·m)為例，氣相F因子的作用變大，液膜波動(dòng)劇烈，破壞了液膜平衡，液膜厚度達(dá)到9.9×10-4m，且隨時(shí)間變化很大，部分液流被氣流吹起形成小液滴.

單位液相流量改變時(shí)，氣相F因子對(duì)速度分量vx和vy的影響基本一致，現(xiàn)只分析vy變化情況，如圖5所示，速度分量vy沿波紋板表面呈周期性變化，波峰處的速度vy>0，說明在波紋板波峰處液膜有向上的流動(dòng)；單位寬度液相流量從1.5×10-4m3/(s·m)增大到4.5×10-4m3/(s·m)時(shí)，隨著單位寬度液相流量的增大，氣相F因子作用下vy峰值有所變小，而vy谷值幾乎不變；當(dāng)Q=6.0×10-4m3/(s·m)時(shí)，即液相流量超過峰值時(shí)，氣相F因子作用下，液膜波動(dòng)劇烈，vy峰值急劇變小.

圖4　液膜分布隨液相流量Q的變化曲線Fig.4　Distribution of liquid surface with different Q

圖5　速度分布Vy隨液相流量Q的變化曲線Fig.5　Distribution of Vy with different Q

3.2對(duì)液膜傳質(zhì)參數(shù)的影響

圖6為單位寬度液相流量改變時(shí)，氣相F因子對(duì)持液量hL及有效界面?zhèn)髻|(zhì)面積與填料比表面積比值ae/ap的影響.單位寬度液相流量Q從1.5×10-4m3/(s·m)增大到4.5×10-4m3/(s·m)時(shí)，隨著液相流量的增加，在氣相F因子作用下，持液量hL隨之增大，但增幅較小；當(dāng)液相流量Q增大到6.0×10-4m3/(s·m)，流量超過臨界值時(shí)，hL增幅較大，而ae/ap的變化情況基本類似；在單位寬度液相流量Q從1.5×10-4m3/(s·m)增大到4.5×10-4m3/(s·m)時(shí)幾乎不變，當(dāng)單位寬度液相流量Q增大到6.0×10-4m3/(s·m)時(shí)，流量超過峰值，氣相F因子作用下，液膜波動(dòng)劇烈，有效傳質(zhì)面積增幅較大.

圖6　不同Q下，hL和ae/ap的變化曲線圖Fig.6　Curve of hL and ae/ap with different Q

4表面張力改變時(shí)氣相F因子對(duì)液膜流動(dòng)的影響

氣相F因子對(duì)液膜流動(dòng)狀態(tài)的影響同樣會(huì)受到液相流量的作用，在氣液兩相體系中，表面張力的作用主要是通過減小界面的面積而使自由能最小化，所以在液相物性中表面張力的作用至關(guān)重要.模擬時(shí)，表面張力σ選取0.02，0.04，0.06，0.08 N/m，將單位寬度液相流量Q設(shè)為3.0×10-4m3/(s·m).

4.1對(duì)液膜流速度場(chǎng)分布的影響

由圖7可得：表面張力σ從0.02 N/m增大至0.04 N/m時(shí)，在氣相F因子作用下液膜厚度有一定增大，但增幅不大，液膜平均厚度為6.0×10-4m左右，此時(shí)液膜達(dá)到臨界值；當(dāng)表面張力σ超過0.04 N/m后，臨界狀態(tài)打破，在氣相F因子的作用下液膜發(fā)生波動(dòng)，氣流吹動(dòng)液流形成小液滴，并且液膜厚度逐漸減?。槐砻鎻埩Ζ以?.06 N/m時(shí)液膜平均厚度為6.59×10-4m，表面張力σ在0.08 N/m時(shí)液膜平均厚度為5.58×10-4m.

圖7　液膜分布隨表面張力σ的變化情況Fig.7　Distribution of liquid surface with different σ

表面張力對(duì)液膜界面速度分量vx和vy的影響基本一致，圖8為界面速度分量vy分布曲線圖，表面張力σ從0.02 N/m增大至0.04 N/m時(shí)，速度分量vy沿波紋板表面呈周期性分布，在氣相F因子作用下vy峰值隨著表面張力增大有所減??；當(dāng)表面張力σ為0.06，0.08 N/m時(shí)，液膜發(fā)生波動(dòng)，但vy主體變化情況幾乎不變，且隨著表面張力增大vy峰值繼續(xù)減小.

圖8　液膜速度分布Vy隨表面張力σ的變化情況Fig.8　Distribution of Vy with different σ

4.2對(duì)液膜傳質(zhì)參數(shù)的影響

圖9為表面張力σ改變時(shí)逆向氣載對(duì)持液量hL及有效界面?zhèn)髻|(zhì)面積與填料比表面積比值ae/ap的影響.從圖9發(fā)現(xiàn)：表面張力σ從0.02 N/m增大至0.04 N/m時(shí)，在氣相F因子作用下持液量hL隨表面張力的增大而增加；表面張力σ超過0.04 N/m后，即表面張力超過臨界值時(shí)，持液量hL隨表面張力的增大而減少；而傳質(zhì)面積與填料比表面積比值ae/ap隨表面張力的增大而增加，但是增長(zhǎng)趨勢(shì)不明顯.

圖9　hL和ae/ap的變化曲線圖Fig.9　Curve of hL and ae/ap

5結(jié)論

通過對(duì)Mellapak 250Y型金屬板波紋填料建立二維模型，分析當(dāng)改變單位寬度液相流量Q和表面張力σ時(shí)，填料表面的液膜厚度、流動(dòng)情況以及氣相F因子對(duì)持液量hL和傳質(zhì)面積與填料比表面積比值ae/ap的影響情況.模擬結(jié)果表明：隨著單位寬度液相流量Q和表面張力σ的變化，填料表面均有一層穩(wěn)定的液膜，但在氣相F因子作用下，液膜流動(dòng)會(huì)發(fā)生變化，并在波峰、波谷處發(fā)生聚集；隨單位寬度液相流量Q的增加，流道變窄，氣相F因子作用增大，液膜波動(dòng)劇烈，液膜厚度隨之增大；速度分量vy波動(dòng)幅度增大，并有液滴形成；hL和ae/ap也隨之增大；而隨著表面張力σ的增加，穩(wěn)定液膜流產(chǎn)生波動(dòng)，形成液滴；速度分量vy主體不變，峰值減??；氣相F因子與表面張力σ共同作用下會(huì)使ae/ap值隨之增大，但hL值卻先增大后減小.同時(shí)，模型采用FLUENT軟件，可以更為直觀的得出液膜在不同液相參數(shù)時(shí)隨時(shí)間的變化情況，這為得出臨界條件、設(shè)置合理的操作參數(shù)提供了便利，同時(shí)FLUENT還可以隨意改變各項(xiàng)參數(shù)，為得出不同狀態(tài)下的試驗(yàn)參數(shù)提供了方便，而且節(jié)約了大量的時(shí)間和資金，也為裝置結(jié)構(gòu)的改進(jìn)提供參考.

參考文獻(xiàn)：

[1]裘俊紅，陳國(guó)標(biāo)，計(jì)建炳.波紋板規(guī)整填料塔液體分布[J].化工學(xué)報(bào)，2003，54(5):646-652.

[2]胡暉，徐世民，李盡鋼.大型填料塔技術(shù)及其工業(yè)應(yīng)用[J].現(xiàn)代化工，2005，25(7):53-55.

[3]陳江波，劉春江.計(jì)算流體力學(xué)方法模擬規(guī)整填料塔內(nèi)流體流動(dòng)行為的研究進(jìn)展[J].石油化工，2007，36(9):961-965.

[4]張兆鑫，趙元虎，徐奔馳.基于CFD法的小型風(fēng)機(jī)非扭曲葉片氣動(dòng)性能分析[J].浙江工業(yè)大學(xué)學(xué)報(bào)，2013，41(1):68-72.

[5]VAN C G. Using computational fluid dynamics to calculate transversal dispersion in a structured packed bed Computers[J]. Computers & chemical engineering,1998,22:767-770.

[6]VAN BATEN J M, KRISHNA R. Liquid-phase mass transfer within KATAPAK-S structures studied using computational fluid dynamics simulations[J]. Catalysis today,2001,69:371-377.

[7]VAN BATEN J M, ELLENBERGER J, KRISHNA R. Radial and axial dispersion of the liquid phase within a KATAPAK-S structure: experiments vs CFD simulations[J]. Chemical engineering science,2001,56:813-821.

[8]VAN BATEN J M, KRISHNA R. Gas and liquid phase mass transfer within KATAPAK-S structures studied using CFD simulations[J]. Chemical engineering science,2002,57:1531-1536.

[9]SZULCZEWAKA B, ZBICINSKI I, GORAK A. Liquid flow on structured packing: CFD simulation and experimental study[J]. Chemical engineering technology,2003,26(5):580-584.

[10]RAYNAL L, BOYER C, BALLAGUET J P. Liquid holdup and pressure drop determination in structured packing with CFD simulations[J]. The canadian journal of chemical engineering,2004,82:871-879.

[11]李相鵬，汪蓉梅.液相物性對(duì)規(guī)整填料內(nèi)氣液兩相流的影響[J].化學(xué)工程，2014，42(1):40-44

[12]李相鵬，李鋒，張超.規(guī)整填料內(nèi)氣相載荷對(duì)液膜流動(dòng)特征的影響[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2015，46(4):1239-1246.

[13]王樹楹.現(xiàn)代填料塔技術(shù)指南[M].北京:中國(guó)石化出版社，1998:9.

[14]ZHAO L Q, CERRO R L. Experimental characterization of viscous film flows over complex surfaces[J]. International journal of multiphase flow,1992,18(4):495-516.

(責(zé)任編輯：陳石平)

Influences of countercurrent gasFfactor on dynamics of liquid film flow over structured packing

LI Xiangpeng, YANG Lijun, LI Feng

(College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310014, China)

Abstract:The influences of F factor in different liquid parameters on the hydrodynamics simulation of gas-liquid two-phase flow in the structured packing were simulated using the CFD numerical techniques. A 2D model based on Navier-Stokes equations and VOF method was developed based on the platform of FLUENT, and the turbulent structure was modeled with RNG K-ε model. The simulations were conducted to study the variations of liquid film hydrodynamic with F factor at different liquid flow rate, surface tension over Mellapak 250Y corrugated packing sheet, and the model was validated by the experimental results presented in literature. The results show that at the different liquid parameters, because of the impact of F factor, the liquid film from a steady flow into droplets; the liquid holdup increases with the Q, decreases with the increase of the σ; When the flow exceeds a critical value, a sudden increase in the ae/ap, and the ae/ap increase with the increase of the σ; as the influence of F factor and liquid parameters, the Vy has a critical value, when the phase parameter is greater than the critical value, the Vy changes in the larger level.

Keywords:F factor; gas-liquid two-phase flow; structured packing; liquid flow rate

收稿日期：2015-09-14

基金項(xiàng)目：浙江省自然科學(xué)基金資助項(xiàng)目(Y20110169)；“十二五”科技支撐計(jì)劃項(xiàng)目(2011BAK06B02)

作者簡(jiǎn)介：李相鵬(1977—)，男，湖南寧鄉(xiāng)人，副教授，主要從事熱能工程、化工過程等數(shù)值模擬與優(yōu)化研究，E-mail：rocleexp@zjut.edu.cn.

中圖分類號(hào)：TQ028.4

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1006-4303(2016)02-0195-07