劉萌萌

(天津美加恒升科技發(fā)展有限公司，天津 300380)

塔器作為分離設(shè)備，是化工、煉油等生產(chǎn)中最重要的設(shè)備之一，分為板式塔和填料塔。板式塔有多種溢流形式。在實(shí)際大規(guī)?；ぱb置中，多溢流塔板因可以承受較大的液相負(fù)荷而得到廣泛應(yīng)用，其中雙溢流塔板應(yīng)用最為廣泛。雙溢流塔板完全對(duì)稱，氣液分配均勻，但是處理量相對(duì)較小，因此很多大型裝置中采用三溢流或四溢流塔板，對(duì)于液量特別大的塔器，需采用六溢流結(jié)構(gòu)。六溢流塔板可以縮短塔板流道長(zhǎng)度，有效減小塔板上的液面落差，從而使塔板操作更加穩(wěn)定。但因六溢流塔板結(jié)構(gòu)復(fù)雜，較難實(shí)現(xiàn)氣液兩相在塔板上的均勻分配，如果設(shè)計(jì)不合理，將導(dǎo)致塔板傳質(zhì)效率降低。張寶樹等[1]對(duì)比了等鼓泡面積法和等通道長(zhǎng)度法在六溢流塔板設(shè)計(jì)中的氣液分配情況，認(rèn)為等通道長(zhǎng)度法與等鼓泡面積法相比更易實(shí)現(xiàn)氣液分配一致，且加工和安裝比較容易；孫蘭義等[2,3]通過(guò)對(duì)比三溢流和四溢流塔板的不同設(shè)計(jì)方法的氣液分配，也認(rèn)為等通道長(zhǎng)度法優(yōu)于等鼓泡面積法。

塔板上流體的流動(dòng)狀況對(duì)塔板效率有著直接的影響，尤其是塔板上液體的返流、滯留都會(huì)顯著降低塔板的分離效率，因此掌握塔板上流體的分布情況對(duì)塔板設(shè)計(jì)有著重要的指導(dǎo)意義。但塔板上氣液流動(dòng)狀態(tài)復(fù)雜，目前的實(shí)驗(yàn)監(jiān)測(cè)手段很難開展較深入的研究，利用計(jì)算流體力學(xué)(computational fluid dynamics, CFD)模擬成為一種行之有效且省時(shí)省力的研究方法。近年來(lái)，針對(duì)篩板、浮閥、垂直篩板等結(jié)構(gòu)的流體力學(xué)研究非常多，但大都專注于傳質(zhì)元件本身性能，多溢流的模擬并不多見，趙丹[4]采用二維模型模擬了四溢流塔板在不同工況下氣液流動(dòng)及對(duì)應(yīng)的傳質(zhì)效率，但并未給出三維模擬結(jié)果。六溢流塔板塔徑較大且結(jié)構(gòu)復(fù)雜，對(duì)計(jì)算機(jī)硬件要求較高，計(jì)算時(shí)間過(guò)長(zhǎng)，故鮮有用CFD方法模擬六溢流塔板氣液流場(chǎng)的報(bào)道。作者以丙烷脫氫裝置中產(chǎn)品分離塔為研究對(duì)象，采用商業(yè)軟件ICEM建立了六溢流塔板的三維模型，并采用CFX 14.5模擬了塔內(nèi)的氣液相分布情況，對(duì)塔板的設(shè)計(jì)具有指導(dǎo)意義。

1 模型的建立

1.1 數(shù)學(xué)模型

常采用的兩相流模型包括歐拉雙流體模型、顆粒軌道模型等，不同模型有著不同的適用體系。當(dāng)分散相體積分率相對(duì)較高且連續(xù)地分布于主相中時(shí)，歐拉雙流體模型比較適用。精餾塔塔板上氣相雖為分散相，但其體積分率較高，塔板上方鼓泡區(qū)的混合物氣液相貫穿，因此作者采用能更好描述塔板上流動(dòng)狀態(tài)的歐拉雙流體模型進(jìn)行模擬，模型基本方程如下：

氣相連續(xù)性方程：

(1)

液相連續(xù)性方程：

(2)

氣相動(dòng)量守恒方程：

?·(αGρGvGvG)=-αG?pG+?·(αGμeff,G(?vG+(?vG)T))-MLG

(3)

液相動(dòng)量守恒方程：

?·(αLρLvLvL)=-αL?pL+?·(αLμeff,L(?vL+(?vL)T))-MLG

(4)

