徐臻偉XU Zhen- 崔政偉 - 袁方洋 - 吳 璨

(1. 江蘇省食品先進制造裝備技術(shù)重點實驗室，江蘇 無錫 214122；2. 江南大學機械工程學院，江蘇 無錫 214122)

靜態(tài)混合器在液體食品工業(yè)上具有廣泛的應用[1]，如果汁、啤酒、葡萄酒及富氧水的制備等，其主要是運用了靜態(tài)混合器對氣相和液相具有良好混合功能的特性。靜態(tài)螺旋切割器是一種基于靜態(tài)混合器發(fā)展演變的在特定使用條件下具有更優(yōu)性能的一種新型氣液混合裝置，與傳統(tǒng)靜態(tài)混合器相比，靜態(tài)螺旋切割器具有更好的混合、反應、分散、傳質(zhì)和傳熱功能[2-4]，并且在特定使用情況下性能更加優(yōu)越。靜態(tài)螺旋切割器作為一種新型的微納米氣泡發(fā)生裝置[5]也具備著良好性能，在現(xiàn)有的微納米氣泡發(fā)生技術(shù)中應用較廣的是溶氣釋氣法和剪切空氣法[6-7]。靜態(tài)螺旋切割器可視為一種基于剪切氣泡的氣泡發(fā)生裝置，其內(nèi)部芯體是由若干薄片進行模塊化堆積組裝而成，具有多變性與易拆分重組性，可以根據(jù)生產(chǎn)的實際需要改變芯體的形狀，這種設計靈活性高，更能優(yōu)化其性能。

在評價靜態(tài)混合器時，壓降是決定其傳質(zhì)傳熱效果的主要參數(shù)之一，通過觀測壓降的大小可以判斷出能耗的高低。目前在對靜態(tài)混合器中氣液兩相流壓降的計算大多采用以試驗為基礎的經(jīng)驗式或半經(jīng)驗式[8]。Muddle等[9]在層流狀態(tài)下對旋流靜態(tài)混合器的流動情況進行了研究，通過數(shù)值模擬得到了壓力降因子公式；Etchells等[10]詳細介紹了Kenics等系列靜態(tài)混合器設備的結(jié)構(gòu)和工作原理，并通過試驗總結(jié)出了壓力降的計算公式；其他學者[11-13]也對靜態(tài)混合器中流場的分析做了相關模擬研究。

本研究在前人對不同種類的靜態(tài)混合器研究方法的基礎上，使用流體力學計算軟件ANSYS Fluent對新型的靜態(tài)螺旋切割器內(nèi)部氧氣-水兩相流動進行了數(shù)值計算，分析不同雷諾數(shù)下靜態(tài)螺旋切割器內(nèi)部流場的流動與壓降規(guī)律，旨在探究其流動和壓降特性，并為工程應用中靜態(tài)螺旋切割器的選用提供參考。

1 靜態(tài)螺旋切割器的結(jié)構(gòu)和工作原理

靜態(tài)螺旋切割器內(nèi)部芯腔由厚度為0.8 mm切割葉片按照螺旋線方程通過軸一片片串起疊加形成，芯體直徑25 mm，總長280 mm，形成空間螺旋面，螺旋表面呈階梯狀。由于切割葉片形狀類似十字，因此又將整個腔體分成4個小的切割空間。切割腔采用變螺距設計，即由流體入口的大螺距過渡到流體出口的小螺距，芯腔是由切割片螺旋疊加形成。圖1為靜態(tài)螺旋切割器整體裝置的結(jié)構(gòu)示意圖，主要由外管道和腔芯兩部分組成。

其中，螺旋線方程：

(1)

式中：

θ——螺旋葉片組件上任一點的扭轉(zhuǎn)角度(弧度)，θ∈[0,2π]；

R——芯體半徑，mm；

l——芯體總長，mm。

氧氣和水在管道內(nèi)在一定壓力作用下從導流錐一側(cè)流入，由于流道是由離散化切割葉片疊加旋轉(zhuǎn)形成，表面呈階梯狀，所以氣液兩相流在一定流速和壓力下經(jīng)過流道時，會受到強力剪切，同時在切割腔里也會形成旋流場，氣液兩相流在切割腔內(nèi)螺旋流動，在螺旋葉片組件的作用下，氧氣大氣泡破碎成小氣泡，從而擴大了傳遞界面，有效提高了氧的傳質(zhì)效率。

1. 連接法蘭 2. 固定環(huán) 3. 空間螺旋切割面 4. 管道 5. 鎖緊螺母 6. 軸 7. 切割腔 8. 導流錐

圖1 靜態(tài)螺旋切割器結(jié)構(gòu)示意圖

Figure 1 Sketch map of the spiral cutting device

2 理論模型及數(shù)值方法

2.1 多相流模型選擇

