陳旭東,胡國玉,2,*,趙騰飛,李忠新,周建平,2,

基于CFD-DEM耦合的核桃殼仁混合物負(fù)壓風(fēng)選仿真研究

陳旭東1,胡國玉1,2,3*,趙騰飛1,李忠新3,周建平1,2,3

1. 新疆大學(xué)機械工程學(xué)院, 新疆 烏魯木齊 830047 2. 新疆維吾爾自治區(qū)農(nóng)牧機器人及智能裝備工程研究中心, 新疆 烏魯木齊 830047 3. 新疆大學(xué)葉城核桃產(chǎn)業(yè)研究院, 新疆 葉城縣 844900

在核桃精深加工過程中，風(fēng)選法可以實現(xiàn)核桃殼仁混合物料高效無損地分離。為了研究不同等級的核桃殼仁顆粒在氣流場中的分選效果，本文設(shè)計一種負(fù)壓垂直氣流風(fēng)選裝置，采用計算流體力學(xué)(CFD)與離散單元法(DEM)耦合的方法對風(fēng)選過程進行仿真模擬，引用清選率與損失率量化風(fēng)選效果。結(jié)果表明：當(dāng)進料量一定，進風(fēng)口風(fēng)速19 m/s時，二分之一殼仁混合物料清選率達(dá)到96.5%，損失率僅為0.41%，分選效果最佳；進風(fēng)口風(fēng)速為17 m/s時，四分之一殼仁混合物料清選率為94.8%，損失率為0.8%，分選效果較為理想；進風(fēng)口風(fēng)速為16 m/s時，八分之一殼仁混合物物料清選率為92%，損失率為1.22%，風(fēng)選效果較為合理。研究表明在垂直風(fēng)道里，核桃殼的流動速度圍繞某一固定值上下波動，核桃仁沉降速度隨著顆粒濃度增加而減小并逐漸趨于穩(wěn)定。

核桃; 負(fù)壓風(fēng)選; 仿真分析

核桃精深加工技術(shù)是制約核桃經(jīng)濟發(fā)展的主要因素，且破殼后殼仁混合物料的分離是核桃精深加工的重要環(huán)節(jié)[1]。核桃經(jīng)過破殼，物料分級，得到二分之一殼仁混合物料、四分之一殼仁混合物料、八分之一殼仁混合物料以及碎料等。設(shè)計殼仁分離裝置，實現(xiàn)對不同等級核桃殼仁混合物料的分離，成為一個亟待解決的問題。

關(guān)于核桃殼仁分離，國內(nèi)董遠(yuǎn)德等[2]設(shè)計了核桃殼仁分離裝置，并進行了水平正交試驗，結(jié)果表明風(fēng)量大小、風(fēng)腔長度和喂入量對分選率有顯著影響。曹成茂等[3]通過控制風(fēng)速均勻性與風(fēng)速精度研究了山核桃物料的風(fēng)選過程，得出風(fēng)速、殼與仁含水率、迎風(fēng)面容量比組合參數(shù)下最優(yōu)的清選率和最低的誤選率。國外Nahal AM等[4]將兩個品種的核桃殼仁物料分為殼、仁、殼仁鑲嵌物三類，通過研究分離時間、進料量和風(fēng)速對殼仁分離效果的影響，分析了不同核桃殼仁混合物的末端速度，最終得出物料最佳分離速度。Jin FH[5]等利用核桃殼仁不同的光學(xué)圖像，采用機器視覺的方法實現(xiàn)殼仁分離，但此種方法不適用大規(guī)模生產(chǎn)，不宜在我國推廣。通過總結(jié)國內(nèi)外學(xué)者對核桃殼仁分離技術(shù)的研究成果以及前往葉城縣核桃產(chǎn)業(yè)基地實地調(diào)研，實現(xiàn)核桃殼仁分離，必須結(jié)合核桃自身物理特性研究，利用核桃殼仁空氣動力學(xué)特性的差異，采用風(fēng)選法可以有效實現(xiàn)分離。

CFD-DEM耦合方法在農(nóng)業(yè)領(lǐng)域已廣泛應(yīng)用。王立軍等[6]利用CFD-DEM耦合方法模擬了不同的風(fēng)速下玉米脫出物的運動過程與規(guī)律。丁賀賀等[7]設(shè)計了垂直風(fēng)選機，采用DEM-CFD耦合方法研究了谷粒-莖稈的分選規(guī)律，分析了風(fēng)選參數(shù)對風(fēng)選效果的影響。

