姚錫文，許開立，張秀敏，徐青偉，李 力，李季碩

(1.東北大學(xué) 資源與土木工程學(xué)院，遼寧 沈陽 110819； 2.深部金屬礦山安全開采教育部重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽 110819； 3.中加深部開采創(chuàng)新研究中心，遼寧 沈陽 110819； 4.河南省永城市安全生產(chǎn)監(jiān)督管理局，河南 商丘 476600)

0 引言

旋風(fēng)除塵器外形看似簡單，但是其內(nèi)部的氣流流動(dòng)較為復(fù)雜，屬于三維、強(qiáng)旋流、各向異性的湍流流動(dòng)，而且其中還存在滯流或回流運(yùn)動(dòng)[1]。對(duì)于生物質(zhì)氣化站而言，旋風(fēng)除塵器屬于第一級(jí)凈化設(shè)備，經(jīng)生物質(zhì)氣化爐產(chǎn)出的高溫粗燃?xì)獠⑽蠢鋮s就直接進(jìn)入到旋風(fēng)除塵器中，這就導(dǎo)致該旋風(fēng)除塵器內(nèi)部氣相流場的溫度較高，通?？梢愿哌_(dá)200℃左右[2]，而且高溫粗燃?xì)庵泻休^多的生物質(zhì)飛灰等顆粒雜質(zhì)。高溫流場的氣流流動(dòng)較常溫流場更為復(fù)雜，這就造成了生物質(zhì)飛灰顆粒在高溫流場中的運(yùn)動(dòng)和分布規(guī)律更加復(fù)雜。

為了提高旋風(fēng)除塵器在高溫下捕集生物質(zhì)飛灰顆粒的效率，必須掌握生物質(zhì)飛灰在高溫旋風(fēng)除塵器內(nèi)部高溫流場的運(yùn)動(dòng)和分布規(guī)律，而研究者們通常會(huì)借助實(shí)驗(yàn)手段對(duì)顆粒的運(yùn)動(dòng)和濃度分布進(jìn)行測試[3-6]。然而，目前高溫除塵理論和技術(shù)的發(fā)展尚不夠成熟，受實(shí)驗(yàn)條件和測試技術(shù)的限制，對(duì)于生物質(zhì)氣化站旋風(fēng)除塵器高溫流場內(nèi)部的顆粒運(yùn)動(dòng)和顆粒濃度分布規(guī)律的測試仍然較為困難。為此，采用數(shù)值模擬手段對(duì)高溫旋風(fēng)除塵器內(nèi)復(fù)雜的強(qiáng)旋湍流與飛灰顆粒的耦合運(yùn)動(dòng)特性進(jìn)行數(shù)值分析，可以較好地解決生物質(zhì)氣化站旋風(fēng)除塵器內(nèi)高溫流場與高飛灰含量等造成其顆粒濃度測量受限的問題，這可為進(jìn)一步提高旋風(fēng)除塵器在高溫下的凈化除灰能力提供理論指導(dǎo)與科學(xué)依據(jù)。

1 生物質(zhì)飛灰顆粒運(yùn)動(dòng)數(shù)學(xué)模型

1.1 顆粒場控制方程

對(duì)于高溫旋風(fēng)除塵器內(nèi)部的生物質(zhì)飛灰顆粒而言，假定其分布函數(shù)f是關(guān)于飛灰顆粒的瞬時(shí)速度us，三維空間位置xs，顆粒溫度Ts，運(yùn)動(dòng)時(shí)間t和顆粒質(zhì)量ms的函數(shù)，則f(us,xs,ms,Ts,t)dusdmsdTs表示時(shí)間為t時(shí)刻在速度間隔(us,us+dus)，質(zhì)量間隔(ms,ms+dms)和溫度間隔(Ts+dTs)內(nèi)單位體積中含有的離散顆粒數(shù)[7]，f輸運(yùn)方程的具體表達(dá)形式如式(1)所示。

(1)

式中：A為離散顆粒的加速度，左側(cè)項(xiàng)可以反應(yīng)顆粒碰撞的傳遞機(jī)制。

