鼓泡塔內(nèi)空氣-醋酸體系流體力學(xué)參數(shù)的CFD-PBM耦合模型數(shù)值模擬

1、鼓泡塔內(nèi)空氣-醋酸體系流體力學(xué)參數(shù)的CFD-PBM耦合模型數(shù)值模擬鼓泡塔反應(yīng)器因操作方便、結(jié)構(gòu)簡(jiǎn)單、傳質(zhì)和傳熱性好等優(yōu)點(diǎn)，被廣泛地應(yīng)用于石油化工、煤化工、環(huán)境工程和食品工程等領(lǐng)域1，如費(fèi)-托合成、環(huán)己烷氧化、甲醇合成和對(duì)苯二甲酸合成等2。近年來(lái)，我國(guó)已經(jīng)成為精對(duì)苯二甲酸（PTA）的生產(chǎn)大國(guó)，其中對(duì)二甲苯（PX）氧化反應(yīng)器是PTA生產(chǎn)過(guò)程中的關(guān)鍵部分，其決定了PTA的產(chǎn)品質(zhì)量和生產(chǎn)速率3-4。反應(yīng)物性與操作條件是影響鼓泡塔內(nèi)流體力學(xué)的主要因素，同時(shí)其氣含率和氣泡尺寸分布是設(shè)計(jì)和分析鼓泡塔反應(yīng)器的重要參數(shù)，因此，研究有機(jī)體系下鼓泡塔流體力學(xué)特性對(duì)反應(yīng)器的設(shè)計(jì)與放大具有重要意義。目前，關(guān)于氣液鼓泡塔

2、數(shù)值模擬的研究仍處于發(fā)展階段5。同時(shí)，雙流體模型中破碎和聚并模型較為復(fù)雜、參數(shù)較多，使得鼓泡塔的數(shù)值模擬成為一個(gè)充滿機(jī)遇與挑戰(zhàn)的研究領(lǐng)域。研究者采用計(jì)算流體力學(xué)（CFD）方法模擬鼓泡塔反應(yīng)器內(nèi)的多尺度流動(dòng)問(wèn)題，模擬結(jié)果的準(zhǔn)確性取決于模型對(duì)塔內(nèi)多相湍流描述是否合理，如相間作用力的設(shè)置、氣泡聚并和破碎的描述以及能量耗散機(jī)制等。Zhang等6-7考察了氣泡聚并、聚并效率、破裂所需的臨界能量、氣泡大小等因素對(duì)CFD-BPM耦合模型數(shù)值計(jì)算的影響；Yu等8采用CFD-PBM模擬來(lái)預(yù)測(cè)分散相液滴尺寸分布，并研究脈沖塔內(nèi)局部液-液流動(dòng)行為；Yang等9-10考察了壓力對(duì)氣泡破碎的影響，并在氣泡破碎模型中加入

3、了密度修正項(xiàng)；Gong等11提出了一種改進(jìn)的模擬湍流引起氣泡聚并的理論模型用來(lái)預(yù)測(cè)臨界聚并速度以及氣泡大小分布；王鐵鋒等12-14的研究認(rèn)為氣泡群曳力系數(shù)與單個(gè)氣泡及氣泡群的運(yùn)動(dòng)速度有關(guān)，考慮了能量和壓力約束條件，并通過(guò)湍流渦體動(dòng)能和氣泡表面能增量的大小來(lái)判斷氣泡破碎的原因；Prince等15的聚并效率模型考慮了液膜排水模型的機(jī)理；Luo等16的聚并效率模型考慮了能量守恒；Lehr等17的聚并效率模型考慮了臨界速度；Luo等18認(rèn)為，當(dāng)湍流渦旋的湍動(dòng)能大于氣泡表面能量增量時(shí)，氣泡發(fā)生破碎，這是一種能量約束的現(xiàn)象；Shi等19結(jié)合這些結(jié)果，通過(guò)引入等效直徑進(jìn)一步修正了氣泡破碎模型；Hinze20

