畢業(yè)設(shè)計(jì)（論文）分離器設(shè)計(jì)

1、氣液分離器設(shè)計(jì)一、概述管柱式氣液旋流分離器是一種帶有傾斜切向入口及氣體、液體出口的垂直管。它依靠旋流離心力實(shí)現(xiàn)氣、液兩相分離，與傳統(tǒng)的重力式分離器相比，具有結(jié)構(gòu)緊湊、重量輕、投資節(jié)省成本等優(yōu)點(diǎn)，是代替?zhèn)鹘y(tǒng)容積式分離器的新型分離裝置。在氣液兩相旋流分析的基礎(chǔ)上，建立了預(yù)測(cè)分離性能的機(jī)理模型，該模型包括了入口分離模型、旋渦模型、氣泡及液滴軌跡模型；依據(jù)機(jī)理模型，提出了管柱式旋流分離器工藝設(shè)計(jì)技術(shù)指標(biāo)和工藝步驟.設(shè)計(jì)根據(jù)管柱式旋流分離器的機(jī)理模型以及設(shè)計(jì)工況，完成了管柱式旋流分離器的結(jié)構(gòu)設(shè)計(jì)、強(qiáng)度分析、理論校核、焊接工藝設(shè)計(jì)以及分離器內(nèi)氣液兩相流的數(shù)值模擬，為工程設(shè)計(jì)和理論設(shè)計(jì)提供一定的理論依據(jù)。2

2、、旋流式分離器的結(jié)構(gòu)及工作原理旋流分離器，是一種利用離心沉降原理將非均相混合物中具有不同密度的相分離的機(jī)械分離設(shè)備。旋流分離器的基本構(gòu)造為一個(gè)分離腔、一到兩個(gè)入口和兩個(gè)出口。分離腔主要有圓柱形、圓錐形、柱-錐形三種基本形式。入口有單入口和多入口幾種，但在實(shí)踐中，一般只有單入口和雙入口兩種。就入口與分離腔的連接形式來分，入口又有切向入口和漸開線入口兩種。出口一般為兩個(gè)，而且多為軸向出口，分布在旋流分離器的兩端?？拷M(jìn)料端的為溢流口，遠(yuǎn)離進(jìn)料端的為底流口。在具有密度差的混合物以一定的方式及速度從入口進(jìn)入旋流分離器后,在離心力場(chǎng)的作用下,密度大的相被甩向四周,并順著壁面向下運(yùn)動(dòng),作為底流排出；密度小

3、的相向中間遷移，并向上運(yùn)動(dòng)，最后作為溢流排出，分離示意圖如圖1。這樣就達(dá)到了分離的目的。旋流分離技術(shù)可用于液液分離、氣液分離、固液分離、氣固分離等。本文設(shè)計(jì)的旋流分離器用于石油鉆井中鉆井液的氣液分離。3、旋流式分離器的優(yōu)缺點(diǎn)在石油化工中裝置中，有各種各樣的分離器，其中以立式重力氣液分離器最為常見，這種氣液分離器具有結(jié)構(gòu)簡(jiǎn)單、操作可靠等持點(diǎn)。立式重力式分離器的主體為一立式圓筒體，多相流一般從該筒體中段進(jìn)入，頂部為氣流出口，底部為液體出口，其結(jié)構(gòu)簡(jiǎn)圖見圖2-2。雖然旋流式氣液分離技術(shù)在石油化工方面的應(yīng)用要晚得多，但與常規(guī)的重力式分離相比較，它具有很多優(yōu)點(diǎn)： = 1 * GB3 分離效率高，由于分離

4、原理的不同使得旋流式分離器具有很高的分離效率； = 2 * GB3 成本低，占用空間較小、維護(hù)費(fèi)用少、能耗低、不需要任何幫助分離的介質(zhì)； = 3 * GB3 安裝靈活方便，旋流器可以任何角度安裝； = 4 * GB3 工作連續(xù)、可靠,操作維護(hù)方便,一旦設(shè)計(jì)、調(diào)試好, 就可自動(dòng)、穩(wěn)定地工作。旋流式氣液分離器有以上優(yōu)點(diǎn)，但也有如下缺點(diǎn): = 1 * GB3 由于旋流器內(nèi)流體的流動(dòng)產(chǎn)生一定的剪切作用，如果參數(shù)設(shè)計(jì)不當(dāng)，容易將液滴(油滴或水滴) 打碎乳化而惡化分離過程； = 2 * GB3 通用性較差。不同的分離要求、不同的處理物料的性質(zhì)往往需要不同結(jié)構(gòu)尺寸或操作條件的旋流器，因此旋流器往往不能互換使

5、用。在欠平衡鉆井中，使用旋流式氣液分離器分離鉆井液中的氣體，能充分發(fā)揮該離器優(yōu)點(diǎn)，同時(shí)又能有效的避免它的缺點(diǎn)。因此，旋流式氣液分離用于分離鉆中的氣體具有廣闊前景。第三章 管柱式氣液旋流分離器29)的四分之一左右。圖3-1 GLCC結(jié)構(gòu)圖 管柱式旋流氣液分離器由入口區(qū)(段)、入口分流區(qū)、漩渦區(qū)、氣泡區(qū)、液滴區(qū)、氣相和液相出口配管等部分組成(圖3-2)。3.2.1 入口區(qū)由于管柱式旋流分離器主要依靠旋流產(chǎn)生的離心力實(shí)現(xiàn)氣液的高效分離,而入口結(jié)構(gòu)決定了分離器的入口氣液相分布及氣液相切向入口速度的大小。因此入口結(jié)構(gòu)是影響管柱式旋流分離器分離特性的關(guān)鍵因素之一。氣液相流速的不同,油、氣兩相或油、氣、水多

6、相流在入口管和噴嘴內(nèi)可能呈現(xiàn)分層流、段塞流、分散氣泡流或環(huán)狀流等多種流型。Kouba的實(shí)驗(yàn)研究表明,采用向下傾斜的入口管,保證入口管流型呈現(xiàn)分層流將在很大程度上改善氣液分離效果、擴(kuò)展管柱式旋流分離器的適用范圍,最佳傾斜角為-27。而傳統(tǒng)分離器采用的入口結(jié)構(gòu)通常為垂直于筒體的結(jié)構(gòu)(目前很多分離器采用的分氣包亦為類似結(jié)構(gòu)),實(shí)驗(yàn)證明采用垂直結(jié)構(gòu)的管柱式旋流分離器,氣液分離效果差,工作范圍大約是傾斜向下的入口結(jié)構(gòu)分離器的一半。入口管向下傾斜,在重力作用下有利于形成分層流,實(shí)現(xiàn)氣液兩相的初步分離。圖3-2 GLCC結(jié)構(gòu)圖3.2.2 入口分流區(qū)入口分流區(qū)即與入口槽連接的筒體部分(如圖3-2)。氣液相經(jīng)入

7、口槽進(jìn)入入口分流區(qū)實(shí)現(xiàn)氣液的初步分離,上部的氣相、下部的液相分別沿筒壁旋轉(zhuǎn)形成旋流場(chǎng)。3.2.3 漩渦區(qū)經(jīng)初步分離的液相以切向速度沿筒壁旋轉(zhuǎn)形成強(qiáng)制旋流。除壁面附近,該旋流可看作剛體轉(zhuǎn)動(dòng),如圖3-3所示,微元控制體以恒定角速度旋轉(zhuǎn)。流動(dòng)是軸對(duì)稱的、各分速度沿軸變化不大,因此假設(shè)微元控制體在軸向、切向上沒有加速度,徑向加速度,應(yīng)用歐拉方程式得徑向和軸向運(yùn)動(dòng)微分方程:圖3-3 旋流控制體受力分析 （1） （2）式中:為含氣泡液相混合物密度;R為筒體半徑,R=d/2?？刂企w壓力微分方程: (3)把方程(1)、(2)代入方程(3),并沿等壓漩渦表面積分()得: (4)整理上式即可得到漩渦方程: (5)