顆粒在運(yùn)動(dòng)過程中需要考慮其受到的阻力、重力、浮力、流體曳力、摩擦力以及顆粒之間碰撞產(chǎn)生的正應(yīng)力等。對(duì)于旋風(fēng)除塵器內(nèi)部的生物質(zhì)飛灰顆粒，其運(yùn)動(dòng)通過積分拉格朗日坐標(biāo)下顆粒作用力微分方程來求解顆粒的軌道，顆粒運(yùn)動(dòng)方程如式(2)所示。

(2)

式中：ug表示氣體的速度矢量；us表示顆粒的速度矢量；β為氣固兩相之間的曳力系數(shù)；g表示重力加速度；αs表示顆粒的體積分?jǐn)?shù)；p為氣體的平均壓力；ρs表示顆粒的密度；τs表示顆粒碰撞正應(yīng)力張量；Fs表示顆粒間的摩擦應(yīng)力張量，只有在離散顆粒濃度較高時(shí)才會(huì)比較重要。

1.2 顆粒場碰撞模型

顆粒之間相互碰撞產(chǎn)生的正應(yīng)力τs實(shí)際上是單個(gè)顆粒附近的多個(gè)顆粒對(duì)其共同作用而產(chǎn)生的合力，計(jì)算時(shí)顆粒正應(yīng)力先由離散顆粒體積插值到網(wǎng)格單元，計(jì)算之后再將其應(yīng)用到離散顆粒場[7]，τs的具體表達(dá)形式如式(3)所示。

(3)

式中：αcp表示最大堆積密度；Ps表示帶有壓力單位的常數(shù)；β為氣固兩相之間的曳力系數(shù)；δ為10-7量級(jí)的小數(shù)，用于消除數(shù)值計(jì)算奇點(diǎn)。

生物質(zhì)飛灰在旋風(fēng)除塵器內(nèi)部湍流運(yùn)動(dòng)過程中，除了受到阻力、重力、浮力、流體曳力、摩擦力以及顆粒碰撞產(chǎn)生的正應(yīng)力之外，其軌跡還受湍流作用影響，而且其湍流強(qiáng)度較大，因此需要把顆粒的湍流擴(kuò)散作用考慮在內(nèi)，故本研究中，飛灰顆粒之間的相互作用采用隨機(jī)軌道模型進(jìn)行計(jì)算。

2 物理模型的建立及邊界條件設(shè)定

2.1 物理模型的建立及其網(wǎng)格生成

利用GAMBIT軟件建立旋風(fēng)除塵器的物理模型，其幾何結(jié)構(gòu)示意如圖1所示。建模過程中，以灰斗中心為坐標(biāo)原點(diǎn)，取向上為正方向，建立坐標(biāo)系，主要分為進(jìn)氣管、環(huán)形空間、筒體空間、錐體空間、灰斗空間、排氣管、灰斗等部分，其中筒體空間和錐體空間合稱為分離空間。

圖1 旋風(fēng)除塵器模型的幾何結(jié)構(gòu)模型Fig.1 Geometry model of cyclone separator

旋風(fēng)除塵器內(nèi)部數(shù)值計(jì)算區(qū)域不規(guī)則，而且流場運(yùn)動(dòng)復(fù)雜，這給網(wǎng)格劃分帶來一定難度，網(wǎng)格劃分的方式將直接影響到FLUENT軟件后續(xù)數(shù)值解算的運(yùn)行效率、準(zhǔn)確性以及迭代結(jié)果的收斂性等情況。本文采用混合網(wǎng)格分塊劃分的模式進(jìn)行網(wǎng)格劃分，對(duì)各子區(qū)域的組成面先采用Tri-Pave面網(wǎng)格進(jìn)行劃分，然后采取Tet/Hybrid四面體網(wǎng)格對(duì)其進(jìn)行體網(wǎng)格劃分。網(wǎng)格劃分結(jié)果表明網(wǎng)格劃分質(zhì)量良好。

2.2 數(shù)值模擬邊界條件及顆粒參數(shù)設(shè)定

