導(dǎo)語：顆粒相濃度數(shù)值研究管理論文一文來源于網(wǎng)友上傳，不代表本站觀點(diǎn)，若需要原創(chuàng)文章可咨詢客服老師，歡迎參考。

摘要：本文采用成熟的湍流和顆粒模型,即基于歐拉系下的k-ε雙方程模型和把顆粒相作為連續(xù)介質(zhì)，部分地考慮顆粒相對流體運(yùn)動影響的無滑移模型，根據(jù)這種百葉窗的使用條件設(shè)定參數(shù)值，模擬出含塵氣流穿過它時的過濾特性。計(jì)算了V型百葉窗過濾器中氣固兩相流動，從而揭示了其內(nèi)部氣相壓力場的分布與顆粒相濃度場的基本特征.并將模擬結(jié)果和試驗(yàn)結(jié)果作了初步的比較和誤差分析，從而有力地闡明了利用CFD技術(shù)在工程設(shè)計(jì)實(shí)施前或產(chǎn)品生產(chǎn)（改造）前預(yù)測方案的實(shí)際效果與產(chǎn)品性能的可行性。

關(guān)鍵詞：氣固兩相流壓力場濃度場過濾特性CFD

0前言

百葉窗式過濾器是慣性除塵器的一種。它是在含塵空氣通道上裝有許多由葉板組成的傾斜的擋灰柵（其截面如圖1所示），含塵氣流在通過擋灰柵時按葉板間的縫隙數(shù)分成多股氣流，每股氣流在通過葉板時會突然改變方向，這時，粉塵顆粒在慣性力的作用下力圖保持原來的流動方向，塵粒掠過擋灰柵的縫隙入口而直接打到后面葉板的表面上，并沿繞過擋灰柵的空氣流動方向相反的方向彈出[1]。本文的研究工作就是根據(jù)這種百葉窗的使用條件設(shè)定參數(shù)值，模擬出含塵氣流穿過它時的過濾特性。

1計(jì)算模型及方法

迄今為止，在工程應(yīng)用中，對氣相流場計(jì)算可行而成熟的模型仍然是氣相作為連續(xù)介質(zhì)基于歐拉系下的k-ε雙方程模型[2]。對顆粒相的計(jì)算模型很多，本文采用的是基于歐拉系下的，把顆粒相作為連續(xù)介質(zhì)，部分地考慮顆粒相對流體運(yùn)動影響的無滑移模型[3]。無滑移模型的優(yōu)點(diǎn)是簡單，可以用較為成熟的處理單相流體的數(shù)值解法來處理兩相流，其缺點(diǎn)是不考慮顆粒相對于流體的速度及溫度滑移（滯后）。但是在本文的研究中，只考慮在室溫情況下流體、顆粒的運(yùn)動情況，不需要考慮其溫度變化，所以可以采用無滑移模型。

本文研究的過濾器是在室溫下使用，內(nèi)部溫度場穩(wěn)定，因而流體能量方程不予考慮，則氣相組的基本方程為：

流體連續(xù)方程：(1)

流體i方向動量方程：

(2)

流體湍能方程：(3)

流體湍能耗散率方程：(4)

根據(jù)無滑移模型的基本假設(shè)可知，顆粒相的動量方程及能量方程均不存在，只需給出其質(zhì)量守恒或連續(xù)方程或擴(kuò)散方程：

(5)

氣相流場計(jì)算采用基于SIMPLE方法的改進(jìn)算法SIMPLEC[4]，該方法具有改善和加快收斂性能。積分區(qū)域網(wǎng)格劃分使用交錯網(wǎng)格，對流項(xiàng)離散使用QUICK格式[5]。由于在歐拉場中，無論是流體的控制方程還是顆粒的控制方程，都具有相同的形式，不同的是源項(xiàng)。因而完全可以將原來求解單相流場的數(shù)值方法推廣到多相流中。

2模擬結(jié)果及分析

2.1氣相壓力場分布

2.1.1壓力分布圖

圖2、圖3、圖4、圖5分別是過濾器迎面風(fēng)速V=1m/s、V=2m/s、V=3m/s、V=4m/s時氣相壓力場的模擬結(jié)果。

圖2v=1m/s時壓力場圖3v=2m/s時壓力場

圖4v=3m/s時壓力場圖5v=4m/s時壓力場

2.1.2壓力場分析

（1）從壓力場可以看出在進(jìn)口處呈現(xiàn)正壓區(qū)，擋板后漩渦區(qū)域比較大，它的存在必然使擋板后形成負(fù)壓區(qū)，擋板前后形成壓差，帶來過濾器阻力損失。且過濾器的壓力損失也是隨著進(jìn)口速度的加大而增加。

（2）漩渦的存在使擋板間氣流繞流通道呈多次折向的弧形狀，這樣的流場分布特點(diǎn)極有利于增強(qiáng)氣固在此分離。

2.2粉塵顆粒的濃度場

2.2.1顆粒相濃度場的基本特征

作者在數(shù)值模擬中選取的模擬顆粒的粒徑為50μm、15μm的實(shí)驗(yàn)用標(biāo)準(zhǔn)滑石粉，其密度為2750Kg/m3。

圖6、圖7采用滑石粉粒徑為50μm，V=3m/s，進(jìn)口粉塵濃度分別為60mg/m3和100mg/m3。

2.2.2濃度場分析