8、 (6)利用漩渦方程積分可得入口分流區(qū)下方液相總體積,應(yīng)用液相體積守恒: (7)式中:指平衡液位高度(見圖3-2)。把方程(5)代入方程(7),積分得到常數(shù),漩渦形狀方程如下:Z= (8)從漩渦方程(8)求解漩渦高度:= (9)3.2.4 氣泡區(qū)和液滴區(qū)分離器上部,氣相為連續(xù)相,液滴分散其中,稱液滴區(qū);下部,液相為連續(xù)相,氣泡分散其中,稱氣泡區(qū)。連續(xù)相做渦旋運(yùn)動(dòng),由于氣液相密度差,分散相粒子(氣泡或液滴)與連續(xù)相間存在滑脫。研究分散相粒子的運(yùn)動(dòng)軌跡,可以分析分離特性。1 .氣泡軌跡分析在漩渦區(qū),較大直徑的氣泡容易被擄獲分出,因此氣泡軌跡的研究區(qū)域是從漩渦底部開始的渦流區(qū)。假設(shè)液相渦流區(qū)小氣泡均

9、勻分布,氣泡的徑向滑脫速度為,軸向滑脫速度為,忽略連續(xù)相-液相的徑向速度,軸向速度即液相折算速度,這樣,氣泡的徑向絕對(duì)速度為 =-,軸向(向下)的絕對(duì)速度為,顯然,實(shí)現(xiàn)氣泡分離的必要條件是壁面處氣泡運(yùn)動(dòng)至中心氣核的時(shí)間應(yīng)少于氣泡隨連續(xù)相液相流出分離器的時(shí)間。2 .液滴軌跡分析 經(jīng)過入口分流區(qū)初步分離后的旋轉(zhuǎn)氣流攜帶液滴進(jìn)入分離器上部的液滴區(qū),與氣泡區(qū)不同的是,由于分散粒子液滴的密度大于連續(xù)相氣相密度,液滴被甩向器壁。假設(shè)徑向液滴滑脫速度為,軸向滑脫速度為,忽略氣相的徑向流速,氣相的軸向速度即折算速度為,這樣,液滴徑向絕對(duì)速度為,軸向(向上)的絕對(duì)速度為,液滴能夠分出的必要條件是分離器中心處的液

10、滴到達(dá)器壁的時(shí)間應(yīng)短于液滴隨旋轉(zhuǎn)氣流流出分離器的時(shí)間。到達(dá)器壁的液滴,在旋轉(zhuǎn)氣流的作用下,將在器壁上形成螺旋狀薄層液流沿器壁向下流動(dòng)分出,完成氣液分離。3 .連續(xù)相旋流特性分離器內(nèi)連續(xù)相旋流場(chǎng)分布對(duì)氣液分離至關(guān)重要。連續(xù)相切向速度、徑向速度和軸向速度構(gòu)成了旋流場(chǎng)。由于粘性耗散和壁面摩擦阻力的作用,切向速度沿軸向逐漸減小,采用旋流強(qiáng)度表征旋流場(chǎng)的衰減特性。旋流強(qiáng)度定義為在某一高度截面上連續(xù)相切向動(dòng)量通量與總軸向動(dòng)量通量的比。 (10)式中:是軸向平均流速,下標(biāo)c表示連續(xù)相。旋流強(qiáng)度沿軸向變化的經(jīng)驗(yàn)關(guān)系式: (11 )式中動(dòng)量通量之比: (12)式中:是連續(xù)相(氣相或液相)入口質(zhì)量流量;是器內(nèi)對(duì)應(yīng)

11、連續(xù)相(氣相或液相)的質(zhì)量流量;指分離器橫截面積;指對(duì)應(yīng)連續(xù)相入口實(shí)際流通面積。根據(jù)實(shí)驗(yàn)與數(shù)值模擬研究,徑向流速通常較切向流速和軸向流速低二個(gè)數(shù)量級(jí),對(duì)氣液分離效率的影響非常小,因此忽略徑向速度。沿旋轉(zhuǎn)半徑方向上,切向速度的變化非常顯著,按照其不同變化規(guī)律,可以將旋流分成兩個(gè)旋轉(zhuǎn)區(qū)域:近壁面處的自由渦旋區(qū)和中心強(qiáng)制渦旋區(qū)。忽略壁面附近的自由渦旋區(qū),將旋流場(chǎng)視為強(qiáng)制渦旋運(yùn)動(dòng),切向流速近似呈線性分布;將軸向流動(dòng)視為勻速運(yùn)動(dòng),軸向流速等于平均流速。對(duì)于工程應(yīng)用,以上的假設(shè)是趨于保守的。切向速度分布: (13)將切向速度表達(dá)式代入旋流強(qiáng)度定義式得: (14)根據(jù)定義是分離器截面位置()的函數(shù),確定了某

12、截面的旋流強(qiáng)度,即可求得壁面處(忽略壁面附近的自由渦旋區(qū))的最大切向速度。 (15)3.3 氣泡、液滴軌跡模型假設(shè):氣泡或液滴分散粒子為球形,在運(yùn)動(dòng)中不變形;粒子間及粒子與器壁間沒有相互作用力;氣泡或液滴的運(yùn)動(dòng)為定常流動(dòng);作用在氣泡或液滴上的各種力的合力為零,分散粒子做勻速運(yùn)動(dòng);分離系統(tǒng)與外界沒有熱交換,視為等溫系統(tǒng);沒有漩渦脫落現(xiàn)象。徑向上,分散粒子受離心力和阻力作用,根據(jù)受力平衡,氣泡徑向滑脫速度: (16)同理,液滴徑向滑脫速度: (17)式中:、為氣泡和液滴直徑;、為氣液連續(xù)相切向速度;、為氣泡和液滴阻力系數(shù)。阻力系數(shù)采用文獻(xiàn)提出的關(guān)系式求解: (18) (19)式中:,軸向上氣泡受重

13、力和阻力的作用,但由于連續(xù)相為油相,柔占度大,滑脫速度小,其流態(tài)一般總處于層流區(qū),因Stokes公式適于求解氣泡的滑脫速度: (20)同樣,液滴在軸向上受重力和阻力作用,根據(jù)受力平衡,液滴軸向滑脫速度: (21)氣泡或液滴所受阻力與氣泡或液滴的合成速度有關(guān),其滑脫速度分別為: (23) (24)根據(jù)上面的氣泡和液滴軌跡分析,時(shí)間間隔中,氣泡或液滴在徑向和軸的位移(見圖3-4)圖3-4 氣泡和液滴軌跡示意圖 消去,即為氣泡和液滴軌跡控制方程,積分得氣泡或液滴的軸向位移, （25） （26）根據(jù)氣泡或液滴實(shí)現(xiàn)分離的必要性,或即為氣泡或液滴區(qū)的最小高度。結(jié)構(gòu)設(shè)計(jì)分析4.1 入口設(shè)計(jì)分析 由于管柱式旋

14、流分離器主要依靠旋流產(chǎn)生的離心力實(shí)現(xiàn)氣液的高效分離,而入口結(jié)構(gòu)決定了分離器的氣液分布及其初始切向入口速度的大小,因此入口結(jié)構(gòu)和尺寸是影響管柱式旋流分離器實(shí)現(xiàn)氣液分離的關(guān)鍵因素。管柱式旋流分離器入口主要由入口管、噴嘴和入口槽3部分組成。氣液相流速的不同,油氣兩相或油氣水多相流在入口管內(nèi)可能呈現(xiàn)分層流、段塞流、分散氣泡流或環(huán)狀流等多種流型。實(shí)驗(yàn)研究表明,采用向下傾斜的入口管,保證入口管流型呈現(xiàn)分層流將在很大程度上改善氣液分離效果,擴(kuò)展管柱式旋流分離器的適用范圍。而傳統(tǒng)分離器采用的入口結(jié)構(gòu)通常為垂直于筒體的結(jié)構(gòu)(目前很多分離器采用的分氣包亦為類似結(jié)構(gòu)),采用垂直結(jié)構(gòu)的管柱式旋流分離器實(shí)驗(yàn)證明,與傾斜