高溫旋風(fēng)除塵器與普通旋風(fēng)除塵器數(shù)值模擬的最大區(qū)別體現(xiàn)在溫度變化對(duì)連續(xù)相流動(dòng)的影響，而這種影響主要是通過影響連續(xù)相的粘度和密度進(jìn)行的。溫度升高，氣體的粘度增大，流動(dòng)中的氣體必將消耗一部分能量用于克服粘性阻力，故耗散的能量隨之增大。然而，氣體的密度隨著溫度的升高而降低，這在一定程度上減少了氣體流動(dòng)中的能量耗散[8]。因此，在數(shù)值模擬過程中，溫度對(duì)連續(xù)相的影響可以用粘度和密度之間的關(guān)系來反映，即通過流體的慣性力與粘性力之比，也就是溫度雷諾數(shù)(ReT)來描述。萬古軍等[8]發(fā)現(xiàn)，溫度T下的雷諾數(shù)隨溫度的升高而逐漸降低，溫度雷諾數(shù)的變化，反映了因溫度變化而使密度和粘度產(chǎn)生的能量耗散變化的關(guān)系。

對(duì)于本文所研究的高溫旋風(fēng)除塵器，其連續(xù)相為高溫氣體，在設(shè)置數(shù)值模擬初始條件和連續(xù)相邊界條件時(shí)，與常溫狀態(tài)下的空氣連續(xù)相不同，這主要體現(xiàn)在湍流動(dòng)力能量、湍流擴(kuò)散比率、湍流強(qiáng)度、溫度雷諾數(shù)、粘度系數(shù)、密度等連續(xù)相參數(shù)的設(shè)置，均與普通旋風(fēng)除塵器數(shù)值模擬的設(shè)置不同，而這些連續(xù)相參數(shù)的變化直接反映到流場的壓力分布、速度變化和顆粒濃度分布中。根據(jù)遼寧省沈陽市周邊某生物質(zhì)氣化站的現(xiàn)場實(shí)際情況可知，最初經(jīng)生物質(zhì)氣化爐產(chǎn)生的粗燃?xì)獾臏囟热愿哌_(dá)200℃左右。因此，本文采用200℃的高溫空氣近似代替生物質(zhì)粗燃?xì)庾鳛檫B續(xù)相進(jìn)行數(shù)值解算。結(jié)合旋風(fēng)除塵器現(xiàn)場實(shí)際及采用的數(shù)值計(jì)算模型與方法，邊界條件及飛灰顆粒主要參數(shù)的設(shè)置如表1所示。

表1 邊界條件主要參數(shù)設(shè)置Table 1 Main setting of boundary conditions

3 數(shù)值模擬結(jié)果分析與討論

3.1 高溫強(qiáng)旋湍流場的壓力分布

旋風(fēng)除塵器三維氣相流場的螺旋運(yùn)動(dòng)造成其壓力分布復(fù)雜，研究其內(nèi)部壓力場的變化規(guī)律對(duì)于理解旋風(fēng)除塵器分離凈化時(shí)的能量損失起到積極作用，對(duì)于提高其凈化效率具有重要的實(shí)際意義。圖2～圖4分別表示旋風(fēng)除塵器內(nèi)不同截面的靜壓分布云圖、動(dòng)壓分布云圖和全壓分布云圖。

圖2 旋風(fēng)除塵器內(nèi)三維流場的靜壓分布云圖(單位：Pa)Fig.2 Static pressure distribution of three-dimensional flow field in cyclone separator (Pa)

圖3 旋風(fēng)除塵器內(nèi)三維流場的動(dòng)壓分布云圖(單位：Pa)Fig.3 Dynamic pressure distribution of three-dimensional flow field in cyclone separator (Pa)

圖4 旋風(fēng)除塵器內(nèi)三維流場的全壓分布云圖(單位：Pa)Fig.4 Total pressure distribution of three-dimensional flow field in cyclone separator (Pa)

從圖2可以看出，沿徑向方向，靜壓分布規(guī)律呈現(xiàn)出明顯的軸對(duì)稱性，其大小由邊緣壁面向中心區(qū)域迅速降低，由此說明，外側(cè)靠近器壁下行流場的靜壓較高，而內(nèi)側(cè)中心區(qū)域上行流場的靜壓較低；另外，沿軸向方向(例如，x=0 m截面內(nèi))，靜壓在分離空間(筒體空間和錐體空間)的變化較小，中心軸線處的靜壓最低且為負(fù)壓，該負(fù)壓區(qū)一直延伸到灰斗；中心區(qū)域負(fù)壓區(qū)的存在揭示了旋風(fēng)除塵器中真空區(qū)的位置及其影響范圍。

