蔡 祿 武曉波 孫治謙 耿 坤 王旱祥 王振波

(1.中國(guó)石油大學(xué)(華東)新能源學(xué)院 2.中國(guó)石油大學(xué)(華東)機(jī)電工程學(xué)院)

0 引 言

在天然氣水合物排水采氣的工程應(yīng)用中，氣液分離系統(tǒng)具有重要作用。氣液旋流分離器作為核心結(jié)構(gòu)，其內(nèi)部流場(chǎng)對(duì)分離效果具有重要影響[1]。在分離過(guò)程中，起初液相以液滴形式存在于氣流中；在高速旋轉(zhuǎn)運(yùn)動(dòng)中，液滴由于受到較大的離心力作用而從氣流中分離，在旋流腔內(nèi)壁面處聚集形成液膜，在重力作用下向下運(yùn)動(dòng)，最終在旋流器排液口排出[2]。

蝸殼式旋流器內(nèi)的流場(chǎng)呈高速旋轉(zhuǎn)流動(dòng)，RSM模型具有較高的計(jì)算精度，適合模擬強(qiáng)旋流動(dòng)[3-5]。謝駿遙等[6]通過(guò)對(duì)比3種湍流模型發(fā)現(xiàn)RSM模型更接近試驗(yàn)結(jié)果。關(guān)于氣液兩相流分離的數(shù)值模擬，前人大多采用多相流模型(Multiphase Model)來(lái)研究液量較多的情況。劉哲等[7]研究了含氣量和湍流擴(kuò)散效應(yīng)對(duì)流場(chǎng)分布的影響；馬粵等[8]采用多相流模型探究液位和處理量對(duì)氣液分離器分離效率的影響。若氣液比很大，液相占比小于10%的情況下，液相以離散的狀態(tài)存在于氣流中，則應(yīng)當(dāng)采用DPM模型[9-10]。吳允苗等[11-12]采用DPM模型來(lái)分析旋流分離器內(nèi)的流場(chǎng)分布；金向紅等[13]采用DPM模型對(duì)旋流器內(nèi)液滴間的碰撞、液滴的破碎和碰壁等過(guò)程進(jìn)行數(shù)值模擬計(jì)算，完善了氣液旋流器的分離機(jī)制和模型；周云龍等[14]分別利用RSM模型和DPM模型對(duì)不同入口高度和寬度的氣液旋流分離器進(jìn)行了數(shù)值模擬。DPM模型將液滴看作離散相顆粒，忽略了液滴的聚并和破碎，進(jìn)而忽略了液膜的形成過(guò)程。這種近似大大簡(jiǎn)化了模擬過(guò)程，但也存在一定誤差。因此，近年來(lái)有學(xué)者不斷開(kāi)發(fā)新的數(shù)學(xué)模型近似模擬氣液比較大的強(qiáng)旋流場(chǎng)。YANG L.L.等[15]采用DPM模型研究了GLCC內(nèi)液滴的聚并和破碎趨勢(shì)，并基于液滴受力分析和流體力學(xué)，建立了液滴遷移模型，準(zhǔn)確預(yù)測(cè)了旋流器的分離性能；岳題等[16-17]采用EWF(Eulerian Wall Film)模型和DPM模型結(jié)合的方法進(jìn)行GLCC內(nèi)氣液兩相數(shù)值模擬，實(shí)現(xiàn)了上部筒體的上行旋流液膜(USLF)流動(dòng)和液滴攜帶的研究。EWF模型能夠?qū)崿F(xiàn)氣液兩相流耦合模擬壁面液膜的分布和流動(dòng)特性[18]，但前人大多針對(duì)GLCC旋流器，很少將該模型應(yīng)用于雙蝸殼氣液旋流分離器中。

為此，筆者采用EWF模型將氣液旋流分離器的流場(chǎng)看作由氣流、液滴和液膜三相構(gòu)成，既能夠模擬實(shí)際工況中較大的氣液比，又能合理考慮液滴的行為，通過(guò)CFD數(shù)值模擬得到雙蝸殼氣液旋流分離器內(nèi)速度分布和壓力分布，并結(jié)合流場(chǎng)分析液膜的形成及運(yùn)動(dòng)特性。所得結(jié)果可為天然氣水合物排水采氣裝置的研發(fā)和應(yīng)用提供理論依據(jù)。