15、向下的入口結(jié)構(gòu)相比,氣液分離效果明顯變差,工作范圍大約減小一半。入口管向下傾斜,在重力作用下有利于形成分層流,實(shí)現(xiàn)氣液兩相的初步分離,同時(shí),向下的傾斜結(jié)構(gòu)使經(jīng)過初步分離的液相在入口下方旋轉(zhuǎn)一圈后形成旋流場(chǎng),避免了對(duì)氣相向分離器上方運(yùn)動(dòng)的阻塞。故入口管采用傾斜入口. 入口管傾角以-27為宜,管長(zhǎng)取1.01.5,入口管直徑的選取應(yīng)保證流型為分層流,由Taitel和Dukler預(yù)測(cè)模型確定,分層流轉(zhuǎn)變?yōu)殚g歇流或環(huán)狀流的判別準(zhǔn)則為: （27）式中： 式中D為入口管直徑,是它的橫截面積,h是無(wú)量綱液位高。迭代求解準(zhǔn)則方程,D和h作為迭代參量,直至準(zhǔn)則方程左端小于1,din即為滿足分層流條件最小入口管直徑

16、。噴嘴是入口段最后一個(gè)影響進(jìn)入分離器氣液相流速分布和入口切向流速大小的因素。通過對(duì)3種不同入口開槽結(jié)構(gòu)(矩形、同心圓形及新月形)的初步實(shí)驗(yàn)發(fā)現(xiàn),同心圓形噴嘴(縮口管)結(jié)構(gòu)的分離特性最差,而矩形結(jié)構(gòu)噴嘴的分離效果最佳,新月形結(jié)構(gòu)噴嘴的分離效果與矩形噴嘴接近,由于矩形槽結(jié)構(gòu)加工困難,推薦采用新月形結(jié)構(gòu)。入口噴嘴面積的選取應(yīng)保證入口液相流速在4.5-6/之間。液相流速過小將難以發(fā)揮旋流離心分離的作用,但液相流速過大將形成過高的漩渦區(qū),在筒體中過早出現(xiàn)氣相夾帶液滴和液相夾帶氣泡現(xiàn)象,影響分離效果。4.1.3雙入口分析 雙傾斜入口將入口流預(yù)分為兩股流動(dòng)：低入口的富液流和高入口的富氣流。雙入口的試驗(yàn)表明中

17、等大小的氣體流量（在入口處段塞流轉(zhuǎn)為分層流）下，氣體帶液率有明顯降低，當(dāng)氣體流量較高時(shí)（在入口處為環(huán)空流），無(wú)多大變化。4.2.1入口位置對(duì)于沒有液位控制的GLCC,將入口段定位于靠近液面的上方是至關(guān)重要的。最新的許多試驗(yàn)都表明,單入口GLCC的最佳液面大約在距離入口下方13L/d處。過低的液面,如距離入口處遠(yuǎn)大于3L/d,會(huì)導(dǎo)致切向入口速度的過度衰減,影響GLCC的性能。如果液面高于入口,氣體會(huì)通過液體而溢出,造成更多液體的攜帶。4.2.2最佳外形比外形比是指GLCC的長(zhǎng)徑比。GLCC的尺寸影響其性能及造價(jià)。對(duì)于一個(gè)給定的直徑,GLCC中入口上方的長(zhǎng)度提供了液流擾動(dòng)的容量,而入口下方的長(zhǎng)度則

18、決定了用于從液體中分離氣泡的存留時(shí)間。另外,離心力和浮力的大小與直徑成反比,切向速度衰減與長(zhǎng)度成正比。由于這一現(xiàn)象的復(fù)雜性,最近才剛剛提出了一套決定最佳外形比的基本標(biāo)準(zhǔn)。4.2.3旋流體錐度針對(duì)反錐型、正錐型和圓柱型旋流體進(jìn)行的研究表明,對(duì)于氣/液分離,圓柱型旋流體要稍優(yōu)于反錐型和正錐型結(jié)構(gòu)。4.3 出口管設(shè)計(jì)分析氣液相出口管線的配置可根據(jù)氣液相流量、配置儀表的要求確定。建議氣相出口流速取330/,液相出口流速取1.212/。容器內(nèi)平衡液位應(yīng)低于入口0 .3,分離氣液相的匯合點(diǎn)一般低于入口點(diǎn)0.30.5,以保持正的靜水壓頭。若分離器配置控制系統(tǒng),匯合點(diǎn)位置可以高于入口。分離器工藝設(shè)計(jì)5.1 計(jì)

19、算分離器直徑考慮分離器上部的氣相分離部分,分離器直徑的選取應(yīng)避免氣相中夾帶液滴,以氣相折算速度表示,即氣相折算速度不能大于氣流中出現(xiàn)液滴時(shí)的臨界速度。氣相臨界速度是: （28）式中We是無(wú)因次Weber數(shù),它決定于液滴的尺寸,這里取值7。對(duì)于分離器入口以下的液相分離部分,應(yīng)充分發(fā)揮離心分離特性,避免液相中夾帶氣泡。研究表明,保持液相入口切向流速和液相流速的比為40時(shí),旋流分離效果最佳。通常切向流速一般取6/,顯然臨界液相流速=0.15/。對(duì)于高氣油比的分離工況,分離器直徑： （29）對(duì)于低氣油比的分離工況,分離器直徑： （30）式中和分別是分離工況下氣體流量和液體流量。根據(jù)入口噴嘴分析,取液相

20、入口切向速度為:=6m/s由上分析知當(dāng)/=40時(shí), 旋流分離效果最佳,故液相流速:已知:日處理液量:900 m33/s, 日處理氣量:90m3/d=0.00104 m3/s,又知本設(shè)計(jì)工況是低氣油比工況,故由公式（30）得分離器的直徑： =2(m)經(jīng)圓整后取 d=300mm5.2分離器的高度計(jì)算以入口為界分離器分為上下兩個(gè)部分。經(jīng)過旋流離心分離,部分液體可能以旋流液膜的形式向上爬升或以液滴的形式隨氣流向上運(yùn)動(dòng)產(chǎn)生氣相夾帶液滴的現(xiàn)象。上半部分的高度應(yīng)足夠高,甚至在極端流型段塞流工況下能夠吸收液相流量的波動(dòng),避免氣相夾帶液滴的現(xiàn)象。參考傳統(tǒng)立式分離器的沉降分離段的處理方法,結(jié)合現(xiàn)場(chǎng)應(yīng)用經(jīng)驗(yàn),管柱式

21、旋流分離器的上半部分高度一般取11.5,根據(jù)要求分離器的日處理量,設(shè)計(jì)半部分的高度為：h=1.00m.因經(jīng)初步分離的液體在入口下方旋轉(zhuǎn)一周后才形成旋流，在這個(gè)過程中形成入口分流區(qū)。在入口分流區(qū)內(nèi)流體作螺旋運(yùn)動(dòng)，所需的時(shí)間t為：t=入口分流區(qū)高度為： h=已知 由公式（9）得 旋渦區(qū)高度： = =下部液相空間的高度選取,應(yīng)保證液相在分離器內(nèi)有足夠的停留時(shí)間使氣泡得以分出進(jìn)入上部氣相空間。根據(jù)實(shí)踐經(jīng)驗(yàn),當(dāng)筒體直徑小于1英寸時(shí),液相空間高度取11.5。對(duì)于大直徑的旋流分離器,根據(jù)建立的氣泡軌跡模型,求解氣泡(氣泡直徑取500)自器壁進(jìn)入中心氣核實(shí)現(xiàn)分離在軸向上穿過的距離,顯然液相空間的高度應(yīng)不小于。

22、因?yàn)榉蛛x器的設(shè)計(jì)直徑是300 mm,故應(yīng)根據(jù)建立的氣泡軌跡模型,求解下部液相空間的高度.氣泡在GLCC內(nèi)穿過的整個(gè)軸向距離計(jì)算公式: 已知 , 由公式（20）得氣泡軸向滑脫速度:氣泡徑向速度分布: （30）式中混合物的密度分布： （31）液相切向速度分布： （32）式中為切向入口速度,m和n是指數(shù),取0.9，是氣泡的直徑。根據(jù)Turton和Levenspie建議阻力系數(shù)為 : （33）雷諾數(shù): （34）因?yàn)榕cr之間存在嵌套關(guān)系，通過與r的關(guān)系式不能直接得出兩者之間的關(guān)系，所以采用試算法確定兩者間的關(guān)系?；舅悸肥牵菏紫仁窃谝粋€(gè)確定的r上給定一給徑向速度v1，計(jì)算出相應(yīng)的 ,，Re(r),v2,