4、研究發(fā)現(xiàn)破碎函數(shù)可以通過(guò)慣性力和表面張力相比來(lái)定義；Kolmogorov在Weber數(shù)的基礎(chǔ)上提出的破碎函數(shù)，認(rèn)為氣泡破碎行為和Weber數(shù)有關(guān)21。因此，上述研究推動(dòng)了氣泡破碎模型與聚并模型在氣液鼓泡塔數(shù)學(xué)模擬的發(fā)展。本文在課題組前期工作的基礎(chǔ)上22-26，在空氣-醋酸體系的鼓泡塔內(nèi)對(duì)曳力模型和聚并模型進(jìn)行修正，通過(guò)CFD-PBM耦合模型進(jìn)行了二維和三維數(shù)值模擬，研究了醋酸濃度對(duì)鼓泡塔內(nèi)流體力學(xué)參數(shù)的影響，并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，驗(yàn)證該模型的可行性。1 數(shù)學(xué)模型1.1 雙流體模型氣液系統(tǒng)采用歐拉-歐拉模型，控制方程如下：質(zhì)量守恒方程tii+iiui=0,i=g,l(1)動(dòng)量守恒方程tiiui+

5、iiuiui=-iP+ieffu+uT+Fi,j+iig(2)式中，i代表氣相g或液相l(xiāng)；、u、t分別表示相含率、密度、速度矢量和時(shí)間；P表示壓力；eff表示流體有效黏度。1.2 湍流方程選擇雷諾時(shí)均法（RANS）中的標(biāo)準(zhǔn)k-模型27，方程具體描述如下：k方程kt+uk=+tkk+Gk-(3)方程t+u=+t+kC1Gk-C2(4)t=lCk2(5)式中，C1=1.44,C2=1.92,C=0.09,k=1.0,= 1.3。1.3 相間作用力模型歐拉-歐拉模型需要相間作用力來(lái)封閉，相間作用力的選擇對(duì)于模擬結(jié)果的精度十分重要。本文模擬選擇了曳力、升力、壁面潤(rùn)滑力和湍流擴(kuò)散力，忽略了虛擬質(zhì)量力的作

6、用。具體表達(dá)式為：Fi,j=FD+FL+FT+FW(6)1.3.1 曳力本文基于Li等28的能量最小多尺度（EMMS）模型，在不同質(zhì)量分?jǐn)?shù)的醋酸-空氣體系中，通過(guò)表面張力函數(shù)f(l/0)修正曳力模型：CDdb=g,smalldsmallCD,smallUg,smallg,small2+g,largedlargeCD,largeUg.largeg,large2g,totalUg2(7)簡(jiǎn)化并優(yōu)化后的關(guān)系式為：CDdb=fl0159-2166Ug+12359Ug2(8)fl0=7.1l0-5l02-1(9)1.3.2 升力氣泡在液相剪切流中上升時(shí)，會(huì)產(chǎn)生與氣泡上升方向垂直的一種徑向力，稱為升力。徑

7、向升力對(duì)于氣含率的徑向分布至關(guān)重要，本文采用文獻(xiàn)29的升力，具體表達(dá)式為：FL=-CLglug-ulul(10)CL=-0.2CTDl2l(11)1.3.3 壁面潤(rùn)滑力壁面潤(rùn)滑力是由靠近壁面處氣液速度梯度引起的，是使氣泡遠(yuǎn)離鼓泡塔壁面的一種力。本文采用文獻(xiàn)30的壁面潤(rùn)滑力，具體表達(dá)式如下：FW=CWLlg(ug-ul)-(ug-ul)nr2nr(12)CWL=CWmax0,1CWD+1-yWCWCdbyWyWCWCdbm-1(13)式中，CWD、CWC為無(wú)量綱常數(shù)；m=1.52。CW為壁面潤(rùn)滑力系數(shù)，是關(guān)于Eo的函數(shù)，具體表達(dá)式如下：CW=0.47Eo1e-0.933Eo+0.1791Eo50

