吴晶晶1，蔡 杰1，高亚东1，席剑飞1，顾中铸1，钟文镇2
(1. 南京师范大学 能源与机械工程学院，江苏 南京 210023；2. 济南大学 机械工程学院，山东 济南 250022)
DOI:10.13732/j.issn.1008-5548.2022.04.010
收稿日期： 2021-07-28， 修回日期：2022-04-02，在线出版时间：2022-06-22。
基金项目：国家自然科学基金项目，编号：51878356；江苏省自然科学基金项目，编号：BK20181385，BK20171477。
第一作者简介：吴晶晶(1997—)，男，硕士研究生，研究方向为颗粒流的CFD-DEM分析。E-mail：976604088@qq.com。
通信作者简介：
顾中铸(1963—)，男，教授，博士，研究方向为细颗粒物的沉降理论与技术。E-mail：guzhongzhu@njnu.edu.cn。
蔡杰(1978—)，男，教授，博士，研究方向为流场特性及其优化。E-mail：caijie@njnu.edu.cn。
摘要：为综合分析不同工况下三元颗粒混合物填充模腔的分离行为，采用计算流体力学-离散单元法(CFD-DEM),模拟高密度、无黏性颗粒混合物填充模具过程，并研究送料速度、颗粒生成方式、送料机构壁面斜度、模具形状与深度以及环境空气等多种因素对其分离行为的影响。结果表明，增加送料机构行进速度，可在行进方向上显著抑制大粒径颗粒的分离行为，降幅高达71%,而在其他方向上对其分离行为的影响较小；降低送料机构壁面斜度、适当加大模具深度均能提高颗粒混合物在模腔内的分布均匀度；分层生成方式产生的颗粒混合物，位于中间层颗粒群的分离行为受到显著抑制，分离指数最小仅为0.8;模腔内的空气会在颗粒填充时发生逸散从而迟滞细颗粒的运动，使得约40%的细颗粒在局部区域集聚，不利于颗粒混合物在模具内的均匀分布。
关键词：三元颗粒混合物；模腔填充；颗粒分离；计算流体力学-离散单元法
Abstract：In order to comprehensively analyze the separation behavior of the ternary particle mixture filling the mold cavity under different working conditions, the computational fluid dynamics-discrete element method(CFD-DEM) was used to simulate the process of filling the mold with a high-density, non-viscous particle mixture. The effects of various factors such as feeding speed, particle generation method, wall inclination of feeding mechanism, mold shape and depth, and ambient air on its separation behavior were studied. The results show that increasing the travel speed of the feeding mechanism can significantly inhibit the separation behavior of large-sized particles in the travel direction, and the drop rate is as high as 71%, while the effect on the separation behavior in other directions is small. Reducing the wall slope of the feeding mechanism and appropriately increasing the depth of the mold can improve the distribution uniformity of the particle mixture in the mold cavity. For the particle mixture produced by the layered generation method, the separation behavior of the particle group located in the middle layer is significantly suppressed, and the minimum separation index is only 0.8. The air in the mold cavity will escape when the particles are filled, which will retard the movement of the fine particles, causing about 40% of the fine particles to accumulate in the local area, which is not conducive to the uniform distribution of the particle mixture in the mold.
Keywords：ternary particle mixture; cavity filling; particle separation; computational fluid dynamics-discrete element method
参考文献（References）：
[1]钟文镇,昃向博,李国平,等.球形玻璃珠填充模腔的流动特性试验[J].中国粉体技术,2019,25(6):32-37.
[2]HUAQIN Y,YUKI M,KAZUYA T,et al.Numerical investigation on the influence of air flow in a die filling process[J].Journal of the Taiwan Institute of Chemical Engineers,2018,90:9-17.
[3]XU Z B,YOSHINAGA S,TSUNAZAWA Y,et al.Numerical investigation of segregation behavior of multi-sized particles during pharmaceutical mini-tablet die filling[J].Journal of Drug Delivery Science and Technology,2020,61:102301.
[4]WIDARTININGSIH P M,MORI Y,TAKABATAKE K,et al.Coarse graining DEM simulations of a powder die-filling system[J].Powder Technology,2020,371:83-95.
[5]TSUJI Y,KAWAGUCHI T,TANAKA T.Discrete particle simulation of twodimensional fluidized bed[J].Powder Technology,1993,77(1):79-87.
[6]HOOMANS B P B,KUIPERS J A M,BRIELS W J,et al.Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidised bed:a hard-sphere approach[J].Chemical Engineering Science,1996,51(1):99-118.
[7]NICOLIN N,WILKE D N,WU C Y,et al.Large-scale GPU based DEM modeling of mixing using irregularly shaped particles[J].Advanced Powder Technology,2018,29(10):2476-2490.
[8]ZHOU M M,WANG S,KUANG S B,et al.CFD-DEM modelling of hydraulic conveying of solid particles in a vertical pipe[J].Powder Technology,2019,354:893-905.
[9]GOU D,AN X,YANG X,et al.CFD-DEM modeling on air impact densification of equal spheres:structure evolution,dynamics and mechanism[J].Powder Technology,2017,322:177-184.
[10]CHEN J L,LI X F,HUAI X L,et al.Numerical study of collection efficiency and heat-transfer characteristics of packed granular filter[J].Particuology,2019,46:75-82.
[11]HERTZ H Z.On the contact of elastic solids[J].Journal Für Die Reine und Angewandte Mathematik (Crelles Journal),1880,92:156-171.
[12]MINDLIN RD,DERESIEWICZ H.Elastic spheres in contact under varying oblique force[J].Journal of Applied Mechanics,1953,20(3):327-344.
[13]YUKI T,YUSUKE S,CHIHARU T,et al.Numerical simulation of industrial die filling using the discrete element method[J].Chemical Engineering Science,2015,138:791-809.
[14]DANIEL M,FERNANDO J M,RAFAEL M.Particle size segregation promoted by powder flow in confined space:the die filling process case[J].Powder Technology,2014,262:215-222.
[15]GUO Y,WU C Y,THORNTON C.The effects of air and particle density difference on segregation of powder mixtures during die filling[J].Chemical Engineering Science,2011,66(4):661-673.
[16]CRISTINA R A,FERNANDO A E,ANA G M,et al.Segregation in the tank of a rotary tablet press machine using experimental and discrete element methods[J].Powder Technology,2018,328:452-469.
[17]GUO Y,KAFUI K D,WU C Y,et al.A coupled DEM/CFD analysis of the effect of air on powder flow during die filling[J].AICHE Journal,2010,55(1):49-62.
[18]RYOICHI F,YUKI S,KAZUNORI K,et al.Size-induced segregation during pharmaceutical particle die filling assessed by response surface methodology using discrete element method[J].Journal of Drug Delivery Science and Technology,2016,35:284-293.
[19]江帆,肖纳,黄春曼.粉末冶金装粉过程的充型均匀性研究[J].特种铸造及有色合金,2017,37(8):820-823.
[20]SAED S R,VIRENDRA M P.Uniformity differentiation analysis of powder deposition characteristics in circular and rectangular shallow dies using feed shoe with square cross-section[J].Particulate Science and Technology,2011,29(3):260-271.