23、然后比較，來確定對(duì)應(yīng)點(diǎn)上氣泡的徑向速度。用此方法把半徑等分成10個(gè)點(diǎn)，通過編程算出其速度。計(jì)算程序如下：#include#includemain() float v1,v2,v3,r,R,R1,p1,p2,p,s,m,D,n1,n2,n3,n4,n5,cd;R=0.15,p1=1.185,p2=1500.00,s=0.05, D=0.0005;scanf(%f,%f,&v1,&r);printf(nv1=%f,r=%f,v1,r);n1=pow(r/R,0.9);p=p1+(p2-p1)*n1;v3=6*n1;R1=p*v1*D/s;n2=pow(R1,0.657);n3=pow(R1,-1.

24、09);cd=24*(1+0.173*n2)/R1+0.413/(1+16300*n3);v2=sqrt(4*(p-p1)*v3*v3*D/(3*p2*r*cd);m=v1-v2;if(m=0) printf(nv1=%f,r=%f,v1,r);else do v1=v1-2*m/3; R1=p*v1*D/s; n4=pow(R1,0.657); n5=pow(R1,-1.09); cd=24*(1+0.173*n4)/(1+16300*n5); v2=sqrt(4*(p-p1)*v3*v3*D/(3*p2*r*cd); while(fabs(v1-v2)1e-6); printf(nv1=%

25、f,r=%f,v1,r); 運(yùn)行結(jié)果如下：r=0.應(yīng)用excel繪出10點(diǎn)的關(guān)系曲線如下圖5-1所示： 圖5-1根據(jù)圖示v1與r的曲線關(guān)系，v1和r之間近似指數(shù)關(guān)系，故設(shè) （35）對(duì)上式兩邊同時(shí)取對(duì)數(shù)得： （36）記 ,b=lgk 則rvr*rlgvrlgv關(guān)于曲線分析圖如下所示：圖5-2由的曲線關(guān)系知近似直線關(guān)系，故采用最小二乘法確定常數(shù)a,b.由方程組： （37） （38）即 （39） （40）把已知代入方程組得：0.086625a+0.825b=0.284776 (41) 0.825a+10b=12.165652 (42)解方程組得：a=15.24,b=-3.13 (43)氣泡在GLCC

26、內(nèi)穿過的整個(gè)軸向距離: GLCC各部分計(jì)算高度總結(jié)如下表：液滴區(qū)（m ）入口分流區(qū)(m)旋渦區(qū)(m)氣泡區(qū)(m)總高(m) = 1 * GB3 徑向方向上旋流中氣泡受力如圖5-3示，由受力可知， 圖5-3 相對(duì)運(yùn)動(dòng)微方程： (47) 式中：為氣泡質(zhì)量，為相對(duì)滑移速度，為氣泡直徑，為半徑r處的旋流速度。假設(shè)旋流為強(qiáng)迫旋渦，即。則上式可寫為： (48)解得： (為常數(shù)) (49)初始條件：時(shí)， (50)設(shè)為時(shí)間常數(shù)，式中第一項(xiàng)，故趨于終端沉降速度：又，故 （t=0時(shí)，）若不考慮旋轉(zhuǎn)時(shí)的能量損失，則 （為切向入口速度）， 解得：一般地，當(dāng)時(shí)，近似認(rèn)為氣泡已遷移到中心，對(duì)應(yīng)的時(shí)間為最小駐留時(shí)間 = 1

27、 * GB3 一般氣泡從邊壁到中心的平均移動(dòng)速度 ，即 = 2 * GB3 分離器的處理量為Q，則分離器內(nèi)液體占據(jù)的最小體積 = 3 * GB3 式中V為旋流器的容積。 = 2 * GB3 豎直方向上計(jì)算在最小駐留時(shí)間內(nèi)混合體運(yùn)動(dòng)的距離時(shí)，可忽略氣泡與液體的相對(duì)滑移，認(rèn)為氣泡隨液體一起在重力的作用下向下運(yùn)動(dòng)，可得在內(nèi)，氣泡下降的距離： = 4 * GB3 即旋流器的最小長(zhǎng)度為。設(shè)計(jì)參數(shù)： 由以上參數(shù)和式 = 1 * GB3 = 2 * GB3 = 3 * GB3 = 4 * GB3 計(jì)算可得：旋流器的半徑。當(dāng)分離器的設(shè)計(jì)直徑時(shí)，相應(yīng)地，圖5-4 平衡液位示意圖 氣液分離器如圖5-3所示，來自井

28、口的氣液混合物進(jìn)入分離器后產(chǎn)生旋流，在離心力的作用下實(shí)現(xiàn)分離，氣體占據(jù)液滴區(qū)并經(jīng)排氣管從氣體出口排出，而液體占據(jù)氣泡區(qū)并經(jīng)液體管從泥漿出口排出。在分離器工作的過程中，可能出現(xiàn)兩種極端情況：（1）氣量特別大，以至氣體壓力p1足以克服排液管線的能量的消耗，使氣體以連續(xù)相從泥漿出口排出；（2）液量特別大，為了克服排液管線的能量消耗，液面上升，有可能出現(xiàn)液相從排氣管線排出。在設(shè)計(jì)設(shè)計(jì)氣液分離器時(shí)，根據(jù)設(shè)計(jì)參數(shù)的要求，要求避免這兩種情況的發(fā)生。這就需要研究氣液相在輸送過程的能量平衡，要求正確的確定液封高度z2和氣液混合物的進(jìn)口高度。5.4.1 氣室壓力p1的確定在管中流動(dòng)的氣體由于溫度和壓力沿管長(zhǎng)變化,

29、其流速和密度也會(huì)有顯著的變化.根據(jù)流體力學(xué)理論穩(wěn)態(tài)流動(dòng)的能量方程有:+ +d=0 (51)通用氣體定律給出: pV=RT或= (52)設(shè)計(jì)參數(shù)所對(duì)應(yīng)的壓力為p、溫度T、體積V,那么: =將上式帶入(2),得 pV=T故管內(nèi)流速等于= (53)將式(3)和式(2)代入式(1),有 pdp ()dp + ()Td=0假設(shè)流動(dòng)為等溫的,積分上式得: p1- 2()ln()=p+() (54)根據(jù)Weymouth提出的公式計(jì)算,即 = (55) 上列各式中:-泥漿密度-氣室內(nèi)絕對(duì)壓力-排氣管出口絕對(duì)壓力-排氣管內(nèi)徑-排氣管長(zhǎng)度M-摩爾質(zhì)量m-氣體質(zhì)量R-氣體常數(shù)T-華氏溫度-出氣管路沿程阻力系數(shù)設(shè)排氣

30、管的直徑DG=20mm,此時(shí)，則= 已知日處理氣量為:90m/d=0.00104 m/s，M=29，R=8.315代入（55）式得 （56）編程計(jì)算程序如下：#include#includemain()double m,n,p,f;clrscr();for(p=1.01;p1.5;p+=0.00001)m=p/1.01;n=log(m);f=p*p-0.267*n-1.192;if(f0.選擇第一斷面為分離器內(nèi)液面,第二斷面為泥漿出口處,建立兩斷面的能量平衡關(guān)系如下:+Z1+=+Z2+(1+) （57） 式中：p-2斷面壓力 Z1-1斷面處液柱高度Z2-2斷面處液柱高度 -1斷面處泥漿流速-2

31、斷面處泥漿流速 D-排液管內(nèi)徑L-排液管長(zhǎng)度 -出液管路沿程阻力系數(shù)由（57）式變形得：Z1-Z2= +(1+) - （58）當(dāng)DL=40mm時(shí)，R代入（58）式得，顯然不滿足條件。當(dāng)DL=80mm時(shí)，R代入（58）式得，為了保證液體的順利排除必須保持一定的壓頭，同時(shí)排液管排出的液體相當(dāng)于自由流，因此要盡可能的使Z1-Z2，故所取DL過大。當(dāng)DL=65時(shí)，R代入（58）式得，當(dāng)DL=50mm時(shí)，R代入（58）式得，由因?qū)嶋H的工況，如果選擇次方案要求的高度太高。故綜上所述，當(dāng)DL=65mm時(shí)比較合理的。故取Z1=1 .9m,Z2,管長(zhǎng)取1,入口直徑的選取是為了保證入口管內(nèi)呈現(xiàn)分層流,入口管直徑根

32、據(jù)Taitel&Duker模型確定。由分流層轉(zhuǎn)變?yōu)閿嗬m(xù)流或環(huán)狀流的判別式: （59）式中： 式中D為入口管直徑，為管內(nèi)橫截面上液面離管底的高度，A為管內(nèi)氣相所占的橫截面積。D和作為迭代量，直至準(zhǔn)則方程左端小于1，D即為滿足分層流條件最小入口直徑。對(duì)判別式變形得: （60）已知 代入上式得: （61）用C語(yǔ)言編程求解最小進(jìn)口管直徑，計(jì)算程序如下#include#includestdio.hmain()float m,k,d,n=0.1; while(n0.04&d0.3&n1) printf(nn=%f,d=%f,n,d);運(yùn)行結(jié)果:n=0.710000 ,n,n,n,n,n,n,n,n,n,5

