氣動噴射數值仿真分析

邵中青，張遼遠,王 碩,呂忠秀

(沈陽理工大學 機械工程學院，遼寧 沈陽 110159)

氣動噴射數值仿真分析

邵中青，張遼遠,王 碩,呂忠秀

(沈陽理工大學 機械工程學院，遼寧 沈陽 110159)

利用CFD軟件Fluent建立數學物理模型，運用歐拉法和拉格朗日方法相結合的計算模型對離散相模型(DPM)進行三維數值模擬。研究當進口處空氣速度為60m/s時，不同直徑的固體顆粒在噴管外部流場的軌跡和橫截面上的分布。仿真得到以下結論：在固體顆粒從噴管噴出后短距離內，固體顆粒大體上呈現圓錐狀噴射，噴射的錐角隨顆粒的減小而逐步減?。活w粒直徑減小時，長軸半徑接近短軸半徑，逐步過渡為圓形。

歐拉法和拉格朗日方法相結合；離散相模型；圓錐狀噴射

固體粉末狀或顆粒狀物料可以借助氣固噴射器輸送裝置進行方便快捷的運輸，主要是通過氣體在管道內高速流動的動能實現運輸功能，其管道布置靈活，輸送的物料不會被污染，與傳統(tǒng)的帶式、鏈式輸送機輸送相比，有不可比擬的優(yōu)點，因而應用領域廣泛[1]。許多國內外學者通過不懈研究，逐漸形成了一套成熟的氣固噴射器結構和流動的理論。目前的氣固噴射器主要有傳統(tǒng)型噴射器和同心型氣固噴射器，熊源泉等[2]也曾對傳統(tǒng)氣固噴射器進行了改進和創(chuàng)新。

目前國內外已有學者對兩相流做了一些研究，主要有兩種模型:一種是歐拉-拉格朗日模型,另一種是歐拉雙流體模型[3]。其主要研究的是管道內部的流場分析，尤其集中在水射流的研究，且逐步趨于成熟，很少研究管道外部的流場分析，特別是對于固體顆粒在外部流場的噴射情況的研究更少。本文主要運用Fuent仿真軟件研究從噴頭噴出后管道外部固體顆粒的分布及軌跡，對氣相場采用歐拉法，用k-ε雙方程進行計算,對固相采用基于拉格朗日法的直接數值模擬法進行計算；對噴播設備進行探索與研究。

1 Fluent的理論依據

1.1 氣相輸送方程

(1)連續(xù)方程

式中：i、j=1、2、3，表示坐標方向，以下相同含義；ρ為流體密度;Uj為氣體速度。

(2)動量方程

式中：

當φ≤0.8時，

當φ>0.8時，

式中：Re=ρφ|Y-U|dP/μ；ξ為顆粒的球形度；dP為顆粒直徑；Re為雷諾數。

1.2 湍流動能方程和湍流動量耗散方程

(1)湍流動能方程(k方程)

由于顆粒對氣相的湍流具有一定的影響，考慮到實際因素，選擇文獻[4]中提出的湍流動能方程:

(2)湍流動量耗散方程(ε方程)

2 建立模型

建模時需滿足：CPU為Inter(R) Celeron(R) E3200 雙核；內存2GB，獨立顯卡，顯存根據主板等的性能配置。為研究固體顆粒從噴管噴出后短距離(2m以內的)的自由噴射情況，建立半徑為3m，長度為20m一個足夠大的外部流域，滿足自由噴射。整體結構是軸對稱結構，采用三維單精度方法計算，建立三維軸對稱計算模型，節(jié)省計算時間。二維模型如圖1所示，三維空間模型如圖2所示。

為研究方便，可通過建立對稱平面(由于尺寸上相差較大，只是畫出大體的示意圖)，研究這一平面上的固體顆粒軌跡，從整體上模擬出固體顆粒的噴射狀態(tài)。AI代表半徑為36mm的平面，作為氣體的入口，邊界類型設置為Velocity-inlet；CD代表半徑為36mm的平面，作為固體顆粒的入口，邊界類型設置為Velocity-inlet；KL代表半徑為3000mm的平面；JKLM為圖中所示較大的圓柱面，作為氣固混合物的出口，邊界類型設置為Outflow。其余平面尺寸大小如表1所示。

