射流鼓泡反應(yīng)器的混合特性

郭天琪，黃正梁，王靖岱，蔣斌波，陽永榮

（浙江大學(xué)化學(xué)工程與生物工程學(xué)院，化學(xué)工程聯(lián)合國家重點實驗室，浙江 杭州 310027）

射流鼓泡反應(yīng)器的混合特性

郭天琪，黃正梁，王靖岱，蔣斌波，陽永榮

（浙江大學(xué)化學(xué)工程與生物工程學(xué)院，化學(xué)工程聯(lián)合國家重點實驗室，浙江 杭州 310027）

射流鼓泡反應(yīng)器以液體射流代替攪拌實現(xiàn)液相混合，具有結(jié)構(gòu)簡單、制造及維護(hù)費用低等諸多優(yōu)點，研究其混合特性對于反應(yīng)器的設(shè)計、優(yōu)化及放大具有重要意義。以空氣-水作為模擬介質(zhì)，采用KCl電解質(zhì)溶液為示蹤劑考察了表觀氣速和射流Reynolds數(shù)的大小對液相宏觀混合時間的影響，并從能量輸入的角度對射流鼓泡反應(yīng)器的混合機(jī)制進(jìn)行分析。研究發(fā)現(xiàn)，在實驗條件下（表觀氣速變化范圍為0.0006～0.0343 m·s-1，射流Reynolds數(shù)的變化范圍為1.75×104～7.00×104），鼓泡的加入使得均相射流反應(yīng)器內(nèi)的液相混合得到改善；隨著表觀氣速增大，液相宏觀混合時間先縮短后延長；當(dāng)氣體輸入功率或液體輸入功率不變時，混合時間隨總輸入功率的增大而縮短。通過對多組實驗數(shù)據(jù)的回歸分析，提出了液相宏觀混合時間與液體輸入功率和氣體輸入功率的經(jīng)驗關(guān)聯(lián)式，計算值與實際值吻合較好。最后基于提出的關(guān)聯(lián)式，發(fā)現(xiàn)當(dāng)總輸入功率一定時，混合時間隨氣體輸入功率的增加先縮短后延長，臨界轉(zhuǎn)變點在氣體輸入功率為總功率的61%處，此時氣液兩相協(xié)同作用最強(qiáng)。

多相反應(yīng)器；混合；氣泡；射流；輸入功率

引　言

射流鼓泡反應(yīng)器是一類帶有液體噴嘴以及氣體發(fā)生鼓泡部件的反應(yīng)器，利用液體噴射代替?zhèn)鹘y(tǒng)攪拌槳，一方面可以降低動設(shè)備的機(jī)械維護(hù)費用，另一方面在液體循環(huán)量較大時，有效地避免了攪拌槳的振動問題。在射流鼓泡反應(yīng)器內(nèi)，液體在泵的輸送下經(jīng)過噴嘴高速噴射進(jìn)入反應(yīng)器內(nèi)，利用高速液體射流的剪切作用來破碎氣泡，實現(xiàn)反應(yīng)器內(nèi)氣液兩相的高效分散混合。眾所周知，反應(yīng)器內(nèi)物料的混合狀態(tài)直接影響傳熱、傳質(zhì)和化學(xué)反應(yīng)。因此，深入研究射流鼓泡反應(yīng)器內(nèi)的混合特性對于反應(yīng)器的設(shè)計、優(yōu)化及放大均具有重要意義。

射流混合的概念廣泛應(yīng)用于均相體系。目前文獻(xiàn)中關(guān)于射流混合的研究可以分為兩個方面，一方面通過實驗研究，考察操作條件以及幾何尺寸，如射流速度、噴嘴直徑、射流角度及位置、液體黏度、反應(yīng)釜直徑及釜底形狀等對混合時間的影響[1-5]；另一方面通過計算流體力學(xué)（CFD）技術(shù)進(jìn)行模擬研究，以期對反應(yīng)器內(nèi)的流體力學(xué)行為及對混合的作用規(guī)律進(jìn)行細(xì)致的刻畫[6-10]。