本文運用CFD-DEM耦合的方法，模擬研究了垂直風(fēng)道內(nèi)負(fù)壓風(fēng)選過程中不同等級核桃殼仁混合物的分離過程，并分析了不同風(fēng)速下殼仁的分離效果，揭示了核桃殼仁在垂直風(fēng)道里的運動規(guī)律，為核桃殼仁分離技術(shù)提供分選方法與理論基礎(chǔ)。

1 風(fēng)選原理

風(fēng)選法是利用物料與雜質(zhì)在分選室內(nèi)懸浮速度的差異實現(xiàn)分離。依據(jù)風(fēng)機工作時分離室內(nèi)風(fēng)壓情況，分正壓吹風(fēng)式分離法與負(fù)壓吸風(fēng)式分離法[8]。本文采用負(fù)壓垂直氣流分離法實現(xiàn)對核桃殼仁混合物料的分離，基本原理是氣流將質(zhì)量較輕的核桃殼吸入到卸料箱中收集，核桃仁依靠自身的重力從分離室內(nèi)流出收集（圖1）。在垂直風(fēng)道里，核桃殼仁顆粒主要受到風(fēng)力與重力的作用。

圖1 核桃殼仁球體化受力圖

在同一流體介質(zhì)中，對于不規(guī)則形狀的核桃混合物料與其當(dāng)量球體的體積和密度保持不變，得到修正的殼仁物料的懸浮速度[9]，如式(1)所示。依據(jù)不同核桃殼仁等級物理參數(shù)，可計算出不同形態(tài)殼仁的懸浮速度，當(dāng)同等級的核桃殼仁在分離室內(nèi)受到的風(fēng)速大于殼的懸浮速度時殼仁完成分離。

式中：υ是料懸浮速度（m/s）；k是球形系數(shù)；d是為當(dāng)量球直徑，m；ρ是物料密度，kg/m3；-是空氣密（kg/m3）；是阻力系數(shù)；是重力加速度（9.8 m/s2）。

2 數(shù)學(xué)模型

2.1 流體動力學(xué)方程

基于CFD-DEM耦合方法，采用標(biāo)準(zhǔn)的湍流模型對核桃殼仁分離過程進行數(shù)值模擬。其流體相連續(xù)方程和動量守恒方程分別為[10]：

式中：ε是氣體體積分?jǐn)?shù)；ρ是氣體密度（kg/m3）；υ氣體流速（m/s）；是時間（s）；是氣體微元體上的的壓力（Pa）；重力加速度（m/s2）；R是動量源項；是氣體動力黏度（Pa·s）。

動力源項R表示在網(wǎng)格單元格內(nèi)氣流所受阻力的總和，其方程式為[11]：

式中：F是顆粒對氣體的阻力作用（N）；V是網(wǎng)格單元的體積（m3）。

由于氣流對核桃殼仁顆粒的主要作用力為曵力，選擇Freestream阻力模型，其計算方程式為：

式中：C是單個顆粒曳力系數(shù)；是顆粒投影面積，m2；是氣體對顆粒的相對流動速度，m/s。

2.2 顆粒接觸碰撞模型

在EDEM與FLUENT耦合過程中，顆粒間的運動行為應(yīng)考慮顆粒與顆粒，顆粒與容器壁面之間相互碰撞過程。本文采用Hertz-Mindlin無粘結(jié)接觸模型，顆粒間碰撞采用軟球面接觸，分離室壁面視為半徑無窮大的球體。顆粒間的法向力F、法向阻尼力F切向力F切向阻尼力F[12]分別為：

式中：是等效彈性模量（Pa）；是等效顆粒半（m）；是法向重疊量（m）；是恢復(fù)系數(shù)；S是法向剛度（N/m）；m是等效質(zhì)量（kg）；u是相對速度的法向速度（m/s）；S是切向剛度（N/m）；是切向重疊量（m）；ul是相對速度的切向速度（m/s）。

2.3 氣固兩相相對運動方程

DEM-CFD耦合過程中，固體相顆粒采用力學(xué)平衡方程進行求解[13]。