表1 模型的尺寸大小 mm

圖1 二維模型

圖2 三維空間模型

3 Fluent下的仿真分析

本文對氣相場采用歐拉法，用三維k-ε雙方程進行計算；對固相采用基于拉格朗日法的直接數值模擬法?？杉俣ㄋ芯康目臻g區(qū)域內的氣體溫度維持穩(wěn)定不變，且忽略固體顆粒與壁面間的相互作用。參照Simple算法，應用離散相模型(DPM)，對氣固兩相進行耦合求解。模型的計算條件見表2所示。

Simple算法計算步驟如下：首先假定一個速度分布，記為u0、v0、w0，以此計算動量離散方程中的系數及常數項；然后假設一個壓力場p0，依次求解動量方程，得u1、v1、w1；再對壓力加以修正，得p1，根據p1改進速度值；利用改進后的速度場求解通過源項物性等與速度場耦合的Φ變量，如果Φ變量不影響流場，則應在速度場收斂后再求解；最后利用改進后的速度場重新計算動量離散方程的系數，并利用改進后的壓力場作為下一層次迭代計算的初值。重復上述步驟，直到收斂。仿真軟件Fluent界面與Simple算法求解參數設置分別如圖3、圖4所示。

表2 本文模型計算條件

圖3 仿真軟件Fluent界面

圖4 Fluent中Simple算法的設置

3.1 同一速度(60m/s)下不同大小的顆粒軌跡

軌跡圖5～7分別是入口空氣速度為60m/s，直徑為5mm、3mm、1mm的固體顆粒在對稱截面上的軌跡分布?？梢钥闯龉腆w顆粒全部噴出所需要時間t(s)分別為1.68、1.64、1.73，時間隨著固體顆粒直徑的增大先減小后增大。

圖5 直徑D=5mm顆粒的軌跡

圖6 直徑D=3mm顆粒的軌跡

圖7 直徑D=1mm顆粒的軌跡

由圖5～7可以看出，固體顆粒全部噴出所需時間t并不隨顆粒的大小單調性的增減。固體顆粒較小或較大時，其運動的平均速度較小，全部噴出所需要的時間t較長；當固體顆粒大小適中時，其運動的平均速度較大，全部噴出所需要的時間t較短。主要原因是：能量的總損失主要包括兩部分—克服重力所造成的能量損失和克服固體顆粒與壁面的摩擦所造成的能量損失。當固體顆粒的直徑較大時，體積較大，重力相對較大，克服重力消耗的能量較多，能量損失較多，因而固體顆粒全部噴出所需要時間相對較長；當直徑很小時，體積較小，重力相對較小，克服重力消耗的能量較小，但此時一些顆粒會懸浮甚至貼附在上壁面上，克服固體顆粒與壁面的摩擦所消耗的能量較多，使能量的總損失較多，固體顆粒全部噴出所需要的時間相對較大；當固體顆粒的直徑適中時，體積剛好，重力適中，克服重力所消耗的能量相對于大顆粒少，同時貼附在上壁面上的固體顆粒數量相對于小顆粒較少，使得能量的總損減小很多，因此固體顆粒運動的平均速度較大，全部噴出所需時間較少。

3.2 同一速度(60m/s)下不同大小的顆粒在橫截面上的分布

研究大小不同的三種固體顆粒從噴嘴噴出后距離噴嘴1000mm、2000mm的橫截面上的分布情況，分布范圍如圖8～10所示(圖像關于Z=0對稱，所給圖像只顯示了Z<0的一部分)。其中D為顆粒直徑，L為橫截面距離噴嘴的距離。

圖8 顆粒直徑D=5mm的截面顆粒分布

圖9 顆粒直徑D=3mm的截面顆粒分布

圖10 顆粒直徑D=1mm的截面顆粒分布

表3 橫截面顆粒各向分布半徑 mm