觀察圖3可以發(fā)現(xiàn)，沿徑向方向，動(dòng)壓云圖總體上也是對(duì)稱分布的，這主要是因?yàn)榍邢蛩俣确植嫉妮S對(duì)稱性所造成的[9]；動(dòng)壓的分布規(guī)律與切向速度的分布規(guī)律較為類似，這是因?yàn)閯?dòng)壓與三維流場的合速度直接關(guān)聯(lián)，而在組成合速度的速度分量(切向速度、軸向速度和徑向速度)中，切向速度對(duì)于合速度產(chǎn)生的影響又起到了主導(dǎo)作用。

對(duì)比分析圖2和圖4發(fā)現(xiàn)，全壓的分布云圖與靜壓較為接近，沿徑向同樣呈現(xiàn)良好的軸對(duì)稱性而且在同一橫截面內(nèi)，隨著半徑減小，全壓逐漸降低，到達(dá)中心軸附近的全壓趨于零。而由伯努利方程[10]可知，全壓的降低程度，實(shí)際上反映了流場可用能量的損失程度，故從能量損失的角度來說，旋風(fēng)除塵器內(nèi)中心區(qū)域的能量損失相對(duì)最為嚴(yán)重。

3.2 生物質(zhì)飛灰顆粒的整體濃度分布

圖5是旋風(fēng)除塵器整體器壁表面的生物質(zhì)飛灰顆粒的質(zhì)量濃度分布云圖。從圖5可以看出，器壁處的顆粒濃度呈現(xiàn)螺旋帶形狀，這與文獻(xiàn)[11]中的實(shí)驗(yàn)結(jié)果吻合較好，之所以會(huì)出現(xiàn)這種螺旋帶分布，是由于生物質(zhì)灰顆粒在自身重力作用和向下的氣流的攜帶作用和離心力綜合作用的結(jié)果。

圖5 旋風(fēng)除塵器整體器壁表面的飛灰顆粒濃度分布Fig.5 Distribution of fly particles concentration on the whole wall surface of cyclone separator

圖6和圖7分別表示旋風(fēng)除塵器內(nèi)部0～180°截面(即x=0縱截面)和90～270°截面(即y=0縱截面)的生物質(zhì)飛灰顆粒的濃度分布云圖。對(duì)比分析圖6和圖7可以發(fā)現(xiàn)，在旋風(fēng)除塵器分離空間的絕大部分區(qū)域，飛灰顆粒濃度較高，呈軸對(duì)稱分布。環(huán)形空間的器壁周圍灰顆粒濃度呈現(xiàn)波浪狀，排塵口處濃度也較高，其主要是由于“顆粒返混”所致，灰斗作為收集裝置顆粒濃度也特別高。

圖6 旋風(fēng)除塵器內(nèi)部0～180°截面的飛灰顆粒濃度分布Fig.6 Distribution of fly ash particles concentration in section of 0～180° in cyclone separator

圖7 旋風(fēng)除塵器內(nèi)部90～270°截面的飛灰顆粒濃度分布Fig.7 Distribution of fly ash particles concentration in section of 90～270° in cyclone separator

圖8 沿著z軸方向不同橫截面的飛灰顆粒濃度分布Fig.8 Distribution of ash concentration in different cross sections along with the z-axis direction

圖9 沿著x軸方向不同縱切面的飛灰顆粒濃度分布Fig.9 Distribution of ash particles concentration in different vertical sections along with the x-axis

圖8為沿z軸方向(即軸向方向)，旋風(fēng)除塵器內(nèi)部不同高度(z=0.96 m，0.90 m，0.80 m，0.60 m，0.40 m，0.20 m和0 m)處的飛灰顆粒濃度分布圖。圖9表示沿x軸方向(即徑向方向)，旋風(fēng)除塵器內(nèi)部不同縱切面(x=-0.06 m，-0.03 m，0 m，0.03 m和0.06 m)的飛灰顆粒濃度分布云圖。綜合分析圖8和圖9，可以更加直觀、明顯地看出旋風(fēng)除塵器中生物質(zhì)飛灰顆粒的濃度分布情況，環(huán)形空間的頂端附近出現(xiàn)“頂灰環(huán)”現(xiàn)象，這與上述分析所得出的規(guī)律保持一致，進(jìn)一步表明了旋風(fēng)除塵器內(nèi)飛灰顆粒濃度場的分布特征。

3.3 燃?xì)膺M(jìn)口速度與分離效率的關(guān)系