33、9n,n,n,n,n,n,n,n,n,n,n,n,故計(jì)算求得最小直徑: D=40mm入口槽選型,根據(jù)分離器工藝設(shè)計(jì)分析入口槽選用新月形,同時(shí)入口噴嘴的截面積的確定保證入口液相流速為6 m/s.綜合考慮本設(shè)計(jì)分離器的結(jié)構(gòu)，最后確定為單進(jìn)口。本設(shè)計(jì)中GLCC工藝尺寸綜合參數(shù)如下表：分離工況GLCC尺寸入口管GLCC回路（m3/s）（m3/s）氣液比工作壓力（Pa）直徑（mm）上部高度（mm）下部高度（mm）直徑（mm）排液管直徑（mm）排氣管直徑（mm）1:103001000406520第六章 強(qiáng)度設(shè)計(jì)1. 設(shè)計(jì)參數(shù)如下表：設(shè)計(jì)壓力（MPa）最大工作壓力 （MPa） 工作溫度容器尺寸(mm)工作介

34、質(zhì)常溫（立式）有毒、易燃首先根據(jù)單層圓筒厚度計(jì)算公式： （62） 式中 計(jì)算厚度，mm 設(shè)計(jì)直徑 焊接接頭系數(shù) 計(jì)算壓力，MPa設(shè)計(jì)溫度下許用應(yīng)力由設(shè)計(jì)參數(shù)，筒體的材料是：16MnR.根據(jù)文獻(xiàn)26查表D1知，在設(shè)計(jì)溫度下16MnR.的許用應(yīng)力，當(dāng)厚度為6-16mm時(shí)，=170MPa,當(dāng)厚度為16-36時(shí)，=163MPa因?yàn)槿萜魉⒔橘|(zhì)易燃，所以壓力容器所有的焊縫采用全焊透結(jié)構(gòu)，根據(jù)文獻(xiàn)26查表4-3知 設(shè)計(jì)壓力：假設(shè)計(jì)算厚度為6-16mm,許用應(yīng)力=170MPa，則=設(shè)計(jì)厚度：對(duì)于16Mn鋼板偏差因而可以取名義厚度但是對(duì)低合金鋼制容器，規(guī)定不包括腐蝕余量的最小厚度應(yīng)力判別式： （63）式中 計(jì)

35、算壓力，MPa 設(shè)計(jì)直徑 設(shè)計(jì)溫度下許用應(yīng)力 -有效厚度已知 則代入上式得故滿足強(qiáng)度要求封頭采用標(biāo)準(zhǔn)橢圓形封頭.材料:16MnR,當(dāng)厚度為616mm時(shí), ,當(dāng)厚度為1636mm時(shí), 1. 橢圓形狀系數(shù): （64）式中:-橢圓長(zhǎng)軸-直邊短軸而對(duì)于標(biāo)準(zhǔn)橢圓形封頭,故=2. 計(jì)算橢圓的厚度: （65）式中： k-橢圓形狀系數(shù)pc-設(shè)計(jì)壓力Di- -圓筒直徑-設(shè)計(jì)溫度下的許用應(yīng)力-焊接接頭系數(shù)假設(shè)計(jì)算壁厚為616mm,則=170MPa,于是有顯然壁厚在假設(shè)范圍內(nèi),則設(shè)計(jì)壁厚對(duì)于16MnR,鋼板負(fù)偏差,故取名義厚度.但對(duì)于低合金鋼制壓力容器規(guī)定不包括腐蝕余量的最小壁厚不小于3mm,加上2mm的腐蝕余量,

36、名義厚度至少應(yīng)取5mm.又由鋼材標(biāo)準(zhǔn)規(guī)格,名義厚度應(yīng)取為 6mm.根據(jù)應(yīng)力強(qiáng)度判別式: （66）式中: pc-設(shè)計(jì)壓力- -圓筒直徑- 有效壁厚-設(shè)計(jì)溫度下的許用應(yīng)力把已知條件代入得: 故強(qiáng)度滿足要求:式中: -設(shè)計(jì)溫度下的許用應(yīng)力 -焊接接頭系數(shù),對(duì)封頭要求100%的全焊縫 -有效壁厚Di- -圓筒直徑k-橢圓形狀系數(shù)代入已知量得: =故封頭滿足最大工作壓力的要求考慮到設(shè)計(jì)壓力較高，介質(zhì)不允許泄漏，擬選用平焊法蘭，凸凹密封面。法蘭材料：16Mn(鍛件)，初擬法蘭的尺寸如下所示D D=100 D=70 D=50 H=10 R=9墊片尺寸：墊片寬度：N=墊片基本密封寬度：墊片有效密封寬度：故墊片

37、的計(jì)算寬度：螺栓選用的材料是35號(hào)鋼，預(yù)緊時(shí)螺栓載荷：操作時(shí)螺栓載荷： 擬采用螺栓4個(gè) ，則，根據(jù)螺紋標(biāo)準(zhǔn)以及管道法蘭螺栓的最小直徑，取螺紋根徑，相當(dāng)于M10的螺栓。實(shí)際螺栓總面積：螺栓平均間距：s=根據(jù)文獻(xiàn)28查表5-10，對(duì)于M12-的螺栓，最小螺栓間距為允許最大螺栓間距：故，即所選螺柱直徑符合安裝和密封要求預(yù)緊時(shí)螺柱載荷：W= M=.m操作時(shí)墊片載荷：、壓力載荷：介質(zhì)靜壓軸向載荷：作用在法蘭端面上的總力矩： 法蘭應(yīng)力及度校核法蘭形狀系數(shù)：根據(jù)文獻(xiàn)28查表（5-28），（5-29），（5-30），（5-31）得：于是 軸向應(yīng)力：MPa徑向應(yīng)力：周向應(yīng)力：MPa擬定法蘭尺寸及選材合適，可安全

38、使用因?yàn)楣艿婪ㄌm的內(nèi)徑符合工稱直徑系列，應(yīng)選用標(biāo)準(zhǔn)管道法蘭。根據(jù)參考文獻(xiàn)24管道法蘭標(biāo)準(zhǔn)選用公稱壓力為1.6MPa帶頸平焊法蘭，材料：16Mn,常溫下允許的工作壓力為1.6MPa.管道法蘭標(biāo)記：管法蘭MFM1.6-22 SH3406-92. 考慮到設(shè)計(jì)壓力較高，介質(zhì)不允許泄漏，擬選用平焊法蘭，凸凹密封面。法蘭材料：16Mn(鍛件)，初擬法蘭的尺寸如下所示D D=150 D=150 D=80 墊片尺寸：墊片寬度：N=墊片基本密封寬度：墊片有效密封寬度： 故墊片的計(jì)算寬度：螺栓選用的材料是35號(hào)鋼，預(yù)緊時(shí)螺栓載荷：操作時(shí)螺栓載荷： 擬采用螺栓4個(gè) ，則，根據(jù)螺紋標(biāo)準(zhǔn)以及管道法蘭螺栓的最小直徑，取螺

39、紋根徑，相當(dāng)于M16的螺栓。實(shí)際螺栓總面積：螺栓平均間距：s=根據(jù)文獻(xiàn)28查表5-10，對(duì)于M16的螺栓，最小螺栓間距為允許最大螺栓間距：故，即所選螺柱直徑符合安裝和密封要求預(yù)緊時(shí)螺柱載荷：W= M=.m操作時(shí)墊片載荷：、壓力載荷：介質(zhì)靜壓軸向載荷：作用在法蘭端面上的總力矩： 法蘭應(yīng)力及度校核法蘭形狀系數(shù)：根據(jù)文獻(xiàn)28查表（5-28），（5-29），（5-30），（5-31）得：于是 軸向應(yīng)力：徑向應(yīng)力： 周向應(yīng)力：MPa擬定法蘭尺寸及選材合適，可安全使用PPa.管道法蘭標(biāo)記：管法蘭MFM1.6-49.5 SH3406-92. 6.4計(jì)算設(shè)備重量載荷 設(shè)備操作重量： （67）-容器殼體和支座的

