蔡文源¹，王利强¹，徐立敏²

（1. 江南大学 江苏省食品先进制造装备技术重点实验室，机械工程学院，江苏 无锡 214122；2. 江苏创新包装科技有限公司，江苏 扬州 225600）

引用格式：

蔡文源，王利强，徐立敏. 基于CFD-DEM的超细碳酸钙螺旋输送仿真分析［J］. 中国粉体技术，2024，30（3）：100-111.

CAI W Y， WANG L Q， XU L M. Simulation analysis of ultrafine calcium carbonate spiral transportation based on CFD-DEM［J］. China Powder Science and Technology，2024，30（3）：100−111.

DOI：10.13732/j.issn.1008-5548.2024.03.009

收稿日期：2024-03-05，修回日期：2024-04-12，上线日期：2024-04-25。

基金项目：中央高校基本科研业务费专项资金项目，编号：JUSRP21115；江苏省食品先进制造装备技术重点实验室自主研究课题资助项目，编号：FMZ202304。

第一作者简介：蔡文源（1999—），男，硕士生，研究方向为包装工艺与机械。E-mail：c2752096870@163. com。

通信作者简介：王利强（1977—），男，教授，博士，江苏省“双创博士”人才计划计划、江苏省“企业创新岗”特聘专家，主要从事食品包装技术研究。E-mail： liqiang-wang@jiangnan. edu. cn

摘要：【目的】 为提高螺旋输送机的输送效率，降低输送机的功耗与磨损，探究在不同进料速率、螺旋轴转速与几何体摩擦系数下，超细碳酸钙在水平变径变距螺旋输送机内的颗粒流动状态、出口质量流量、输送机功耗与磨损分布。【方法】使用计算流体动力学（computational fluid dynamics， CFD）与离散单元法（discrete element method， DEM）双向耦合数值模拟的方法，对螺旋输送机在不同转速下的质量流率进行分析对比，验证数值模型的正确性。【结果】 摩擦系数对颗粒的运动有较大影响，颗粒流的轴向速度峰值和质量流率峰值随着摩擦系数的增加先增大再减小；随着下料速度和摩擦系数的增大，输送机功率明显增大，且摩擦系数在高进料速度与低转速的情况下对功耗的影响相对于低进料速度和高转速更加明显；磨损较严重的区域集中在下料口处的螺旋轴与螺旋叶片的边缘处。【结论】 简单增大或减小摩擦系数并不能提高颗粒的轴向速度和质量流量，而是存在一个局部最优参数组合；适当地提高转速能够减小颗粒密实度与颗粒停留时间，从而减小输送机的功耗与几何体磨损。

关键词：超细碳酸钙；计算流体动力学；离散单元法；螺旋输送机；颗粒流动

Abstract

Objective As an important inorganic non-metallic powder material， ultrafine calcium carbonate is widely used in plastics， rubber， coatings， and paper making， etc. Due to its extremely small particle size and adhesive nature， there are difficulties in its conveying process leading to lower conveying efficiency and increased equipment wear. Therefore， it requires improved design and operation of conveying equipment. Screw conveyor， as an efficient and widely used solid material conveying equipment， is characterized by its simple structure， low maintenance costs， and ability to convey materials in horizontal， inclined， or vertical directions. With the advancement of computer simulation technology， the discrete element method （DEM） technique has become an important tool to study the particle flow and interactions in screw conveyor conveying process. Current research on CFD-DEM coupling method for screw conveyors typically focuses on operational parameters for large particle materials. However， they overlook investigations into the role of geometric friction coefficients in the ultrafine powder conveying. In order to improve the efficiency of screw conveyors and reduce power consumption and wear， the study was conducted to examine the particle flow state， outlet mass flow rate， conveyor power consumption， and wear distribution of ultrafine calcium carbonate in a horizontal screw conveyor with variable diameter and pitch. This study provided theoretical and technical support for effective transport of ultrafine powders， as well as new perspectives and methodological foundation for engineering applications and scientific research in related fields.