射流現(xiàn)象用于氣液混合反應(yīng)主要有以下兩種反應(yīng)器形式，一種是噴射式氣液反應(yīng)器，另一種是射流鼓泡反應(yīng)器。前者多用來代替鼓泡塔[11-12]，而后者則多用來代替攪拌反應(yīng)器。其中噴射式氣液反應(yīng)器又根據(jù)射流噴嘴在反應(yīng)器上下位置的不同，分為下噴式自吸反應(yīng)器和上噴式環(huán)流反應(yīng)器。在傳統(tǒng)的鼓泡塔內(nèi)，氣泡尺寸決定氣液接觸面積，影響氣液傳質(zhì)效果，而其受向下流動的液體和向上流動的氣體影響，這使得鼓泡塔內(nèi)的氣液兩相行為在一定程度上受到限制[13]。因此提出采用噴射式氣液反應(yīng)器來代替鼓泡塔以期得到更好的氣液傳質(zhì)效果[11-12]。在噴射式氣液反應(yīng)器內(nèi)，液體射流不僅可以對氣泡產(chǎn)生剪切作用，形成較小的氣泡，同時還會與小氣泡接觸進(jìn)行兩相反應(yīng)。反應(yīng)器內(nèi)氣、液兩相一同從噴嘴處進(jìn)料，當(dāng)具有一定壓力的液體經(jīng)反應(yīng)器上部的噴嘴噴射時，會產(chǎn)生很高的液體流速，噴嘴周圍由此形成低壓區(qū)，反應(yīng)原料被吸入吸氣室，然后在吸氣室下方的混合室內(nèi)充分混合。噴嘴主要用于將氣體卷吸入液體中，形成小氣泡。與之不同的是，在射流鼓泡反應(yīng)器中，氣、液兩相分別進(jìn)料，噴嘴既用于攪拌液體，又可以使液體形成循環(huán)，分散氣體。此外，噴射式氣液反應(yīng)器中下噴式的噴嘴在液面上方，上噴式的噴嘴在反應(yīng)器底部全部淹沒在液相中，二者均是氣液順流；而射流鼓泡反應(yīng)器中噴嘴伸入到液面下方，氣液兩相形成逆流。

與噴射式氣液反應(yīng)器相比，射流鼓泡反應(yīng)器的研究相對較少。Amiri等[14]首次對射流鼓泡反應(yīng)器的混合性能進(jìn)行了研究，得到射流可以有效改善液相混合以及越接近氣液兩相相遇位置混合時間越短的結(jié)論，并根據(jù)不同噴嘴Reynolds數(shù)下，隨表觀氣速變化得到的混合時間的4種不同變化趨勢，總結(jié)出4類氣液兩相流動結(jié)構(gòu)。但文獻(xiàn)只依據(jù)流動結(jié)構(gòu)圖對實驗現(xiàn)象進(jìn)行解釋，并未對結(jié)果進(jìn)行深入分析，挖掘產(chǎn)生不同流動結(jié)構(gòu)的原因。

在氣液兩相體系中，液相混合是影響反應(yīng)器性能的主要因素之一[15]。良好的液相混合可以防止組分在氣液相界面間的過度積累或缺失，進(jìn)而增強(qiáng)質(zhì)量傳遞的驅(qū)動力，加速有效組分的轉(zhuǎn)移。因此，可以通過強(qiáng)化液相混合和增加小氣泡的數(shù)量來增大氣液間的相界面積，進(jìn)而從根本上強(qiáng)化傳質(zhì)。此外，在氣液兩相體系中，傳質(zhì)效果受氣體在液膜中溶解度限制，而液相混合則可減小液膜厚度，進(jìn)一步改善傳質(zhì)效果[16]。

本文針對射流鼓泡反應(yīng)器的液相混合特性這一重要參數(shù)開展實驗研究，采用電解質(zhì)示蹤法測量液相宏觀混合時間，系統(tǒng)考察了表觀氣速對液相宏觀混合時間的作用規(guī)律，并從能量輸入的角度對宏觀混合特性進(jìn)行分析，提出液相宏觀混合時間的經(jīng)驗關(guān)聯(lián)式，最終希望為工業(yè)射流鼓泡反應(yīng)器的開發(fā)提供理論指導(dǎo)。

1　實驗部分

1.1實驗裝置及流程

圖1為冷模實驗系統(tǒng)流程，包括反應(yīng)器、氣路系統(tǒng)、水路系統(tǒng)、檢測儀器以及計算機(jī)數(shù)據(jù)采集系統(tǒng)5個部分。反應(yīng)器由有機(jī)玻璃制成，由筒體、半球形下封頭、液體注入管、氣體分布環(huán)、擋板等部分組成。反應(yīng)器筒體和半球形封頭的內(nèi)徑為 380 mm，反應(yīng)器筒體高度為1000 mm，反應(yīng)器內(nèi)裝有氣體分布環(huán)、液體噴嘴以及擋板。其中，氣體分布環(huán)上方布有均勻分布的出氣孔，下方布有少量排液口，液體噴嘴為縮徑式，具體結(jié)構(gòu)如圖2所示，氣體分布環(huán)距反應(yīng)器底部垂直距離為0.15 m，液體噴嘴出口距氣體分布環(huán)垂直距離為 0.4 m。反應(yīng)器設(shè)有進(jìn)氣口（A）的公稱直徑為31 mm，排液口（B）的公稱直徑為63 mm，取樣口（C）四周對稱分布，公稱直徑為30 mm。