生物質(zhì)飛灰的形狀不規(guī)則，而且比表面積較大，在旋風(fēng)除塵器離心力作用下，飛灰顆粒很容易發(fā)生漂移或翻轉(zhuǎn)，導(dǎo)致其運(yùn)動(dòng)復(fù)雜，這對(duì)高溫旋風(fēng)除塵器的分離效率將產(chǎn)生較大影響。因此，為了掌握高溫下燃?xì)膺M(jìn)口速度與旋風(fēng)除塵器分離效率之間的關(guān)系，在其他條件全部相同的情況下，通過設(shè)置不同的進(jìn)口風(fēng)速來分析旋風(fēng)除塵器進(jìn)、出口處飛灰顆粒的質(zhì)量流量，高溫燃?xì)獾娜肟谒俣确謩e設(shè)置為v=12 m/s，14 m/s，16 m/s，18 m/s，20 m/s，22 m/s，24 m/s和26 m/s，并由此計(jì)算旋風(fēng)除塵器在不同進(jìn)口風(fēng)速下的分離效率，進(jìn)而研究進(jìn)口風(fēng)速與分離效率的關(guān)系，具體結(jié)果如表2所示。由表2可知，旋風(fēng)除塵器的分離效率，隨著燃?xì)膺M(jìn)口速度的增加呈現(xiàn)先增大后減小的趨勢；燃?xì)鉁囟葹?00℃、進(jìn)口速度為24 m/s時(shí)，旋風(fēng)除塵器的分離效率最大，約為89.53%。

為了證明本文數(shù)值模擬研究的準(zhǔn)確性與可靠性，將本文數(shù)值模擬結(jié)果與文獻(xiàn)[12]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比研究，該實(shí)驗(yàn)測得5條不同溫度和進(jìn)口風(fēng)速下旋風(fēng)除塵器的分離效率曲線，如圖10所示。

從圖10可以看出，當(dāng)旋風(fēng)除塵器的進(jìn)口溫度不變時(shí)，隨著進(jìn)口風(fēng)速的增加，每條分離效率曲線都存在一個(gè)最高點(diǎn)，該點(diǎn)對(duì)應(yīng)的進(jìn)口風(fēng)速為最佳進(jìn)口風(fēng)速，而且分離效率以最佳進(jìn)口風(fēng)速為界，隨著進(jìn)口風(fēng)速的增加而先增大后減小，這與本文在研究燃?xì)膺M(jìn)口速度與分離效率的關(guān)系時(shí)得出的結(jié)論相符，進(jìn)而驗(yàn)證了本文數(shù)值模擬的可靠性與準(zhǔn)確性。

表2 旋風(fēng)除塵器在不同進(jìn)口風(fēng)速下的分離效率Table 2 Separation efficiencies of cyclone separatorunder different inlet velocities

圖10 不同溫度下分離效率與進(jìn)口風(fēng)速的關(guān)系Fig.10 The relationship between separation efficiency and inlet velocity at different temperature

4 結(jié)論

1) 沿徑向方向，旋風(fēng)除塵器高溫強(qiáng)旋湍流場的靜壓、動(dòng)壓和全壓分布均呈現(xiàn)良好的軸對(duì)稱性；中心區(qū)的靜壓最低且為負(fù)壓，負(fù)壓區(qū)的存在揭示了內(nèi)部真空區(qū)的位置；從能量損失的角度來說，旋風(fēng)除塵器內(nèi)部中心區(qū)的能量損失最為嚴(yán)重。

2) 高溫旋風(fēng)除塵器內(nèi)部生物質(zhì)飛灰顆粒的整體濃度分布呈螺旋狀分布；環(huán)形空間和分離空間的絕大部分區(qū)域，飛灰的顆粒濃度較高，呈軸對(duì)稱分布規(guī)律；環(huán)形空間頂端出現(xiàn)“頂灰環(huán)”現(xiàn)象，而器壁附近則呈現(xiàn)波浪狀；排塵口處濃度同樣較高，其主要是由于“顆粒返混”現(xiàn)象所致。