40、重量-容器內(nèi)部構(gòu)件的重量-容器保溫材料的重量-平臺(tái)、扶梯的重量-操作時(shí)容器內(nèi)物料的重量-人孔、接管、法蘭等附件的重量(2)計(jì)算分離器各段載荷，所取截面如圖1所示圖6-1各段載荷計(jì)算如下表所示各段載荷0-11-22-3（kgf）計(jì)算風(fēng)載荷: （68） 以1-2段為例計(jì)算風(fēng)載荷：-設(shè)備所在地的基本風(fēng)載荷，取烏魯木齊的基本風(fēng)載荷進(jìn)行計(jì)算=60-高度變化系數(shù)，由2表16-6得-1-2段分離器高度，=2m-空氣動(dòng)力系數(shù)，對(duì)圓筒形設(shè)備-風(fēng)振系數(shù)，-計(jì)算段頂截面距地面高度系數(shù), 由2表16-4得 -設(shè)備基本振型自振周期變化系數(shù)因?yàn)楸驹O(shè)計(jì)氣液分離器是等直徑厚壁設(shè)備，故由2表16-5得=1，-直立設(shè)備有效直徑，

41、-容器各段外直徑，mm-容器各段保溫層厚度，mm-籠式扶梯當(dāng)量寬度，mm-操作平臺(tái)當(dāng)量寬度，mm=m=同理可求出分離器各段風(fēng)載荷于下表所示計(jì)算各段載荷0-11-22-31（kgf）602平臺(tái)數(shù)0110(2)計(jì)算風(fēng)彎矩容器任意計(jì)算截面I-I的風(fēng)彎矩計(jì)算公式: （69）分離器各截面的風(fēng)彎矩計(jì)算結(jié)果如下表:截面彎矩0-0 1-12-2任意截面的地震彎矩: （70）等直徑、等壁厚設(shè)備任意截面計(jì)算截面I-I地震彎矩: （71）底部截面的地震彎矩: （72）(當(dāng)時(shí)，視設(shè)備為柔性結(jié)構(gòu)，則需考慮高振型的影響，在進(jìn)行穩(wěn)定或其他驗(yàn)算時(shí)，所取的地震彎矩值應(yīng)為上列計(jì)算值的1.25倍)-地震影響系數(shù)，這里取最大地震時(shí)影

42、響系數(shù)Q0 -自計(jì)算截面以上的操作載荷-塔頂至第 I段(含第i段)的高度,mh-計(jì)算截面距地面的高度,m-計(jì)算截面以上集中重量QX的作用點(diǎn)距地面的高度,m- 度處,集中載荷引起的水平地震載荷,N各截面的地震載荷計(jì)算如下表:截面各截面的地震彎矩0-0，故分離器為柔性構(gòu)件1-12-2-偏心彎矩因?yàn)闆]有偏心載荷，故設(shè)計(jì)壓力引起的軸向應(yīng)力: （73）操作或者非操作時(shí)重力引起的軸向應(yīng)力: （74）彎矩引起的軸向引力: （75）式中： -設(shè)計(jì)壓力-分離器直徑-分離器各段壁厚-任意計(jì)算截面處以上筒體承受的操作或非操作時(shí)的重量載荷，kgf中較大者各種載荷引起的軸向引力如下表:設(shè)計(jì)壓力引起的軸向力 操作重量引起

43、的軸向力 16MPa MPa最大彎矩引起的軸向力 MPaMPaMPa最大組合軸向拉力，出現(xiàn)正常操作的情況下： （76）強(qiáng)度條件： （77）-焊接系數(shù)由以上分析知筒體最大組合軸向拉力，發(fā)生在正常操作時(shí)的1-1截面MPa=1701.00=170MPa 故筒體的強(qiáng)度滿足要求最大組合軸向壓應(yīng)力，出現(xiàn)在停車情況下： （78）穩(wěn)定條件：中的較小值 （79）計(jì)算B值。 MPaMPa故截面2-2滿足穩(wěn)定性要求MPaMPa故截面1-1滿足穩(wěn)定性要求故截面0-0滿足穩(wěn)定性要求綜上所述,分離器滿足穩(wěn)定性要求第七章 焊接工藝設(shè)計(jì)1.筒體與筒體以及筒體與封頭焊縫焊接接頭型式和尺寸選用UG24(HG20583-1998)

44、.2.根據(jù)GB/T14957-947選用焊絲的牌號(hào)H10Mn;根據(jù)GB12470-70選用焊劑的牌號(hào)HJ431型號(hào)HJ401-H08A. 3.焊接采用埋弧焊,對(duì)焊縫進(jìn)行100%的射線探傷檢測(cè),要求符合GB3325-87中的級(jí)為合格. 筒體與接管的焊接接頭型式和尺寸選用G2 GB20583-1998,其示結(jié)構(gòu)如圖2所示.根據(jù)GB/T518-95低合金鋼電焊條選用焊條牌號(hào)J502，型號(hào)E5003。焊接采用手工電弧焊,對(duì)焊縫進(jìn)行100%的射線探傷檢測(cè),要求符合GB3325-87中的級(jí)為合格. 帶頸平焊法蘭與接管焊接接頭尺寸選用F6(JB47004703-92),其結(jié)構(gòu)如圖4所示。2.根據(jù)GB/T51

45、8-95低合金鋼電焊條選用焊條牌號(hào)J507，型號(hào)E5015。3.焊后對(duì)焊縫進(jìn)行100%的射線探傷檢測(cè),要求符合GB 3325-87中的級(jí)為合格.第八章 氣液兩相流場(chǎng)的數(shù)值模擬8.1 數(shù)值計(jì)算方法簡(jiǎn)介 計(jì)算流體力學(xué)作為流體力學(xué)研究中的一門新興分支，正在工業(yè)和科研領(lǐng)域內(nèi)發(fā)揮越來越重要的作用。將CFD工具運(yùn)用到分離機(jī)械的研究中，也成為工程技術(shù)人員改進(jìn)設(shè)計(jì)、提高效率的有效手段，是CFD應(yīng)用的前沿。一些成熟的算法，模型也以商業(yè)軟件的形式出現(xiàn)在工程及科研領(lǐng)域。相比研究單位自行開發(fā)的計(jì)算程序，商業(yè)計(jì)算軟件一般具有以下特點(diǎn)： 通用性廣。由于商業(yè)軟件面向的用戶對(duì)象廣泛，處理的實(shí)際問題多種多樣，因此其覆蓋的應(yīng)用范

46、圍要盡可能廣。計(jì)算穩(wěn)定性好。多數(shù)軟件經(jīng)過不同研究領(lǐng)域內(nèi)的算例測(cè)試，對(duì)不同類型的問題具有較好的適應(yīng)能力。使用方便，商業(yè)軟件都提供了比較友好的用戶界面，方便用戶的使用。一般商業(yè)軟件也存在一些明顯的不足，例如：算法相對(duì)陳舊，不能緊跟CFD研究領(lǐng)域內(nèi)的最新成果；與不同行業(yè)內(nèi)的實(shí)際要求存在一定的距離，難以將各研究單位已有的研究成果結(jié)合到商業(yè)軟件中。這在一定程度上限制了商業(yè)軟件在工程實(shí)際中的應(yīng)用。本文運(yùn)用Fluent軟件對(duì)離心是分離器的內(nèi)流場(chǎng)進(jìn)行分析計(jì)算，F(xiàn)luent公司是享譽(yù)世界的最大計(jì)算流體力學(xué)（CFD）軟件供應(yīng)商，F(xiàn)luent軟件能夠精確地模擬無(wú)粘流、層流、湍流、化學(xué)反應(yīng)、多相流等復(fù)雜的餓流動(dòng)現(xiàn)象。