相體積分率方程：

αG+αL=1

(5)

壓力限制方程：

PG=PL=P

(6)

上式中αG和αL分別表示氣相體積分率和液相體積分率，SLG代表氣液兩相間質(zhì)量傳遞源項(xiàng)，ρG和ρL分別為氣相和液相的密度，vG和vL分別為氣相和液相的速度，μeff,G和μeff,L分別為氣相和液相的有效黏度，PG和PL分別為氣相和液相的壓力。MLG是兩相間動(dòng)量傳遞源項(xiàng)，它直接反映了氣液兩相相間相互作用的動(dòng)量傳遞情況。塔板上氣液兩相的接觸狀態(tài)，多認(rèn)為是氣體以氣泡群形式穿過(guò)塔板上的液層，兩相間的動(dòng)量傳遞主要通過(guò)兩相間的相互作用力來(lái)完成，兩相間作用力主要有曳力、升力和虛擬質(zhì)量力[5]等，其中由于升力和虛擬質(zhì)量力對(duì)塔板模擬結(jié)果影響不大，模擬中通常被忽略，曳力系數(shù)CD常采用Fischer 和Quarini[6]提出的常數(shù)0.44。

塔板模擬常用的湍流模型有標(biāo)準(zhǔn)k-ε模型、Wilcox k-ω 模型、Shear Stress Transport (SST)模型，經(jīng)過(guò)以往模擬與實(shí)驗(yàn)的對(duì)比，本模擬中湍流模型選用精度較高的SST模型來(lái)封閉求解氣液兩相的方程。流體力學(xué)方程以及傳質(zhì)方程離散時(shí)，采用高階離散格式，湍流方程采用一階迎風(fēng)格式，并采用時(shí)間推進(jìn)的全隱式多方程耦合求解線性方程組。目標(biāo)均方根殘差為1.0×10-4。

1.2 幾何建模

以國(guó)內(nèi)某大型丙烷脫氫裝置產(chǎn)品分離塔為對(duì)象建立模型，塔徑為9 200 mm，塔板間距485 mm，采用本公司自有軟件對(duì)該塔進(jìn)行設(shè)計(jì)，塔板上布置梯形固閥，建立模型如圖1。模型主要模擬中間一層塔板的流體流動(dòng)情況，包括兩個(gè)邊降液管和兩個(gè)偏心降液管，降液管按常規(guī)設(shè)計(jì)進(jìn)行建模，塔板上僅布置固閥，未考慮特殊結(jié)構(gòu)。

圖1 六溢流塔板三維建模

因模型較大，考慮到計(jì)算機(jī)硬件及模擬時(shí)間的限制，只針對(duì)1/4塔進(jìn)行模擬，如圖2所示。

圖2 物理模型及邊界條件

分區(qū)域?qū)δＰ瓦M(jìn)行網(wǎng)格劃分，因塔板上方閥孔位置為主要傳質(zhì)區(qū)，因此對(duì)塔板附近的網(wǎng)格進(jìn)行加密劃分，如圖3所示，總網(wǎng)格數(shù)568萬(wàn)。

圖3 模型網(wǎng)格劃分

如圖2所示物理模型及邊界條件，上層降液管作為液相入口，設(shè)置為質(zhì)量入口；氣相入口為下層塔板開孔區(qū)，設(shè)置為速度入口。入口的操作條件按閥孔氣速進(jìn)行輸入，假定入口處速度分布均勻，且垂直于入口平面。頂部對(duì)應(yīng)塔板區(qū)域?yàn)闅庀喑隹冢撞肯聦咏狄汗艹隹跒橐合喑隹?，出口均設(shè)置為壓力出口，假設(shè)流動(dòng)已充分發(fā)展，氣相出口壓力為 0 kPa(表壓)，液相出口壓力為 150 Pa(表壓)。塔壁、塔板及固閥壁面采用無(wú)滑移壁面條件，近壁區(qū)處理方法選用SST模型默認(rèn)的Automatic近壁處理方法。液相出口以上填充一定高度液柱作為初始條件，起到液封作用，以縮短收斂時(shí)間，塔內(nèi)其他區(qū)域均為全氣相，模型主要參數(shù)見表1。

表1 模型主要參數(shù)

2 模擬結(jié)果分析

2.1 六溢流塔盤濃度場(chǎng)分布