圖1　冷模實驗流程Fig.1　Schematic diagram of experimental setup

冷模實驗裝置的結(jié)構(gòu)參數(shù)為自行設(shè)計確定，具體尺寸如表1所示。其中，di為筒體和半球形封頭內(nèi)徑，h為筒體高度，da為噴嘴縮徑前內(nèi)徑，db為噴嘴縮徑后內(nèi)徑，α為收縮角度，l為收縮部分長度，hr為縮徑后圓柱段長度，dg為氣體分布環(huán)直徑，ds為分布環(huán)內(nèi)徑，dk為分布環(huán)開孔直徑，n為開孔數(shù)。

表1　實驗裝置結(jié)構(gòu)參數(shù)Table 1　Parameters of experimental apparatus

冷模實驗在常溫、常壓下進(jìn)行，以空氣和水作液體從反應(yīng)器底部排出，在循環(huán)泵的泵送作用下經(jīng)流量計計量后通過噴嘴噴射進(jìn)入反應(yīng)器，構(gòu)成液體循環(huán)。實驗過程中，通入反應(yīng)器內(nèi)的氣體流量的變化范圍為0.25～14 m3·h-1，以反應(yīng)器截面積計算，表觀氣速的變化范圍為0.0006～0.0343 m·s-1，該操作區(qū)間可包含氣泡的鼓泡態(tài)和聚并態(tài)，液體射流Reynolds數(shù)的變化范圍為1.75×104～7.00×104。射流Reynolds數(shù)Rej的計算如式（1）所示。

式中，dj為射流出口直徑，uj為射流出口速度，ρ為射流液體密度，η為射流液體黏度。

圖2　液體噴嘴及氣體分布環(huán)結(jié)構(gòu)Fig.2　Structure diagram of liquid jet and gas distribution ring

1.2實驗測試方法及數(shù)據(jù)處理

采用電解質(zhì)示蹤法對液相宏觀混合時間進(jìn)行檢測。實驗開始前，向反應(yīng)器內(nèi)注入自來水約64 L，液面高度約 0.625 m，打開離心泵及風(fēng)機(jī)，調(diào)整閥門開度，使得氣液流量分別達(dá)到規(guī)定值。測量過程中，將一定體積(50 ml)的飽和KCl溶液通過注入裝置，經(jīng)由反應(yīng)器液體噴嘴脈沖加入。為消除氣泡對檢測造成的影響，在反應(yīng)器下部檢測點安裝特殊設(shè)計的引流裝置[17]，將塔內(nèi)液體引出，分離氣液兩相。采用位于筒體側(cè)壁底部的電導(dǎo)率儀（上海雷磁公司，DDSJ-308A）測量液體電導(dǎo)率，并將檢測結(jié)果輸入計算機(jī)內(nèi)保存。每個數(shù)據(jù)點實驗次數(shù)不少于5次，混合時間以電導(dǎo)率值的均勻度達(dá)到99%所用的時間為標(biāo)準(zhǔn)[14]。

圖3為檢測過程中液體電導(dǎo)率隨時間變化的典型結(jié)果。為了便于比較及后續(xù)的分析，采用如式（2）所示的方法對電導(dǎo)率儀的測量結(jié)果進(jìn)行量綱1化處理。

圖3　液體電導(dǎo)率隨時間的變化Fig.3　Variation of conductivity with time in mixing time

式中，κ（t）為t時刻的電導(dǎo)率，κi為實驗初始時的電導(dǎo)率，κ∞為實驗結(jié)束時達(dá)到穩(wěn)定的電導(dǎo)率。采用式（2）對圖3中的數(shù)據(jù)進(jìn)行處理，結(jié)果如圖4所示。當(dāng)示蹤劑濃度達(dá)到平衡濃度的 99%～101%之間時，該時刻與示蹤劑注入的時間間隔即為混合時間tM。

圖4　量綱1處理后電導(dǎo)率cj隨時間的變化Fig.4　Variation of dimensionless cjwith time in mixing time

2　實驗結(jié)果與討論