式中：m是顆粒的質(zhì)量（kg）；I是顆粒的轉(zhuǎn)動慣量（kg·m2）；V是顆粒的速度（m/s）；ω是顆粒的角速度（rad／s）；是重力加速度（m/s2）；是吹風(fēng)傾角；是受風(fēng)面積（m2）；ρ是氣體密度（kg／m3）；υ是氣體流速（m/s）；F是法向力（Ｎ）；F是切向力（Ｎ）；T是切向力矩（Ｎ·ｍ）；T是摩擦力矩（Ｎ·ｍ）。

3 風(fēng)選室及物料幾何模型

3.1 負(fù)壓風(fēng)選室?guī)缀文Ｐ团c原理

本文模型采用垂直氣流負(fù)壓風(fēng)選裝置，如圖2所示，運用Solidworks軟件建立，其特征在于：進料口為120 mm×400 mm，保證物料是單層進入分離腔，與分離腔的夾角為45°；進風(fēng)口的幾何尺寸為400 mm×120 mm，距離進料口560 mm處；分離腔的幾何尺寸為1100 mm×400 mm×120 mm，上端設(shè)計物料加速流動區(qū)域，確保殼與仁完全分開后核桃殼加速進入卸料箱內(nèi)；卸料箱上端面為820 mm×400 mm的矩形，下端面為400 mm×300 mm的矩形，高度為1200 mm；出風(fēng)口直徑為Φ250 mm的圓形敞口，卸料箱下端面安裝有星型卸料器，確保工作時產(chǎn)生負(fù)壓環(huán)境。

1.進風(fēng)口 2.進料口 3.分離腔 4.加速流動區(qū) 5.出風(fēng)口 6.卸料箱 7.星型閉風(fēng)器

核桃殼仁混合物料從進料口落入分離腔內(nèi)，由于核桃殼與仁之間懸浮速度的差異，核桃仁所受的浮力小于自身重力，從分離腔內(nèi)落下，核桃殼所受浮力大于自身重力，在風(fēng)力的作用下進入加入流動區(qū)，流入卸料箱，由于卸料箱空間面積的急劇擴大，所受浮力減小，落入卸料箱底部，由星型卸料器收集。上述就是整個負(fù)壓風(fēng)選裝置具體工作原理。

3.2 核桃殼仁混合物料模型

核桃破殼后產(chǎn)生的物料經(jīng)過分級得到不同規(guī)格的殼仁混合物，對核桃殼仁從棱長、長徑、厚度三個方向進行測量，測量尺寸如表1所示。依據(jù)尺寸將核桃殼仁顆粒分為四個等級，分別是1/2殼仁混合物、1/4殼仁混合物、1/8殼仁混合物、1/16殼仁混合物等。1/16殼仁懸浮速度差異較小，使用風(fēng)選法效果不佳[14]，本文不對其進行研究分析。

圖1 不同等級核桃殼仁平均尺寸

依據(jù)核桃殼仁的實際尺寸，首先利用Solidworks軟件進行顆粒輪廓建模，然后將模型導(dǎo)入EDEM軟件，采用微元球形顆粒填充，得到如圖3的幾何模型。

(a) 二分之一殼仁 (b) 四分之一殼仁 (c) 八分之一殼仁

4 物料力學(xué)特性及仿真參數(shù)設(shè)置

4.1 物料力學(xué)特性及接觸系數(shù)

EDEM軟件中需要設(shè)置各個材料的物理特性參數(shù)以及材料間接觸系數(shù)，根據(jù)曹成茂等人前期的研究成果[3]，各個材料的物理特性參數(shù)以及材料間接觸系數(shù)如表2和表3所示。

表2 材料物理參數(shù)

圖3 材料接觸系數(shù)

4.2 仿真參數(shù)設(shè)置

EDEM與FLUENT耦合過程是一個瞬態(tài)雙向數(shù)據(jù)傳遞的過程，DEM方法能夠準(zhǔn)確地分析顆粒的形狀、級配等因素對顆粒的運動的情況，結(jié)合計算CFD方法，能夠準(zhǔn)確地分析顆粒相與流體相間的相互作用。其中EDEM與FLUENT耦合過程如下：

（1）首先在Solidworks軟件中對負(fù)壓風(fēng)選裝置進行三維建模，然后將該模型導(dǎo)入ANSYS Workbench的Geometry模塊中，利用流體填充功能對內(nèi)部氣流場進行建模，使用mesh功能對流體域進行網(wǎng)格劃分，啟動自適應(yīng)尺寸調(diào)整，設(shè)置分辨率為5，插入補丁適形法，采用四面體網(wǎng)格，并對進風(fēng)口與出風(fēng)口進行定義，如圖4所示。

