1. 引言
中药配方颗粒粒径较小,极易发生扬尘,这对中药颗粒配药机的颗粒传输和汇集提出了较高要求。气力输送设计易于操作,可以简单地改变管道配置来调整输送路径,在保障物料纯净度方面具有显著优势[1]。基于气力输送设备优点,研发了多品种颗粒类物料的汇集装置。由于在试验条件下不能对颗粒和气流进行有效监控,本文选用CFD-DEM耦合来仿真模拟颗粒在汇集装置中的汇集过程,该方法在农业和工业等[2] [3]工程领域已被证明是有效的。王晓晨[4]采用CFD-DEM耦合的方式,对油茶果负压吸附系统开展了数值模拟与试验研究,得到了各因素的影响大小和最优的工况,试验与数值模拟结果基本一致。Zhou [5]采用CFD-DEM研究了旋风除尘器中微米级颗粒的速度和压力分布,并将仿真模拟与试验进行对比,二者数据基本吻合。张凯[6]采用CFD-DEM对气吸式排种器进行了仿真研究,得到了排种器的最佳工况。计算流体动力学和离散元耦合模型可以跟踪单个颗粒的运动和受力,考虑颗粒形状、材料属性等因素,可以更加直观观察到颗粒的复杂流动过程,更适合于气固两项流的研究。
本文以颗粒汇集装置为研究对象,考察汇集装置在不同出口速度下的颗粒流态的变化,寻找与汇集装置想匹配的速度区间,结合试验验证,证明了颗粒汇集装置对颗粒传输和汇集的可行性,并对比了仿真模拟与汇集试验结果的一致性,为中药配方颗粒汇集工艺优化提供参考。
2. 数学模型
2.1. 几何模型
汇集装置采用多支路管道与旋风分离器结合的结构形式,实现多品种中药颗粒的传输与汇集。颗粒汇集装置几何模型如图1所示,汇集装置左右两侧共设置了24个管道入口,可以多路并行对颗粒的汇集,工作时,颗粒依靠自身重力和旋风分离器上方风机提供的负压吸入管道及旋风分离器中,并在旋风分离其中完成气固分离,气体从旋风分离器的排气口排出,中药颗粒下落到料仓,几何结构参数如表1和表2。
Figure 1. Geometric model of particle collecting device
图1. 颗粒汇集装置几何模型
Table 1. Particle collecting device geometry parameters (mm)
表1. 颗粒汇集装置几何结构参数(mm)
Ta |
Tb |
Tc |
Im |
Io |
Ii |
Ic |
I |
Oc |
Oi |
965 |
150 |
300 |
1558 |
350 |
500 |
400 |
400 |
18 |
25 |
Table 2. Cyclone geometry parameters (D = 100 mm)
表2. 旋风分离器几何结构参数(D = 100 mm)
a/D |
b/D |
De/D |
h/D |
Ht/D |
S/D |
Be/D |
Li/D |
Hc/D |
Dc/D |
0.4 |
0.24 |
0.415 |
1.2 |
3 |
0.5 |
0.375 |
1 |
1 |
0.5 |
2.2. 数学模型
2.2.1. 固相控制方程
在CFD-DEM中,固相被视为离散项,粒子在任何时间t的平移和旋转运动都可以用牛顿定律来描述:
(1)
(2)
式中:
,
,
和
分别是粒子i的质量,惯性矩,平移和旋转速度。作用在固相上的力是气固相互作用力
,粒子i和j之间的粒子之间的力包括,接触力
,粘性阻尼力
和重力
。粒子间力在与粒子i接触的
粒子上求和。切向力产生力矩
,并导致粒子i旋转,因为粒子间力作用在粒子i和j之间的接触点,而不是粒子中心。
是与第i个粒子旋转相反的滚动摩擦转矩。
2.2.2. 固相控制方程
流体被视为连续相,其微观尺度运动由局部平均Navier-stokes方程求解,质量和动量守恒由如下式:
(3)
(4)
式中:
为孔隙率;
为平均流体速度和波动流体速度;
为体积流体–颗粒相互作用力;
和g分别为时间,流体密度,压力,粘性应力张量和重力加速度。
是CFD单元中的网格数量,
是网格体积;
是湍流引起的雷诺应力项。
2.3. 网格划分
使用fluent meshing中对颗粒汇集装置几何模型进行三维非结构六面体网格划分,图2显示了颗粒汇集装置CFD网格,入口打开数量为10的条件下,边界层数为3层,网格最大歪斜律为0.3,达到计算要求。对网格无关性进行验证,在出口流速10 m/s,分别设置了三组不同计算节点的网格,分别为282,926,363,630,487,625个网格数量,稳态模式迭代稳定。计算压降差异小于5%,这表明计算结果与网格大小无关,见表3。因此最终计算采用节点数量为363,630的网格。
Figure 2. Particle collecting device CFD grid
图2. 颗粒汇集装置CFD网格
Table 3. Grid independence verification
表3. 网格无关性验证
网格数目 |
压降(pa) |
最大计算误差/% |
282,926 |
660.8 |
|
363,630 |
689.7 |
4.19% |
487,625 |
662.7 |
|
2.4. 边界条件及计算设置
在CFD-DEM耦合计算中,采用SIMPLEC压力–速度耦合算法,标准壁面,RSM湍流模型,插值算法PRESTO,求解动量方程为QUICK,湍流动能和耗散曲率为二阶迎风格式,RSM为一阶迎风。对照中国药典常用处方的中成药配伍,主要药味在10味以下的应占70%以上,在颗粒汇集装置中选取汇集装置左右边各五个入口,共10个多支路管道入口。入口设为压力入口,出口设为速度出口,壁面处为标准无滑移边界条件,入口和出口处采用5%的湍流强度,fluent时间步长为1 × 10−4 s。在模拟时,DEM中管道入口的颗粒由颗粒工厂随机产生,DEM时间步长为1 × 10−6 s,模拟仿真均采用双向耦合。共进行了5种工况的模拟,当汇集装置对颗粒进行传输汇集时,排气口出口速度分别8 m/s,10 m/s,12 m/s,当颗粒完全进入旋风分离器进行颗粒汇集时,出口速度为6 m/s,3 m/s,根据常见内科疾病中药配方颗粒处方精选[7],据统计颗粒剂处方一次最多服用均在20 g以内,设置每个入口颗粒质量流率和给料量分别为40 g/s和2 g,本文采用Furukawa R [8]文献中对制药颗粒的参数标定数据。
3. 模拟结果和分析
3.1. 不同出口速度下的颗粒流
图3显示了3 s时不同出口速度下汇集装置的颗粒流态(a = 8 m/s,b = 12 m/s),在锥体区域和料仓区域建立grid bin group统计颗粒质量。如图3(a)所示,当出口速度为8 m/s时,有颗粒在管道内发生堆积,管道中的气体流速较小,不足以使颗粒完全进入旋风分离器中进行气固分离;同时可以观察,管道内残留的颗粒位置有所不同,这是由于开口位置不同,造成的气路管道内气体流量的不匹配;如图3(b)所示,当出口速度为12 m/s时,颗粒完全进入旋风分离器中,汇集过程中有部分颗粒从汇集装置排气口逃逸,分离效率下降。随着出口速度逐渐增加,锥体部分的颗粒向上运动,且圆周旋转的颗粒逐渐增加。
Figure 3. Particle flow patterns of particle collecting devices at different exit speeds
图3. 不同出口速度下颗粒汇集装置颗粒流态
3.2. 颗粒流
颗粒经多支路管道进入旋风分离器后,在离心力的作用下沿着分离器的内壁面向下流动[9],图4显示了出口速度为10 m/s时颗粒汇集装置内部颗粒流态,经过3 s的模拟时间,在当前的CFD-DEM模型中,颗粒在管道无残留,并且颗粒没有从出口逃逸的现象发生,在颗粒进入分离空间后,旋风分离器锥体区域产生一些“停滞”颗粒做水平圆周旋转,没有到达底部,出现这种现象是因为圆锥部分的速度分布剖面值较高,对旋风分离器底部料仓的颗粒收集产生负面影响[10],颗粒将悬浮在锥形部分,不会从旋风分离器中出来,根据这些观测结果,在颗粒经多支路管道进入旋风分离器后,需要设计的一个合适的速度来收集旋流器底部的颗粒。
Figure 4. Particle flow pattern in particle collecting device
图4. 颗粒汇集装置内部颗粒流态
3.3. 不同出口速度下旋风分离器内部颗粒收集
图5显示了颗粒进入旋风分离器后,不同出口速度下的颗粒收集状态(a = 6 m/s,b = 3 m/s)。如图5(a)所示,当出口速度为6 m/s时,仍有部分颗粒悬浮在旋风分离器的锥体区域,如图5(b)所示,当出口速度下降到3 m/s时,颗粒完全收集在料仓中,达到宏观稳定状态。因此,当颗粒通过管道完全进入旋风分离器进行汇集时,出口速度至少要下降到3 m/s及以下。
Figure 5. Particle collection state inside cyclone separator at different exit speeds
图5. 不同出口速度下旋风分离器内部颗粒收集状态
4. 汇集装置试验与结果对比
在试验时,气流是通过安装在旋风分离器及过滤器上的负压风机产生并引入旋风分离器及管道的,而在模拟过程中,风机和过滤器被移除,并且仅设置其出口风速的影响以节省计算资源。根据机械结构和仿真模拟结果,利用多支路管道和旋风分离器结合的方式组成的汇集装置可以实现多种颗粒的传输与汇集,因此基于汇集装置的结构,加工并组装了一套汇集装置试验样机,图6显示了颗粒汇集装置试验样机。
Figure 6. Experimental prototype of particle collecting device
图6. 颗粒汇集装置试验样机
试验所用的中药颗粒选用华润三九医药股份有限公司生产的漏芦配方颗粒,试验时,打开与仿真模拟位置相同的十个管道支路入口并连接颗粒给料器,图7显示了不同出口速度下颗粒汇集装置试验,当汇集装置出口速度为8 m/s时,有部分颗粒残存在旋风分离器入口处无法进入旋风分离器中,而同一条件下仿真值也为8 m/s时,同样有部分颗粒残留在旋风分离器入口处;当出口速度为9.2 m/s时,可以观察到颗粒在旋风分离器入口处没有残留,并有部分颗粒在锥体区域做水平圆周运动。当旋风分离器出口速度由9.2 m/s下降到4 m/s后,颗粒完全落入汇集装置下方的料仓中,而同条件下的仿真模拟值为3 m/s。分别进行了八组分离效率试验,见表4,由表4可知汇集装置平均分离效率达98%,而仿真模拟出口速度为10 m/s时,分离效率为100%,在分离分离效率方面,二者误差为2%,根据中国药典第四部通则0104 [11]颗粒剂6.0 g以上,装量差异限度±5%的规定,该颗粒汇集装置符合设计要求。对比试验与仿真的方案,对误差产生的原因进行归纳分析,结果如下。
Figure 7. Experiments of particle collecting devices at different exit speeds
图7. 不同出口速度下颗粒汇集装置试验
Table 4. The separation efficiency of the granules of the gulu formula (Outlet velocity 9.2 m/s)
表4. 漏芦配方颗粒分离效率(出口速度9.2 m/s)
实验组数 |
给料器给料量(g) |
料仓收集量(g) |
分离效率% |
平均值(%) |
1 |
20.00 |
19.70 |
98.5 |
98 |
2 |
20.00 |
19.42 |
97.1 |
3 |
20.00 |
19.44 |
97.2 |
4 |
20.00 |
19.50 |
97.5 |
5 |
20.00 |
19.56 |
97.8 |
6 |
20.00 |
19.75 |
98.7 |
7 |
20.00 |
19.83 |
99.1 |
8 |
20.00 |
19.62 |
98.1 |
1) 试验过程中的颗粒粒径不统一,漏芦配方颗粒的粒径范围为200~1500 um,仿真过程中粒径统一设置为800 um,并且仿真使用的球形模型与实际的中药配方颗粒形状存在差异。
2) 在试验中,多叉路管道采用PTEE四氟圆管,而仿真模拟采用的是不锈钢圆管参数,两圆管在力学性能和表面粗糙度方面的差异也会扩大试验误差。
5. 结论
利用CFD-DEM的耦合方法,对颗粒在汇集装置中的颗粒运动过程进行仿真,当出口气体流速为8 m/s时,部分颗粒在管路发生堆积,无法进入旋风分离器进行气固分离;当出口速度为10 m/s时,颗粒经过气路管道完全进入旋风分离器中;当出口达一定速度对旋风分离器内的颗粒进行汇集时,旋风分离器锥体区域产生一些“停滞”颗粒做水平圆周旋转,随着出口速度逐渐增加,锥体部分的颗粒向上运动,且圆周旋转的颗粒逐渐增加,颗粒进入旋风分离器进行气固分离时,出口速度逐渐减小到3 m/s及以下时,锥体区域的颗粒下落到料仓,达到稳定状态。多支路管道和旋风分离器的结构形式能高效地完成多种颗粒的传输和汇集,试验数据与模拟结果之间有很好的一致性。
基金项目
2023辽宁省教育厅面上项目(揭榜挂帅服务地方):“大健康”背景下智能中药免煎颗粒收集及给药装置的研发(JYTMS20230395)。