2.1鼓泡的影響

為探究氣相鼓泡對均相射流混合的影響，在射流反應(yīng)器內(nèi)通入氣體，對比了表觀氣速ug分別為0 和0.0061 m·s-1時，液相宏觀混合時間隨液體射流Reynolds數(shù)的變化情況，實驗結(jié)果如圖5所示。從圖中的結(jié)果可以看出，有鼓泡時反應(yīng)器內(nèi)的液相宏觀混合時間小于無鼓泡時的液相宏觀混合時間；當(dāng)射流Reynolds數(shù)為7×104時，兩條曲線之間的差距最小，此時有鼓泡時的液相宏觀混合時間是無鼓泡時的70%；而當(dāng)射流Reynolds數(shù)為1.75×104時，有鼓泡時的液相宏觀混合時間只有無鼓泡時的60%。由此可見，鼓泡對均相體系內(nèi)的混合情況有明顯的改善作用，這主要是由于氣泡引起的反應(yīng)器內(nèi)液體的局部湍動強(qiáng)化了液相混合。進(jìn)一步觀察可以發(fā)現(xiàn)，有無鼓泡時液相宏觀混合時間的差值隨著液體射流Reynolds數(shù)的增大而減小，分析其原因是由于射流Reynolds數(shù)較大時，反應(yīng)器內(nèi)液體主要受射流循環(huán)混合作用，表觀氣速一定的氣泡所引起的局部液相湍動對液相混合的強(qiáng)化程度相對降低。

圖5　有無鼓泡時液相宏觀混合時間隨射流Reynolds數(shù)的變化Fig.5　Effect of bubble injection（ug=0.0061 m·s-1）on mixing time with various jet Reynolds numbers

圖6　射流反應(yīng)器內(nèi)實際流動情況Fig.6　Real flow conditions in jet reactor

結(jié)合圖6實驗過程中所拍攝的照片及射流理論等對實驗現(xiàn)象進(jìn)一步分析，未加入鼓泡時，射流以一定的速度離開噴嘴出口后，與周圍靜止的流體形成速度不連續(xù)的間斷面，而速度間斷面是不穩(wěn)定的，射流引起的摻混相應(yīng)產(chǎn)生了對射流的阻力，使射流流速降低，能量減弱[18]。通入氣體后，環(huán)狀分布管鼓出的氣泡在上升過程中，受到反應(yīng)器中心豎直向下的液體射流的剪切分散作用，較好地分布在反應(yīng)器內(nèi)，引起周圍液體的局部湍動，在氣泡的帶動下向上流動，到達(dá)液面附近后在反應(yīng)器壁面處形成向下流動的液體，構(gòu)成液相流動循環(huán)。相比未加入鼓泡時，氣泡引起的液體湍動使得反應(yīng)器壁面處及液面附近的液體的流動速度增加，死區(qū)及低速區(qū)大大減少。綜上可知，加入鼓泡后，反應(yīng)器內(nèi)氣泡帶動液體局部湍動，形成液面附近的流動循環(huán)，液相混合在很大程度上得到了改善。

2.2表觀氣速的影響

圖7為不同射流Reynolds數(shù)下表觀氣速與液相宏觀混合時間的關(guān)系，其中橫坐標(biāo)為表觀氣速，縱坐標(biāo)是液相宏觀混合時間，不同曲線代表不同的射流Reynolds數(shù)。從圖中可以看出，當(dāng)射流Reynolds數(shù)一定時，液相宏觀混合時間隨表觀氣速的增大先縮短后延長，在不同的射流Reynolds數(shù)下均表現(xiàn)出類似的規(guī)律。當(dāng)射流Reynolds數(shù)較?。≧ej=1.75×104）時，隨著表觀氣速進(jìn)一步增大，液相混合時間的增幅較大。進(jìn)一步分析，當(dāng)射流Reynolds數(shù)較?。≧ej=1.75×104）時，極小值點對應(yīng)的表觀氣速約為0.007 m·s-1，而當(dāng)射流Reynolds數(shù)較大時，極小值點對應(yīng)的表觀氣速約為 0.017 m·s-1。也就是說，在不同射流Reynolds數(shù)下，液相宏觀混合時間的臨界轉(zhuǎn)變氣速是不同的。在低射流Reynolds數(shù)下，反應(yīng)器內(nèi)為氣相鼓泡控制，因此出現(xiàn)了與鼓泡塔（射流Reynolds數(shù)為0）中類似的研究結(jié)果[19-20]。文獻(xiàn)表明，鼓泡塔在低表觀氣速下(ug<0.01 m·s-1)，液相宏觀混合時間隨表觀氣速增加急劇減少；在較高的表觀氣速下(ug>0.01 m·s-1)，氣速影響趨緩，甚至出現(xiàn)了隨氣速增加混合時間延長的結(jié)果。但在高射流Reynolds數(shù)下，表觀氣速增大所形成的并聚的大氣泡可以被能量較高的射流破碎和剪切成為小氣泡，而小氣泡有利于液相湍動及混合，因此液相宏觀混合時間的臨界轉(zhuǎn)變氣速增大。下文將從能量輸入的角度對此做更深入的分析。