圖4 流體域網(wǎng)格劃分

（2）EDEM中導(dǎo)入上述圖2的簡易模型，然后將上述圖3建立的同等級的顆粒模型導(dǎo)入Bulk Material中，并將上述表2和表3中的材料物理參數(shù)及材料接觸系數(shù)輸入到EDEM中，設(shè)置顆粒工廠為390 mm×110 mm的虛擬矩形平面，能同時產(chǎn)生核桃殼與核桃仁兩種物料，由于核桃是不規(guī)則的球體，且大小不一，同等級的核桃殼仁也有一定的尺寸差距，為了更加接近實際尺寸，對輸入的模型在體積設(shè)置上采用用戶自定義的方法，物料的尺寸在物料模型尺寸0.7～1.2倍之間浮動。綜合考慮核桃殼與仁的質(zhì)量占比與計算機的處理能力，設(shè)置物料進料量為540 kg/h，殼仁占比為1:2，核桃殼的生成速率為0.05 kg/s，核桃仁的生成速率為0.1 kg/s。EDEM中設(shè)置瑞麗時間步長為1e-5 s，每0.01 s保存一次數(shù)據(jù)，每個數(shù)據(jù)點迭代1×103次。

（3）選擇耦合接口為DDPM類型接口，考慮顆粒在流場中的所占體積分?jǐn)?shù)，在FLUENT中讀取耦合文件，采用k-epsilon模型。設(shè)定進風(fēng)口邊界條件為速度入口(velocity-inlet)，速度為仿真模擬值。出風(fēng)口與外界大氣壓相連，設(shè)定其邊界條件為壓力出口(pressure-outlet)，其他保持默認(rèn)。顆粒由入料口流入，設(shè)定邊界條件為速度入口，速度值為0 m/s。FLUENT中每100步自動保存數(shù)據(jù)，設(shè)置雷利時間步長為瑞麗時間步長的50倍，即5 e-4 s，總步長為10000步，總的仿真時間為5 s。

5 不同等級殼仁混合物仿真結(jié)果與分析

為了更直觀地研究不同等級的核桃殼仁混合物的風(fēng)選效果，采用清選率和損失率作為衡量標(biāo)準(zhǔn)[15]。如圖5所示，在核桃仁出口劃分出清選率1的統(tǒng)計區(qū)域Ⅰ，統(tǒng)計該區(qū)域仿真結(jié)束時物料總質(zhì)量和核桃仁的質(zhì)量m1；在核桃殼出口劃分出損失率2的統(tǒng)計區(qū)域Ⅱ，統(tǒng)計=0～5 s內(nèi)流出的核桃仁的質(zhì)量m2，則：

圖 5 核桃殼仁顆粒分離瞬態(tài)圖

圖6 分離室氣流場流速云圖

運用FLUENT軟件得到分離室對稱截面的內(nèi)部氣流場流速云圖，由圖6可知，氣流在分離腔內(nèi)分布比較均勻，在氣流加速流動區(qū)，由于空間面積的縮小，氣流流速明顯高于進風(fēng)口流速，有利于核桃殼加速流動，避免在彎形風(fēng)道里出現(xiàn)堆積；在卸料箱內(nèi)空間面積迅速增大，氣流沿著側(cè)壁流動，有利于核桃的迅速沉降，卸料箱中心區(qū)域的氣流流速接近0 m/s，實現(xiàn)卸料功能；在卸料箱底部風(fēng)速逐漸減小，氣流速度小于核桃殼的懸浮速度，核桃殼從底部沉降完成卸料。

5.1 二分之一殼仁混合物單因素仿真實驗

5.1.1 風(fēng)量風(fēng)速對分選效果影響為探究不同的風(fēng)速對殼仁分離效果的影響，對二分之一殼仁混合物在風(fēng)速18 m/s、19 m/s、20 m/s時分離過程進行模擬。如圖7所示，當(dāng)風(fēng)速由18 m/s增加到20 m/s時，二分之一仁的清選率由93.28%增加到98.52%，損失率由0.2%增加到0.65%。這是由于隨著風(fēng)速的增大，物料在分離室內(nèi)經(jīng)過風(fēng)力作用，一些懸浮速度相近的殼與仁被吹入卸料箱內(nèi)，導(dǎo)致清選率與損失率上升的原因。由于二分之一殼仁之間懸浮速度差異明顯，在提高風(fēng)速的情況下，能夠較好的完成風(fēng)選效果。綜合考慮，風(fēng)速為19 m/s時，清選效果最優(yōu)。

