基于EMMS模型垂直向上密相懸浮氣力輸送過程數(shù)值模擬研究

潘剛

摘 要：基于EMMS模型，建立垂直向上密相懸浮輸送過程的數(shù)學(xué)模型。對輸送管道內(nèi)固體顆粒的速度、體積分?jǐn)?shù)以及管道內(nèi)壓力變化進(jìn)行了研究。研究結(jié)果表明：固體顆粒在管道發(fā)生團(tuán)聚現(xiàn)象；在高度方向上，固體顆粒速度呈現(xiàn)出先減小后趨于穩(wěn)定的變化趨勢，固相顆粒平均體積分?jǐn)?shù)呈現(xiàn)出先增加后平穩(wěn)的趨勢；在高度h=3.5～8 m范圍內(nèi)氣力輸送過程趨于穩(wěn)定。

關(guān)鍵詞：EMMS模型 氣力輸送 數(shù)值模擬

中圖分類號：TH411 文獻(xiàn)標(biāo)識碼：A 文章編號：1674-098X（2015）09（c）-0009-02

Numerical Simulation of Vertical Upward Dense Phase Suspension Pneumatic Conveying Based on EMMS Model

Pan Gang

（Inner Mongolia University of Science and Technology， School of energy and environment，Baotou Inner Mongolia，014010，China）

Abstract：The mathematical model of vertical upward dense phase suspended transportation process is established based on EMMS model.The velocity and volume fraction of solid particles in the conveying pipeline and the pressure changes in the pipeline are studied.Research results show that agglomeration of solid particles in the pipeline.The velocity of solid particles decreases firstly and then tends to be stable in height direction.The average volume fraction of the solid particles increased first and then tends to be stable in height direction.The pneumatic conveying process is stable in the height h=3.5～8 m range.

Key Words：EMMS model；Pneumatic conveying；Numerical simulation

氣力輸送顆粒物料是工業(yè)中常見的輸送方式，是復(fù)雜的氣-固兩相流動過程[1]。輸送管道內(nèi)氣固兩相的流動狀態(tài)影響氣力輸送過程的安全性和穩(wěn)定性，相關(guān)研究者對氣力輸送過程進(jìn)行了研究。東南大學(xué)[2]通過實驗分析干煤粉在管道中發(fā)生兩相流動的過程，并且利用其特性研究工業(yè)問題。而原田幸夫等人[3]通過管道內(nèi)氣力輸送測速實驗研究，得到影響固相懸浮速度提高的原因很大一部分是由于壓降的增加。針對垂直放置的管道內(nèi)大粒徑固相顆粒數(shù)量少的情形，Tsuji等人[4]利用Lagrange的研究方法進(jìn)行了數(shù)值模擬計算。Ottjies等人[5]對發(fā)生在水平放置的管內(nèi)固相顆粒的運動進(jìn)行了數(shù)值模擬，并研究了固相顆粒與管子壁的碰撞。浙江大學(xué)[6]研究了在管道內(nèi)發(fā)生稠密兩相流動的情形，并建立了包含有連續(xù)介質(zhì)以及顆粒軌道模型的模型。以上研究通過實驗手段，采用雙流體模型的數(shù)值模擬方法對氣力輸送過程中氣固兩相的速度分布，固相顆粒體積分?jǐn)?shù)進(jìn)行了研究。但采用EMMS模型對氣力輸送過程的研究較少，該文基于EMMS模型，對垂直向上密相懸浮氣力輸送過程進(jìn)行了數(shù)值模擬研究。

1 物理模型

該文所研究的密相懸浮氣力輸送過程是發(fā)生在垂直的管道中，管長為8 m，管徑為8 cm，如圖1所示。固體顆粒密度為1020 kg/m3，粒徑為1 mm；固相顆粒入口速度為6.5 m/s，氣體入口速度為17.5 m/s。采用二維軸對稱模型研究氣力輸送過程的氣固兩相流動過程。

2 數(shù)學(xué)模型

采用雙流體模型模擬密相懸浮氣力輸送過程。氣固兩相之間的作用力主要考慮相間曳力，曳力模型采基于EMMS（Energy- Minimization Multi-Scale）的曳力模型，由中科院所開發(fā)的EMMS軟件中的曳力計算模塊（EMMS/Matrix）進(jìn)行計算。EMMS/Matrix曳力模型是基于EMMS思想[8]，EMMS模型對非均勻結(jié)構(gòu)和氣固相間作用進(jìn)行了多尺度分解，用結(jié)構(gòu)參數(shù)描述各個尺度內(nèi)氣固兩相的相互作用。模擬采用FLUENT6.3.26作為求解器。模擬計算的時間步長為0.0005 s。

3 計算結(jié)果及分析

3.1 固相顆粒平均速度

