顧宗幸 李浙昆 蔡培良 王勝枝 楊劍鋒 陸菲

摘 要：為研究離心式分離腔體內(nèi)部氣流流動軌跡以及梗絲顆粒的運動規(guī)律，進一步深入研究梗絲分離效果，基于ANSYS建立流體域模型，利用流體仿真軟件Fluent，采用RNG湍流模型模擬腔體內(nèi)氣流軌跡，拉格朗日離散相模型（DPM）模擬顆粒相運動規(guī)律，得到腔體內(nèi)氣流軌跡、壓力和顆粒相運動軌跡云圖，并在入口處設(shè)定不同大小的氣流速度，當入射速度為6m/s時，氣流流動軌跡發(fā)生改變，達到8m/s時氣流軌跡比較理想，且顆粒運動軌跡正確。此時梗絲在離心分離腔體內(nèi)可以得到有效分離，梗中含絲率降低到2.55%左右，煙絲利用率進一步提高。

關(guān)鍵詞：離心腔體;顆粒相;數(shù)值模擬;入射速度

DOI：10. 11907/rjdk. 191359 開放科學(xué)（資源服務(wù)）標識碼（OSID）：

中圖分類號：TP319文獻標識碼：A 文章編號：1672-7800（2019）008-0159-03

Simulation of Flow Field in Centrifugal Separation Chamber Based on Fluent

GU Zong-xing1，LI Zhe-kun1，CAI Pei-liang2，WANG Sheng-zhi2，YANG Jian-feng2，LU Fei1

（1. Faculty of Mechanical and Electrical Engineering， Kunming University of Science and Technology， Kunming 650500， China;

2. Qujing Cigarette Factory， HongyunHonghe Tobacco （Group） Co.， Ltd.， Qujing 655001， China）

Abstract：In order to analyze the flow trajectory of the airflow inside the centrifugal separation chamber and the movement law of the stem-splitting particles， the separation effect of the stems was further studied. Based on ANSYS， the fluid domain model was established. Using the fluid simulation software Fluent， the RNG turbulence model was used to simulate the airflow trajectory in the cavity， and the Lagrangian discrete phase model （DPM） was used to simulate the motion law of the particle phase， then the air velocity trajectory ， pressure in cavity and motion trajectory of particle phase were obtained. The airflow velocity of different sizes was set at the entrance. When the incident velocity was 6m/s， the airflow trajectory changed. When the velocity reached 8m/s， the airflow trajectory was ideal and the particle motion trajectory was correct. The stem wire can be effectively separated in the centrifugal separation chamber， and the silk content in the stem was reduced to about 2.55%， which improved the utilization rate of the tobacco.

Key Words：centrifugal cavity;particle phase;numerical simulation;incidence speed

基金項目：云南中煙工業(yè)有限責(zé)任公司科技項目（2017GY05）

作者簡介：顧宗幸（1993-），男，昆明理工大學(xué)機電工程學(xué)院碩士研究生，研究方向為數(shù)字化設(shè)計與制造;李浙昆（1957-），男，博士，昆明理工大學(xué)機電工程學(xué)院教授、博士生導(dǎo)師，研究方向為機電系統(tǒng)設(shè)計及控制技術(shù)等。本文通訊作者：李浙昆。

0 引言

目前，大多數(shù)煙草企業(yè)采用真空浮選的方式對成品煙絲中多余梗簽進行分離[1-3]，將風(fēng)送進入梗簽分離裝置進行分離后，還會有一大部分結(jié)團煙絲被分離進梗簽收集箱內(nèi)[4-5]，為了使其在梗簽中分離、提高煙絲利用率，許多煙草企業(yè)對此進行了研究。劉程[6]在PROTOS70的基礎(chǔ)上設(shè)計一套新型梗絲分選裝置，該裝置加長梗絲提升通道，使梗絲一邊上升一邊左右擺動，使塊狀梗絲松散，調(diào)節(jié)風(fēng)選氣壓使梗絲分離;廖艷培[7]針對YJ17卷煙機設(shè)計一套梗絲分離3次風(fēng)選裝置，利用梗絲懸浮速度差異通過調(diào)壓裝置、吹風(fēng)裝置使梗絲分離;劉瀾波[8]設(shè)計了一套煙梗打散型漂浮室裝置，將其安裝在“之”字形漂浮室垂直方向上，使煙絲與煙梗徹底分離，提高梗絲二次分離效果;張皓天[9]通過模擬仿真方形和圓形懸浮腔，得到的方形懸浮腔分離效果較好;蔡培良[10]對梗簽二次分離裝置工作原理、主要結(jié)構(gòu)參數(shù)以及梗絲受力進行了分析研究，但并未對整個離心分離腔體內(nèi)的氣流以及梗絲顆粒運動軌跡進行模擬仿真。本文通過Fluent軟件對離心分離腔體內(nèi)的氣流及顆粒相運行軌跡進行模擬仿真，同時研究不同入射速度下腔體內(nèi)速度變化，找出最優(yōu)入射速度，進一步探討梗絲經(jīng)過離心分離腔體后梗中含絲率。

