選粉機三維動態(tài)流場數(shù)值研究

豆海建，曾榮，唐清華，柴星騰，聶文海，申占民

選粉機屬于粉磨系統(tǒng)的核心裝備，其性能指標(biāo)的好壞直接影響粉磨系統(tǒng)的產(chǎn)量、質(zhì)量、粉磨電耗及穩(wěn)定性。選粉機的性能指標(biāo)又直接決定于其內(nèi)部的流場結(jié)構(gòu)，因此國內(nèi)外均對選粉機流場開展了較多的理論研究[1～9]，這些研究方法同本文不同：

（1）采用多重參照系（Multiple Reference Frame model，簡寫MRF）方法處理轉(zhuǎn)子的轉(zhuǎn)動過程

多重參照系主要用于計算域內(nèi)不同子域以不同的速度旋轉(zhuǎn)或移動的流體計算，由于在各自的計算子域內(nèi)采用其對應(yīng)的參照系，相對速度梯度小，降低數(shù)值求解誤差。但多重參照系沒有考慮旋轉(zhuǎn)或移動子域?qū)ζ湎噜徸佑虻南鄬\動，只是對旋轉(zhuǎn)或移動參照系下的控制方程中添加由于相對運動產(chǎn)生的動量、質(zhì)量、能量及它標(biāo)量源項，即旋轉(zhuǎn)或移動子域的計算網(wǎng)格是固定的，這相當(dāng)于將旋轉(zhuǎn)或運動子域中的動部件固定在某一特定位置，研究流體相對于該特定位置的流動情況。對于選粉機流場研究而言，采用多重參照系解決選粉機轉(zhuǎn)子的轉(zhuǎn)動問題就是將選粉機轉(zhuǎn)子固定在相對于靜子的某一特定位置，然后人為讓轉(zhuǎn)子域的流體以給定的速度轉(zhuǎn)動，研究其相對于轉(zhuǎn)子流動情況，即轉(zhuǎn)子內(nèi)流體的流動是主動因素，這背離了選粉機轉(zhuǎn)子域流體的流動是由于轉(zhuǎn)子葉片轉(zhuǎn)動而被動流動的實際物理過程。

（2）采用定常流動計算方法

選粉機轉(zhuǎn)子動葉片在轉(zhuǎn)動一周的過程中，周期性地經(jīng)過相鄰兩靜葉片之間扇形區(qū)域，在旋轉(zhuǎn)經(jīng)過每扇形區(qū)域過程中，由于動葉片相對于靜葉片的位置在不停地變化造成動靜葉片間的流體空間也在不停地變化，因此即便轉(zhuǎn)子以恒定的轉(zhuǎn)速轉(zhuǎn)動，選粉機內(nèi)的流場仍為典型的非定常流動。

針對國內(nèi)外目前對選粉機流場研究存在的不足之處，本文采用滑移網(wǎng)格技術(shù)和非定常流動求解方法研究選粉機三維動態(tài)流場，不僅證實了轉(zhuǎn)子葉片間及轉(zhuǎn)子內(nèi)流體的流動是由于轉(zhuǎn)子轉(zhuǎn)動引起的被動流動，而且能夠同時研究固體顆粒在選粉機空間和時間上的運動行為，更接近選粉機內(nèi)部氣固流動實際過程。

1 幾何物理模型

本文以TRM56.4輥磨選粉機為研究對象，以降低選粉機阻力為結(jié)構(gòu)優(yōu)化目的，采用三維建模軟件構(gòu)建選粉機流體域?qū)嶓w幾何模型，為采用滑移網(wǎng)格技術(shù)，以選粉機動-靜葉片1/2間隙柱面位置及轉(zhuǎn)子上面端面位置對選粉機幾何模型進行分割（圖1）。將生成的三維幾何模型導(dǎo)出parasolid格式文件（*.x_t），并導(dǎo)入網(wǎng)格生成軟件Gambit2.4.6，采用網(wǎng)格加密和網(wǎng)格過渡技術(shù)，生成圖2所示的計算網(wǎng)格。為研究幾種常見動葉片的阻力特性，對TRM、LV、Leosche及L型動葉片結(jié)構(gòu)并基于相同的動葉片間隙風(fēng)速（動力相似）工藝計算分別進行幾何建模和計算，動葉片結(jié)構(gòu)形式及尺寸見圖3。

表1 選粉機計算邊界參數(shù)

2 計算邊界條件

選粉機計算邊界參數(shù)見表1。