固相顆粒平均速度沿高度方向的變化規(guī)律如圖2所示。固相顆粒平均速度的變化規(guī)律呈先下降后穩(wěn)定的變化趨勢。在管道高度方向h=0～3.5 m區(qū)間內(nèi)，固相顆粒平均速度曲線呈線性下降的趨勢，隨著管道高度的增加，平均速度Us越小，從6 m/s下降為2 m/s；即在氣固混合物沿管道上升的過程中，固相顆粒的平均速度是逐漸減小，顆粒團(tuán)聚現(xiàn)象的發(fā)生導(dǎo)致固相顆粒平均速度的減小。在高度方向h=3.5～8 m區(qū)間內(nèi)，固相顆粒平均速度雖然有微小波動，但大致呈穩(wěn)定的狀態(tài)，固相顆粒速度約為2 m/s。即隨著管道高度的增加，固相顆粒平均速度的變化越來越平緩，這意味著氣固混合物在管道中的運動逐漸趨于穩(wěn)定。

3.2 固相顆粒的平均體積分?jǐn)?shù)

圖3為固體顆粒體積分?jǐn)?shù)的瞬時值和時間平均值，在管道內(nèi)呈不均勻分布。從圖3（a）可知固體顆粒局部體積分布較高，表明固體顆粒在管道內(nèi)發(fā)生了團(tuán)聚；從圖3（b）可知固體顆粒體積分?jǐn)?shù)的時間平均值，表現(xiàn)為中間的顆粒體積分?jǐn)?shù)較高，而邊壁的固體顆粒體積分?jǐn)?shù)較低。

圖4為固體顆粒平均體積分?jǐn)?shù)在高度方向上的變化規(guī)律。在高度方向上，顆粒體積分?jǐn)?shù)表現(xiàn)出先增加，再減小，并逐漸趨于穩(wěn)定的變化趨勢，如圖4所示。在管道高度方向h=0～3.5 m區(qū)間內(nèi)，固相顆粒速度逐漸減小，顆粒出現(xiàn)聚團(tuán)，導(dǎo)致固相顆粒體積分?jǐn)?shù)逐漸增加，從0.04逐漸增加到0.1。而在管道高度方向h=3.5～8 m區(qū)間內(nèi)，固相顆粒平均體積分?jǐn)?shù)曲線有微小先下降后平穩(wěn)的趨勢，即從0.1逐＼漸下降為0.08，并趨于穩(wěn)定。即隨著管道高度的增加，固相顆粒平均速度越來越平緩，使固相顆粒平均體積分?jǐn)?shù)越來越平緩，管道內(nèi)固體顆粒的流動趨于穩(wěn)定。

3.3 管道內(nèi)壓力的變化

壓力隨著管道高度方向的變化規(guī)律如圖5所示，管道內(nèi)的壓力沿著管道高度方向呈逐漸降低的趨勢。在管道高度方向h=0～3.5 m的這個區(qū)間內(nèi)，隨著管道高度的升高，氣固混合物的絕對壓力下降的較為緩慢，在高度方向h=3.5～8 m區(qū)間內(nèi)，氣固混合物的絕對壓力呈線性趨勢下降。表明在管道中由于氣固運動過程中能量的耗散，壓力始終是呈下降趨勢。

4 結(jié)語

建立了垂直向上密相懸浮氣力輸送過程的數(shù)學(xué)模型，氣-固之間的曳力采用EMMS模型，對氣力輸送過程的氣固兩相流動進(jìn)行了研究，得到主要結(jié)論如下所述。

（1）垂直向上密相懸浮氣力輸送管道內(nèi)存在顆粒團(tuán)聚現(xiàn)象，導(dǎo)致了固體顆粒平均速度在高度方向呈現(xiàn)出先減小，隨后趨于穩(wěn)定的變化趨勢。

（2）固相顆粒平均體積分?jǐn)?shù)呈現(xiàn)出先增加后平穩(wěn)的趨勢。即隨著管道高度的增加，固相顆粒平均速度越來越平緩，使固相顆粒平均體積分?jǐn)?shù)趨于平緩，管道內(nèi)固體顆粒的流動趨于穩(wěn)定。在高度h=3.5～8 m范圍內(nèi)氣力輸送過程趨于穩(wěn)定。

（3）由于輸送過程中能量的耗散，管道內(nèi)的壓力呈逐漸下降趨勢。

參考文獻(xiàn)

[1] 謝灼利.密相懸浮氣力輸送過程及其數(shù)值模擬研究[D].北京化工大學(xué)，2001：3-4.

[2] 許傳龍.氣固兩相流顆粒荷電及流動參數(shù)檢測方法研究[D].南京：東南大學(xué)，2006.

[3] 原田幸夫，成井浩，島田克彥，等.利群的速度[J].日本機(jī)械學(xué)會論文集（第二部），1965，30（210）：231.

[4] TSUJI Y.Numerical simulation of gas-solid two-phase flow in a vertical pipe[J].Gas-Solid Flows ASME，1991，121（2）：123-128.

[5] OTTJIES J A.Digital simulation of pneumatic transport[J].Chemical Engineering Science，1978，33（9）：783-786.

[6] 馬銀亮.高濃度氣固兩相流的數(shù)值模擬研究[D].浙江大學(xué)，2001.

[7] 李靜海.顆粒流體復(fù)雜系統(tǒng)的多尺度模擬[M].北京：科學(xué)出版社，2005.