8、.00599Eo-0.01875Eo330.17933Eo(14)1.3.4 湍流擴(kuò)散力湍流擴(kuò)散力是由液相湍流旋渦引起的，此力使得徑向氣含率分布更均勻。本文采用文獻(xiàn)31的湍流擴(kuò)散力，具體表達(dá)式為：FT=CTDlklg(15)式中，kl為湍流強(qiáng)度；CTD為湍流擴(kuò)散力系數(shù)，取值范圍0.11，本文中取值為1。1.4 群體平衡模型群體平衡模型（PBM）主要是描述反應(yīng)器中顆?；驓馀荽笮》植嫉囊环N方法。在鼓泡塔中，氣泡主要發(fā)生破碎和聚并，因此該狀況下群體平衡模型表示如下：nV,tt時(shí)間項(xiàng)+nV,tubV,t對(duì)流項(xiàng)+vdVbdtnV,t傳質(zhì)和壓力變化的影響=120VnV-V,tnV,tcV-V,VdV-聚并

9、引起的源項(xiàng)0nV,tnV,tnV,VdV聚并引起的匯項(xiàng)+0V,VbVnV,tdV破碎引起的源項(xiàng)-bVnV,t破碎引起的匯項(xiàng)(16)式中，V為母氣泡體積；V為子氣泡體積；n(V,t)為體積為V的氣泡數(shù)密度函數(shù)；c(V-V,V) 為氣泡聚并速率；（V,V）為體積V的氣泡破裂成體積V的子氣泡分布函數(shù)。1.4.1 聚并模型氣泡聚并速率可表示為：cdi,dj=cdi,djPcdi,dj(17)氣泡間碰撞頻率為：cdi,dj=241/3di+dj2di2/3+dj2/31/2(18)本文基于文獻(xiàn)16的聚并模型，在空氣-醋酸體系中，通過(guò)引入表觀氣速和表面張力修正項(xiàng)對(duì)聚并系數(shù)進(jìn)行了修正，如表1所示。表1聚并系

10、數(shù)的修正Table 1Correction of the coalescence coefficientUg/(m/s)Concentration/%(mass)/0CeexpsimError/%0.047501.441.200.1200.113-6.190.071501.441.150.1520.1583.950.094501.441.100.1850.2022.650.047601.361.150.1390.130-6.920.071601.361.100.1740.1771.720.094601.361.050.2120.2214.250.047701.291.100.1430.141-

11、1.400.071701.291.050.1820.1872.750.094701.291.000.2200.2325.450.047801.201.050.1400.134-6.830.071801.201.000.1740.1824.600.094801.200.950.2070.2269.170.047901.111.000.1210.116-4.130.071901.110.950.1600.1653.130.094901.110.900.1950.2086.670.0471001.000.950.1170.108-7.700.0711001.000.900.1530.1551.310

12、.0941001.000.850.1930.1993.11新窗口打開|下載CSV修正的聚并系數(shù)Ce具體表達(dá)式如下：Ce=0.9Ug-0.0200.62(19)式中，0為25下標(biāo)準(zhǔn)濃度醋酸的表面張力;為實(shí)際醋酸溶液的表面張力。1.4.2 破碎模型本文主要采用Luo等18的氣泡破碎模型進(jìn)行數(shù)值模擬，該模型具體表達(dá)式為：bfv,d=0.92381-gd21/3min11+211/3exp-12cfl2/311/3d5/3d(20)式中，的取值為2.047；min為湍流渦和母氣泡最小尺寸比，即min=min/d，其中min為Kolmogorov最小湍流渦尺寸的11.4倍；cf為氣泡表面能增加量。2 實(shí)