47、應(yīng)用領(lǐng)域包括：航空航天、汽車設(shè)計(jì)、生物醫(yī)藥、化學(xué)處理、石油天然氣，發(fā)電系統(tǒng)、電子半導(dǎo)體、蝸輪設(shè)計(jì)、HVAC、玻璃加工等。FLUENT具有精度高，收斂快，穩(wěn)定性好等特點(diǎn)，同時(shí)可通過添加擁護(hù)自定義的函數(shù)（UDF）解決實(shí)際具體問題。Gambit是前置處理器，能針對(duì)極其復(fù)雜的幾何外行生成三維四面體，六面體的非結(jié)構(gòu)化網(wǎng)絡(luò)及混合網(wǎng)絡(luò)。該模塊還具有方便的網(wǎng)絡(luò)檢查功能，對(duì)網(wǎng)絡(luò)單元體積、扭曲率、長(zhǎng)細(xì)比等影響收斂和穩(wěn)定的參數(shù)進(jìn)行統(tǒng)計(jì)并生成報(bào)告。8.1.1 控制方程對(duì)于所有的流動(dòng)問題，F(xiàn)LUENT需要求解質(zhì)量和動(dòng)量守恒方程。對(duì)于熱傳導(dǎo)或可壓縮流動(dòng)，需要解能量守恒的附加方程。對(duì)于包括組分混合和反應(yīng)的流動(dòng)，需要解組分守

48、恒方程或者使用PDF模型來解混合分?jǐn)?shù)的守恒方程。當(dāng)流動(dòng)是湍流是，還要解附加的輸運(yùn)方程。FLUENT可以在慣性坐標(biāo)系（無(wú)加速坐標(biāo)系）和具有加速度的參考坐標(biāo)系（旋轉(zhuǎn)坐標(biāo)系）中建立流動(dòng)模型。在旋轉(zhuǎn)坐標(biāo)系中，通過建立一個(gè)與旋轉(zhuǎn)設(shè)備一起運(yùn)動(dòng)的相對(duì)坐標(biāo)系來建摸，近似認(rèn)為流體旋轉(zhuǎn)角速度為常數(shù)，旋轉(zhuǎn)邊界相對(duì)與參考系靜止。在旋轉(zhuǎn)坐標(biāo)系中，絕對(duì)速度v或相對(duì)速度v的關(guān)系如下： （80）其中 角速度向量（即旋轉(zhuǎn)坐標(biāo)系的角速度）； 旋轉(zhuǎn)坐標(biāo)系中的位置向量 旋轉(zhuǎn)坐標(biāo)系中的質(zhì)量守恒方程（連續(xù)性方程）： （81）其中源項(xiàng)上假如到連續(xù)性相的第二相質(zhì)量（比方說由于液滴的蒸發(fā)，質(zhì)量發(fā)生變化），源項(xiàng)也可以是自定義源項(xiàng)。慣性坐標(biāo)系中的

49、動(dòng)量方程： （82）其中p是靜壓，是應(yīng)力張量，和分別是中立體積力和外部體積力（如離散相相互作用的產(chǎn)生的升力）。包含了其它的模型相關(guān)源項(xiàng)，如多孔介質(zhì)和自定義源項(xiàng)。應(yīng)力張量 （83） 其中：分子粘性； I單位張量。旋轉(zhuǎn)坐標(biāo)系中的動(dòng)量方程為： （84）由于在FLUENT中忽略了 項(xiàng)，因此不能用動(dòng)量方程的相對(duì)速度表達(dá)式準(zhǔn)確的計(jì)算隨時(shí)間變化的角速度。湍流是由大小不同尺度的旋渦組成，對(duì)時(shí)間和空間都是非線性的隨機(jī)運(yùn)動(dòng)。它最本質(zhì)的特征是“湍動(dòng)”，既隨機(jī)的脈動(dòng)。湍流流場(chǎng)是無(wú)數(shù)不同尺寸渦旋相互摻混的流動(dòng)場(chǎng)。湍流的出現(xiàn)影響著整個(gè)流場(chǎng)的速度、壓力、溫度和物質(zhì)濃度的分布。目前，關(guān)于紊流的數(shù)值計(jì)算可分為細(xì)觀模擬和統(tǒng)觀模擬

50、。完全模擬、大渦流模擬屬于細(xì)觀模擬。完全模擬（Direct Numberical Simulation DNS）在湍流尺度的網(wǎng)格尺寸內(nèi)求解N S方程而不使用任何湍流模型。該方法必須采用很少的時(shí)間和空間步長(zhǎng)。如文獻(xiàn)估算，長(zhǎng)。如文獻(xiàn)估算，對(duì)一個(gè)渦旋進(jìn)行數(shù)值計(jì)算，至少要設(shè)置十個(gè)節(jié)點(diǎn)，這樣對(duì)于一個(gè)小尺寸范圍內(nèi)的紊流運(yùn)動(dòng)要在1cm 的流場(chǎng)中布置10個(gè)節(jié)點(diǎn)。顯然，完全模擬在短期尚無(wú)法用于求解實(shí)際工程中的復(fù)雜湍流問題。大渦流模擬（Large Eddy Simulation LES）是在大渦尺度的網(wǎng)格尺寸內(nèi)求結(jié)NS方程，由于計(jì)算量仍很大，只能模擬一些簡(jiǎn)單流動(dòng)，如彎道等，目前也不能直接用于工程。而基于求解雷諾時(shí)

51、均方程的模擬，即統(tǒng)觀模擬，利用某些假設(shè)將雷諾時(shí)均方程中高階未知關(guān)聯(lián)項(xiàng)或者時(shí)均量來表達(dá)，從而使雷諾時(shí)均方程封閉。因此統(tǒng)觀模擬當(dāng)前成為解決工程實(shí)際問題的有效手段。所謂湍流模型理論就是以雷諾平均運(yùn)動(dòng)方程為基礎(chǔ)，依靠理論與經(jīng)驗(yàn)引進(jìn)一些模擬假設(shè)，建立一組描述湍流平均值的封閉方程組的計(jì)算方法。該模型的平均行為，應(yīng)與實(shí)際的湍流統(tǒng)計(jì)平均行為基本一致。 1925年P(guān)randtl提出的動(dòng)量傳遞理論以及后來提出的自由剪切層模型、泰勒的渦量傳遞理論的馮-卡門的相似性理論等一系列半經(jīng)驗(yàn)理論，其基本思想都是建立關(guān)于雷諾應(yīng)力的模型假設(shè)，使雷諾平均運(yùn)動(dòng)方程得以封閉。由于只考慮了一階湍流統(tǒng)計(jì)量的動(dòng)力學(xué)微分方程，即平均運(yùn)動(dòng)方程，

52、僅引人附加的代數(shù)關(guān)系而沒有引進(jìn)任何高階統(tǒng)計(jì)量的微分方程，因此屬于零方程模型。這種模型對(duì)運(yùn)動(dòng)的預(yù)測(cè)性很差，不適用于有回流的復(fù)雜運(yùn)動(dòng)。Kolmogorov和Prandtl提出的單方程模型雖較半經(jīng)驗(yàn)理論有所改進(jìn)，且Cd和Cu值也較容易確定，但L值的確定并不比混合長(zhǎng)度的確定容易。因而單方程模型同樣只適用于簡(jiǎn)單流動(dòng)，不適用于帶回流的復(fù)雜流動(dòng)。在所有的雙方程模型中，K-雙方程模型的應(yīng)用最為普遍，先后由周培源（1945）、Harlow-Nakayama(1968) 、Jones-Launder(1972)提出。在進(jìn)一步簡(jiǎn)化的模型中，人們放棄給雷諾應(yīng)力建立方程的想法，將它們直接用推廣的Bossinesq渦粘性

53、模式來表示。大量的預(yù)報(bào)及實(shí)驗(yàn)結(jié)果對(duì)照表明，K-模型可以成功或基本成功用于以下幾種情形：無(wú)浮力平面射流、平壁邊界層、管流、通道流或噴管內(nèi)流動(dòng)、無(wú)旋渦及弱旋的二維和三維流動(dòng)。由于它采用了同向性湍流疏運(yùn)的假設(shè)，故它不適應(yīng)具有非同向性湍流輸運(yùn)的強(qiáng)旋流。在離心分離器內(nèi)流場(chǎng)中，因切速度遠(yuǎn)遠(yuǎn)大于徑向和軸向分量，因而，通常采用雷諾應(yīng)力模型（RSM）、代數(shù)應(yīng)力模型（ASM）或RNG模型來代替K-模型。雷諾應(yīng)力模型RSM對(duì)雷諾應(yīng)力及通量項(xiàng)采用微分方程直接求解，具有很大通用性，這一模型的優(yōu)點(diǎn)在于可準(zhǔn)確地考慮各向異性效應(yīng)。代數(shù)應(yīng)力模型（ASM）是雷諾應(yīng)力模型在一定條件下的簡(jiǎn)化表達(dá)式，表達(dá)式形式隨簡(jiǎn)化條件而異，由應(yīng)力代