對兩相流中各相的數(shù)值模擬通常有2種處理方法[14]：歐拉方法和拉格朗日方法。在Fluent中，共有3種不同的歐拉多相流模型可供選擇。本試驗中氣液混合的次相體積分數(shù)為6%，氧氣經(jīng)過切割破碎后變成小量級氣泡分散在連續(xù)的液相中，在靜態(tài)螺旋切割器中充滿著單相或兩相混合物，由于含氣率較低，所以可以忽略兩相之間的速度差和溫度差以及氣體的可壓縮性對混合流體總體積的影響，即氣液兩相在流動空間內(nèi)具有連續(xù)的動量分布，因此選用流體體積(VOF)模型[15]。

2.2 控制方程

螺旋切割器氣液兩相流體介質(zhì)為氧氣和水，在模擬中選用VOF模型,因為內(nèi)部流道為變螺距的螺旋流道，氣液兩相流在流道內(nèi)發(fā)生旋轉(zhuǎn)流動，從入口端到出口端流線彎曲變化程度較大，標準k-ε模型對于該類流場流動問題的分析并不合適[16]，因此對湍流的數(shù)值計算模擬選用Realizablek-ε模型。

2.3 網(wǎng)格劃分與無關性驗證

由于切割葉片將芯腔劃分為4個旋轉(zhuǎn)對稱的獨立流道，為了簡化模型，只取其中一個流道進行仿真。在進行模擬之前，使用Workbench中的Mesh模塊對螺旋切割器結(jié)構(gòu)進行網(wǎng)格劃分。整體采用四面體非結(jié)構(gòu)化網(wǎng)格結(jié)構(gòu)，網(wǎng)格邊長為0.4 mm，網(wǎng)格節(jié)點數(shù)量為240萬個，同時劃分網(wǎng)格尺寸分別為0.20，0.24，0.30，0.34 mm，并分別計算這5種網(wǎng)格在其他條件一致的情況下壓力的變化(表1)。

表1 網(wǎng)格無關性驗證Table 1 Verification of grid independence

由表1可看出，隨著網(wǎng)格尺寸的增大，進出口壓力損失的相對變化率均小于5%，認為已經(jīng)達到網(wǎng)格無關性要求。

2.4 物性參數(shù)

在不考慮溫度壓力變化的情況下，氣液兩相流體的物理參數(shù)設定為：操作溫度20 ℃，操作壓力為標準大氣壓，氧氣密度ρ0=1.299 kg/m3，氧氣黏度μ0=1.919×10-5Pa·s；水密度ρ1=998.2 kg/m3，水黏度μ1=1.003×10-3Pa·s。

2.5 初始條件與邊界條件

進口邊界條件為速度入口條件，水的流量為0.7 m3/h，氧氣的流量為0.8 L/min，折算速度為0.72 m/s，出口邊界條件采用自由出流邊界條件。壁面采用無滑移邊界條件，采用標準壁面函數(shù)法處理邊界湍流。

2.6 離散格式、算法及控制

本模擬采用Fluent軟件,采用SIMPLE算法耦合壓力場和速度場。求解器中壓力項采用PRESTO！算法，主要參數(shù)動量、體積分數(shù)和湍動能均采用二階迎風離散格式，其他為默認值。

3 結(jié)果與討論

3.1 螺旋切割器內(nèi)的速度分布

圖2為在水流量0.7 m3/h，氧氣流量0.8 L/min條件下氣液兩相的速度跡線圖，可看出氣液兩相在流道內(nèi)螺旋流動的軌跡，圖3為切割器經(jīng)向等值切面的速度云圖。圖4為距離進口端240 mm處的截面圖，在上面取線1和線2，分別得到在這兩條線上的速度數(shù)值見圖5、6。

在進入靜態(tài)螺旋切割器時，氧氣速度高于水的速度，兩相進入管道后，在速度差的影響下，隨著氧氣和水的混合，氣相速度降低，液相速度增加，由于黏性阻力的作用，流體在管道中靠近管道壁面的速度會低于管道內(nèi)部的速度。隨著靜態(tài)螺旋切割器從入口到出口的螺距越來越小，在管道中產(chǎn)生的離心力將會越來越大。在離心力和壁面對流體阻流的共同作用下，兩相流軸向的速度在經(jīng)向截面上的分布會變得不均勻，由圖5可看出，在內(nèi)側(cè)壁面附近的速度下降，而靠近外側(cè)壁面的速度增加，其軸向速度最大值逐漸向流道的外側(cè)壁面靠近。

由圖3可看出，速度在沿著切割器旋轉(zhuǎn)方向上隨著螺距的減小而變得越來越大，由圖6可得出，速度在內(nèi)側(cè)螺旋面附近越來越大，其原因是伴隨著螺距變小，流體在流道中受到外側(cè)螺旋面的阻力越來越大，將流體向內(nèi)側(cè)螺旋面擠壓，從而使得速度增加。