13、驗(yàn)裝置和模擬設(shè)置2.1 實(shí)驗(yàn)部分實(shí)驗(yàn)在玻璃鼓泡塔內(nèi)進(jìn)行，實(shí)驗(yàn)裝置和測(cè)量方法詳見(jiàn)文獻(xiàn)32，塔內(nèi)直徑D=0.15 m，塔高H=2.2 m，采用差壓變送器、光纖探針、ERT技術(shù)手段測(cè)量實(shí)驗(yàn)數(shù)據(jù)。同樣塔底部用以儲(chǔ)存氣體的氣室以及多孔氣體分布板均采用不銹鋼材料制作完成，氣體分布板上共有2 mm19的開孔，開孔率約為0.338%。2.2 物理特性實(shí)驗(yàn)采用符合標(biāo)準(zhǔn)的工業(yè)冰醋酸（100%的醋酸），并根據(jù)其密度配制不同濃度的醋酸溶液，對(duì)所配溶液進(jìn)行密度、表面張力、黏度的測(cè)量，不同濃度醋酸溶液的物性如圖1所示。圖1圖1不同質(zhì)量分?jǐn)?shù)的醋酸物理性質(zhì)Fig.1Physical properties of acetic

14、acid with different mass fractions2.3 數(shù)值模擬條件的設(shè)置數(shù)值模擬使用ANSYS FLUENT軟件作為計(jì)算平臺(tái)，氣相為常溫常壓下的空氣，其作為離散相；液相為不同質(zhì)量分?jǐn)?shù)的醋酸溶液，其作為連續(xù)相。對(duì)于軟件程序，設(shè)置速度為鼓泡塔塔底入口邊界條件，壓力為塔頂出口邊界條件。鼓泡塔塔高為2.2 m，初始靜止液面高度為1.1 m，實(shí)驗(yàn)值在塔高為0.86 m處進(jìn)行采集。模擬的時(shí)間步長(zhǎng)固定為 0.002 s，并認(rèn)為在80 s內(nèi)達(dá)到了準(zhǔn)穩(wěn)態(tài)。3 空氣-醋酸體系二維數(shù)值模擬3.1 網(wǎng)格無(wú)關(guān)性驗(yàn)證圖2描述了表觀氣速為0.094 m/s、系統(tǒng)壓力為101325 Pa下，采用網(wǎng)格數(shù)為

15、2200、3080和4590個(gè)三種網(wǎng)格驗(yàn)證徑向氣含率分布的計(jì)算精度。從結(jié)果看出，徑向氣含率隨著網(wǎng)格數(shù)的增加基本不發(fā)生變化，因此，在綜合考慮模擬的計(jì)算精度和計(jì)算時(shí)間的情況下，本模擬中最終選擇網(wǎng)格數(shù)為3080個(gè)的Grid網(wǎng)格。圖2圖2不同網(wǎng)格質(zhì)量對(duì)徑向氣含率的影響Fig.2Influence of different mesh quality on radial gas holdup3.2 不同曳力模型對(duì)徑向氣含率的影響在80%醋酸的條件下，將S-N、Tomiyama單氣泡曳力模型以及本文修正的曳力模型等模擬的徑向氣含率結(jié)果進(jìn)行對(duì)比分析，如圖3所示。單氣泡直徑的設(shè)定采用光纖探針實(shí)驗(yàn)測(cè)得的平均值，在

16、0.071和0.094 m/s的表觀氣速下，測(cè)得的氣泡直徑分別為7.0和6.5 mm。從圖中可以看出，徑向氣含率隨表觀氣速的增加而增加，從塔中心到塔壁處以拋物線形式分布，越接近塔壁處，徑向氣含率越小。而且可以發(fā)現(xiàn)壁面處的徑向氣含率模擬結(jié)果明顯低于光纖探針測(cè)得的實(shí)驗(yàn)值，這是因?yàn)楫?dāng)鼓泡塔的塔徑較小時(shí)，具有較強(qiáng)的邊壁效應(yīng)，導(dǎo)致塔壁面處預(yù)測(cè)值偏低。圖3圖3不同曳力模型和表觀氣速下徑向氣含率對(duì)比Fig.3Comparison of radial gas holdup under different drag force models and superficial gas velocities通過(guò)對(duì)三種

