生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)數(shù)值模擬

張玉1，翟明1，張波1，董芃1，朱群益1，鄒洪順2

(1.哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院，黑龍江哈爾濱150001；

2.沈陽前進(jìn)鍋爐廠，遼寧沈陽110141)

摘要：本文運(yùn)用Fluent軟件，分析了生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)。研究結(jié)果表明，RSM模型最適合生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)的模擬。與RSM模型相比，RNGk-ε模型相差最大10%，標(biāo)準(zhǔn)k-ε模型和k-ω模型相差最大50%，DES模型則在整體趨勢(shì)上不符。雙對(duì)稱入口改善了流場(chǎng)的軸對(duì)稱性，斜切入口能夠極大地減小旋風(fēng)熱解爐內(nèi)的局部渦流。綜合考慮短路流、切向速度、軸向速度、停留時(shí)間，旋風(fēng)熱解爐入口采用雙對(duì)稱30°傾角、排氣管不插入旋風(fēng)筒的形式。

關(guān)鍵詞：旋風(fēng)熱解爐；氣相等溫流場(chǎng)；湍流模型；雙對(duì)稱入口結(jié)構(gòu)；入口截面角；數(shù)值模擬

中圖分類號(hào)：TK229.6`+5

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1002-6339 (2015) 05-0413-07

Abstract：This paper analyzes the gas-phase isothermal flow field of a biomass cyclone pyrolysis burner with the Fluent software. The results show that RSM is the best suitable turbulence model for simulating the gas-phase isothermal flow field of a biomass cyclone pyrolysis burner. Compared with the RSM, the maximum difference of RNG k-εis 10%, the maximum difference of standard k-ε and k-ω is 50%. DES does not match the overall trends above. The form of double symmetrical inlet improves the axial symmetry of the biomass cyclone pyrolysis burner. The inlet with oblique angles greatly reduces the internal local vortex. The short volume is very small when the inlet section angle is over 25°. When the inlet section angle is 25°, the tangential velocity reaches maximum, when the inlet section angle is 45°, the axial velocity reaches maximum. Considering the short-circuit flow, tangential velocity, axial velocity and particle retention time synthetically, a double symmertical inlet with the inlet section angle of 30° and an exhaust pipe without stretching into the cyclone are adopted for the biomass cyclone pyrolysis burner.

收稿日期2015-01-08修訂稿日期2015-05-07

作者簡介：張玉(1985~)，女，博士研究生，研究方向?yàn)樯镔|(zhì)熱解氣化。

Numerical Simulation for the Gas-phase Isothermal Flow Field of a Biomass Cyclone Pyrolysis Burner

ZHANG Yu1,ZHAI Ming1,ZHANG Bo1,DONG Peng1,ZHU Qun-yi1,ZOU Hong-shun2

(1.Harbin Institute of Technology,School of Energy Science and Engineering,

Harbin 150001,China;2.Shenyang Qianjin Boiler Factory,Shenyang 110141,China)

Key words：biomass cyclone pyrolysis burner; gas-phase isothermal flow field; turbulence model; double symmetrical inlet; inlet section angle; numerical simulation

0引言

國內(nèi)外關(guān)于生物質(zhì)熱解氣化研究及生物質(zhì)氣化技術(shù)的應(yīng)用工藝有很多報(bào)導(dǎo)[1-4]，但眾多的氣化工藝過程仍有焦油產(chǎn)生。由于焦油粘度較大，其不僅會(huì)導(dǎo)致產(chǎn)氣率降低，還會(huì)堵塞和腐蝕設(shè)備，嚴(yán)重時(shí)使系統(tǒng)無法正常運(yùn)行。對(duì)于焦油的處理，工業(yè)上大都采用水洗法，產(chǎn)生的洗焦廢水會(huì)造成二次污染且目前洗焦廢水的生化處理工藝仍不成熟。這些問題都極大地影響了熱解氣體的后續(xù)利用及推廣使用。

在生物質(zhì)熱解氣化過程中，解決焦油最徹底的方法之一是利用高溫將焦油熱解成沸點(diǎn)低于常溫的氣體[5-6]。本文的研究對(duì)象生物質(zhì)旋風(fēng)熱解爐[7]主要是通過高溫裂解法來降低燃?xì)庵械慕褂秃?。運(yùn)行時(shí)，氣體燃料高速燃燒器產(chǎn)生高溫厭氧煙氣以高速射流的形式切向噴入生物質(zhì)旋風(fēng)熱解室中，形成強(qiáng)烈的旋風(fēng)渦旋，同時(shí)與給料倉送來的生物質(zhì)物料迅速混合、加熱并充分的熱解，最后生成各種生物質(zhì)氣，并經(jīng)排出管排出。由于此旋風(fēng)熱解爐是在高溫下運(yùn)行，且內(nèi)部等溫流場(chǎng)為強(qiáng)旋轉(zhuǎn)和高強(qiáng)度湍流運(yùn)動(dòng)，流動(dòng)規(guī)律非常的復(fù)雜，很難通過實(shí)驗(yàn)研究出其內(nèi)部流動(dòng)狀況，因此，本文應(yīng)用Fluent軟件對(duì)生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)進(jìn)行模擬，分析各種湍流模型的適用性，及結(jié)構(gòu)參數(shù)對(duì)生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)的影響，進(jìn)而為生物質(zhì)旋風(fēng)熱解爐的結(jié)構(gòu)優(yōu)化設(shè)計(jì)奠定基礎(chǔ)。

1物理建模

1.1　旋風(fēng)熱解爐的基本結(jié)構(gòu)

生物質(zhì)旋風(fēng)熱解爐如圖1所示，具體運(yùn)行過程為：利用熱解后產(chǎn)生的1 600℃以上的生物質(zhì)氣余熱的一部分，加熱加壓空氣達(dá)到800℃以上，與燃料氣在高速燃燒器中進(jìn)行燃燒，產(chǎn)生的1 600℃以上高溫高速低氧煙氣射入熱解爐，將生物質(zhì)物料充分熱解。

圖1　生物質(zhì)高溫旋風(fēng)熱解爐結(jié)構(gòu)圖

1.2　網(wǎng)格劃分

利用Gambit建立幾何模型和網(wǎng)格劃分。筒體直徑1 800 mm，進(jìn)氣口尺寸250 mm×300 mm，模型以環(huán)形空間的頂蓋中心為坐標(biāo)原點(diǎn)，沿軸向?yàn)閆軸，向下為正。采用結(jié)構(gòu)化網(wǎng)格，如圖2所示，最小網(wǎng)格尺寸為35 mm，網(wǎng)格節(jié)點(diǎn)數(shù)為419 641，網(wǎng)格數(shù)為399 071。

圖2　網(wǎng)格結(jié)構(gòu)

1.3　網(wǎng)格無關(guān)性驗(yàn)證

取入口為1 600℃的不可壓縮理想氣體，入口速度為80 m/s，入口當(dāng)量直徑為272.7 mm，出口邊界為Outflow。圖3為采用標(biāo)準(zhǔn)k-ε模型4種網(wǎng)格尺寸在截面z=2 000 mm的計(jì)算數(shù)據(jù)?？梢钥闯觯W(wǎng)格尺寸為35 mm和網(wǎng)格尺寸為25 mm、20 mm的模型參數(shù)相差很小。從計(jì)算的經(jīng)濟(jì)角度出發(fā)，選用網(wǎng)格尺寸為35 mm的模型。

2湍流模型分析

在生物質(zhì)旋風(fēng)熱解爐內(nèi)，切向速度占主導(dǎo)地位，它帶動(dòng)氣體做高速旋轉(zhuǎn)運(yùn)動(dòng)。軸向速度相對(duì)較小，但它是導(dǎo)致各種二次流的主要原因，不利于生物質(zhì)燃?xì)獾臒峤?。相?duì)于軸向速度和切向速度，徑向速度非常小(0~5 m·s)；壓降也是考查旋風(fēng)熱解爐工作性能的一個(gè)重要指標(biāo)[8]。因此，從切向速度、軸向速度、靜壓力及動(dòng)力壓力4個(gè)不同方面來考察旋風(fēng)熱解爐不同湍流模型下氣相等溫流場(chǎng)狀況。

圖4為旋風(fēng)熱解爐在不同湍流模型下計(jì)算結(jié)果對(duì)比。除了DES模型以外，其他模型模擬結(jié)果趨勢(shì)相同，其中RSM和RNGk-ε模型幾乎重合。切向速度分布曲線呈明顯的駝峰形狀，速度場(chǎng)為典型的蘭金渦[9]。軸向速度包括兩部分，即外部的下行區(qū)和內(nèi)部的上行區(qū)，兩區(qū)的分界點(diǎn)處軸向速度為0[10]，說明旋風(fēng)熱解爐中存在兩股旋流，再由切向速度分布(切向速度方向相同)可以得出上行流和下行流旋轉(zhuǎn)方向相同。爐內(nèi)靜壓力沿徑向由外向內(nèi)遞減，靠近中心軸線處的靜壓最低，靠近氣化爐器壁的靜壓最大[11]。動(dòng)壓在強(qiáng)制渦和準(zhǔn)自由渦的分界面處最大。在強(qiáng)制渦區(qū)，沿氣化爐半徑，動(dòng)壓逐漸減小，并在軸線處趨于零；在準(zhǔn)自由渦內(nèi)，沿氣化爐半徑，動(dòng)壓逐漸減小。

圖4　不同湍流模型計(jì)算結(jié)果對(duì)比(z=2 500 mm)

在生物質(zhì)高溫旋風(fēng)熱解爐氣相等溫流場(chǎng)中，因切向速度分量遠(yuǎn)大于軸向和徑向速度分量，且RSM模型直接求解各獨(dú)立雷諾應(yīng)力分量的輸運(yùn)方程，考慮了各向異性湍流輸運(yùn)的特性，因而更加適合于三維強(qiáng)旋流和高強(qiáng)度的旋風(fēng)熱解爐氣相等溫流場(chǎng)的計(jì)算。標(biāo)準(zhǔn)k-ε模型有較高的穩(wěn)定性、經(jīng)濟(jì)性和計(jì)算精度，但它是基于各向同性湍流輸運(yùn)假設(shè)，適合高雷諾數(shù)湍流，不適合強(qiáng)旋流等各向異性較強(qiáng)的流場(chǎng)，因而與RSM模型相比相差較大。RNGk-ε模型與標(biāo)準(zhǔn)k-ε模型一樣，也是二方程湍流模型，它更一般、更基本，可以計(jì)算低雷諾數(shù)湍流，考慮到旋轉(zhuǎn)效應(yīng)，模擬中選擇了Swirl Dominated Flow選項(xiàng)，對(duì)強(qiáng)旋流氣相流場(chǎng)的旋風(fēng)熱解爐計(jì)算精度有很好的改進(jìn)效果，因而與RSM模型相比相差較小。k-ω模型中考慮了低雷諾數(shù)影響、可壓縮性影響和剪切流影響，更適用于尾跡流動(dòng)、混合層和射流情況，因而與RSM模型相比相差較大。DES是最精細(xì)的湍流模型，精度高，但是需要的網(wǎng)格數(shù)量大，要求的計(jì)算量、內(nèi)存需要都非常大，由于限于設(shè)備計(jì)算能力，與RSM模型相比偏差非常大。

與RSM相比, 標(biāo)準(zhǔn)k-ε模型相差最大50%，RNGk-ε模型相差最大10%，k-ω模型相差最大45%，DES模型則在整體趨勢(shì)上不符。如果考慮工作效率，RNGk-ε模型也較適合于旋流占優(yōu)的旋風(fēng)熱解爐氣相等溫流場(chǎng)模擬，但標(biāo)準(zhǔn)k-ε模型、k-ω模型和DES模型則不適合旋風(fēng)熱解爐氣相等溫流場(chǎng)的模擬。因此，后續(xù)模擬選用RSM湍流模型。

3結(jié)構(gòu)參數(shù)分析

與旋風(fēng)分離器相似，對(duì)于常規(guī)的切向單進(jìn)口旋風(fēng)熱解爐，在爐頂板附近，排氣管外壁與筒體之間，由于徑向速度與軸向速度的存在，將形成一個(gè)流動(dòng)緩慢的邊界層，邊界層中的靜壓變化小于筒體中的靜壓變化，因此在頂板外側(cè)沿筒體的流體在靜壓的作用下進(jìn)入邊界層內(nèi)，遇到排氣管的外壁面轉(zhuǎn)彎向下，從其出口進(jìn)入排氣管，形成局部渦流，最后隨中心上升氣流逸出排氣管[12-13]，影響生物質(zhì)氣體在旋風(fēng)熱解爐中的充分熱解。因此，在旋風(fēng)熱解爐的結(jié)構(gòu)設(shè)計(jì)中，有效避免這種短路流所造成的不利影響是優(yōu)化的關(guān)鍵。

3.1　入口形式

為了避免短路流，采用了雙對(duì)稱即四進(jìn)口入口結(jié)構(gòu)形式。對(duì)于單進(jìn)口的旋風(fēng)熱解爐，入口邊界條件為1 600℃的不可壓縮理想氣體，入口速度為320 m/s，入口當(dāng)量直徑為272.7 mm，入口雷諾數(shù)為15.671×105，出口邊界條件為Outflow。雙對(duì)稱入口速度為80 m/s，入口雷諾數(shù)為3.918×105，其它條件不變。

(1)入口結(jié)構(gòu)對(duì)短路流的影響

圖5為旋風(fēng)熱解爐不同入口結(jié)構(gòu)形式速度云圖。

圖5　不同入口結(jié)構(gòu)形式速度云圖

從圖中可以看出，無論單入口，還是雙對(duì)稱入口，在靠近排氣管外壁都存在短路流現(xiàn)象。由于單進(jìn)口不對(duì)稱，使得氣流進(jìn)入旋風(fēng)熱解爐時(shí)更加紊亂，渦核容易出現(xiàn)變形。采用雙對(duì)稱入口后，氣體在旋風(fēng)熱解爐內(nèi)為多點(diǎn)對(duì)稱進(jìn)氣，使速度分布呈現(xiàn)明顯的對(duì)稱性，減小了流場(chǎng)內(nèi)部的渦流，氣體攜氣量減少。在相同的氣體流量下，由于進(jìn)口氣流速度減為原來的1/4，阻力降低，流場(chǎng)更穩(wěn)定[14]。

(2)入口結(jié)構(gòu)對(duì)切向速度的影響

圖6所示為x=0截面上z=1 000 mm時(shí)不同入口結(jié)構(gòu)的切向速度分布。切向速度分布曲線呈明顯的駝峰形狀，速度場(chǎng)為典型的蘭金渦。入口改成雙對(duì)稱結(jié)構(gòu)后，切向速度軸對(duì)稱性明顯提高，而且旋渦的中心與幾何模型的中心相重合；在相同的入口氣量下，雙對(duì)稱結(jié)構(gòu)流場(chǎng)的切向速度明顯降低，且速度分布更均勻。

圖6　Z=1 000 mm時(shí)不同入口結(jié)構(gòu)切向速度分布圖

(3)入口結(jié)構(gòu)對(duì)軸向速度的影響

圖7所示為x=0截面上z=1 000 mm時(shí)不同入口結(jié)構(gòu)的軸向速度分布。軸向速度構(gòu)成了旋風(fēng)熱解爐的外層下行、內(nèi)層上行的氣體雙層旋轉(zhuǎn)流動(dòng)結(jié)構(gòu)，兩層氣流的分界點(diǎn)處軸向速度為0。采用雙對(duì)稱入口結(jié)構(gòu)后，軸向速度對(duì)稱性明顯改善。

圖7　z=1 000 mm時(shí)不同入口結(jié)構(gòu)軸向速度分布圖

(4)入口結(jié)構(gòu)對(duì)總壓的影響

圖8所示為x=0截面上z=1 000 mm時(shí)不同入口結(jié)構(gòu)的總壓分布?？倝毫ρ貜较蛴赏庀騼?nèi)遞減，靠近軸線處的總壓最低，靠近壁面的總壓最大。在相同的流量下，雙對(duì)稱進(jìn)口的壓力損失比單進(jìn)口的明顯減小，這是由于采用雙對(duì)稱進(jìn)口后，入口速度降低了3/4。改進(jìn)后總壓在2 000～3 000 Pa之間，非常穩(wěn)定，且更加節(jié)能。

圖8　Z=1 000 mm時(shí)不同入口結(jié)構(gòu)總壓分布圖

3.2　入口截面角

氣體進(jìn)入常規(guī)的直切式進(jìn)口旋風(fēng)熱解爐后，必定要經(jīng)過排氣芯管外壁和筒體內(nèi)壁之間，不可避免會(huì)存在“短路流”。如果具有一定的入口截面角，會(huì)使生物質(zhì)氣體在重力和慣性力的作用下直接進(jìn)入芯管以下，不僅能使生物質(zhì)充分熱解，還能更大程度地減小芯管底部“短路流”。同時(shí)，旋風(fēng)熱解爐是在高溫下工作，為了防止磨損，排氣管不伸入旋風(fēng)筒中。

(1)不同入口截面角速度云圖

圖9所示是不同入口截面角在x=0截面上速度云圖。當(dāng)進(jìn)口角度為0°~15°時(shí)，在環(huán)形空間頂蓋處確實(shí)存在一定的短路流，當(dāng)進(jìn)口角度為15°~30°時(shí)，短路流量逐漸減少，進(jìn)口角度為40°~50°時(shí)，幾乎沒有短路流，且進(jìn)口周邊出現(xiàn)速度為零的區(qū)域。因此，氣體從進(jìn)氣口進(jìn)入旋風(fēng)筒后，如果入口角度安排的適當(dāng)，可減少短路流，使生物質(zhì)能夠在爐中充分熱解。

圖9　不同入口截面角速度云圖

(2)不同入口截面角對(duì)切向速度的影響

圖10所示為入口風(fēng)速為80 m/s時(shí)截面z=1 000 mm和z=2 000 mm處切向速度分布。不論入口是水平還是具有一定截面角，爐內(nèi)切向速度均呈現(xiàn)典型的蘭金渦結(jié)構(gòu)：外旋渦邊界速度為零，速度梯度較大，內(nèi)渦旋中速度梯度較小，內(nèi)外渦旋交界面上切向速度值最大。入口具有一定的截面角后，切向速度在數(shù)值上均有所降低，主要是由于傾斜速度一部分用來旋轉(zhuǎn)，一部分用來向下流動(dòng)，但切向速度仍非常大，可以保持一定的旋轉(zhuǎn)速度，使生物質(zhì)氣充分混合。隨著傾斜角度的增加，切向速度先減小后增加再減小。入口截面角為25°時(shí)的切向速度值最大。

圖10　不同入口截面角切向速度分布

(3)不同入口截面角對(duì)軸向速度的影響

圖11所示為入口風(fēng)速為80 m/s時(shí)截面z=1 000 mm和z=2 000 mm處切向速度分布。不論入口是水平還是具有一定截面角，爐內(nèi)流動(dòng)軸向速度沿徑向可分為上行流和下行流?？拷诿嫣幍南陆盗髁孔钅芊从尺M(jìn)入爐內(nèi)的氣體流量，隨著深度的增加，軸向速度在減小，但是下降流量在逐漸增加，表明進(jìn)入爐內(nèi)的氣體大部分向下流動(dòng)。隨著入口截面角的增加，軸向速度先減小后增加再減小，當(dāng)入口截面角為45°時(shí)，軸向速度達(dá)到最大。

(4)平均停留時(shí)間

根據(jù)Rosin等人提出來的停留時(shí)間模型，可以計(jì)算出傾斜角度為25°時(shí)，氣體平均停留時(shí)間為0.471 77 s；當(dāng)傾斜角度為30°時(shí)，氣體平均停留時(shí)間為0.493 s；當(dāng)傾斜角度為40°時(shí)，氣體平均停留時(shí)間為0.553 s；當(dāng)傾斜角度為45°時(shí)，氣體平均停留時(shí)間為0.595 s。

綜上分析，取傾斜角度大于20°時(shí)，短路流量非常少。從切向速度分布來看，角度為25°~45°之間時(shí)切向速度較大，25°最大，其次是30°，相差不到10 m/s，45°與最大值相差15~20 m/s，為使生物質(zhì)充分混合角度取25°~30°合適。從軸向速度來看，45°時(shí)軸向速度最大，其它角度下的軸向速度相差不大。從停留時(shí)間來看，各角度差別不大。綜合考慮，采用30°可以滿足各種條件，因此，旋風(fēng)熱解爐入口采用30°的傾斜角度。

4結(jié)論

本文對(duì)生物質(zhì)旋風(fēng)熱解爐氣相等溫流場(chǎng)進(jìn)行了數(shù)值模擬，得到的結(jié)論如下：

(1)RSM模型最適合旋風(fēng)熱解爐流場(chǎng)模擬，如果考慮計(jì)算效率，RNGk-ε模型也較適合。與RSM模型相比，RNGk-ε模型相差最大10%，標(biāo)準(zhǔn)k-ε模型和k-ω模型相差最大50%，DES模型則在整體趨勢(shì)上不符。

(2)雙對(duì)稱入口改善了流場(chǎng)的軸對(duì)稱性，切向速度、軸向速度及總壓均有所降低；斜切入口能夠極大地減小旋風(fēng)熱解爐內(nèi)的局部渦流：傾斜角度大于25°時(shí)，短路量非常少；傾斜角度為25°時(shí)切向速率達(dá)到峰值，傾斜角度為45°時(shí)軸向速率達(dá)到峰值。

(3)綜合考慮短路流、切向速度、軸向速度、停留時(shí)間等因素，旋風(fēng)熱解爐入口采用雙對(duì)稱30°傾角、排氣管不插入旋風(fēng)筒的形式。

參考文獻(xiàn)

[1]Prabir Basu.Biomass Gasification,Pyrolysis and Torrefaction,Second Edition: Practical Design and Theory[M].Academic Press,New York,2013.

[2]Erik Dahlquist.Technologies for Converting Biomass to Useful Energy:Combustion,Gasification,Pyrolysis,Torrefaction and Fermentation[M].CRC Press,London,2013.

[3]韋姝，魏海坤，朱婷婷，等.考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究[J].電網(wǎng)與清潔能源，2015，31(3)：93-123.

[4]馮語晴，楊建華，張菁，等.我國石巖氣發(fā)電的經(jīng)驗(yàn)借鑒及發(fā)展前景分析[J].電網(wǎng)與清潔能源，2015,31(6):104-110.

[5]Hamaoka T,Inoue S,Uno S,et al.Effect of Woody Biomass Components On Air-steam Gasification[J].Biomass & Bioenergy.2005,28(1):69-76.

[6]Pian C C P,Yoshikawa K.Development of a High-Temperature Air-Blown Gasification Systerm[J].Bioresource Technol.2001,79(3):231-241.

[7]張玉，董芃，翟明，等.兩段式生物質(zhì)旋風(fēng)高溫?zé)峤鈿饣癄t:中國.201110094281.8[P].2011-10-09.

[8]孫蘭義，崔銘偉，李軍，等.不同結(jié)構(gòu)循環(huán)旋風(fēng)分離器流場(chǎng)的數(shù)值分析[J].石油化工設(shè)備,2009(5):32-36.

[9]韓婕，劉阿龍，彭東輝，等.旋風(fēng)分離器兩相流動(dòng)數(shù)值模擬研究進(jìn)展[J].天然氣化工,2012(5):55-61.

[10]萬古軍，魏耀東，時(shí)銘顯.高溫條件下旋風(fēng)分離器內(nèi)氣相流場(chǎng)的數(shù)值模擬[J].過程工程學(xué)報(bào),2007(5):871-876.

[11]趙宏強(qiáng)，蔣海華，謝武裝.基于大渦模擬的旋風(fēng)分離器內(nèi)流場(chǎng)數(shù)值模擬研究[J].環(huán)境工程學(xué)報(bào),2009(4):759-763.

[12]陳啟東，左志全.不同側(cè)向入口旋風(fēng)分離器流場(chǎng)數(shù)值分析[J].中國工程科學(xué),2014(2):58-67.

[13]Alex C.Hoffmann,Louis E.Stein.Gas Cyclones and Swirl Tubes:Principles,Design,and Operation[M].Springer,Berlin,2014.

[14]吳彩金，馬正飛，韓虹.斜切雙入口旋風(fēng)分離器流場(chǎng)的數(shù)值模擬[J].煤礦機(jī)械,2010(9):61-63.