圖7 風(fēng)速對二分之一殼仁清選效果的影響

5.1.2 殼仁顆粒在分離腔內(nèi)運動速度由于分離室在方向?qū)ΨQ，在方向氣流對物料的速度影響較小，并且物料主要在方向運動，因此本文僅對殼仁顆粒的方向速度進行分析。

在進風(fēng)口不同風(fēng)速條件下，核桃殼仁顆粒在方向上的速度隨時間的變化曲線如圖8所示。在=0～5時間內(nèi)，二分之一殼在不同的風(fēng)速下呈現(xiàn)出不同的運動速度，總體在0.9 m/s～2.1 m/s的區(qū)間內(nèi)波動，這是因為二分之一殼是不規(guī)則的非球形顆粒，在流場中的運動姿態(tài)隨著流體的運動發(fā)生變化，導(dǎo)致物料的迎風(fēng)面積發(fā)生改變，從而使得物料的懸浮也在變化。

圖8 二分之一殼仁運動速度隨時間變化規(guī)律

Fig.8 Velocities of half shell and kernel moving with increase of time

當(dāng)風(fēng)速為18 m/s時，二分之一殼的運動速度在1.1 m/s上下波動，二分之一仁的沉降速度在0.9 m/s上下波動。當(dāng)風(fēng)速為19 m/s時，殼的運動速度在1.4 m/s上下波動，仁的運動速度在=0～5 s時間內(nèi)呈緩慢下降趨勢，從0.8 m/s降到0.5 m/s。當(dāng)風(fēng)速為20 m/s，核桃殼的運動速度在1.7 m/s上下波動，核桃仁的運動速度在=0～5 s時間內(nèi)從0.7 m/s下降到0.4 m/s。這是由于隨著時間變化，分離腔內(nèi)物料占流體區(qū)域的體積分?jǐn)?shù)增大，顆粒濃度隨之增大，仁與仁之間空隙減小，氣流作用在核桃仁顆粒上的實際風(fēng)速變大，導(dǎo)致核桃仁沉降速度變小。

5.1.3 殼仁顆粒在分離腔內(nèi)運動軌跡為研究不同等級殼仁顆粒在分離腔內(nèi)的分離過程，選取部分顆粒運動軌跡圖進行規(guī)律探究。

二分之一殼仁顆粒運動軌跡如圖9所示，在不同的風(fēng)速下，殼仁分離區(qū)域明顯不同。當(dāng)風(fēng)速為18 m/s時，殼仁分離區(qū)與核桃仁沉降區(qū)重合度高，在分離腔內(nèi)殼的運動速度最高達(dá)到3 m/s，仁的最低沉降速度為0.9 m/s；當(dāng)風(fēng)速為19 m/s時，分離區(qū)與沉降區(qū)重合部分減少，有少量核桃殼的運動軌跡出現(xiàn)在沉降區(qū)內(nèi)，殼的運動速度最高達(dá)到3.51 m/s，仁的最低沉降速度為0.59 m/s；當(dāng)風(fēng)速為20 m/s時，分離區(qū)與沉降區(qū)分布清晰，殼仁運動軌跡之間無交叉現(xiàn)象，殼的運動速度最高達(dá)到4 m/s，仁最低沉降速度僅為0.12 m/s。

圖9 二分之一殼仁運動軌跡

5.2 四分之一殼仁混合物單因素仿真實驗

5.2.1 風(fēng)量風(fēng)速對分選效果影響對四分之一殼仁混合物在風(fēng)速16 m/s、17 m/s、18 m/s時分離過程進行模擬。如圖10所示，當(dāng)風(fēng)速由16 m/s增加到18 m/s時，四分之一仁的清選率由88.45%增加到96.5%，損失率由0.54%增加到1.12%。隨著風(fēng)速的增加，四分之一仁的清選率增加顯著，但伴隨著損失率也明顯增加。綜合考慮，風(fēng)速為17 m/s時，清選率與損失率相對合理，清選效果比較理想。

圖10 風(fēng)速對四分之一殼仁清選效果的影響