結合表3和圖8～10，通過對比分析可以發(fā)現：顆粒的分布整體近似呈現出橢圓形，長軸沿Z方向，短軸沿Y方向；距離噴嘴的距離L值相同時，k1

此外，氣動能量對顆粒噴射的影響可通過氣體速度和顆粒質量流量反映出來。氣體速度一定情況下，只考慮質量流量的影響。分別選取固體顆粒質量流量(kg/s)為0.5、1、3，顆粒半徑(mm)為5、3、1進行研究。通過仿真分析可知，質量流量為1kg/s時的三種顆粒的分布半徑大致為0.5kg/s時的2倍，但質量流量為3kg/s時的分布半徑基本與0.5kg/s相同，且像1mm這樣的小顆粒噴射效果較差；同時發(fā)現質量流量對噴射錐角的影響很小。隨著顆粒質量流量的增加，顆粒分布半徑先呈現線增長，之后減小。表明在一定范圍內較大的質量流量有利于顆粒獲得較大的噴射覆蓋面積，但超出這個范圍，反而會適得其反。二者的具體關系尚需進一步深入研究。

4 結論

(1)在固體顆粒從噴管噴出后短距離內，固體顆粒大體上呈現圓錐狀噴射，噴射的錐角隨固體顆粒的增大而增大。

(2)固體顆粒全部噴出所需時間并不隨固體顆粒的大小單調性的增減，固體顆粒較小或較大時，所需時間相對較長，固體顆粒運動的平均速度較大；當固體顆粒大小適中時，所需時間較小，固體顆粒運動的平均速度較大。

(3)顆粒的分布整體近似呈現出橢圓形，當顆粒直徑減小時，長軸半徑接近短軸半徑，逐步過渡為圓形；對直徑相同的顆粒，顆粒分布的集中性隨離噴嘴的距離的增大變得相對均勻；顆粒越大，其分布的范圍就越大。

[1]盧孟龍，賈明生.直立螺旋——氣固噴射器小型給粉裝置[J].硫磷設計與粉體工程，2007，(2)：24-27.

[2]熊源泉，袁竹林，章名耀.加壓條件下氣固噴射器輸送特性的三維數值模擬[J].化工學報，2004，55(10)：1638-1643.

[3]鮑仲輔,杜群貴,謝愷.基于歐拉雙流體模型對氣固噴射器三維數值仿真[J].流體機械，2010，38(5)：24-28.

[4]Crowe C T.On models for turbulence modulation in fluid-particle flows[J].Int，J，Multiphase Flow，2000，26(5)：719-727.

[5]Bertodano M L，Lahey R T，Jones O C.Development of a K-ε model for bubbly two-phase flow[J].Trans，of the ASME，J.Fluids Eng，1994，116(1)：128-134.

[7]Elghobashi S E，Abou-Arab T W.A two-equation turbulence model for two-phase flows[J].Phys Fluids，1983，26(4)：931-938.

PneumaticJetNumericalSimulation

SHAO Zhongqing,ZHANG Liaoyuan,WANG Shuo,LV Zhongxiu

(Shenyang Ligong University,Shenyang 110159,China)

The simulation idea of CFD software Fluent was adopted,and the mathematical physical model was established,and by the calculation model combining Euler method and Lagrange method,three-dimensional numerical simulation was applied to discrete phase model (DPM).Different diameter of solid particles on the outside of the nozzle flow field on the particle trajectories and the distribution of solid particles on cross section mainly were studied,under the circumstance of the inlet air speed of 60 m/s.Through the simulation the following conclusion was got:within short distance after solid particles ejected from nozzle,solid particles generally presented the conical spray,spray cone angle increased gradually with the reduction of particle diameter decreasing;with the particle diameter decreasing,long axis radius is closer to the short axis radius,gradually transforming to the circular.

combining Euler method and Lagrange method;discrete phase model;conical spray

2013-07-10

邵中青(1986—)，男，碩士研究生；通訊作者：張遼遠(1962—)，男，教授，博士，研究方向：精密及特種加工.

1003-1251(2014)01-0081-06

O355

A

趙麗琴)