Methods In this study， we adopted the CFD-DEM coupling method， using FLUENT， a fluid dynamics software， and EDEM， a discrete element software， to explore the gas-solid two-phase flow characteristics of ultrafine calcium carbonate powder during horizontal conveying with variable diameter and pitch of the spiral. It also explored the effects of different operational and geometric parameters on the conveying process. In CFD-DEM gas-solid coupling， the gas phase was treated as a continuous phase， governed by continuity and momentum conservation equations. RNG k-ε model was used in the turbulence model， which was more suited for the complex turbulent flow inside the screw conveyor. The particles were regarded as a discrete phase， and their motion was described by the Newtonian kinetic equations. The Archard wear model was used to simulate the wear process on the spiral blade surfaces correlating the amount of material worn on the surface of the spiral blade to the friction work exerted by the frictional action of the particles on that surface. The simulation model of the horizontal screw conveyor with variable diameter and pitch was constructed in this paper and its specific structure was shown in Fig. 1 and 2. A tetrahedral unstructured mesh with high adaptability was chosen to be applied for meshing in this study， as shown in Fig. 3. Parameters of the simulation experiment were detailed in Tab. 1.

Results Simulated mass flow rates at different rotational speeds closely aligned with experimental results， showing a small margin of error. This high level of accuracy confirmed the model's reliability for simulating the screw conveying of ultrafine calcium carbonate. In the horizontal screw conveyor with variable diameter and pitch， particles uniformly flowed towards the outlet along the screw shaft， demonstrating stable conveying performance. The coefficient of friction between particles and geometry significantly influenced particle behavior in various conveying sections. However， optimal conveying relied not solely on higher or lower friction coefficients but instead on finding an appropriate balance. Variations in peak mass flow rates were linked to particle axial velocity， which was influenced by friction coefficients. Excessive friction impeded particle flow， while insufficient friction reduced particle forward momentum， leading to equipment clogging and reduced mass flow rates. Power consumption analysis at a constant feed rate of 2 kg/s revealed decreasing consumption with increasing rotational speed， with the impact of the friction coefficient diminishing at higher speeds. Wear simulation results indicated that higher rotational speeds reduced material residence time on blade surfaces， thus reducing wear interactions. Conversely， increased friction coefficients significantly elevated wear rates due to heightened sliding or rolling resistance between material and blade surfaces. This abrasive interaction intensified wear rates under conditions of higher friction.

Conclusion In a screw conveyor， materials moved axially and in a circular flow within each section. The coefficient of friction between the particles and the geometry had a significant effect on the particle flow state. Changes in the friction coefficient did not necessarily result in increased axial velocity of the particles or the mass flow rate of the conveyor. Instead， there existed a locally optimal combination of parameters that could be achieved. As the discharge speed and friction coefficient increased， the power consumption of the conveyor also rose. An appropriate increase in rotational speed could reduce power consumption. However， the influence of the friction coefficient on power consumption was more pronounced at high feed speeds and low rotational speeds， compared to low feed speeds and high rotational speeds. The more serious wear area concentrated at the edges of the spiral shaft and spiral blades near the lower feed opening. For any given rotational speed， higher coefficient of friction between the particles and the geometry with greater feed speed resulted in more severe wear of the spiral blades and spiral shaft. Conversely， higher rotational speeds of the spiral shaft led to lower overall wear rates. The study provides theoretical and technical support for efficient conveying of ultrafine powders， and offers new perspectives and methodological basis for its engineering applications and scientific research in related fields.

Keywords：ultrafine calcium carbonate； Computational fluid dynamics； Discrete element method； Screw conveyor； Particle flow

参考文献（References）

［1］TRAN V T， BUI N T， BUI T A. Application of EDEM simulation for calculating and optimizing a closed coal fly ash screw conveyor［J］. Applied Sciences，2023，13（22）：12169.

［2］TROFIMOV A D， IVANOV A A， ZYUZIN M V， et al. Porous inorganic carriers based on silica， calcium carbonate and calcium phosphate for controlled/modulated drug delivery： fresh outlook and future perspectives［J］. Pharmaceutics，2018，10（4）：167-168.

［3］ALOBAID F， ALMOHAMMED N， FARID M M， et al. Progress in CFD simulations of fluidized beds for chemical and energy process engineering［J］. Progress in Energy and Combustion Science，2022，91：100-930.

［4］XIE C H， ZHAO Y Z. Investigation of the ball wear in a planetary mill by DEM simulation［J］. Powder Technology，2022，398：117-057.