圖2 螺旋切割器內(nèi)單流道的速度跡線圖Figure 2 Path lines of single flow channel's velocity

圖3 螺旋切割器內(nèi)單流道經(jīng)向截面速度沿軸線變化云圖Figure 3 The contours of single flow channel's velocity along the axis of spiral cutter

圖4 z=240 mm處的截面圖Figure 4 The sectional view of z=240 mm

圖5 線1上的速度數(shù)值圖Figure 5 The velocity magnitude of line 1

圖6 線2上的速度數(shù)值圖Figure 6 The velocity magnitude of line 2

3.2 螺旋切割器內(nèi)的壓降特性

雷諾數(shù)是流體力學中的一個無量綱數(shù)，其公式為：

(2)

式中：

Re——雷諾數(shù)；

ρ——流體密度，kg/m3；

u——流體流速，m/s；

L——特征長度，m；

μ——動力黏度，Pa·s。

本試驗所算雷諾數(shù)中黏度、密度為水的參數(shù)。

對于常見管道流動，壓力降有2種表示方法，一種表示為達西公式[17]：

(3)

式中：

△Psm——混合管的壓降，Pa；

λ——摩擦因子；

l——特征長度，m；

d——管道直徑，m；

ρm——流體密度，kg/m3；

um——流動速度，m/s。

另一種表示為Z因子表達法：

Z=△Psm/△P0，

(4)

式中：

Z——流體經(jīng)過混合器所產(chǎn)生的壓降與水平圓直管壓降的比值；

△Psm——混合器壓降，Pa；

△P0——水平圓直管壓降，Pa。

圖7所示以螺旋切割器為研究對象，得到切割器軸截面平均壓力沿軸線的變化結(jié)果，截面壓力云圖見圖8。

在保持氧氣含量不變的情況下，通過改變?nèi)肟谒俣葋砀淖兝字Z數(shù)，等距離檢測7個截面上的平均壓力，不同雷諾數(shù)下切割器內(nèi)各截面壓降沿軸線方向的變化見圖9。

為研究螺旋切割器內(nèi)氣液兩相流壓力的分布情況，在Re=5 800時考察了在z=0，20，40，60，80，100，120，140，160，180，200，220，240，260，280 mm經(jīng)向等值切面的壓力，見圖7。從圖7可以看出，裝置內(nèi)壓力從入口到出口逐漸減小，根據(jù)伯努利方程可知，此時流體的靜壓轉(zhuǎn)化為動壓，表現(xiàn)為速度增加，靜壓下降。在螺旋管道中，螺旋葉片組件對氣液兩相流體會產(chǎn)生強制擾流作用，使氣液兩相在管道流動過程中旋轉(zhuǎn)流動，湍流程度加強，從而致使裝置內(nèi)部出現(xiàn)明顯的壓力降。由圖8可看出，當氣液兩相流發(fā)生旋轉(zhuǎn)流動時，其徑向截面上也存在一定的壓力梯度，且越接近出口的截面壓力梯度越大。

在研究不同Re下壓降變化時，考察了z=0，40，80，120，160，200，240，280 mm經(jīng)向等值切面的壓力。由圖9可看出，隨著Re的增大，螺旋切割器內(nèi)的總壓降增加，在低Re的管道中，壓降變化較為舒緩，在高Re的管道中，靠近出入口的壓降變化較明顯。

3.3 氣液兩相流在螺旋切割器與水平圓直管中的壓降對比

為了進一步研究靜態(tài)螺旋切割器中壓降的變化，分別計算了不同雷諾數(shù)下的Z因子。圖10顯示了當氣體體積分數(shù)約為6%時，Z因子在不同雷諾數(shù)下的結(jié)果。通過比較，可以看出，氣液兩相流通過靜態(tài)螺旋切割器的壓降比水平圓管要大得多。Z因子隨雷諾數(shù)的增大而增大。在含氣率范圍內(nèi)，當Re=800～6 600時，Z因子的增大幅度較低，約為150；當Re=6 600～13 000時，Z因子增大幅度較大，達到220左右；當Re=13 000～26 000時，Z因子增大幅度減緩，此時螺旋切割器內(nèi)壓降是水平圓直管的260倍左右。

Re=5 800圖7 螺旋切割器內(nèi)截面平均壓力沿軸線變化圖Figure 7 Average pressure along the axis of spiral cutter

圖8 螺旋切割器內(nèi)單流道截面壓力沿軸線變化云圖Figure 8 The contours of single flow channel's pressure along the axis of spiral cutter

圖9 不同Re下螺旋切割器內(nèi)截面的壓降沿軸線變化圖Figure 9 Pressure distribution along the axis of spiral cutter under different Re