對實驗結(jié)果進(jìn)一步分析，當(dāng)表觀氣速較小時，射流可以到達(dá)氣體分布環(huán)上方，對反應(yīng)器內(nèi)的氣泡起到向器壁推動和分散作用，可以減少單一射流混合時反應(yīng)器內(nèi)的死區(qū)。同時，由于射流形成的液體循環(huán)，使得氣泡引起的液體局部湍動更加劇烈，射流和鼓泡相互協(xié)同影響，形成了小規(guī)模的動量和能量交換。隨著表觀氣速的增大，氣泡增多，氣體引起的湍動作用得以加劇，同時射流和氣泡相互碰撞引起的液體的湍動和攪動更加劇烈，因此混合時間減少。當(dāng)表觀氣速進(jìn)一步增加的時候，氣泡開始發(fā)生并聚，射流可破碎大氣泡，但是引起液體湍動的小氣泡未較之前增多，因此混合時間基本保持不變。當(dāng)表觀氣速再增加時，如果射流Reynolds數(shù)相對不足以破碎并聚的大氣泡，而大氣泡引起的液相循環(huán)作用又強(qiáng)于局部的湍流脈動，則會出現(xiàn)混合時間稍有增長，混合效果變差的情況。由此可知，表觀氣速不同，氣泡的數(shù)量、大小及狀態(tài)均不同，進(jìn)而影響反應(yīng)器內(nèi)液體流動狀態(tài)，對液相混合構(gòu)成影響。

圖7　不同射流Reynolds數(shù)下表觀氣速對液相宏觀混合時間的影響Fig.7　Effect of aeration rate on mixing time with different jet velocities

2.3能量輸入分析

在傳統(tǒng)的鼓泡反應(yīng)器中，反應(yīng)器內(nèi)液相混合主要為湍動混合機(jī)制控制，氣泡上升卷吸液體促進(jìn)了液相的混合。而在射流鼓泡反應(yīng)器中，液相混合則存在湍動混合和循環(huán)混合兩種控制機(jī)制，一方面射流會破碎氣泡、卷吸液體，影響氣泡和釜內(nèi)液體耦合控制的湍動混合強(qiáng)度；另一方面射流同時促使反應(yīng)器內(nèi)形成較規(guī)則的液體循環(huán)，影響射流和釜內(nèi)液體耦合控制的循環(huán)混合強(qiáng)度。兩種混合機(jī)制共同決定了射流鼓泡反應(yīng)器中的混合狀態(tài)。本文將從能量輸入的角度對射流鼓泡反應(yīng)器的混合機(jī)制進(jìn)行定量分析。

射流鼓泡反應(yīng)器中，無論是從氣體分布器進(jìn)入反應(yīng)器的氣體，還是從液體噴嘴進(jìn)入反應(yīng)器的液體，其實質(zhì)均是向反應(yīng)器內(nèi)輸入一定的能量。反應(yīng)器能量的輸入是控制反應(yīng)器內(nèi)液相混合、氣液傳質(zhì)等的重要參數(shù)[21]。在射流鼓泡反應(yīng)器中，由于氣體和液體兩相逆流進(jìn)入反應(yīng)器，其流動方向及對釜內(nèi)液體混合的控制機(jī)制均不同，氣泡和射流兩種作用機(jī)制并非單一疊加關(guān)系。而這兩種不同輸入能量間的相互協(xié)同作用會對反應(yīng)器的宏觀混合產(chǎn)生重要的影響。

氣液射流反應(yīng)器中氣體和液體的能量輸入可分別用式（3）和式（4）進(jìn)行計算[21]。

式中，PG和PL分別為氣體輸入功率和液體噴射輸入功率，W；DN為噴嘴直徑，m；VG和VLN分別為氣體和液體的體積流量，m3·s-1；ρG和ρL分別為氣相和液相的密度，kg·m-3；M為空氣的相對分子質(zhì)量；u0為氣體的表觀氣速，m·s-1；β為最大靜壓頭與反應(yīng)器頂部壓力（pT）之比，其大小可根據(jù)式（5）計算[21]