1 離心腔體數(shù)值仿真

1.1 幾何建模與網(wǎng)格劃分

離心分離腔體利用離心分離的原理[11]，在正壓吹風(fēng)的作用下，梗絲混合物進入腔體內(nèi)，由于腔體結(jié)構(gòu)的改變使氣流發(fā)生改變，進而使重量較輕（小于氣流的升力）的煙絲、重量較重（大于氣流升力）的梗簽達到分離。其結(jié)構(gòu)模型如圖1所示。

Fluent主要用于流場內(nèi)的數(shù)值模擬[12]，流體域模型建好后利用Workbench中的Mesh模塊對已經(jīng)建好的流域進行前處理工作，首先選擇右側(cè)底部圓形面作為入口，標為inlet，然后選擇左側(cè)上部矩形平面為出口1，標為outlet1，左側(cè)底部圓形平面為出口2，標為 outlet2，其余面定義為壁面，標記為wall，如圖2所示;然后將物理場選項設(shè)置為CFD，求解器設(shè)置選項為FLUENT，Relevance Center選項設(shè)置為fine，Smoothing設(shè)置為high，即將其劃分為分布較密、質(zhì)量較高的網(wǎng)格[13-14]。點擊 Generate Mesh 按鈕可生成相應(yīng)網(wǎng)格。

圖1 離心分離腔體? ? ? ? ? ? ? ? ? ? ? 圖2 進出口及壁面定義

1.2 計算模型構(gòu)建

氣相與顆粒相分離過程屬于兩相流，在腔體內(nèi)部發(fā)生。固體顆粒主要在內(nèi)部氣流場的作用下運動，在分析顆粒相分離前對離心分離裝置內(nèi)的氣流場研究很重要，有助于進一步研究離心分離裝置性能。

離心腔內(nèi)部流體遵循質(zhì)量守恒和動量守恒規(guī)律。本模型數(shù)值仿真中，利用拉格朗日方程計算固相顆粒運動軌跡，研究對象采用RNG k-ε湍流模型，該模型與標準k-ε湍流相似，但RNG k-ε湍流模型應(yīng)用范圍廣、結(jié)果可靠[15]。這是因為RNG k-ε方程中出現(xiàn)了新函數(shù)項，得到的湍動能和耗散率方程為[16-19]：

[？（pk）？τ+？（pkui）？xi=？？xiakueff+？k？xj+Gk+pε？（pk）？τ+]

[？（pkui）？xi=？？xjazueff+？k？xj+？igεkGk-？2gpg2k] （1）

式（1）中，[ak]、[az]、[？2g]、[η]和[β]為經(jīng)驗常數(shù)，均采用Fluent默認數(shù)值。

采用DPM模型跟蹤固相顆粒運動，固相顆粒作用力平衡方程在笛卡爾坐標系下的形式為：

[dupdt=FD（μ-μP）+gx（pp-p）px+Fx]? ? （2）

式（2）中，[FD（μ-μP）]為顆粒單位質(zhì)量曳力，[μ]為液相速度，[up]為顆粒速度，[μP]為流體動力黏度，[p]為流體密度，[pp]為顆粒密度，[gx]為顆粒受到的其它加速度，[Fx]為其它作用力。

1.3 材料物性、邊界條件與初始條件

（1）離心分離腔體的選材為1Cr18Ni9Ti，而Fluent默認固體材料為鋁，修改其相應(yīng)屬性。離心分離腔仿真中，在腔體入口通入氣體;固體顆粒材料設(shè)置為wood，直徑取0.335～3.815mm。

（2）固體顆粒在氣體推動作用下進入離心腔體內(nèi)，由于結(jié)構(gòu)原因?qū)崿F(xiàn)分離。因此設(shè)定氣體及固體顆粒速度相同，分別取5～8m/s之間速度進行模擬仿真。

（3）當模擬梗簽顆粒軌跡時，選擇面射流源將固體顆粒從進口面均勻射入腔體內(nèi)，同時將集梗斗底面設(shè)為捕集面（trap），而集絲斗底面設(shè)為反射面（reflect）;當模擬煙絲顆粒軌跡時，將集梗斗底面設(shè)為反射面（reflect），而集絲斗底面設(shè)為捕集面（trap），出口永遠為逃逸面（escape）。