54、數(shù)表達(dá)式加上及方程構(gòu)成，因此又稱為擴(kuò)展的K-模型或者2個(gè)半方程模型。應(yīng)用該模型可避開求解雷諾應(yīng)力方程所面臨的十分復(fù)雜的計(jì)算工作，又較好地預(yù)報(bào)了流動(dòng)的各向異性的特點(diǎn)。RNG K-模型是一種修正的K-模型，在文獻(xiàn)3334中有較詳細(xì)的討論。它從原始的基本方程推導(dǎo)而來，其中使用了所謂的“Kolmogorov”數(shù)學(xué)技巧。Kolmogorov定律：E（k）=k 。RNG模型是一個(gè)更一般，更基本的模型，尤其對(duì)強(qiáng)旋流流場(chǎng)及高曲率流線的離心分離器有著很好的改進(jìn)效果。在FLUENT中有多種湍流模型可選擇，包括常用的Spalart-Allmaras單方程模型、標(biāo)準(zhǔn)K-模型、和K-o模型、雷諾應(yīng)力模型以及大旋渦模型。

55、選擇合適的湍流模型是數(shù)值計(jì)算中及其重要的環(huán)節(jié)，不同湍流模型的使用會(huì)產(chǎn)生不同的流場(chǎng)細(xì)節(jié)。本文應(yīng)用RNGK-模型預(yù)測(cè)離心式分離器的內(nèi)流場(chǎng)。多相流的數(shù)值計(jì)算方法有兩種：歐拉拉格朗日（Euler-Lagrange）方法和歐拉歐拉（Euler- Euler）方法。 在Fluentzhong中的拉格朗日離散項(xiàng)模型遵循歐拉拉格朗日方法。流體相為連續(xù)相，直接求解時(shí)均N-S方程，而離散相通過計(jì)算流場(chǎng)中大量的粒子，氣泡或是液滴的運(yùn)動(dòng)得到。離散相和流體相之間有動(dòng)量、質(zhì)量和能量的交換。該模型假設(shè)離散相（第二相）的體積比率很低。粒子或液滴運(yùn)行軌跡的計(jì)算是獨(dú)立的，被安排在流相計(jì)算的指定間隙完成。這樣處理能較好的符合噴霧干

56、燥，煤和液體燃料燃燒，以及一些粒子負(fù)載的流動(dòng)情況，但是不適合用于液-液混合物，流化床和其它第二相體積率不容忽略的情形。歐拉模型在歐拉歐拉方法中，不同的相被處理成互相貫穿的連續(xù)介質(zhì)。由于一種相所占的體積不能被其他相占有，故引入相體積分?jǐn)?shù)（phasic volume fraction）。體積分?jǐn)?shù)是時(shí)間和空間的連續(xù)函數(shù)，各相的體積分?jǐn)?shù)之和等于1。從各相的守恒方程可以推導(dǎo)出一組方程，這些方程對(duì)于所有的相都具有類似的形式。從實(shí)驗(yàn)得到的數(shù)據(jù)可以建立一些特定的關(guān)系，從而使上述方程封閉，另外，對(duì)于小顆粒流（granular flows），可以通過應(yīng)用分子運(yùn)動(dòng)論的理論使方程封閉。在FLUENT中，有三種歐拉歐拉

57、多相流模型：流體體積模型（VOF），混合（Mixture）模型以及歐拉（Euler）模型?；旌夏Ｐ秃蜌W拉模型主要用于模擬相間的混合和分離。（1）流體體積模型（VOF） VOF模型通過求解單獨(dú)的動(dòng)量方程和處理穿過區(qū)域的每一相的體積分?jǐn)?shù)來模擬兩種或三種不能混合的流體。典型的應(yīng)用包括預(yù)測(cè)流體中大氣泡的運(yùn)動(dòng)和氣液界面的穩(wěn)態(tài)和瞬態(tài)變化。（2）混合模型（Mixture） 混合模型求解混合物的動(dòng)量方程，并通過相對(duì)速度描述離散相。混合模型是歐拉模型在幾種情況下的很好替代。當(dāng)顆粒相廣泛分布或界面規(guī)律未知時(shí)，完善的多相流模型是不切實(shí)際的。在這種情況下，混合模型能取得較好的結(jié)果。 混合模型的應(yīng)用包括低負(fù)載的粒子負(fù)載

58、流，氣泡流，沉降，以及旋風(fēng)分離等多相流。混合模型也可用于沒有相對(duì)速度的均勻多相流。（3）歐拉模型（Eulerian） 歐拉模型是Fluent中最復(fù)雜的多相流模型。它建立含有n個(gè)動(dòng)量方程和連續(xù)方程的方程組來求解各相。壓力項(xiàng)和各界面交換系數(shù)耦合在一起，耦合方式依賴于所含相的情況。顆粒流（液固）與非顆粒流（液液）的處理是不同的。對(duì)于顆粒流，可應(yīng)用分子運(yùn)動(dòng)理論求得流動(dòng)特性。不同相之間的動(dòng)量交換也依賴于混合物的類別。該模型的應(yīng)用包括氣泡柱，上浮，顆粒懸浮，以及流化床。 解決多相流問題，先選擇最能符合實(shí)際的流體模式。然后根據(jù)不同的模式，選擇恰當(dāng)?shù)亩嘞嗔髂Ｐ汀Ｈ缭跉庖夯蛞阂憾嘞嗔髦?，大致有泡狀流，彈狀流，?/p>

59、塞流及自由液面流四種模式。對(duì)于離心分離器，主要是分離在來連續(xù)液相中分散的小氣泡（含氣率為10%左右），氣泡的運(yùn)動(dòng)一般認(rèn)為是泡狀流，可選用混合模型或歐拉模型。以下為混合模型和歐拉模型的求解策略。 啟動(dòng)混合模型求解，采用0.2或更小的滑流速度欠松弛因子開始計(jì)算。如果解顯示出好的收斂趨勢(shì)，可逐漸增加欠松弛因子。初始計(jì)算是可以不求解題解分?jǐn)?shù)和滑流速度方程，當(dāng)計(jì)算收斂后，在求解這些方程。 啟動(dòng)歐拉模型求解，為了提高收斂性，在求解多相流模型前可以先獲得初始解再繼續(xù)計(jì)算。初始解可通過混合模型得到。由于歐拉模型對(duì)內(nèi)存及收斂要求較高，在下面的數(shù)值計(jì)算中，多相流模型選用混合模型。8.1.4 數(shù)值計(jì)算方法目前關(guān)于數(shù)

60、值計(jì)算的方法有：有限差分法、有限元法、邊界元法等。在CFD中，有限差分法占主導(dǎo)地位，且最為成熟，目前已經(jīng)發(fā)展了多種收斂性好、精度高的離散格式，較常用的有：Taylor（泰勒）展開法、有限容積法。與有限差分法相比，有限單元法在流體力學(xué)應(yīng)用還遠(yuǎn)遠(yuǎn)不夠成熟，易出現(xiàn)數(shù)值溢出、數(shù)值發(fā)散等問題。Flent采用有限容積離散動(dòng)量方程，速度和壓力耦合采用SIMPLE、SIMPLEC及PISO算法。對(duì)流項(xiàng)差分格式納入了一階迎風(fēng)、中心差分及QUICK等格式。在多種情況下，中心格式的計(jì)算穩(wěn)定性好于迎風(fēng)格式。代數(shù)方程可以采用多重網(wǎng)格及最小殘差法。對(duì)于旋轉(zhuǎn)坐標(biāo)系中的流動(dòng)問題，求解過程經(jīng)常變得不穩(wěn)定。這是因?yàn)楫?dāng)旋轉(zhuǎn)項(xiàng)的影響