17、曳力模型的模擬結(jié)果進(jìn)行對(duì)比，結(jié)果表明，修正后的曳力模型，直接以CD/db的形式計(jì)算曳力的大小，因此不會(huì)受到氣泡尺寸大小的影響，與光纖探針測(cè)得的實(shí)驗(yàn)值吻合較好；而S-N單氣泡曳力模型的預(yù)測(cè)結(jié)果偏低，Tomiyama單氣泡曳力模型的預(yù)測(cè)結(jié)果明顯偏高，與王鈺等33通過(guò)基于EMMS方法的鼓泡塔反應(yīng)器CFD模擬結(jié)果相一致。3.3 軸向高度對(duì)氣含率的影響圖4為醋酸濃度為100%，塔軸向高度為0.76 m和0.86 m處的徑向氣含率分布。從圖中可以看出，徑向氣含率和軸向液速呈現(xiàn)出塔中心大，邊壁減小的趨勢(shì)。在0.76 m和0.86 m處徑向氣含率分布基本重合，說(shuō)明在鼓泡塔的這兩個(gè)軸向高度處已經(jīng)為充分發(fā)展階段，

18、與文獻(xiàn)34中列管型鼓泡塔中流動(dòng)發(fā)展規(guī)律相一致。圖4圖4不同軸向高度下的徑向氣含率分布Fig.4Radial gas holdup and axial liquid velocity distribution at different column heights3.4 空氣-醋酸體系徑向氣含率圖5為空氣-醋酸體系下徑向氣含率分布。從圖中可以看出，在塔高0.86 m處，模型的模擬結(jié)果與光纖探針測(cè)得的實(shí)驗(yàn)值基本吻合。在相同濃度下，徑向氣含率隨表觀氣速的增大而增大。且徑向氣含率從塔中心到塔壁處呈現(xiàn)減小趨勢(shì)，在塔壁處徑向氣含率減小幅度較大。圖5圖5空氣-醋酸體系下徑向氣含率分布Fig.5Radial

19、gas holdup distribution with air-acetic acid system3.5 空氣-醋酸體系平均氣含率圖6為不同濃度和表觀氣速下醋酸平均氣含率分布。從圖中可以看出，實(shí)驗(yàn)值和模擬結(jié)果在10%的誤差范圍內(nèi)基本吻合。醋酸濃度在50%100%的范圍內(nèi)時(shí)，平均氣含率先增加后減小，濃度在70%80%之間時(shí)，平均氣含率存在最大值。醋酸濃度在50%80%范圍內(nèi)時(shí)，表面張力呈現(xiàn)減小的趨勢(shì)，表面張力減小使得塔內(nèi)氣泡穩(wěn)定性減弱，大氣泡破碎成小氣泡，氣含率增加。醋酸濃度在80%100%范圍內(nèi)時(shí)，表面張力減小，氣泡破碎的概率增大。但同時(shí)發(fā)現(xiàn)，在此醋酸濃度范圍內(nèi)，黏度在13 mPas的低黏

20、度范圍內(nèi)是逐漸減小的。Ruzicka等35研究認(rèn)為黏度在03 mPas時(shí)，氣含率隨黏度的增加而增加。因此，在此濃度范圍內(nèi)，黏度減小導(dǎo)致氣含率減小。黃娟等36通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn)，氣含率在醋酸濃度為60%80%的范圍內(nèi)存在最大值，也驗(yàn)證了本文的模擬結(jié)果。圖6圖6不同醋酸濃度和表觀氣速下平均氣含率分布Fig.6Distribution of average gas holdup under different acetic acid concentrations and superficial gas velocities3.6 空氣-醋酸體系軸向液速圖7為不同濃度和表觀氣速下軸向液速分布。由圖可以看

21、出，塔內(nèi)軸向液速沿塔中心向上運(yùn)動(dòng)，且軸向液速?gòu)乃行牡剿谔幹饾u減小，在塔壁處軸向液速為負(fù)值，說(shuō)明塔壁面處軸向液速沿壁面向下運(yùn)動(dòng)，而且塔中心軸向液速都在0.20.4 m/s范圍內(nèi)。在相同的醋酸濃度下，軸向液速隨著塔內(nèi)表觀氣速的增加而增加，這是因?yàn)楸碛^氣速越大，向上運(yùn)動(dòng)的氣體對(duì)塔中心的液體施加的軸向向上的作用力越大，因此軸向液速越大。在相同的表觀氣速下，軸向液速隨著醋酸濃度的增大基本不變。Yan等24研究認(rèn)為，在不同的表面張力下，塔內(nèi)的軸向液速基本不變。Krishna等37通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn)，液體黏度對(duì)塔內(nèi)軸向液速的影響作用可以忽略。圖7圖7空氣-醋酸體系下軸向液速分布Fig.7Axial liq