5.2.2 殼仁顆粒在分離腔內(nèi)運動速度由圖11可知，在=0～5 s時間內(nèi)，四分之一殼在不同的風(fēng)速下總體在0.75 m/s～2.0 m/s的區(qū)間內(nèi)波動，四分之一仁在分離腔內(nèi)隨著顆粒濃度的增大，沉降速度逐漸減小。在風(fēng)速為16 m/s時，核桃殼的運動速度在1.0 m/s上下波動，核桃仁的運動速度在=0～5 s時間內(nèi)從0.9 m/s下降到0.5 m/s。當(dāng)風(fēng)速為17 m/s時，核桃殼的運動速度在1.25 m/s上下波動，核桃仁的沉降速度在=0～5 s時間內(nèi)呈緩慢下降趨勢，從0.7 m/s降到0.4 m/s。當(dāng)風(fēng)速為18 m/s，核桃殼的運動速度在1.5 m/s上下波動，核桃仁的運動速度在0～5s時間內(nèi)從0.5 m/s下降到0.2 m/s。隨著殼仁顆粒尺寸變小，殼仁質(zhì)量減小，分離所需風(fēng)速也相應(yīng)減小，但減小速率較為緩慢，這主要與核桃殼的迎風(fēng)面積有關(guān)。綜合考慮，運用風(fēng)選法對核桃殼仁分離的過程中，不僅要考慮核桃仁的清選率與損失率，也要分析殼仁顆粒在分離腔內(nèi)的運動速度，避免核桃仁因風(fēng)速過大在分離室內(nèi)下降速度較慢或者呈懸浮狀態(tài)。

圖11 四分之一殼仁運動速度隨時間變化曲線

Fig.11 Velocities of quarter shell and kernel moving with increase of time

5.2.3 殼仁顆粒在分離腔內(nèi)運動軌跡四分之一殼仁顆粒在不同風(fēng)速下運動軌跡如圖12所示，當(dāng)風(fēng)速為16 m/s時，殼仁分離區(qū)域與沉降區(qū)域重合度大，不利于殼仁顆粒的分離，分離腔內(nèi)殼的最高運動速度為3.15 m/s，仁的最低沉降速度為0.89 m/s；當(dāng)風(fēng)速為17 m/s時，分離區(qū)域與沉降區(qū)域重合度明顯減小，存在少部分核桃殼顆粒出現(xiàn)在沉降區(qū)域內(nèi)，殼最高運動速度達(dá)到3.96 m/s，仁最低沉降速度為0.52 m/s；當(dāng)風(fēng)速為18 m/s時，分離區(qū)域與沉降區(qū)域重合度較小，分離腔內(nèi)殼的運動速度達(dá)到4.57 m/s，仁的沉降速度為0.11 m/s。這主要是因為物料從進料口流入，到達(dá)分離腔時獲得一定的速度，在不同的風(fēng)速作用下，核桃殼速度減為0時的位移不同，導(dǎo)致分離區(qū)域位置與面積不同。

圖12 四分之一殼仁運動軌跡

5.3 八分之一殼仁混合物單因素仿真實驗

5.3.1 風(fēng)量風(fēng)速對分選效果影響對八分之一殼仁混合物在風(fēng)速15 m/s、16 m/s、17 m/s時分離過程進行模擬。如圖13所示，當(dāng)風(fēng)速由15 m/s增加到17 m/s時，八分之一仁的清選率由87.9%增加到96.75%，損失率由0.67%增加到1.12%。隨著風(fēng)速的增大，損失率較清選率增幅較快。這是由于核桃殼仁尺寸越小，殼與仁之間懸浮速度差縮小，分離難度增加。綜合考慮，風(fēng)速為16 m/s時，八分之一殼仁混合物清選效果比較合理。

圖13 風(fēng)速對八分之一殼仁清選效果的影響

5.3.2 殼仁顆粒在分離腔內(nèi)運動速度由圖14可知，在=0～5 s時間內(nèi)，八分之一殼在不同的風(fēng)速下總體在0.5 m/s～1.75 m/s的區(qū)間內(nèi)波動，八分之一仁在分離腔內(nèi)隨著顆粒濃度的逐漸增大，沉降速度緩慢減小并趨于穩(wěn)定。在=2 s以后，不同風(fēng)速下殼與仁的運動速度逐漸趨于穩(wěn)定，與二分之一殼仁和四分之一殼仁相比，八分之一殼仁在氣流場內(nèi)速度上下波動范圍較小，這是因為八分之一殼仁在流場中運動時迎風(fēng)面積的變化范圍小，對殼仁的懸浮速度影響較小。