圖4和圖5分別為豎直和水平截面上液相體積分率分布云圖。從水平截面云圖可以看出，在邊降液一側(cè)，由于弓形區(qū)較長(zhǎng)，液層堆積嚴(yán)重，相對(duì)較厚，因此導(dǎo)致圖4所示的該區(qū)域出現(xiàn)較為明顯的漏液現(xiàn)象。隨著塔器直徑越來(lái)越大，塔內(nèi)件設(shè)計(jì)中的放大效應(yīng)主要體現(xiàn)在氣液相分布均勻性方面，通過(guò)流體力學(xué)計(jì)算可以較直觀地發(fā)現(xiàn)其中的問(wèn)題。對(duì)于以上問(wèn)題，在設(shè)計(jì)中需采用輔堰及柵欄堰對(duì)該處降液管進(jìn)行特殊設(shè)計(jì)，以實(shí)現(xiàn)液相在塔板上的均勻分布，減少液相滯留區(qū)的存在，進(jìn)而更好地發(fā)揮固閥的傳質(zhì)效率。輔堰在多溢流塔板設(shè)計(jì)中較為常見，多用于邊降液管，在六溢流塔板中，由于塔徑較大，除邊降液管外，三降液的偏心降液管也需要考慮輔堰。柵欄堰利用凸起的柵欄阻擋部分液相通道，縮短了出口堰的實(shí)際有效堰長(zhǎng)，多用于中心降液管或偏心降液管靠近塔中間的部分，可解決降液管不同位置溢流強(qiáng)度不一致的問(wèn)題。除弓形區(qū)存在液相滯留問(wèn)題外，可以看出塔盤上液相體積分離分布比較均勻，證明我們所采用的計(jì)算方法可以使六溢流塔板各鼓泡區(qū)上氣液相分布趨于平均。

圖4 豎直截面上液相體積分率分布云圖

圖5 水平截面上液相體積分率分布云圖

2.2 六溢流塔盤速度場(chǎng)分布

圖6為氣相速度矢量圖，可以看出氣相在塔板上分布比較均勻，閥孔氣速與計(jì)算值接近，約1 m·s-1。由于降液管傾斜段下方已排有閥孔，加之傾斜段的導(dǎo)流，該區(qū)域氣相流動(dòng)到上層塔板后會(huì)在相應(yīng)位置出現(xiàn)氣相集中區(qū)(最高氣速1.5～2 m·s-1)，因此該區(qū)域液層會(huì)相對(duì)較薄。工程設(shè)計(jì)時(shí)，降液管下方需采用專用鼓泡內(nèi)件以實(shí)現(xiàn)對(duì)氣相的導(dǎo)流并減小液體沖擊漏液。

圖6 豎直截面上氣相速度分布矢量圖

2.3 壓力分布云圖

圖7為塔內(nèi)豎直截面上壓力分布云圖，可以看出，模擬所得該層塔板壓降約280～320 Pa。工程設(shè)計(jì)中，塔板壓降分為干板壓降和液層阻力兩部分，計(jì)算模擬壓降較低的原因可能是弓形區(qū)及第一排閥漏液導(dǎo)致模擬液層比計(jì)算值稍低，因此液層阻力也比計(jì)算值小一些。

圖7 豎直截面上壓力分布云圖

3 結(jié)論

作者以實(shí)際項(xiàng)目中的大型板式塔為對(duì)象進(jìn)行建模，模擬了六溢流塔板上的氣液分布狀態(tài)，可以直觀地看到弓形區(qū)液相的滯留以及氣相集中區(qū)的存在。針對(duì)大型板式塔六溢流流場(chǎng)的模擬結(jié)果，對(duì)塔板的設(shè)計(jì)進(jìn)行了優(yōu)化：(1)弓形區(qū)塔板進(jìn)行導(dǎo)流排布且降液管設(shè)計(jì)輔堰和柵欄堰結(jié)構(gòu)，以實(shí)現(xiàn)液相的平均分配并減小弓形區(qū)出現(xiàn)液相滯留區(qū)的可能；(2)降液管下方需采用專用鼓泡內(nèi)件以實(shí)現(xiàn)對(duì)氣相的導(dǎo)流并減小沖擊漏液。該塔目前已順利開車且運(yùn)行狀態(tài)良好。由此可見，在進(jìn)行大型塔器塔內(nèi)件的設(shè)計(jì)時(shí)，采用CFD模擬作為輔助工具，通過(guò)直觀的流場(chǎng)數(shù)據(jù)分析塔板設(shè)計(jì)的合理性，使塔板的大型化設(shè)計(jì)更加高效和準(zhǔn)確。