22、uid velocity distribution with air-acetic acid system3.7 空氣-醋酸體系氣泡直徑分布圖8為不同濃度和表觀氣速下醋酸的氣泡直徑分布。由圖8可看出，在相同的醋酸濃度下，隨著表觀氣速的增加，徑向氣泡直徑分布也隨著增大。這是因?yàn)殡S著塔內(nèi)表觀氣速增大，增加了塔內(nèi)氣泡的聚并，小氣泡數(shù)量減少，大氣泡數(shù)量增多，徑向氣泡直徑分布增大。在相同的表觀氣速下，徑向氣泡直徑分布受到液相醋酸黏度和表面張力的影響，黏度和表面張力的變化導(dǎo)致氣泡發(fā)生聚并和破碎。醋酸濃度在50%100%的范圍內(nèi)，隨著黏度增加，徑向氣泡直徑分布相應(yīng)增加，當(dāng)醋酸濃度達(dá)到80%時(shí)，黏度減小，徑向

23、氣泡直徑分布相應(yīng)減小；而在此醋酸濃度范圍內(nèi)，表面張力在一直減小，徑向氣泡直徑分布也相應(yīng)減小。因此，醋酸濃度為100%時(shí)，其黏度和表面張力都為最小值，此時(shí)氣泡直徑分布最小。圖8圖8不同醋酸濃度和表觀氣速下徑向氣泡直徑分布Fig.8Radial bubble diameter distribution under different acetic acid concentrations and superficial gas velocities3.8 空氣-醋酸體系氣泡數(shù)密度分布圖9為不同濃度和表觀氣速下氣泡數(shù)密度分布。由圖9可看出，氣泡在直徑為7 mm左右所占有的比例最大，這也與圖8光纖探針測(cè)

24、得的氣泡直徑值分布相符合。在相同表觀氣速下，醋酸濃度為50%100%范圍內(nèi)，氣泡數(shù)密度變化程度不大。在醋酸濃度為50%時(shí)，氣泡尺寸分布較窄，而醋酸濃度為70%和80%時(shí)，氣泡尺寸分布較寬，此時(shí)黏度較大。而在醋酸濃度為100%時(shí)，表面張力和黏度等物理性質(zhì)達(dá)到最小值，此時(shí)氣泡尺寸分布最窄。因?yàn)殡S著醋酸濃度的增加，液相醋酸黏度先增大后減小，表面張力一直減小。液相黏度增加，氣泡發(fā)生聚并概率增加，小氣泡聚并成大氣泡，氣泡尺寸分布變寬；表面張力減小，大氣泡變得不穩(wěn)定，破碎成小氣泡，塔內(nèi)小氣泡數(shù)量增多增加了氣泡在塔內(nèi)的停留時(shí)間，氣含率增加，氣泡分布范圍變窄。在相同的醋酸濃度下，表觀氣速的增加使得進(jìn)入塔內(nèi)的氣

25、體流量增加，氣泡數(shù)密度增加。圖9圖9氣泡數(shù)密度分布（z=0.86 m）Fig.9Bubble number density distribution (z=0.86 m)4 醋酸體系三維數(shù)值模擬4.1 三維網(wǎng)格劃分本部分為不同濃度醋酸三維數(shù)值模擬，采用三維非結(jié)構(gòu)性網(wǎng)格的劃分，三維網(wǎng)格的網(wǎng)格數(shù)為93856個(gè)。與二維數(shù)值模擬相比，三維數(shù)值模擬可以得到塔截面的氣含率分布狀況、塔內(nèi)瞬時(shí)氣含率分布以及塔內(nèi)流型的變化情況。因此，進(jìn)行三維數(shù)值模擬對(duì)了解塔內(nèi)氣液相的瞬時(shí)變化情況具有一定的指導(dǎo)意義。4.2 三維徑向截面氣含率分布通過(guò)ERT實(shí)驗(yàn)測(cè)得的云圖和三維數(shù)值模擬的云圖進(jìn)行了對(duì)比分析，圖10為50%80%醋酸