當(dāng)風(fēng)速為15 m/s時，核桃殼的運動速度在0.75 m/s上下波動，核桃仁的運動速度在=0～5 s時間內(nèi)從0.9 m/s下降到0.5 m/s并逐漸趨于穩(wěn)定。當(dāng)風(fēng)速為16 m/s時，核桃殼的運動速度在0.9 m/s上下波動，核桃仁的運動速度在=0～5 s時間內(nèi)從0.8 m/s降到0.4 m/s并趨于穩(wěn)定。當(dāng)風(fēng)速為17 m/s，核桃殼的運動速度在1.1 m/s上下波動，核桃仁的運動速度在=0～5 s時間內(nèi)從0.7 m/s下降到0.2 m/s，核桃殼與核桃仁的相對運動速度差值較大，核桃仁沉降速度緩慢，接近懸浮狀態(tài)，這種情況會造成顆粒在分離腔內(nèi)大量堆積而無法完成分離。

圖14 八分之一殼仁運動速度隨時間變化規(guī)律

5.3.3 殼仁顆粒在分離腔內(nèi)的運動軌跡八分之一殼仁顆粒在不同風(fēng)速下運動軌跡如圖15所示，當(dāng)風(fēng)速為15 m/s時，分離區(qū)與沉降區(qū)重合度高，殼仁運動軌跡相互交錯，分離腔內(nèi)殼的最高運動速度為3.42 m/s，仁的最低沉降速度為0.53 m/s；當(dāng)風(fēng)速為16 m/s時，重合區(qū)域減少，殼仁運動軌跡分布合理，分離腔內(nèi)殼的最高運動速度為3.61 m/s，仁的最低沉降速度為0.28 m/s；當(dāng)風(fēng)速為17 m/s時，有少量的核桃殼運動軌跡出現(xiàn)在沉降區(qū)內(nèi)，分離腔內(nèi)殼的最高速度達(dá)到4.3 m/s，核桃仁最低沉降速度僅為0.16 m/s。

圖15 八分之一殼仁運動軌跡

6 結(jié)論

（1）利用CFD-DEM耦合方法模擬研究了核桃殼仁顆粒在垂直風(fēng)道里的運動及分離情況，基于EDEM建立了不同規(guī)格的核桃殼仁顆粒模型，準(zhǔn)確反映了殼仁顆粒的物理特性與空氣動力學(xué)特性，證明該方法可以準(zhǔn)確模擬核桃殼仁顆粒風(fēng)選效果，揭示核桃殼仁顆粒在流場中的運動本質(zhì)。該模擬方法可對核桃殼仁分離裝置的設(shè)計與優(yōu)化提供理論參考；

（2）本文設(shè)計的核桃殼仁風(fēng)選室模型，當(dāng)物料進料量一定，進風(fēng)口風(fēng)速為19 m/s時，二分之一殼仁混合物料清選率達(dá)到96.5%，損失率僅為0.41%，分選效果最佳；進風(fēng)口風(fēng)速為17 m/s時，四分之一殼仁混合物料清選率為94.8%，損失率為0.8%，分選效果較為理想；進風(fēng)口風(fēng)速為16 m/s時，八分之一殼仁混合物料清選率為92%，損失率為1.2%，風(fēng)選效果較為合理；

（3）隨著進風(fēng)口風(fēng)速增加，不同等級的殼仁顆粒清選率增加，損失率也隨之增加；隨著核桃殼仁顆粒尺寸的減小，殼仁顆粒懸浮速度重合度增加，分離難度變大；

（4）核桃殼在垂直分離腔內(nèi)的運動姿態(tài)時刻發(fā)生變化，且在不同的風(fēng)速下運動速度圍繞不同固定值上下波動，隨著核桃殼尺寸減小，波動幅度隨之減小；核桃仁在分離腔內(nèi)隨著顆粒濃度的增加，仁與仁之間空隙減小，氣流作用在顆粒上的實際風(fēng)速變大，沉降速度逐漸減小并趨于穩(wěn)定。

[1] 牛浩.核桃殼仁分離裝置設(shè)計及試驗研究[D].阿拉爾:塔里木大學(xué),2020

[2] 董遠(yuǎn)德,張學(xué)軍,史建新.核桃殼仁風(fēng)選機的設(shè)計與試驗研究[J].農(nóng)業(yè)機械,2012(24):156-158

[3] 曹成茂,羅坤,彭美樂,等.山核桃物料風(fēng)選機理與風(fēng)選性能研究[J].農(nóng)業(yè)機械學(xué)報,2019,50(9):105-112