1 物理模型

本文擬對(duì)雙蝸殼氣液旋流分離器進(jìn)行數(shù)值模擬。建立雙蝸殼氣液旋流分離器的幾何模型，如圖1所示，并采用ICEM進(jìn)行網(wǎng)格劃分；分別采用50萬(wàn)、60萬(wàn)、70萬(wàn)、80萬(wàn)和90萬(wàn)的網(wǎng)格進(jìn)行入口-溢流壓降驗(yàn)證，結(jié)果如圖2所示。

圖1 雙蝸殼旋流分離器模型結(jié)構(gòu)Fig.1 Model structure of double volute cyclone separator

由圖2可知，當(dāng)網(wǎng)格數(shù)量達(dá)到70萬(wàn)即可得到準(zhǔn)確的計(jì)算結(jié)果。

雙蝸殼氣液旋流分離器具有良好的分離性能。蝸殼式入口結(jié)構(gòu)具有造旋作用，使流體在蝸殼內(nèi)進(jìn)行預(yù)分離，因此蝸殼式旋流器具有壓降低的特點(diǎn)[19-20]。內(nèi)部流場(chǎng)的切向速度分布與入口結(jié)構(gòu)密切相關(guān)，在0°～180°區(qū)間蝸殼流道空間逐漸縮小，切向速度升高；在180°～360°區(qū)間流道面積不變，一部分氣流進(jìn)入柱段，氣量減少，切向速度降低。單蝸殼式旋流器的流場(chǎng)呈現(xiàn)明顯非對(duì)稱性，徑向速度由于數(shù)值較小，對(duì)偏心距的變化很敏感，可能造成離心或數(shù)值變大的假象。雙蝸殼式旋流分離器的入口由2個(gè)單蝸殼相對(duì)排列組合而成，雙蝸殼入口的徑向力大大減小，流場(chǎng)呈現(xiàn)良好的對(duì)稱性[21]。氣液兩相流體從雙蝸殼入口進(jìn)入旋流器后，向下旋轉(zhuǎn)運(yùn)動(dòng)至柱段，排氣心管下端伸入柱段上方，有效削弱了短路流。柱段下方連接漸縮段，漸縮段直徑逐漸減小，便于收集分離出來(lái)的液相。旋流腔直徑100 mm，高485 mm。

2 數(shù)學(xué)模型及邊界條件

本文采用針對(duì)液膜開(kāi)發(fā)的EWF壁面液膜模型與DPM離散相模型結(jié)合的方法進(jìn)行計(jì)算。不考慮能量傳遞，液膜流動(dòng)的質(zhì)量方程為：

(1)

液膜流動(dòng)的動(dòng)量方程為[16]：

(2)

(3)

(4)

設(shè)置連續(xù)相流體為密度1.225 kg/m3的空氣，離散相為密度998.2 kg/m3的液態(tài)水。

邊界條件設(shè)置：雙入口設(shè)置為velocity-inlet；溢流口與大氣連通，因此設(shè)置為pressure-outlet；底流口封閉用于收集液相，設(shè)置為wall；其余壁面均設(shè)置為wall。

3 數(shù)值模擬結(jié)果及分析

3.1 數(shù)值模擬準(zhǔn)確性驗(yàn)證

為驗(yàn)證數(shù)值模擬的準(zhǔn)確性，將模擬得到的分離效率曲線與前人試驗(yàn)結(jié)果進(jìn)行對(duì)比。王建軍等[19]所研究的雙蝸殼旋風(fēng)管與本文研究對(duì)象結(jié)構(gòu)相似、尺寸接近，將其試驗(yàn)結(jié)果和本文的數(shù)值模擬結(jié)果進(jìn)行對(duì)比(見(jiàn)圖3)，以驗(yàn)證本文模型的正確性。

圖3 數(shù)值模擬準(zhǔn)確性驗(yàn)證Fig.3 Verification of numerical simulation accuracy