式中，L為泡沫液位，m；εG為氣含率；g為重力加速度。εG可通過測量靜液層高h(yuǎn)和通氣時泡沫液位L算出

圖8為不同氣體輸入功率下液相宏觀混合時間與液體輸入功率的關(guān)系。其中橫坐標(biāo)為液體輸入功率，縱坐標(biāo)為混合時間，不同曲線代表不同的氣體輸入功率，實驗過程中通過改變氣體和液體的流速來改變其輸入功率。從圖中可以發(fā)現(xiàn)，在一定的氣體輸入功率下，液體輸入的功率越大，反應(yīng)器的液相宏觀混合時間越短。這是因為在固定了一種能量輸入的情況下，增加另一種能量輸入的大小，反應(yīng)器內(nèi)的總能量輸入增大。如當(dāng)反應(yīng)器內(nèi)的液體能量輸入增大時，射流帶動的液體循環(huán)混合更加劇烈，反應(yīng)器內(nèi)的混合狀況得以改善，相應(yīng)地，液相宏觀混合時間縮短。即在此情況下，反應(yīng)器的能量輸入越高，反應(yīng)器內(nèi)的混合情況越好，液相宏觀混合時間越短。但需要注意的是，此結(jié)論適用于固定一種能量輸入的大小，僅考慮氣體或液體單因素對宏觀混合時間的影響。當(dāng)二者的輸入能量均發(fā)生變化時，情況會更加復(fù)雜。

圖8　不同氣體功率下液體輸入功率對液相宏觀混合時間的影響Fig.8　Effect of liquid power input on mixing time with different gas power input

圖9為實驗過程中混合時間隨反應(yīng)器總輸入功率的變化。圖中的點代表相應(yīng)總輸入功率下所對應(yīng)的混合時間。從圖中可以看出一個基本的變化規(guī)律，即總輸入功率越大，宏觀混合時間越短。

圖9　總輸入功率對液相宏觀混合時間的影響Fig.9　Effect of total power input on mixing time

2.4液相宏觀混合時間經(jīng)驗關(guān)聯(lián)式

從前文的研究結(jié)果可以發(fā)現(xiàn)，通過改變氣速（表觀氣速范圍為0.0006～0.0343 m·s-1）或射流Reynolds數(shù)（射流 Reynolds數(shù)范圍為 1.75×104～7.0×104）可以改變反應(yīng)器的能量輸入，進(jìn)而在一定程度上影響反應(yīng)器的液相宏觀混合時間。在分析了氣體及液體能量輸入對液相宏觀混合時間的影響規(guī)律的前提下，本文提出了液相宏觀混合時間關(guān)聯(lián)式

擬合得到a=85.87, b= -0.12, c= -0.19。實驗值和計算值的比較如圖10所示，實驗值和計算值的相對誤差基本控制在10%以內(nèi)，模型的準(zhǔn)確性較高。

為了定量研究氣體和液體輸入功率對液相宏觀混合時間的影響，即氣泡和射流的協(xié)同作用，在總輸入功率P一定的前提下對式（7）求導(dǎo)

令式（8）為零，則 PG=[0.19/(0.19+0.12)]P= 0.61P。當(dāng)PG＜0.61P時，式（8）小于0，當(dāng) PG>0.61P時，式（8）大于0。也就是說，在總輸入功率一定的前提下，氣體輸入功率不足總輸入功率的61%時，氣體輸入功率越大則混合時間越短，而當(dāng)其大小超過總輸入功率的61%時，氣體輸入功率越大混合時間也越長，PG=0.61P為混合時間的極小值點。

由前文的分析可知，射流鼓泡反應(yīng)器內(nèi)存在兩種型式的混合，即由氣體鼓泡引起的湍動混合以及由液體射流形成的循環(huán)混合。湍動混合強(qiáng)度越大，混亂度越高；循環(huán)混合強(qiáng)度越大，則液體湍動趨于規(guī)則有序，混亂度降低，混合效果較湍動混合差。因此，在氣液總輸入功率一定的條件下，當(dāng)氣體輸入功率較小時，增加氣體輸入功率可以改善反應(yīng)器內(nèi)的湍動混合程度，縮短液相宏觀混合時間；但當(dāng)氣體輸入功率進(jìn)一步增大時，勢必會使得液體輸入功率減小，此時液體射流能量過低，無法有效破碎由于較高的氣體輸入功率所形成的并聚的大氣泡，反應(yīng)器內(nèi)混合效果變差。

