撞擊流反應(yīng)器的流場測量及數(shù)值模擬研究

孫懷宇， 王 宇

(沈陽化工大學(xué)化學(xué)工程學(xué)院，遼寧沈陽110142)

撞擊流是使兩股流體沿軸相向流動，并在中點處的撞擊區(qū)域產(chǎn)生撞擊.撞擊流技術(shù)最早應(yīng)追溯到1953年Koppers-Totsek粉煤氣化爐的研發(fā)和應(yīng)用［1］.20世紀(jì)60年代初原蘇聯(lián)的 Elperin等學(xué)者率先開展了這一系統(tǒng)的研究［2］.20世紀(jì)70年代末，以色列Tamir［3］及其領(lǐng)導(dǎo)的研究組從設(shè)備形式及應(yīng)用上進(jìn)行了較深入系統(tǒng)的研究.我國學(xué)者從20世紀(jì)90年代開始對撞擊流進(jìn)行研究.伍沅設(shè)計了用于制備超細(xì)粉體的浸沒循環(huán)撞擊流反應(yīng)器［4］.

由流體的不同，撞擊流分為以氣相為連續(xù)相的撞擊流及以液相為連續(xù)相的撞擊流.近年來，撞擊流領(lǐng)域的研究明顯轉(zhuǎn)向以液體為連續(xù)相.這種研究重點的轉(zhuǎn)移是由于液體的密度和黏度分別比氣體大3個和2個數(shù)量級，致使兩股相向流體撞擊時發(fā)生強(qiáng)烈的動量傳遞和流團(tuán)間相互剪切、擠壓作用，還在撞擊區(qū)產(chǎn)生波動，其結(jié)果將有效地促進(jìn)微觀混合［5］.

本文主要研究撞擊流反應(yīng)器內(nèi)流體的壓力分布及流速分布，以便更好地了解流動情況，為反應(yīng)器的放大設(shè)計及應(yīng)用提供理論依據(jù).同時應(yīng)用FLUENT軟件對流場進(jìn)行數(shù)值模擬，更深入地解流場情況.

1 實驗部分

1.1 反應(yīng)器及實驗方法

使用有機(jī)玻璃制作撞擊流反應(yīng)器，設(shè)備尺寸為20 cm×10 cm×5 cm，上端開有10 cm×2 cm的長方形溢流口，兩導(dǎo)流管為Φ0.75 cm×0.1 cm，管口間距為2 cm.定義以撞擊面中心為坐標(biāo)中心的三維坐標(biāo)系，x軸與導(dǎo)流管軸線重合，y軸為垂直方向，y-z平面稱為撞擊面，x-z平面稱為水平面，x-y平面稱為垂直面.實驗使用水為工質(zhì).實驗流程如圖1所示.

圖1 設(shè)備流程圖Fig.1 Equipment flow chart

測量時用閥門調(diào)節(jié)兩管流量.壓力探頭垂直插入液面下指定位置，使用雙液倒U形管壓差計(水和柴油做指示劑)對撞擊流反應(yīng)器內(nèi)壓力進(jìn)行測量.用熱電阻作為測溫探頭，使用PCI-1710數(shù)據(jù)采集卡建立數(shù)據(jù)采集系統(tǒng)，以1 kHz的采樣率對探頭處溫度進(jìn)行測量.在流動方向上放置2個間隔一定距離的探頭，并在管出口處迅速加入一定量的熱水，通過分析兩探頭測量得到溫度脈沖峰的時間差，可得到兩探頭間的流體速度.

1.2 實驗結(jié)果與分析

圖2和圖3分別是撞擊速度u=0.4 m/s時x-z面的壓力分布和流速分布.

圖2 u=0.4m/s時x-z面壓力分布Fig.2 Pressure distribution on x-z plane when u=0.4 m/s

從圖2可以看出:壓力分布關(guān)于撞擊面基本對稱，在撞擊面的中心周圍存在一個高壓力區(qū)域，撞擊面中心的壓力最大，這是因為在兩股流體撞擊過程中部分動壓頭轉(zhuǎn)換為靜壓頭.而在此區(qū)域外的壓力較低.

圖3 u=0.4m/s時x-z面流速分布Fig.3 Velocity distribution on x-z plane when u=0.4 m/s

從圖3可以看出:速度分布關(guān)于撞擊面基本對稱.撞擊面中心周圍存在一個低速區(qū)域，撞擊面中心的速度基本為0.在此低速區(qū)域周圍存在一個錐形的流速較高的區(qū)域.從流體力學(xué)角度分析，在低速區(qū)域周圍的流動速度較高，而兩側(cè)的壓力差又較大，所以在此區(qū)域流動的流體容易產(chǎn)生旋渦，并產(chǎn)生較強(qiáng)的混合及反應(yīng)促進(jìn)作用.

2 使用FLUENT進(jìn)行流體力學(xué)數(shù)值模擬