（1）從四個模擬結(jié)果可以看出，在過濾器的進(jìn)口處粉塵濃度基本成平行水平線；在擋板的夾角處出現(xiàn)粉塵高濃度區(qū)，并且粉塵濃度的最大值出現(xiàn)在迎風(fēng)面的擋板夾角處。

（2）由于在擋板夾角區(qū)域存在渦流區(qū)，使得粉塵顆粒受到慣性作用易于沉積在擋板上。并且在迎風(fēng)面的擋板處出現(xiàn)氣流流速最小值，使得顆粒濃度最大值出現(xiàn)在此處。

（3）從擋板積灰效果圖可以看出擋板夾角內(nèi)壁面的攔截效果最明顯，但在夾角外壁面還是有一部分顆粒沉積。

3模擬結(jié)果與試驗(yàn)結(jié)果比較及分析

3.1過濾器壓力損失比較

把模擬和試驗(yàn)所得的數(shù)據(jù)通過最小二乘法擬合成曲線，如圖10所示。

圖10壓力損失曲線比較圖

可以看出兩者的壓力損失都是隨著過濾風(fēng)速的增加而增加，這個基本的變化趨勢是正確的，并且曲線形狀有一定的相似度，都是按開口向上的拋物線變化。但模擬預(yù)測值比試驗(yàn)測試值偏大。

3.2過濾效率比較

試驗(yàn)測試時采用中位徑為15μm試驗(yàn)用，的標(biāo)準(zhǔn)滑石粉，3m/s的過濾風(fēng)速。

圖11過濾效率曲線比較圖（d=15μm）

3.3誤差分析

經(jīng)過分析，作者認(rèn)為模擬值與實(shí)測值之間的誤差來源主要有以下幾個方面：

（1）本文中顆粒相選用的是無滑移模型（單流體模型），不考慮其相對于流體的速度及溫度滑移（滯后），與實(shí)際差別較大（預(yù)報(bào)的過程的發(fā)展大大超前于實(shí)際），因而通過該模型模擬出的結(jié)果必然與實(shí)際值之間存在誤差。

（2）障礙物固壁邊界上的速度、紊流動能、紊流動能耗散率都近似取為零，與實(shí)際情況相比存在著誤差。

（3）試驗(yàn)平臺也存在著一定的系統(tǒng)誤差，測量、讀數(shù)、計(jì)算及修正過程中又存在著隨機(jī)誤差，這都將會導(dǎo)致模擬值與測試值之間的誤差加大。

4結(jié)論

（1）本文結(jié)合過濾器的實(shí)際使用情況，建立適合的數(shù)學(xué)模型，采用無滑移模型研究了V型百葉窗過濾器內(nèi)特定的流動問題，其氣相壓力場和顆粒相濃度場分布基本上都是符合流場定性分析的結(jié)果和實(shí)際運(yùn)行情況的特征。

（2）利用數(shù)值模擬結(jié)果，計(jì)算出過濾器的壓力損失，并與試驗(yàn)測試值進(jìn)行比較和分析，就可以對模擬結(jié)果和試驗(yàn)測試值做雙向的檢驗(yàn)，以期分別對這兩種方法提出更為合理和有效的改進(jìn)方案。

（3）在工程設(shè)計(jì)及產(chǎn)品開發(fā)中應(yīng)用CFD技術(shù)，可改變傳統(tǒng)的設(shè)計(jì)程序，由于CFD軟件可以相對準(zhǔn)確地給出流體流動的細(xì)節(jié)，如速度場、壓力場、溫度場或濃度場分布的時變特性（不定常特性），因而不僅可以準(zhǔn)確預(yù)測流體產(chǎn)品的整體性能，而且很容易從流場的分析中發(fā)現(xiàn)產(chǎn)品或工程設(shè)計(jì)中的問題，據(jù)此提出改進(jìn)方案[6]。在本文中所研究的V型百葉窗過濾器，如對其擋板的形狀、排列順序、擋板的夾角進(jìn)行一系列的變化，模擬出它在不同情況下的壓力損失、過濾效率、氣粒分離特性，就可為該系列過濾器的設(shè)計(jì)，改進(jìn)提供強(qiáng)有力的技術(shù)支持。

符號說明：

xj張量坐標(biāo)k紊流動能

ρ空氣密度Gk紊能產(chǎn)生項(xiàng)

v空氣流速ε紊能耗散項(xiàng)

p空氣壓力f混合物分?jǐn)?shù)

μe有效粘性系數(shù)s顆粒源項(xiàng)

g重力加速度

模型中各通用常數(shù)據(jù)經(jīng)驗(yàn)可取為：

c1=1.44c2=1.92σk=1σε=1.33σp=0.09

參考文獻(xiàn)：

[1]吳杰，楊海霞，王淑玲.石化技術(shù),2002,8(1):45～48

[2]湯廣發(fā),呂文瑚，王漢青.室內(nèi)氣流數(shù)值模擬與模型實(shí)驗(yàn).長沙：湖南大學(xué)出版社，1989

[3]周力行.湍流兩相流動與燃燒的數(shù)值模擬.北京：清華大學(xué)出版社，1991

[4]S.V.帕坦卡.傳熱與流體流動的數(shù)值計(jì)算.北京：科學(xué)出版社，1984.

[5]陶文銓.計(jì)算流體力學(xué)與傳熱學(xué).北京：中國建筑工業(yè)出版社，1991.

[6]龔光彩.暖通空調(diào)，1999,29(6):25～27