3) 旋風(fēng)除塵器分離效率隨燃?xì)膺M(jìn)口速度的不斷增加，呈現(xiàn)先增大后減小的趨勢，燃?xì)鉁囟葹?00℃時(shí)，分離效率在速度為24 m/s時(shí)最高，可達(dá)89.53%。模擬結(jié)果與相關(guān)實(shí)驗(yàn)數(shù)據(jù)吻合較好，驗(yàn)證了其可靠性。旋風(fēng)除塵器前端可通過安設(shè)高溫旋風(fēng)風(fēng)筒來提高其分離能力。

4) 旋風(fēng)除塵器作為生物質(zhì)粗燃?xì)鈨艋闹饕O(shè)備，在實(shí)際生產(chǎn)中，應(yīng)根據(jù)進(jìn)入旋風(fēng)除塵器粗燃?xì)獾膶?shí)際溫度來調(diào)節(jié)風(fēng)機(jī)參數(shù)，以使燃?xì)膺M(jìn)口速度調(diào)到最佳，進(jìn)而使分離效率達(dá)到最優(yōu)，提高旋風(fēng)除塵器在實(shí)際生產(chǎn)中的凈化能力。

[1] 董玉平，董磊，強(qiáng)寧，等.旋風(fēng)分離器內(nèi)生物質(zhì)焦油湍流特性的數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26(9)：171-175，385.

DONG Yuping，DONG Lei，QIANG Ning，et al. Numerical simulation of biomass gas and tar torrential flow characteristics in cyclone separator [J]. Transactions of the CSAE，2010，26(9)：171-175，385.

[2] 張秀敏. 基于風(fēng)險(xiǎn)分析的生物質(zhì)氣化及凈化模擬研究[D]. 沈陽：東北大學(xué)，2017.

[3] 吳飛雪，董守平，時(shí)銘顯.激光粒子成像技術(shù)測定旋風(fēng)分離器內(nèi)顆粒濃度場的實(shí)驗(yàn)研究[J].石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2000，24(6)：72-76.

WU Feixue，DONG Shouping，SHI Mingxian. Experimental study on particle concentration distribution in a cyclone by particle image technology [J]. Journal of the University of Petroleum, China (Edition of Natural Science)，2000，24(6)：72-76.

[4] 李永軍，易維明，何芳，等.旋風(fēng)除塵器對(duì)生物炭粉的除塵特性[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28(3)：60-64.

LI Yongjun，YI Weiming，HE Fang，et al. Removal effect of biochar powder from biomass pyrolysis by clone dust remover [J]. Transactions of the CSAE，2012，28(3)：60-64.

[5] BOHNET M. Influence of the gas temperature on the separation efficiency of aerocyclones[J]. Chemical Engineering & Processing Process Intensification, 1995, 34(3):151-156.

[6] 陳曦，葛少成，葛斐，等.針對(duì)礦山細(xì)粉塵的新型微霧旋風(fēng)除塵系統(tǒng)[J].中國安全科學(xué)學(xué)報(bào)，2014，24(9)：122-127.

CHEN Xi，GE Shaocheng，GE Fei，et al. Study on novel micro-fog & cyclonic system for removing particulate [J]. China Safety Science Journal，2014，24(9)：122-127.

[7] SNIDER D M, CLARK S M, O'ROURKE P J. Eulerian-Lagrangian method for three-dimensional thermal reacting flow with application to coal gasifiers[J]. Chemical Engineering Science, 2011, 66(6): 1285-1295.

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

WAN Gujun，WEI Yaodong，SHI Mingxian. Numerical simulation of flow field in a cyclone separator under high temperature [J]. The Chinese Journal of Process Engineering，2007，7(5)：871-876.

[9] 欒一剛.軸流旋風(fēng)分離器數(shù)值模擬與實(shí)驗(yàn)研究[D].哈爾濱：哈爾濱工程大學(xué)，2011.

[10] 吳頌平，劉趙淼. 計(jì)算流體力學(xué)基礎(chǔ)及其應(yīng)用[M]. 北京：機(jī)械工業(yè)出版社，2015.

[11] 盧小珍.柴油機(jī)排氣旋風(fēng)分離器中的流場及微粒分離規(guī)律的研究[D].北京：北京交通大學(xué)，2006.

[12] 李文琦，陳建義.旋風(fēng)分離器高溫性能試驗(yàn)研究[J].中國石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2006，30(3)：97-100，105.

LI Wenqi，CHEN Jianyi. Experimental research of cyclone performance at high temperature [J]. Journal of China University of Petroleum，2006，30(3)：97-100，105.