使用FLUENT軟件建立反應(yīng)器模型［6］，并對撞擊流反應(yīng)器的流場進(jìn)行三維數(shù)值模擬.幾何模型和網(wǎng)格模型使用前處理器GAMBIT完成，求解器選為FLUENT 5/6，利用Tgrid程序?qū)ㄈ肟诠艿姆磻?yīng)器劃分為四面體網(wǎng)格.設(shè)置兩入水管出口為velocity-inlet邊界條件，頂面長方形溢流口為outflow邊界條件.為簡化計算，假設(shè)撞擊流反應(yīng)器內(nèi)流體的時均運動為穩(wěn)定流動;流體為連續(xù)、不可壓縮牛頓流體;反應(yīng)器與外界無熱量交換;不考慮重力影響;壁面為無滑移光滑界面.選用標(biāo)準(zhǔn)k-ε雙方程作為流體力學(xué)模型，此模型可對湍流發(fā)展非常充分的流體進(jìn)行計算.通過模擬計算，得到指定流速下反應(yīng)器中的壓力場及速度場.

3 實驗與模擬結(jié)果比較分析

圖4和圖5為撞擊速度u=0.4 m/s、中心水平面x=0時沿z軸壓力和速度的實驗結(jié)果與模擬結(jié)果比較.

圖4 u=0.4 m/s時z軸壓力模擬值與實驗值的比較Fig.4 Simulation compared with the experimental values of pressure on z-axis when u=0.4 m/s

在入口流速相同情況下，數(shù)值模擬與實驗所得z軸上壓力分布和速度分布基本一致.表明使用數(shù)值模擬軟件對此反應(yīng)器中的流場進(jìn)行模擬可行.但實驗值與模擬值存在一定的差別，說明使用流體力學(xué)模型在模擬撞擊流時還有一定的誤差，這方面有待于在流體力學(xué)模型方面進(jìn)行進(jìn)一步的研究及改進(jìn).數(shù)值模擬得到的在u=0.4 m/s時x-y平面和y-z平面的壓力分布和速度分布如圖6～圖9所示.

圖5 u=0.4 m/s時z軸速度模擬值與實驗值的比較Fig.5 Simulation compared with the experimental values of velocity on z-axis when u=0.4 m/s

圖6 u=0.4 m/s時x-y面壓力分布Fig.6 Pressure distribution on x-y plane when u=0.4 m/s

圖7 u=0.4 m/s時y-z面壓力分布Fig.7 Pressure distribution on y-z plane when u=0.4 m/s

圖8 u=0.4 m/s時x-y面流速分布Fig.8 Velocity distribution on x-y plane when u=0.4 m/s

圖9 u=0.4 m/s時y-z面流速分布Fig.9 Velocity distribution on y-z plane when u=0.4 m/s

從圖6～圖9可進(jìn)一步證實，在撞擊面的中心有一個高壓力低流速區(qū)域，區(qū)域關(guān)于撞擊面對稱，基本為軸對稱.流體在進(jìn)行撞擊時，繞這個區(qū)域進(jìn)行流動，在區(qū)域的周圍產(chǎn)生強(qiáng)烈的湍流，這個強(qiáng)烈的湍流區(qū)域?qū)旌线^程及反應(yīng)過程有很強(qiáng)的促進(jìn)作用.

4 結(jié)論

(1)通過實驗得到撞擊流反應(yīng)器內(nèi)的壓力分布及速度分布.證明在撞擊面中心部分存在一個撞擊區(qū)，此區(qū)域關(guān)于撞擊面對稱，基本上軸對稱.區(qū)域內(nèi)壓力較高，流速很低，圍繞此區(qū)域有一個強(qiáng)烈湍動的區(qū)域，會促進(jìn)混合及反應(yīng)的進(jìn)行.

(2)使用FLUENT對反應(yīng)器建模，使用k-ε方程進(jìn)行流體力學(xué)模擬計算，得到的模擬結(jié)果與實驗結(jié)果基本一致.

［1］ Tamir A.撞擊流反應(yīng)器——原理和應(yīng)用［M］.伍沅，譯.北京:化學(xué)工業(yè)出版社，1996:172－173.

［2］ Elperin I T.Heat and Mass Transfer in Opposing Currents［J］.J Engng Physics，1961(6):62－68.

［3］ 伍沅.撞擊流性質(zhì)及其應(yīng)用［J］.化工進(jìn)展，2001，20(11):8－13.

［4］ 伍沅.浸沒循環(huán)撞擊流反應(yīng)器:中國，CNZL002303264［P］.2000－07－06.

［5］ 伍沅.撞擊流中連續(xù)相研究重點的轉(zhuǎn)移［J］.化工進(jìn)展，2003，22(10):1066－1071.

［6］ 于勇.Fluent入門與進(jìn)階教程［M］.北京:北京理工大學(xué)出版社，2008:1－2.