（4）求解，選用穩(wěn)態(tài)求解：利用二階迎風(fēng)格式離散的方法計算連續(xù)性方程、動量方程、湍動能方程和湍動能耗方程，殘差均下降至1×10-4;采用欠松弛迭代方法避免非線性問題的發(fā)散[20-21]。

2 仿真結(jié)果及分析

2.1 腔體內(nèi)速度分析

為了獲得最佳工作風(fēng)速使梗絲得到有效分離，利用數(shù)值模擬的方法可以直接觀察到整個離心腔體內(nèi)氣相流場的速度分布規(guī)律。取不同大小的速度，得到腔體內(nèi)氣流軌跡及速度變化，如圖3所示。

從圖3中可以看出，進口風(fēng)速在5m/s時氣流軌跡主要流向左下通道，在左上通道流出比較稀疏，這樣當梗絲進入腔體后，大部分梗絲會從左下通道流出，無法實現(xiàn)梗絲進一步分離。當進口速度進一步增大到6m/s時，氣流軌跡發(fā)生變化，主要從左上通道流出，左下通道氣流變得稀疏。進一步增大進口速度至6.5、7.0、7.5、8.0m/s時，發(fā)現(xiàn)氣流也主要從左上通道流出，左下通道很少，相互比較發(fā)現(xiàn)當進口速度為8m/s時氣流軌跡較好，在分離口處不會有漩渦出現(xiàn)，且在分口處速度也降低到4.1m/s左右，這是通道截面積增大引起的，如此會使煙絲與梗簽受力發(fā)生變化，煙絲不能克服離心力，繼續(xù)在氣流作用下從左上通道流出，而梗簽本身比煙絲沉，由于重力大于離心力便會下落，從左下通道流出，從而實現(xiàn)梗簽與煙絲分離。

（a）5.0m/s? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（b）5.5m/s

（c）6.0m/s? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （d）6.5m/s

（e）7.0m/s? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （f）8.0m/s

圖3 不同進口速度下腔體內(nèi)速度軌跡

2.2 腔體內(nèi)壓力分析

經(jīng)過后處理得到整個離心分離腔內(nèi)及z=0.05m時的壓力云圖，在腔體右半部分壓力逐漸減小，左半部分壓力基本保持不變，如圖4所示。

圖4 離心分離腔體內(nèi)壓力云圖

2.3 腔體內(nèi)顆粒軌跡

梗絲顆粒運動軌跡跟蹤是進行梗絲分離裝置性能研究的重要環(huán)節(jié)。由于該分離裝置結(jié)構(gòu)設(shè)計相對簡單，顆粒在分離裝置內(nèi)的運動情況不是很復(fù)雜，考慮除受到氣流曳力外，還受到顆粒與顆粒、顆粒與壁面之間的碰撞。本文主要跟蹤顆粒粒徑在0.335～3.815mm梗簽顆粒與煙絲顆粒的軌跡，在入口處釋放一定數(shù)量的顆粒，找到顆粒運動規(guī)律，如圖5-6所示。

圖5 煙絲顆粒軌跡? ? ? ? ? ? ? ? ? ? ? ? 圖6 梗簽顆粒軌跡

由圖5、圖6可以看出，煙絲及梗簽顆粒受湍流氣流影響，其隨機性較大，相互之間與壁面之間均會發(fā)生碰撞，但煙絲及梗簽顆粒的相對運動正確，在圖5中可以看出大部分煙絲在氣流作用下直接沿外壁面流出，只有極少量流入左下通道。而在圖6中梗簽在到達分口處時，由于氣流速度變小以及自身重密度大，會有向下運動趨勢，梗簽從左下通道流出，不過也有一些細碎的梗簽從左上通道流出，但數(shù)量極少，這是不可避免的。

為了研究分離效率，通過多次模擬在入口處加入不同數(shù)量的煙絲及梗簽顆粒，得到梗簽顆粒平均捕集效97.83%，煙絲顆粒捕集效率為76.14%。若按照原始梗中含絲率為20.27%[～]22.56%計算最后收集到的梗中含絲率，設(shè)梗絲混合物的數(shù)量為[γ]個，則按照公式（3）可計算出經(jīng)過離心分離腔體后梗中含絲率：

[γ（1-97.83%）（20.27%～22.56%）γ1-97.83%77.44%～79.73%+γ76.14%（20.27%～22.56%）][=（2.51%～2.59%）] （3）

從而利用該分離裝置使煙絲及梗簽得到進一步分離，梗中含絲率由原先的20.27%[～]22.56%降低到2.51%～2.59%，可極大提高經(jīng)濟效益。

4 結(jié)語