Voxy=0.06圖10 Z因子隨雷諾數(shù)的變化Figure 10 Z factor-Re curve

4 結(jié)論

本研究采用ANSYS Fluent對一種具有特殊螺旋結(jié)構(gòu)的靜態(tài)螺旋切割器內(nèi)氣液兩相流動進行仿真計算，分析該設備的內(nèi)部流場的流動與壓降規(guī)律，為評價該設備的性能提供依據(jù)。螺旋切割器內(nèi)軸線方向速度分布顯示，該設備對氣液兩相流在管道中的流動混合具有促進作用，且這樣的促進作用是由螺距的變化所決定的，由管道入口的大螺距過渡到出口的小螺距時促進效果逐漸加強。當雷諾數(shù)>13 200時，兩相流在管道中壓降顯著增大，說明此時管道對兩相流的流動作用開始加強。通過計算兩相流動的Z因子，發(fā)現(xiàn)Z因子隨Re的增加而增大；在Re處于15 000左右時，Z因子的增

幅減緩，此時螺旋切割器內(nèi)壓降是水平圓直管的225倍左右。本研究探究了螺旋切割器內(nèi)氣液兩相流動特性，著重分析了雷諾數(shù)對壓降的影響規(guī)律，但是未得出螺旋切割器內(nèi)壓力降的預測公式，這會在進一步的工作中體現(xiàn)。

[1] 孟建軍. 靜態(tài)混合器在啤酒及其他液體食品工業(yè)中的應用[J]. 中外酒業(yè)·啤酒科技, 2012(3): 43-47.

[2] 劉濤. SK型靜態(tài)混合器與噴射混合器的CFD模擬研究[D]. 青島: 青島科技大學, 2009: 3-4.

[3] 王修鋼, 郭瓦力, 吳劍華. 靜態(tài)混合器中液液分散的實驗及CFD模擬[J]. 化工學報, 2012(63): 767-768.

[4] 孟輝波, 禹言芳, 王豐, 等. 水平與豎直布置方式下靜態(tài)混合器內(nèi)瞬態(tài)壓力波動混沌特征[J]. 北京化工大學學報: 自然科學版, 2013(40): 93-93.

[5] 鄧超, 楊麗, 陳海軍, 等. 微納米氣泡發(fā)生裝置及其應用的研究進展[J]. 石油化工, 2014(43): 1 206-1 207.

[6] 劉季霖. 微納米氣泡發(fā)生裝置研究[D]. 杭州: 浙江大學, 2012: 4-6.

[7] 李恒震. 微納米氣泡特性及其在地下水修復中的應用[D]. 北京: 清華大學, 2014: 10-13.

[8] 居迎軍, 劉曉光, 高永亮, 等. 多相管流摩阻計算方法綜述[J]. 中國石油和化工, 2008(10): 55-58.

[9] MUDDE R F, PIJPEN C, BEUGELS R. Simulation of the laminar flow in a primix static mixer[C]//ASME 2002 Pressure Vessels and Piping Conference. New York: [s.n.], 2002, 448(2): 67-73.

[10] ETCHELLS A W, MEYER C F. Handbook of industrial mixing：Science and practice[M]. Hoboken: John Wiley Sons Inc, 2004: 426-431.

[11] 翟俊霞, 陳嘉南, 涂善東. 含靜態(tài)混合元件管內(nèi)流體速度測量與實驗系統(tǒng)[J]. 石油機械, 2003, 31(12): 1-3.

[12] REGNER M, ?STERGREN K, TRG?RDH C. Effects of geometry and flow rate on secondary flow and the mixing process in static mixers: A numerical study[J]. Chemical Engineering Science, 2006, 61(18): 6 133-6 141.

[13] REGNER M, ?STERGREN K, TRG?RDH C. Influence of viscosity ratio on the mixing process in a static mixer: Numerical study[J]. Industrial & Engineering Chemistry Research, 2008, 47(9): 3 030-3 036.

[14] 劉德新. 精餾塔板氣液兩相流體力學和傳質(zhì) CFD 模擬與新塔板的開發(fā)[D]. 天津: 天津大學, 2008: 23-28.

[15] 張春梅, 孔祥玲, 王宗勇, 等. SMV靜態(tài)混合器內(nèi)氣液兩相流壓降的研究[J]. 當代化工, 2011, 40(5): 490-493.

[16] 楊春玲. 新型離氏螺旋曝氣設備氣液兩相流數(shù)值模擬及實驗研究[D]. 北京: 北京化工大學, 2011: 27-29.

[17] 孫西歡, 閻慶紱, 武鵬林, 等. 圓管螺旋流起旋器結(jié)構(gòu)參數(shù)與阻力研究[J]. 流體機械, 2000, 28(10): 7-9.