3 數(shù)值求解

本文采用CFD商用軟件Fluent6.3非定常求解器，湍流模型采用RNG[9]湍流模型，速度-壓力耦合方程采用PISO算法[10]求解，差分格式除壓力采用PRESTO!格式[11]外，其余均采用 Third-Order MUSCL 格式[12]，滑移網(wǎng)格求解設(shè)置如下：

（1）定義轉(zhuǎn)子流體域為rotor，類型為Fluid，將rotor內(nèi)轉(zhuǎn)子葉片、筋板、轉(zhuǎn)子軸等所有隨轉(zhuǎn)子一起轉(zhuǎn)動的幾何結(jié)構(gòu)面統(tǒng)一定義為固壁rotorwall，類型為WALL。

（2）定義圖2所示滑移網(wǎng)格界面為INTERFACE類型邊界，輸出各流體域網(wǎng)格，并逐個導(dǎo)入CFD求解器，分別將Slide1-1與Slide1-2、Slide2-1與Slide2-2、Slide3-1與Slide3-2定義 Slide1、Slide2、Slide3三個滑移網(wǎng)格副；其次，將轉(zhuǎn)子流體域rotor設(shè)置為動網(wǎng)格（Moving mesh），設(shè)置自轉(zhuǎn)revolving運動類型，并根據(jù)轉(zhuǎn)子轉(zhuǎn)動方向按右手定則定義轉(zhuǎn)軸矢量坐標(biāo)，依據(jù)轉(zhuǎn)子設(shè)計或工作轉(zhuǎn)速計算并設(shè)置rotor轉(zhuǎn)速（rad/s）。

（3）將轉(zhuǎn)子rotor-wall壁面設(shè)置為跟隨相鄰網(wǎng)格運動，將相對運動速度值設(shè)置為0m/s，實現(xiàn)轉(zhuǎn)子葉片及內(nèi)部筋板、轉(zhuǎn)軸按照流體域rotor設(shè)定的轉(zhuǎn)速轉(zhuǎn)動。

（4）按表1所示的選粉機喂料量及喂料粒度分布定義DPM多相流模型，設(shè)置選粉機入口為顆粒相注入口，入口類型為面Surface，入射方向沿氣流方向。

（5）求解器采用非定常求解器，并按照轉(zhuǎn)子轉(zhuǎn)速α（rad/s）、單個轉(zhuǎn)子葉片掃過相鄰靜子葉片間隙所經(jīng)過的角度θ（rad），計算數(shù)值求解最大時間步長τmax、τinitial：

其中n≥2整數(shù)，視實際求解過程的穩(wěn)定性取值，n值越大求解越穩(wěn)定，但計算收斂速度越慢。

（6）初始化求解程序，設(shè)置求解殘差及選粉機進、出口壓力監(jiān)控曲線，開始迭代計算，求解過程中視殘差曲線的穩(wěn)定性調(diào)整計算時間步長τ（τinitial≤τ≤τmax），待選粉機進、出口壓力監(jiān)控曲線水平或長時間規(guī)律波動、選粉機入口與成品出口和粗粉回料出口質(zhì)量流速平衡，即可視為求解收斂。

4 計算結(jié)果與討論

表2 幾種常見動葉片結(jié)構(gòu)選粉機計算阻力

圖3所示四種不同的動葉片結(jié)構(gòu)在57r/min、35r/min兩種轉(zhuǎn)子轉(zhuǎn)速工況下選粉機阻力計算數(shù)據(jù)見表2。各方案采用相同結(jié)構(gòu)靜葉片，且動葉片數(shù)目基于動葉片間相同的靜態(tài)間隙風(fēng)速條件計算得到，在相同的轉(zhuǎn)子轉(zhuǎn)速、選粉機風(fēng)量和喂料量條件下，表2所示的選粉機方案具有動力相似的條件，因而計算得到的選粉機阻力具有可比性，即計算得到的各方案選粉機阻力主要是由于動葉片結(jié)構(gòu)不同造成的。由表2可見，在兩種不同的轉(zhuǎn)子轉(zhuǎn)速條件下，LV動葉片方案具有最小的選粉機阻力；此外，在35r/min轉(zhuǎn)子轉(zhuǎn)速下，LV動葉片同其他方案的阻力差最大613.4Pa，最小345.2Pa；在57r/min條件下，阻力差最大 957.7Pa，最小 424.8Pa，這說明轉(zhuǎn)子速度越高，其低阻優(yōu)勢越明顯。