本文采用離散相模型對離心分離腔體內(nèi)的流場及顆粒軌跡進行仿真分析，進口通入的風(fēng)速大小不同，當進口風(fēng)速在5.0～6.0m/s時氣流軌跡主要流向左下通道，而達到6.5m/s時氣流軌跡發(fā)生改變，主要流向左上通道，當進口氣流速度為8m/s時達到最佳。同時在入口處入射一定量的煙絲和梗簽顆粒，大部分煙絲會從左上通道流出，多數(shù)梗簽會直接落入左下通道流出，煙絲及梗簽通過該分離裝置后可得到有效分離，梗中含絲率大幅降低。通過該數(shù)值模擬分析可為后續(xù)梗絲分離裝置研究與開發(fā)提供支持。

參考文獻：

[1] 曾靜，李斌，馮志斌，等. 卷煙機剔除梗簽物中含絲量的檢測[J]. 煙草科技，2012（8）：5-7，11.

[2] BUFFO M M， CORREA L J， ESPERANCA M N， et al. Influence of dual-impeller type and configuration on oxygen transfer， power consumption， and shear rate in a stirred tank bioreactor[J].? Biochemical Engineering Journal，2016，114：130-139.

[3] 陳鍇，周義民. ZJ 116卷接機組梗簽含絲率分析與改進[J]. 機械管理開發(fā)，2018（3）：13-14+68.

[4] 白雪，蔣曉偉，蔡培良，等. ZJ1卷煙機在線梗簽分離系統(tǒng)的設(shè)計[J]. 包裝工程，2018，（21）：32-36.

[5] 王勝枝，肖偉忠，蔡培良，等. 卷煙生產(chǎn)中梗簽二次分離裝置的設(shè)計[J]. 機械制造，2018，（11）：64-66+70.

[6] 劉程. PROTOS70卷煙機新型梗絲分選裝置的設(shè)計[J]. 裝備制造技術(shù)，2013，（4）：80-81.

[7] 廖艷培， 陳志堅. YJ17型卷煙機梗絲分離三次風(fēng)選裝置的設(shè)計[J]. 裝備制造技術(shù)，2016（3）：244-246.

[8] 劉瀾波，韓艷軍，陳恩，等.? PROTOS70 卷接機組煙梗打散型漂浮室的設(shè)計[J]. 煙草科技，2016，49（5）：90-93.

[9] 張皓天，李浙昆，蔡培良，等. 外置式卷接機組梗簽分離裝置的設(shè)計[J]. 煙草科技，2018，51（10）：97-103.

[10] 蔡培良，顧宗幸，王勝枝，等. 梗簽二次分離裝置離心分離部分的結(jié)構(gòu)分析研究[J]. 價值工程，2018，37（28）：164-166.

[11] 楊德武，湯明東，馮穎，張賢彬. 離心三相分離過程數(shù)值模擬研究[J]. 過濾與分離，2007（2）：21-24.

[12] NICOLET C. Hydro-acoustic modelling and numerical simulation of unsteady operation of hydroelectric systems[C]. Lausanne： Ecole Polytech-nique Federale de Lausanne，2007.

[13] 張樹女，馮立斌，徐連奎. 基于FLUENT的長短葉片水輪機數(shù)值仿真[J]. 軟件導(dǎo)刊，2018，17（5）： 121-123.

[14] 趙雪寧. 基于FLUENT軟件的pVTt法標準容器流場仿真研究[J]. 軟件導(dǎo)刊，2018，17（2）： 154-157.

[15] 石鑫，向陽，文利雄，等. 基于離散相模型的旋轉(zhuǎn)填充床內(nèi)的流場分析[J]. 高校化學(xué)工程學(xué)報，2012（3）：388-394.

[16] TSAY J，F(xiàn)OX R D，OZKAN H E，et al. Evaluation of a pneumatic-shield spraying system by CFD simulation[J]. Transactions of the ASAE，2002，45（1）：47-54.

[17] DELELE M A，MOOR A D E，SONCK B，et al. Modelling and validation of the air flow generated by a cross flow air sprayer as affected by travel speed and fan speed[J]. Biosystems Engineering，2005，92（2）：165-174.

[18] ADRIAN R J. Particle-image techniques for experimental fluid mechanics[J]. Annual Review Fluid Mechanics，1991（23）：261-304.

[19] ADRIAN R J.Twenty years of particle image velocity[J]. Experiments in Fluids，2005，39（2）：159-169.

[20] 范強，李永杰. 旋轉(zhuǎn)旋流離心機轉(zhuǎn)鼓內(nèi)流體動力學(xué)研究[J]. 機械研究與應(yīng)用，2009（4）：47-49.

[21] 孫國祥，汪小旵，丁為民，等. 基于CFD離散相模型霧滴沉積特性的模擬分析[J]. 農(nóng)業(yè)工程學(xué)報，2012（6）：13-19.

（責(zé)任編輯：江 艷）