當(dāng)氣體輸入功率占總輸入功率61%時液相宏觀混合時間最短，說明反應(yīng)器內(nèi)氣體鼓泡和液體射流相互協(xié)同作用最強(qiáng)。

圖10　液相宏觀混合時間經(jīng)驗公式計算值與實驗值的比較Fig.10　Comparison of liquid mixing time between calculated data and experimental data

3　結(jié)　論

（1）應(yīng)用電解質(zhì)示蹤法測定了射流鼓泡反應(yīng)器內(nèi)的液相宏觀混合時間。實驗結(jié)果表明，鼓泡的加入很大程度上改善了均相射流反應(yīng)器內(nèi)液體的混合情況；液相宏觀混合時間隨著表觀氣速的增大呈先縮短后延長的趨勢，進(jìn)一步增大表觀氣速，當(dāng)射流Reynolds數(shù)過小不足以破碎大氣泡時，混合時間稍有增長。

（2）從能量輸入的角度探討鼓泡注入的氣體輸入功率以及射流注入的液體輸入功率與液相宏觀混合時間的關(guān)系，分析發(fā)現(xiàn)當(dāng)氣體輸入功率或液體輸入功率不變時，混合時間隨總輸入功率的增大而縮短；當(dāng)總輸入功率一定時，混合時間隨氣體輸入功率的增加先縮短后延長，臨界轉(zhuǎn)變點在氣體輸入功率為總功率的61%處，此時氣體鼓泡和液體射流協(xié)同作用最強(qiáng)。

（3）在實驗條件下（以反應(yīng)器截面積計，表觀氣速的變化范圍為 0.0006～0.0343 m·s-1，射流Reynolds數(shù)的變化范圍為1.75×104～7.0×104），回歸得到射流鼓泡反應(yīng)器內(nèi)液相宏觀混合時間的經(jīng)驗關(guān)聯(lián)式為：

符號說明

L——泡沫液位，m

M——空氣的相對分子質(zhì)量

PG——氣體流動輸入功率，W

PL——液體流動輸入功率，W

pT——反應(yīng)器頂部壓力，Pa

tM——混合時間，s

ug——表觀氣速，m·s-1

uj——射流出口液速，m·s-1

VG——氣體體積流速，m3·s-1

VLN——液體體積流速，m3·s-1

εG——氣含率

η ——黏度，Pa·s

κ——電導(dǎo)率，mS·cm-1

ρG——氣相的密度，kg·m-3

ρL——液相的密度，kg·m-3

References

[1] Meng Huibo(孟輝波), Wang Yanfen(王艷芬), Yu Yanfang(禹言芳), Wang Wei(王偉), Wang Feng(王豐), Wu Jianhua(吳劍華). Research progress of mixing time in jet mixing equipment [J]. Chemical Industry and Engineering Process(化工進(jìn)展), 2012, 31 (12): 2615-2625.

[2] Fossett H, Prosser L E. The application of free jets to the mixing of fluids in bulk [J]. Proceedings of the Institution of Mechanical Engineers, 1949, 160 (1): 224-232.

[3] Fox E A , Gex V E. Single-phase blending of liquids[J]. AIChE Journal, 1956, 2 (4): 539-544.

[4] Lane A G C, Rice P. Comparative assessment of the performance of the three designs for liquid jet mixing [J]. Industrial & Engineering Chemistry Process Design and Development, 1982, 21 (4): 650-653.

[5] Simon M, Fonade C. Experimental study of mixing performances using steady and unsteady jets [J]. The Canadian Journal of Chemical Engineering, 1993, 71 (4): 507-513.

[6] Brooker L. Mixing with the jet set [J]. Chemical Engineer-London,1993, (550): 16.

[7] Ranade V V. Towards better mixing protocols by designing spatially periodic flows: the case of a jet mixer [J]. Chemical Engineering Science, 1996, 51 (11): 2637-2642.

[8] Giorges A T, Forney L J, Wang X. Numerical study of multi-jet mixing [J]. Chemical Engineering Research and Design, 2001, 79 (5): 515-522.

[9] Jayanti S. Hydrodynamics of jet mixing in vessels [J]. Chemical Engineering Science, 2001, 56 (1): 193-210.

[10] Rahimi M, Parvareh A. Experimental and CFD investigation on mixing by a jet in a semi-industrial stirred tank [J]. Chemical Engineering Journal, 2005, 115 (1/2): 85-92.