[4] Nahal AM, Arabhosseini A, Kianmehr MH. Separation of shelled walnut particles using pneumatic method [J]. International Journal of Agricultural and Biological Engineering, 2013,6(3):88-93

[5] Jin FH, Qin L, Jiang L,. Novel separation method of black walnut meat from shell using invariant features and a supervised self-organizing map [J]. Journal of Food Engineering, 2008,88(1):78-85

[6] 王立軍,李洋,梁昌,等.貫流風(fēng)篩清選裝置內(nèi)玉米脫出物運動規(guī)律研究[J].農(nóng)業(yè)機械學(xué)報,2015,46(9):122-127

[7] 丁賀賀,石加聯(lián),馬學(xué)東,等.基于DEM-CFD耦合的谷物垂直氣流風(fēng)選仿真試驗[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報,2019,54(6):190-197

[8] 朱占江,李忠新,楊莉玲,等.核桃殼仁分離技術(shù)研究進展[J].食品工業(yè),2014,35(2):216-219

[9] 馬豪.核桃殼仁分離機的設(shè)計及試驗研究[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2011

[10] 蔣恩臣,孫占峰,潘志,等.基于CFD-DEM的收獲機分離室內(nèi)谷物運動模擬與試驗[J].農(nóng)業(yè)機械學(xué)報,2014,45(4):117-122

[11] 胡國明.顆粒系統(tǒng)的離散元素法分析仿真[M].武漢:武漢理工大學(xué)出版社,2010

[12] 王國強.離散單元法及其在EDEM上的實踐[M].西安:西北工業(yè)大學(xué)出版社,2010

[13] 楊航,劉芳,李加祥,等.基于CFD-DEM耦合的氣吸式三七播種器工作參數(shù)仿真分析[J].中國農(nóng)機化學(xué)報,2019,40(1):22-25

[14] 王維,王亞妮,盧軍黨,等.核桃破殼后各級物料懸浮速度的試驗研究[J].農(nóng)業(yè)技術(shù)與裝備,2020(11):5-6,8

[15] 陳廣慧,李軍,騰杰,等.基于EDEM-Fluent耦合的風(fēng)篩式清選裝置仿真及分析[J].中國農(nóng)機化學(xué)報,2020,41(5):69-73

Simulation Study on Negative Pressure Winnowing of Walnut Kernel Mixture Based on CFD-DEM Coupling

CHEN Xu-dong1, HU Guo-yu1,2,3*, ZHAO Teng-fei1, LI Zhong-xin3, ZHOU Jian-ping1,2,3

1./,830047,2.,830047,3.,844900,

A method of winnowing was proposed to achieve the efficient and non-destructive separation of the walnut shell and kernel mixed material in the deep processing of walnut. In order to study the cleaning effect of different scales of walnut shell and kernel mixed particles in the airflow field, a negative pressure and vertical airflow ventilation device was designed. The winding process was simulated by the method of coupling computational fluid dynamics (CFD) and discrete element method (DEM), and the air separation effect was quantified by the cleaning rate and loss rate.The results showed that with the constant of feed quantity,when the wind speed of the inlet was 19 m/s, the half walnut shell and kernel mixed material had the best cleaning effect that the cleaning rate was 96.5% and the loss rate was only 0.41%; when the air inlet wind speed was 17 m/s, the 1/4 walnut shell and kernel mixed material had more ideal cleaning effect that the cleaning rate was 94.8% and the loss rate was 0.8%; when the air inlet wind speed was 16 m/s, the 1/8 walnut shell and kernel mixed material had more reasonable cleaning effect that the cleaning rate was 92% and the loss rate was 1.22%. It showed that in the vertical air duct, the flow velocity of walnut shell fluctuated around a fixed value, and the sedimentation velocity of walnut kernel decreased with the increase of particle concentration and gradually tended to be stable.

Walnut; negative pressure air separation; simulation analysis

S233.5

A

1000-2324(2021)06-1017-011

2021-08-06

2021-08-21

葉城縣農(nóng)產(chǎn)品銷售“雙線九進”和滬喀品牌推廣項目(KSHSY20190901)

陳旭東(1996-),男,碩士研究生,主要從事智能農(nóng)業(yè)機械. E-mail:2393800511@qq.com

通訊作者:Author for correspondence. E-mail:xjhuguoyu@126.com