由圖3可知，試驗(yàn)結(jié)果和數(shù)值模擬結(jié)果差距最大點(diǎn)位于氣相流量220 m3/h處，最大計(jì)算誤差為11.9%。誤差產(chǎn)生的原因：本文所用物理模型為DPM+EWF模型，考慮了液膜的形成但沒(méi)有研究液滴的破碎；并且采用理想液滴粒徑分布，而試驗(yàn)過(guò)程中無(wú)法控制液滴粒徑呈規(guī)律性分布；再者，由于旋流器內(nèi)流場(chǎng)的不穩(wěn)定性，壓降可能在一定范圍內(nèi)波動(dòng)。所以試驗(yàn)結(jié)果和模擬結(jié)果存在誤差合理。觀察曲線變化趨勢(shì)基本相同，且誤差在可接受范圍內(nèi)，認(rèn)為數(shù)值模擬采用的模型和計(jì)算方法合理。

3.2 旋流器內(nèi)速度場(chǎng)分布

以氣液旋流分離器柱段中心軸線方向?yàn)閦軸，溢流口方向?yàn)檎较颍x取如圖4所示的7個(gè)(柱段截面z=228、172、116 mm，橫截面柱段與漸縮段交界面z=83 mm，漸縮段截面z=53和23 mm，底流口截面z=0)進(jìn)行速度場(chǎng)分析。所選截面涵蓋了排氣心管插入深度段、柱段和漸縮段，具有較強(qiáng)的代表性。

圖4 流體域橫截面位置分布Fig.4 Cross section of fluid domain

圖5 不同截面處的速度分布Fig.5 Velocity distribution at different cross sections

圖5為不同截面處的速度分布圖。圖5a和圖5b分別為柱段和漸縮段截面處切向速度分布圖。切向速度是產(chǎn)生離心力和離心加速度的主要?jiǎng)恿Γ趦上喾蛛x過(guò)程中起到重要作用。柱段上半部分(z=228和172 mm 2個(gè)截面)受排氣管底部的影響，半徑為0.03～0.04 m處排氣管兩側(cè)速度驟降。在不受排氣管影響的柱段下半部分和漸縮段，各個(gè)截面處切向速度沿直徑分布呈“M”形的軸對(duì)稱分布，且位置越靠近漸縮段，峰值速度越穩(wěn)定，因此流動(dòng)越穩(wěn)定。由于切向速度越大越有利于分離，因此柱段和漸縮段是分離的關(guān)鍵區(qū)域。

圖5c和圖5d分別為柱段和漸縮段截面處軸向速度沿直徑方向分布圖。軸向速度以指向溢流口方向?yàn)檎较?。靠近幾何中心處軸向速度為正，靠近邊壁處軸向速度為負(fù)。內(nèi)旋流為上行流，方向?yàn)檎瑥囊缌骺谂懦?；外旋流為下行流，方向?yàn)樨?fù)，最終匯聚于排液口。兩者交界形成的零速包絡(luò)面較規(guī)整，其直徑與排氣心管直徑較接近。

圖5e和圖5f分別為柱段和漸縮段截面處徑向速度沿直徑方向分布圖。徑向速度相對(duì)于其他2個(gè)分速度數(shù)值最小。從壁面到軸心沿半徑方向呈現(xiàn)先增大后減小的趨勢(shì)，柱段中下段和漸縮段徑向速度為負(fù)，指向軸心方向，有利于連續(xù)相(氣相)的聚集和排出。柱段上段的徑向速度指向邊壁方向，能夠攜帶液滴向排氣心管壁面運(yùn)動(dòng)，有利于液相的二次分離。

3.3 旋流器內(nèi)壓力場(chǎng)分布

本文模擬了不同氣相流量下旋流器內(nèi)的壓力分布情況。圖6為入口與溢流口壓力之差(溢流口壓降)隨入口氣相流量的變化曲線。由圖6可以看出，壓力差隨著入口氣量增大而增大，與前人研究結(jié)論相符，即壓降近似與入口氣速的二次方成正比。

圖6 溢流口壓降隨入口氣相流量的變化曲線Fig.6 Pressure drop at overflow outlet vs.gas flow rate at inlet