[11] Evans G M. A study of a plunging jet bubble column[D]. New Castle: University of Newcastle, 1990.

[12] Blenke H. Biochemical loop reactors [J]. Biotechnology, 1985, 2: 465-517.

[13] Evans G M, Jameson G J. Hydrodynamics of a plunging liquid jet bubble-column [J]. Chemical Engineering Research and Design, 1995, 73(6): 679-684.

[14] Amiri T Y, Moghaddas J S, Moghaddas Y. A jet mixing study in two phase gas-liquid systems [J]. Chemical Engineering Research and Design, 2011, 89(3A): 352-366.

[15] Paul E L, Atiemo-Obeng V, Kresta S M. Handbook of Industrial Mixing: Science and Practice[M]. New York: John Wiley & Sons, 2004.

[16] Zlokarnik M. Stirring[M]. Wiley Online Library, 1988.

[17] Ma Yuelong(馬躍龍), Huang Juan(黃娟), Shen Chunyin(沈春銀), Dai Gance(戴干策). Liquid phase mixing time in bubble columns [J]. Journal of East China University of Science and Technology(華東理工大學(xué)學(xué)報), 2010, 36 (2): 165-172.

[18] Tian Yanli(田艷麗).Numerical simulation on flow field of rotory jet mixing system and nozzle structure optimization[D]. Hangzhou: Zhejiang University, 2008.

[19] Miron A S, Garcia M C C, Camacho F G, Grima E M, Chisti Y. Mixing in bubble column and airlift reactors [J]. Chemical Engineering Research and Design, 2004, 82 (10): 1367-1374.

[20] Haque M W, Nigam K, Joshi J B. Optimum gas sparger design for bubble columns with a low height-to-diameter ratio [J]. Chemical Engineering Journal, 1986, 33(2): 63-69.

[21] Zhang Yuanxing(張元興), Xu Xueshu(許學(xué)書). Biochemical Reactor Engineering(生物反應(yīng)器工程)[M]. Shanghai: East China University of Science and Technology Press, 2001.

Mixing characteristics in jet bubbling reactor

GUO Tianqi, HUANG Zhengliang, WANG Jingdai, JIANG Binbo, YANG Yongrong
(State Key Laboratory of Chemical Engineering, College of Chemical and Biological Engineering, Hangzhou 310027, Zhejiang, China)

The jet bubbling reactor uses liquid jet to achieve the liquid mixing instead of mechanical stirring, which brings several advantages such as simple structure and low cost of maintenance and manufacturing. The research of its mixing characteristics plays a significant role in the design, optimization and scaling up of the reactor. Based on the air-water system, the electrolyte tracer (KCl solution) method was applied to investigate the influences of gas velocity and jet Reynolds number on the liquid mixing time with the cold model experimental apparatus. The mixing mechanism in jet bubbling reactor had also been analyzed from the perspective of power input. The results showed that within the experimental range (ugfrom 0.0006 to 0.0343 m·s-1, Rejfrom 1.75×104to 7.00×104), the introduction of gas bubbling strengthened the liquid mixing conditions. With the increase of superficial gas velocity, the liquid mixing time decreased at first and then increased. When the gas or liquid power input kept constant, the mixing time decreased with the increase of the total power input. Through the regression analysis of all the experimental data, relationship between liquid mixing time, and liquid and gas power inputs hadbeen built up. An empirical correlation was proposed, and the calculated value was fitted well with the experimental data. Based on the obtained equation, the liquid mixing time was found to decrease at first and then increase with the increase of the gas power input if the total power input was remained constant. The transition point was around where gas input power occupied 61% of the total input power. At this point, the synergistic effect was the strongest.

date: 2015-04-07.

HUANG Zhengliang, huangzhengl@zju.edu. cn

supported by the National Natural Science Foundation of China (21406194), the Natural Science Foundation of Zhejiang Province (LQ13B060002, R14B060003) and the National Basic Research Program of China (2012CB720500).

multi-phase reactor; mixing; bubble; jet; power input

10.11949/j.issn.0438-1157.20150435

TQ 027.1

A

0438—1157（2015）11—4438—08

2015-04-07收到初稿，2015-05-25收到修改稿。

聯(lián)系人：黃正梁。第一作者：郭天琪（1990—），女，碩士。

國家自然科學(xué)基金項目(21406194)；浙江省自然科學(xué)基金項目（LQ13B060002，R14B060003）；國家重點基礎(chǔ)研究發(fā)展計劃項目（2012CB720500）。