26、的平均氣含率云圖分布。從圖中可以看出，三維數(shù)值模擬的云圖和ERT測(cè)得的云圖變化規(guī)律基本一致。為了可視化不同濃度醋酸在不同表觀氣速下的變化規(guī)律，使用統(tǒng)一的圖例（00.40）進(jìn)行說(shuō)明。圖中中心處黃色區(qū)域代表較高的氣含率，而邊壁藍(lán)色區(qū)域代表較低的氣含率。隨著表觀氣速的增加，圖像中的黃色區(qū)域逐漸增加，藍(lán)色區(qū)域逐漸減小。說(shuō)明在相同醋酸濃度下，氣含率隨著表觀氣速的增加而增加。同時(shí)，從圖中可以看出，當(dāng)醋酸濃度為70時(shí)，塔中心氣泡量相對(duì)較大，與黃娟等36實(shí)驗(yàn)結(jié)果一致。圖10圖10不同表觀氣速下不同濃度醋酸橫截面氣含率云圖分布（z= 0.86 m）Fig.10Cross-section contours of

27、gas holdup distribution of different mass fractions of acetic acid under different superficial gas velocities（z= 0.86 m）4.3 三維柱體及軸向截面氣含率分布從三維柱體圖11可以看出，氣含率隨表觀氣速的增加而增加，且塔中心氣含率較高，塔壁面處氣含率較低。三維柱體圖可以更全面地看出塔中心和塔壁面氣含率沿塔軸向高度的變化趨勢(shì)。也可以看出氣泡在上升過(guò)程中的流動(dòng)形態(tài)，即以S形螺旋式上升。圖11圖1180%醋酸三維柱體氣含率分布Fig.11Three-dimensional cylind

28、er gas holdup distribution map of 80% acetic acid solution4.4 不同時(shí)間三維軸向截面氣含率分布在表觀氣速為0.094 m/s的情況下，研究了130 s時(shí)間軸向截面的氣含率變化，圖12為濃度80%的醋酸在不同時(shí)間三維軸向截面氣含率云圖分布。從圖中可以看出，隨著模擬時(shí)間的增加，初始液面從1.1 m上升到1 s時(shí)的1.2 m，最后穩(wěn)定在1.4 m。氣泡流動(dòng)從一開始沿塔高呈現(xiàn)出向上緩沖狀態(tài)，到5 s時(shí)垂直貫穿整個(gè)液面高度，且呈現(xiàn)出中心對(duì)稱分布，達(dá)到25 s后呈現(xiàn)為S形上升的趨勢(shì)。而且從整體上來(lái)看，軸向瞬時(shí)氣含率分布呈現(xiàn)出塔中心高、塔壁面處低的

29、分布情況，這也與塔內(nèi)實(shí)際的流動(dòng)情況相符。圖12圖1280%醋酸軸向截面氣含率云圖分布Fig.12Axial-section contours of gas holdup distribution of 80% acetic acid solution4.5 不同軸向高度處三維徑向截面氣含率分布圖13為濃度80%的醋酸三維軸向截面氣含率充分發(fā)展?fàn)顟B(tài)下的云圖分布。從圖中可以看出，在表觀氣速為0.094 m/s的條件下，不同軸向高度的氣含率呈現(xiàn)中心高、邊壁低的趨勢(shì)，而且較高的氣含率基本分布在塔中心，軸向高度的變化對(duì)徑向截面氣含率的大小分布影響不大。圖13圖13不同濃度醋酸不同軸向高度處徑向截面氣含率云圖分布Fig