由于流場(chǎng)呈對(duì)稱分布，壓力場(chǎng)也具有良好的對(duì)稱性。圖7為氣相流量Q=440 m3/h時(shí)流場(chǎng)壓力分布云圖。由圖7可知，從雙蝸殼入口到排氣心管處整體壓力不斷降低，入口段和柱段旋流腔之間存在一次明顯的壓力降低。這是由于雙蝸殼入口的結(jié)構(gòu)形式對(duì)兩相流體起到預(yù)分離作用。氣流從入口段進(jìn)入旋流腔，流道突然增大，流速降低，壓力減小。柱段和漸縮段沿直徑方向，從外旋流向內(nèi)旋流壓力減小，結(jié)合此處的速度分布得知，柱段和漸縮段是兩相流體分離的關(guān)鍵部位。

圖7 Q=440 m3/h時(shí)旋流器內(nèi)壓力分布云圖Fig.7 Cloud chart for pressure distribution in cyclone at Q=440 m3/h

3.4 液膜的形成及運(yùn)動(dòng)

EWF模型可以用來(lái)預(yù)測(cè)壁面上液膜的產(chǎn)生和流動(dòng)，適用于液膜厚度較薄的情況，因此可用于旋流器內(nèi)液膜的計(jì)算，計(jì)算結(jié)果如圖8～圖16所示。

圖8 液膜厚度隨時(shí)間的積累圖Fig.8 Accumulation of liquid film thickness with time

圖9 液膜速度隨時(shí)間的變化圖Fig.9 Variation of liquid film velocity with time

結(jié)合圖8和圖9可以看出，液膜的積累形式與旋流場(chǎng)的流動(dòng)軌跡相吻合，氣流中攜帶的液滴在高速旋轉(zhuǎn)運(yùn)動(dòng)的過(guò)程中不斷被分離出來(lái)，一邊向旋流腔內(nèi)壁運(yùn)動(dòng)；一邊在重力的作用下向下運(yùn)動(dòng)，在柱段的中下部積累形成液膜。液膜厚度較大的位置，即液膜初步形成的位置具有較大的液膜速度。較大的液膜速度有利于液膜向下遷移，在2.5 s時(shí)均勻分布于漸縮段。隨著液膜的質(zhì)量不斷增加，漸縮段的液膜不斷積累，在3.0 s時(shí)沿旋轉(zhuǎn)軌跡增大液膜厚度。此時(shí)液膜的速度較大，切向速度和軸向速度占主導(dǎo)，有利于液膜在底流口收集。

圖10 z=144 mm處截面液膜厚度隨時(shí)間的變化曲線Fig.10 Variation of liquid film thickness at the z=144 mm cross-section with time

圖11 z=43 mm處截面液膜厚度隨時(shí)間的變化曲線Fig.11 Variation of liquid film thickness at the z=43 mm cross-section with time

結(jié)合圖10可以看出：在液膜形成初期，在柱段中下部，90°、135°和270°、315°對(duì)稱位置開(kāi)始出現(xiàn)液膜；隨著液膜的不斷積累，在45°和225°對(duì)稱位置液膜厚度不斷增加；增加到一定程度后，液膜在重力的作用下向下運(yùn)動(dòng)。當(dāng)液膜運(yùn)動(dòng)到漸縮段，如圖11所示，前2.0 s液膜厚度較小且均勻，鋪滿整個(gè)漸縮段；2.5 s后液膜厚度突然增加并維持在0.06～0.08 mm。

分析圖12得到，柱段的液膜厚度分布仍與氣流的旋轉(zhuǎn)流動(dòng)軌跡相吻合，每一橫截面處液膜厚度存在峰值，除峰值以外的其余位置液膜的厚度很小。液膜厚度在柱段的中下部出現(xiàn)峰值，即圖12中z=116 mm處，位于排氣管插入深度以下、柱段與漸縮段交界處以上。觀察圖13可得，漸縮段的液膜厚度整體波動(dòng)小，分布均勻。結(jié)合圖8可知，漸縮段中下部存在液膜厚度較大的2段區(qū)域，厚度約為0.2 mm，其余區(qū)域液膜厚度較小，穩(wěn)定在0.05～0.10 mm。由于漸縮段的截面直徑逐漸減小，越靠近排液口位置的液膜厚度起伏越小，液膜平均厚度越大，越有利于液相的收集。