根據(jù)圖3所示選粉機動葉片結(jié)構(gòu)，LV動葉片區(qū)別于其他動葉片結(jié)構(gòu)的共同點在于其外側(cè)擋板，該擋板能有效克服動葉片間氣流的旋轉(zhuǎn)離心力阻礙來自靜葉片的氣流進入動葉片間隙（圖4）；如無該擋板，動葉片間的旋轉(zhuǎn)氣流的離心力就在一定程度上阻礙靜葉片氣流進入動葉片間隙（圖5），從而增加阻力。Leosche動葉片雖有外側(cè)擋板，但由于該擋板同葉片直段存在42°的夾角，在隨其內(nèi)側(cè)轉(zhuǎn)子高速旋轉(zhuǎn)的過程中，該擋板會直接對來自靜葉片的氣流產(chǎn)生外推作用，此外，42°夾角相對于LV的垂直角，其平衡動葉片間氣流離心作用的效果也受到削弱，這必然導(dǎo)致Leoshce動葉片方案阻力最高，且轉(zhuǎn)子速度越高，高阻現(xiàn)象越明顯。

LV動葉片方案計算阻力低，同工業(yè)應(yīng)用實踐中LV選粉機相比TRM、Leosche、L型輥磨選粉機阻力低的現(xiàn)象吻合，這說明本文所采用的選粉機流場數(shù)值研究方法合理，基本能夠正確反映工業(yè)生產(chǎn)實際。

5 結(jié)論

（1）采用滑移網(wǎng)格技術(shù)成功解決了選粉機流場研究中轉(zhuǎn)子轉(zhuǎn)動問題，改善了目前國內(nèi)外普遍采用的MRF方法處理轉(zhuǎn)子轉(zhuǎn)動過程研究的不足，提高了選粉機流場研究的精度。

（2）LV動葉片結(jié)構(gòu)相對于TRM、Leosche、L型選粉機動葉片方案阻力性能最低，且低阻優(yōu)勢隨轉(zhuǎn)子速度的增加而增加。

[1]黃億輝,姜大志.轉(zhuǎn)子底盤結(jié)構(gòu)對O-SEPA選粉機流場特性的影響[J].裝備制造技術(shù),2011,(2):10～12.

[2]童聰,李雙躍,綦海軍,任朝富.立磨選粉機葉片結(jié)構(gòu)對分級區(qū)速度場影響分析[J].化工進展,2012,Vol.31(4):778～783.

[3]童聰,李雙躍,綦海軍,任朝富.立磨選粉機葉片參數(shù)的分析與優(yōu)化設(shè)計[J].過程工程學(xué)報,2012,Vol.12(1):14～18.

[4]綦海軍,李雙躍,任朝富,李庭婷.立磨選粉機導(dǎo)流圈的數(shù)值模擬與分析[J].浙江工業(yè)大學(xué)學(xué)報，2012，Vol.40(1)：70～74.

[5]陳杰來，姜大志，黃億輝.O-Sepa選粉機轉(zhuǎn)子結(jié)構(gòu)對流場特性的影響[J].中國粉體技術(shù)[J]，2011，Vol.17(6)：38～41.

[6]黃億輝,姜大志.O-Sepa選粉機內(nèi)三維流場的數(shù)值模擬[J].江南大學(xué)學(xué)報(自然科學(xué)版),2010,Vol.9(6):666～670.

[7]范增曉,姜大志,黃億輝.O-SEPA選粉機內(nèi)部流場數(shù)值模擬的方法選擇[J].江南大學(xué)學(xué)報(自然科學(xué)版),2011,Vol.10(5):568～572.

[8]Karunakumari L，Eswaraiah C，Jayanti S，et al.Experimental and numerical study of a rotating wheel air classifier［J］. AIChE Journal，2005，vol.51(3):780.

[9]D.Choudhury.Introduction to the Renormalization Group Method and Turbulence Modeling[G].Fluent Inc.Technical Memorandum TM-107,1993.

[10]R.I.Issa.Solution of Implicitly Discretized Fluid Flow Equations by Operator Splitting[J].J.Comput.Phys.,62:40-65,1986.

[11]S.V.Patankar.Numerical Heat Transfer and Fluid Flow[J].Hemisphere,Washington,DC,1980.

[12]B.Van Leer.Toward the Ultimate Concervative Difference Scheme.IV.A Second Order Sequel to Godunov's Method[J].Journal of Computational Physics,32:101-136,1979.■