圖12 計(jì)算穩(wěn)定后柱段不同截面處液膜厚度Fig.12 Liquid film thickness at different cross sections of column section after calculation stabilization

圖13 計(jì)算穩(wěn)定后漸縮段不同截面處液膜厚度Fig.13 Liquid film thickness at different cross sections of reducing section after calculation stabilization

由圖14可知：壁面處液膜的軸向速度在10-2m/s數(shù)量級(jí)，在重力作用下液膜以較低的速度穩(wěn)定下移，最終可在排液口收集。排氣心管處液膜的運(yùn)動(dòng)方向?yàn)樯仙\(yùn)動(dòng)；經(jīng)過(guò)旋流分離后，未能即時(shí)分離的液相在排氣管內(nèi)壁繼續(xù)旋轉(zhuǎn)運(yùn)動(dòng)并生成厚度極小的液膜，液膜厚度如圖15所示。排氣心管內(nèi)壁面的液膜起初在管中部沿旋轉(zhuǎn)軌跡形成，并逐漸向排氣管上部擴(kuò)展。此時(shí)在溢流口收集到的液膜厚度約0.006 mm，這是由于液膜在達(dá)到較大厚度后可能滑落或者被再次吹起。溢流口液膜越薄對(duì)分離越有益。

圖14 壁面和排氣心管處液膜的軸向速度分布云圖Fig.14 Cloud chart for axial velocity distribution of liquid film at wall and exhaust pipe

圖16 柱段液膜厚度和分離效率隨入口流量的變化曲線Fig.16 Variation of column section liquid film thickness and separation efficiency with inlet flow rate

由圖16可以看出：隨著氣相流量的增加，液膜厚度不斷增大；且入口氣相流量越大，液膜厚度增大速度逐漸緩慢。分析可知，當(dāng)氣相流量較小時(shí)，隨著流體速度的增大，液相的轉(zhuǎn)速增大，液滴的碰撞概率增加，因此聚并的概率增大，液滴不斷被甩向壁面后被液膜捕捉，液膜厚度不斷增加。液膜厚度的形成有利于分離效率的提升[22]，入口氣相流量越大，即入口氣速越大，液滴越容易聚并形成液膜，分離效果越好。當(dāng)氣量增大到一定程度后，圖16中所示氣量大于400 m3/h后，液膜厚度增加速度緩慢。分析原因，其一此時(shí)液膜厚度已經(jīng)增大到0.08 mm，液膜運(yùn)動(dòng)加劇，向底流口方向液膜厚度較小的位置移動(dòng)；其二隨著氣速的增加，已經(jīng)形成的液膜表面可能會(huì)因氣速過(guò)大而被再次吹起，因此液膜厚度沒(méi)有顯著增大。

4 結(jié) 論

(1)本文所得速度場(chǎng)分布情況與前人研究結(jié)果相符。切向速度在三相速度場(chǎng)中占主導(dǎo)，切向速度越大，越有利于分離。柱段上半部分受排氣心管底部的影響，排氣管近壁面處的切向速度驟減；柱段下半部分切向速度呈“M”形的軸對(duì)稱分布。軸向速度呈外圍下行流，中心上行流的分布趨勢(shì)。漸縮段三相速度均較穩(wěn)定，數(shù)值相對(duì)較小且呈現(xiàn)良好的規(guī)律性。徑向速度以向心運(yùn)動(dòng)為主，數(shù)值最小。

(2)入口-溢流壓降隨入口流量(氣速)的增大而增大，近似與入口氣速的二次方成正比。壓力最低處位于排氣心管底部中心處。柱段和漸縮段壓力呈軸對(duì)稱分布，外旋流壓力高，內(nèi)旋流壓力低。

(3)柱段和漸縮段流場(chǎng)分布較穩(wěn)定，是兩相分離的關(guān)鍵部位，液膜也在此段形成。柱段中下部液膜初步形成的軌跡與旋轉(zhuǎn)流動(dòng)軌跡重合，逐漸運(yùn)動(dòng)至漸縮段后，液膜厚度逐漸趨于均勻。因此液膜厚度最大值出現(xiàn)在柱段中下部，而漸縮段的液膜厚度平均值較大。液膜初步形成的位置具有較大的液膜運(yùn)動(dòng)速度。