留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于CFD-EDM的自动投饵饲料颗粒气力输送数值模拟

胡昱 黄小华 陶启友 袁太平 王绍敏 刘海阳

引用本文:
Citation:

基于CFD-EDM的自动投饵饲料颗粒气力输送数值模拟

    作者简介: 胡 昱(1982—),男,助理研究员,从事渔业设施工程技术与装备研究。E-mail: scshuyu@163.com;
  • 中图分类号: S 969.31+4

Simulation of automatic pneumatic conveying process of automatic feeding system based on computational fluid dynamics-discrete element method

  • CLC number: S 969.31+4

  • 摘要: 针对深水网箱自动投饵系统颗粒饲料气力输送时易阻塞与破损的问题,为揭示自动投饵系统中颗粒饲料的气力输送运动特性,开展了管道颗粒饲料气固两相流的数值模拟研究。基于气固两相流理论,分别建立了计算流体模型与离散元模型,通过计算流体力学与离散单元法(computational fluid dynamics-discrete element method, CFD-EDM)耦合求解,对饲料颗粒从气力输送管道初始阶段到稳定阶段的运动过程进行了分析,得到颗粒从初始状态到运动稳定阶段的颗粒位置分布情况。对不同位置管道内压力进行了分析对比,得到从管道入口到颗粒运动稳定状态时刻的压降。模拟仿真结果为自动投饵装备的模拟仿真与优化设计提供了参考,使其可以对饲料颗粒运动有更好的导向性,更好地降低弯管处颗粒碰撞的能量损失。
  • 图 1  自动投饵系统装备原理结构图

    Figure 1.  Principle and structure diagram of automatic feeding system

    图 2  多球面饲料离散元模型

    Figure 2.  Discrete element model of multi sphere feed

    图 3  自动投饵设备结构示意图

    Figure 3.  Schematic diagram of auto-feeding equipment

    图 4  饲料管道内输送过程

    Figure 4.  Transportation process of pellet feed pipeline

    图 5  颗粒饲料动力统计区域

    Figure 5.  Dynamic statistics area of pellet feed

    图 6  不同阶段饲料颗粒动能

    Figure 6.  Kinetic energy of feed particles at different stages

    图 7  饲料颗粒不同运动段形态分布

    Figure 7.  Morphological distribution of feed particles at different motion stages

    图 8  颗粒运动不同距离处管道截面总压云

    Figure 8.  Total pressure cloud of pipe cross section at different distance of particle movement

    表 1  EDEM中物料属性及接触参数

    Table 1.  Material properties and contact parameters in EDEM

    名称
    name
    参数
    parameter
    数值
    volume
    物料 material 泊松比 0.4
    剪切模量/Pa 1×108
    密度/kg·m–3 390
    壳体 shell 泊松比 0.3
    剪切模量/Pa 7.0×1010
    密度/kg·m–3 7 800
    物料-物料 materiel-materiel 恢复系数 0.2
    静摩擦系数 0.6
    动摩擦系数 0.01
    物料-壳体 materiel-shell 恢复系数 0.15
    静摩擦系数 0.3
    动摩擦系数 0.01
    下载: 导出CSV
  • [1] 郭根喜, 陶启友, 黄小华, 等. 深水网箱养殖装备技术前沿进展[J]. 中国农业科技导报, 2011, 13(5): 44-49. doi: 10.3969/j.issn.1008-0864.2011.05.07
    [2] 庄保陆, 郭根喜. 水产养殖自动投饵装备研究进展及其应用[J]. 南方水产, 2008, 4(4): 67-72. doi: 10.3969/j.issn.2095-0780.2008.04.013
    [3] 胡昱, 郭根喜, 汤涛林, 等. 基于MCGS的深水网箱自动投饵远程控制系统的设计[J]. 渔业科学进展, 2010, 31(6): 110-115. doi: 10.3969/j.issn.1000-7075.2010.06.016
    [4] 胡昱, 郭根喜, 黄小华, 等. 基于PLC的深水网箱自动投饵系统[J]. 南方水产科学, 2011, 7(4): 61-68. doi: 10.3969/j.issn.2095-0780.2011.04.0010
    [5] 宋协法, 路士森. 深水网箱投饵机设计与试验研究[J]. 中国海洋大学学报(自然科学版), 2006, 36(3): 405-409.
    [6] TURMELLE C A, SWIFT M R, CELIKKOL B, et al. Design of a 20-ton capacity finfish aquaculture feeding buoy[C]. 2006 Oceans MTS, Boston, 2006: 1-6.
    [7] TSUNODA T, KITAZAWA D, KINOSHITA T, et al. Concept of an offshore aquaculture system with an automated feeding platform[C]. Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering, Estoril, 2008: 527-534.
    [8] De MORAES M S, MUIÑOS TORNEIROS D L, da SILVA ROSA L, et al. Experimental quantification of the head loss coefficient K for fittings and semi-industrial pipe cross section solid concentration profile in pneumatic conveying of polypropylene pellets in dilute phase[J]. Powder Technol, 2017, 310(4): 250-263.
    [9] VASQUEZ N, FAN Y, JACOB K. Effect of material elasticity and friction on particle dynamics during dilute pneumatic conveying[J]. Powder Technol, 2016, 303(12): 90-99.
    [10] BHATTARAI S, OH J H, EUH S H, et al. Simulation study for pneumatic conveying drying of sawdust for pellet production[J]. Drying Technol, 2014, 32(6): 1142-1156.
    [11] 祝先胜. 气力输送管内气固两相流动的数值模拟[D]. 上海: 华东理工大学, 2015: 6-8.
    [12] 王明旭, 秦超, 李永祥, 等. 气力输送过程中粮食颗粒的输送特性研究[J]. 农机化研究, 2014(9): 18-22. doi: 10.3969/j.issn.1003-188X.2014.09.004
    [13] 卢洲, 刘雪东, 潘兵. 基于CFD-DEM方法的柱状颗粒在弯管中输送过程的数值模拟[J]. 中国粉体技术, 2011, 17(5): 65-69. doi: 10.3969/j.issn.1008-5548.2011.05.017
    [14] OEHME M, AAS T S, SØRENSEN M, et al. Feed pellet distribution in a sea cage using pneumatic feeding system with rotor spreader[J]. Aquacult Engin, 2012, 51(11): 44-52.
    [15] HALSTENSEN M, IHUNEGBO F N, RATNAYAKE C. Online acoustic chemometric monitoring of fish feed pellet velocity in a pneumatic conveying system[J]. Powder Technol, 2014, 263(9): 104-111.
    [16] TURIDSYNNØVE A, MAIKE O, METTE S. Analysis of pellet degradation of extruded high energy fish feeds with different physical qualities in a pneumatic feeding system[J]. Aquacult Engin, 2011, 44(1): 25-34. doi: 10.1016/j.aquaeng.2010.11.002
    [17] 杜俊, 胡国明, 方自强, 等. 稀相气力输送CFD-DEM仿真中两种模型的比较[J]. 海军工程大学学报, 2015, 27(2): 57-60.
    [18] 丁岩峰, 孙奉昌, 张春霞. 射流喷射器用于干渣气力输送的实验研究[J]. 中国粉体技术, 2014, 20(1): 75-78. doi: 10.3969/j.issn.1008-5548.2014.01.018
    [19] 刘晓斌, 段广彬, 刘宗明, 等. 水平弯管内气固两相流的试验研究与数值模拟[J]. 济南大学学报(自然科学版), 2012, 26(3): 230-235.
    [20] 李志华, 刘凯, 李州, 等. 炭黑在水平管道中气力输送的数值模拟与分析[J]. 起重运输机械, 2012(1): 46-50. doi: 10.3969/j.issn.1001-0785.2012.01.014
    [21] 杨洪雷, 靳世平, 刘小康. 煤粉气力输送过程中采用不同输送管径的数值模拟研究[J]. 工业加热, 2015, 44(4): 62-64. doi: 10.3969/j.issn.1002-1639.2015.04.020
    [22] 李丹阳, 刘姝, 王晓宁. 颗粒粒径对气力输送影响的数值模拟与分析[J]. 中国粉体技术, 2015, 21(6): 8-11.
  • [1] 庄保陆郭根喜 . 水产养殖自动投饵装备研究进展与应用. 南方水产科学, 2008, 4(4): 67-72.
    [2] 胡昱郭根喜黄小华陶启友 . 深水网箱养殖自动投饵计量装置设计. 南方水产科学, 2014, 10(2): 80-85. doi: doi:10.3969/j.issn.2095-0780.2014.02.012
    [3] 胡昱郭根喜黄小华陶启友张小明 古恒光 . PLC在深水网箱自动投饵系统中的应用. 南方水产科学, 2011, 7(4): 61-68. doi: 10.3969/j.issn.2095-0780.2011.04.010
    [4] 郭根喜庄保陆王良运胡昱黄小华陶启友 . 基于PLC的远程气力输送自动投饵控制系统的设计与实现. 南方水产科学, 2008, 4(6): 7-16.
    [5] 刘朝阳王印庚孙晓庆 . 颗粒饲料携带细菌与大菱鲆疾病发生的相关性. 南方水产科学, 2009, 5(4): 13-21. doi: 10.3969/j.issn.1673-2227.2009.04.003
    [6] 胡昱郭根喜黄小华陶启友张小明 . 基于ANSYS的水下洗网机流场模拟仿真. 南方水产科学, 2010, 6(1): 7-11. doi: 10.3969/j.issn.1673-2227.2010.01.002
    [7] 李敏张汉华朱长波 . 环境因子与饲料对鲍生长影响的研究进展. 南方水产科学, 2007, 3(2): 76-80.
    [8] 王旭李娜耿安静杨慧王富华 , . 高效液相色谱-电感耦合等离子体质谱联用技术测定海产品中5种形态砷. 南方水产科学, 2013, 9(2): 50-56. doi: 10.3969/j.issn.2095-0780.2013.02.009
    [9] 王信超孙敬敬范美华武梅石戈王日昕廖智 . 厚壳贻贝血细胞颗粒的蛋白质组学分析. 南方水产科学, 2012, 8(2): 7-14. doi: 10.3969/j.issn.2095-0780.2012.02.002
    [10] 李坚明甘晖冯广朋黄凯龚竹林 . 饲料脂肪含量与奥尼罗非鱼幼鱼肝脏形态结构特征的相关性. 南方水产科学, 2008, 4(5): 37-43.
    [11] 许晓娟区又君李加儿 . 延迟投饵对卵形鲳鲹早期仔鱼阶段摄食、成活及生长的影响. 南方水产科学, 2010, 6(1): 37-41. doi: 10.3969/j.issn.1673-2227.2010.01.007
    [12] 陈培基李刘冬黄珂 . 自动定氮仪测定水产品中的挥发性盐基氮. 南方水产科学, 2005, 1(6): 50-54.
    [13] 徐文其 . 小包装鱿鱼丝自动精确称量分析. 南方水产科学, 2011, 7(6): 72-76. doi: 10.3969/j.issn.2095-0780.2011.06.012
    [14] 粟丽朱长波陈素文 . 混养罗非鱼对凡纳滨对虾养殖围隔中悬浮颗粒物和弧菌数量的影响. 南方水产科学, 2015, 11(6): 19-26. doi: 10.3969/j.issn.2095-0780.2015.06.003
    [15] 吴凡蒋明文华刘伟田娟喻丽娟陆星 . 饲料碳水化合物与脂肪比例对吉富罗非鱼成鱼生长、体成分和血清生化指标的影响. 南方水产科学, 2019, 15(4): 53-60. doi: 10.12131/20190047
    [16] 许霄霄刘伟文华蒋明吴凡 . 高糖饲料对吉富罗非鱼生长性能、饲料利用和糖脂代谢的影响. 南方水产科学, 2017, 13(5): 94-102. doi: 10.3969/j.issn.2095-0780.2017.05.013
    [17] 林黑着叶乐陈燕军周竞晖李卓佳吴开畅 . 速大肥在凡纳滨对虾饲料中的应用. 南方水产科学, 2005, 1(6): 63-65.
    [18] 李海涛胡云华王银东 . 酸化剂对奥尼罗非鱼生长性能和饲料利用的影响. 南方水产科学, 2009, 5(5): 67-71. doi: 10.3969/j.issn.1673-2227.2009.05.012
    [19] 黄忠林黑着牛津吕国敏陈旭陈明强 . 肌醇对卵形鲳鲹生长、饲料利用和血液指标的影响. 南方水产科学, 2011, 7(3): 39-44. doi: 10.3969/j.issn.2095.0780.2011.03.006
    [20] 齐野孙向军于刚梁拥军宋立伟白东清乔秀亭 . 饲料碳水化合物水平对宝石鲈餐后代谢的影响. 南方水产科学, 2010, 6(2): 59-65. doi: 10.3969/j.issn.1673-2227.2010.02.010
  • 加载中
图(8)表(1)
计量
  • 文章访问数:  956
  • HTML全文浏览量:  515
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-28
  • 录用日期:  2019-02-21
  • 网络出版日期:  2019-03-21
  • 刊出日期:  2019-06-01

基于CFD-EDM的自动投饵饲料颗粒气力输送数值模拟

    作者简介:胡 昱(1982—),男,助理研究员,从事渔业设施工程技术与装备研究。E-mail: scshuyu@163.com
  • 中国水产科学研究院南海水产研究所,广东省网箱工程技术研究中心,农业农村部南海渔业资源开发利用重点实验室,广东 广州 510300

摘要: 针对深水网箱自动投饵系统颗粒饲料气力输送时易阻塞与破损的问题,为揭示自动投饵系统中颗粒饲料的气力输送运动特性,开展了管道颗粒饲料气固两相流的数值模拟研究。基于气固两相流理论,分别建立了计算流体模型与离散元模型,通过计算流体力学与离散单元法(computational fluid dynamics-discrete element method, CFD-EDM)耦合求解,对饲料颗粒从气力输送管道初始阶段到稳定阶段的运动过程进行了分析,得到颗粒从初始状态到运动稳定阶段的颗粒位置分布情况。对不同位置管道内压力进行了分析对比,得到从管道入口到颗粒运动稳定状态时刻的压降。模拟仿真结果为自动投饵装备的模拟仿真与优化设计提供了参考,使其可以对饲料颗粒运动有更好的导向性,更好地降低弯管处颗粒碰撞的能量损失。

English Abstract

  • 网箱养殖是我国海水鱼类养殖的主要方式,发展深水网箱养殖是我国海洋开发的重要战略和设施养殖的主要趋势[1-2]。深水网箱鱼类养殖过程主要采用自动投饵系统进行投饵,在自动投饵系统中,采用管道气力输送颗粒饲料。颗粒饲料在输送过程中,常常会由于输送工艺参数选取不合理导致饲料颗粒堆积阻塞管道或饲料破损等问题[3-7]

    目前对气固两相流的管道气力输送的研究主要针对粮食、橡胶等颗粒在直管弯管的运动特性模拟[8-10],祝先胜[11]模拟了三维90°弯管结构内的气固两相流动过程,获得了不同工况下的颗粒运动轨迹、速度、碰撞作用力、碰撞频率及流体压力、速度分布等。王明旭等[12]以小麦颗粒为研究对象,建立了以气流速度、小麦颗粒在输送过程中的体积百分含量为优化参数,流场内部压力为目标的小麦颗粒输送特性优化模型。卢洲等[13]采用计算流体力学和离散单元法耦合模拟弯管内的柱状颗粒气力输送过程,对不同弯径比的90°弯管内柱状颗粒的运动状态、碰撞特性、破碎原因及相关的力学特性进行了研究。针对气力投饵管道输送的研究主要是通过现场实测比较输送工艺参数对投饵的影响。Oehme等[14]研究了深水网箱鲑鱼养殖气力投饵系统中,不同气力输送速度和不同的颗粒粒径,颗粒饲料旋转撒料装置的抛洒面积分布状况。Halstensen等[15]提出了一种鱼饲料颗粒速度的无损监测方法,给出了颗粒速度和空气流量的部分最小二乘验证结果。Turidsynnøve等[16]利用3种不同物理品质三文鱼颗粒饲料,在不同风速及进料速率下,研究气力输送时的颗粒破碎情况,发现在高空速和低进料速率下,颗粒破碎率最高。国内外众多学者就气固两相流建立起了许多数学模型,并将这些模型应用于输送对象广泛的颗粒气力输送中,对气力输送下的颗粒运动特性进行了一些研究,而这些气力输送的数值模拟主要以二维气固两相流体模型为主,这与实际气力输送过程尚存在不少偏差[17-22]

    本文基于气固两相流理论,利用计算流体力学和离散单元法进行耦合模拟,对气力输送下颗粒饲料在投饵系统管道初始阶段到稳定阶段运动过程中的流动状态与运动特性进行研究,并分析不同阶段饲料颗粒动能、管内的压降、物料所受的曳力。最后根据仿真分析结果,对设备结构改进提出了指导意见。

    • 气体的流动规律用气固两相耦合的N-S方程描述。其中气相连续性方程:

      $ \displaystyle\frac{{\partial ({\alpha _g}{\rho _g}\,{u_{gj}})}}{{\partial {x_j}}} = 0 $

      气相动量方程:

      $ \begin{split} \displaystyle\frac{{\partial ({\alpha _g}{\rho _g}\,{u_{gi}}{u_{gj}})}}{{\partial {x_j}}} =& \displaystyle\frac{\partial }{{\partial {x_j}}}\left[{\mu _e}(\displaystyle\frac{{\partial {u_{gi}}}}{{\partial {x_j}}} +\displaystyle\frac{{\partial {u_{gj}}}}{{\partial {x_i}}})\right] + {\alpha _g}{\rho _g}{g_i} - \\ &{\alpha _g}\displaystyle\frac{{\partial p}}{{\partial {x_i}}} + {\beta _i}({u_{pi}} - {u_{gi}}) \end{split}$

      固相连续性方程:

      $ \frac{{\partial ({\alpha _p}{\rho _p}\,{u_{pj}})}}{{\partial {x_j}}} = 0 $

      固相动量方程:

      $ \begin{aligned} \frac{{\partial ({\alpha _p}{\rho _p}\,{u_{pi}}{u_{pj}})}}{{\partial {x_j}}} = \frac{\partial }{{\partial {x_j}}}\left[{\mu _p}(\frac{{\partial {u_{pi}}}}{{\partial {x_j}}} + \frac{{\partial {u_{pj}}}}{{\partial {x_i}}}) - {p_p}{\delta _{ij}}\right] +\\ {\alpha _p}{\rho _p}{g_i} - {\alpha _p}\frac{{\partial p}}{{\partial {x_i}}} + {\beta _i}({u_{pi}} - {u_{gi}})\qquad\quad \end{aligned}$

      $ {\alpha _g}$'/>为气相体积分数,$ {\alpha _p}$'/>为固相体积分数,$ {\alpha _g}$'/>+ $ {\alpha _p}$'/>=1,

      $ {u_{gi}}$$ {u_{gj}}$$ {u_{pi}}$$ {u_{pj}}$分别为气相、固相在$ i{\text{、}}j$方向上的速度分量;$ {\rho _g}$$ {\rho _p}$分别是气相、固相的密度,$ {g_i}$为重力加速度在$ i$'/>方向上的分量;$ {\mu _e}$$ {\mu _p}$为气相、固相在有效粘度;$ {\beta _i}$为气固两相交换系数在$ i$'/>方向上的分量。

    • 采用Hertz接触理论,根据Mindlin切向接触力求解方法,建立Hertz-Mindlin无滑动接触模型。其中:

      $ {\text{颗粒间法向力,}} {F_n} = {E^*}\sqrt {{R^*}} \delta _n^{\frac{3}{2}} $

      $ {\text{法向阻尼力,}} F_n^d = - 2\sqrt {\frac{5}{6}} \beta \sqrt {{S_n}{m^*}} v_n^{rel} $

      $ {\text{切向阻尼力,}} F_n^d = - 2\sqrt {\frac{5}{6}} \beta \sqrt {{S_t}{m^*}} v_t^{rel} $

      $ \beta = \frac{{\ln {\rm e}}}{{\sqrt {{{\ln }^2}{\rm e} + {\pi ^2}} }} $

      $ {\text{法向刚度,}} \ {S_n} = 2{E^*}\sqrt {{R^*}\alpha } $

      $ {\text{切向刚度,}} {S_t} = 8{G^*}\sqrt {{R^*}\alpha } $

      $ {m^*} = \frac{{{m_1}{m_2}}}{{{m_1} + {m_2}}} $

      其中$ {R^*}$为等效粒子半径,$ {E^*}$为等效弹性模量,$ {m^*}$为等效质量,$ {G^*}$为等效剪切模量,$ \delta $为切向重叠量,$ v_n^{rel}$为相对速度的法向分量值,$ v_t^{rel}$为切向相对速度。

    • 自动投饵系统装备主要由罗茨风机、旋转供料器、料仓、水平管道、S形管道与输料管道等构成(图1)。罗茨风机产生的空气流进入管道,料仓的颗粒饲料通过旋转供料器排入输送管道,空气流混合颗粒饲料后,进入S形输料管,由出料管输送到输料管道。

      图  1  自动投饵系统装备原理结构图

      Figure 1.  Principle and structure diagram of automatic feeding system

    • 实际饲料颗粒为膨化颗粒,形状可近似为高6 mm、底面直径7 mm的圆柱体,其尺寸见图2-a。为了更好地使离散元模型与真实饲料形状接近,在EDEM中利用particle工具创建颗粒模型,采用18个直径为3.5 mm的球形颗粒聚合成饲料离散元模型(图2-bc)。物料属性及接触参数见表1

      图  2  多球面饲料离散元模型

      Figure 2.  Discrete element model of multi sphere feed

      名称
      name
      参数
      parameter
      数值
      volume
      物料 material 泊松比 0.4
      剪切模量/Pa 1×108
      密度/kg·m–3 390
      壳体 shell 泊松比 0.3
      剪切模量/Pa 7.0×1010
      密度/kg·m–3 7 800
      物料-物料 materiel-materiel 恢复系数 0.2
      静摩擦系数 0.6
      动摩擦系数 0.01
      物料-壳体 materiel-shell 恢复系数 0.15
      静摩擦系数 0.3
      动摩擦系数 0.01

      表 1  EDEM中物料属性及接触参数

      Table 1.  Material properties and contact parameters in EDEM

      为了提高仿真计算速度,节约计算机资源,对自动投饵颗粒饲料气力输送设备结构进行简化,只保留料仓、供料器、管道等部件。模型长2 919 mm,宽800 mm,高1 574 mm,简化后见图3

      图  3  自动投饵设备结构示意图

      Figure 3.  Schematic diagram of auto-feeding equipment

    • 仿真过程中,忽略供料器与输送管道衔接部位的流场区域,只对输送管道进行网格划分。为了准确增加流场内网格节点和捕捉数据变化,采用Ansys Workbench中meshing对管道内部区域进行非结构四面体网格划分,共产生116 675个网格,节点数42 206个。导入fluent后经检查无负体积网格产生。

      供料器转速为23 r·min–1,EDEM中仿真时间步长设置为5e–6 s,为瑞利时间步的35.3%,共进行4 s。仿真开始前利用颗粒静态生成法在料仓内产生了4万个多球面饲料颗粒。

      在fluent中选择k-$ \varepsilon $湍流模型,入口边界速度为28 m·s–1,湍流强度为5%,湍流黏度静速率10。速度压力耦合求解器选择SIMPLE,采用二阶迎风离散格式。

      耦合算法选择欧拉-欧拉双向模型,fluent时间步长为5e–4 s,为EDEM的100倍,共8 000个时间步长。

    • 饲料颗粒从供料腔进入输送管道后,由于主要受到气流的曳力作用开始做加速运动。饲料颗粒在气流入口段速度较小,基本维持在10 m·s–1以下,随着运动距离的不断加大,颗粒也获得较高速度(图4)。图中管道末端饲料颗粒速度达到最大(26 m·s–1),并且出现悬浮现象,这与实际现象一致。该设备由于供料器是圆柱形壳体,落料口并没有直接与供料器壳体两端面连接,仿真过程中发现供料器靠近两端面处有颗粒滞留,导致无法充分供料。在后续结构改进时,要对供料器壳体做进一步优化,使其可以对颗粒运动有良好的导向性。

      图  4  饲料管道内输送过程

      Figure 4.  Transportation process of pellet feed pipeline

      供料器接口处至管道末端位置水平距离为2 780 mm,对此段长度按695 mm等分成4段,对管道内颗粒平均动能进行统计。颗粒在气流推动下获得的动能沿运动方向不断增加(图5)。第1段区域内颗粒处于刚开始运动阶段,平均动能较小,基本为0.002 5 J;第2段区域颗粒速度维持在13 m·s–1左右,处于加速阶段,平均动能为0.02 J,相比第一阶段获得运动较多;第3区域中管道有距离较短的弧形弯道,颗粒经过此段后速度继续增加,平均动能为0.025 J,同时该阶段较前2个阶段平均动能波动性增加;第4区域中,颗粒运动速度达到最大(26 m·s–1),平均动能接近0.03 J,颗粒运动趋于稳定状态。对比4个阶段可以发现,颗粒在第1、第2阶段处于初始运动状态,气流要对饲料颗粒做功,使其发生运动,所以这2个阶段较第3、第4阶段耗能更大,气流压降也更大。颗粒经过加速后进入第3、第4阶段,速度逐渐达到最大值,并维持稳定状态,但由于颗粒空隙率增加,动能波动性加强(图6)。

      图  5  颗粒饲料动力统计区域

      Figure 5.  Dynamic statistics area of pellet feed

      图  6  不同阶段饲料颗粒动能

      Figure 6.  Kinetic energy of feed particles at different stages

      饲料颗粒具有不同的速度时在管道内分布的形态也不一样。颗粒从供料器进入输送管道后,由于喂入量较大,颗粒加速缓慢,颗粒与颗粒之间空隙率低,颗粒成整体移动方式前进(图7-a)。在颗粒整体前进一段距离后,由于存在颗粒间碰撞作用力、颗粒与壁面碰撞作用力、场内漩涡等较复杂情况,导致颗粒间空隙率增加,单个颗粒速度增加明显,团聚移动方式减弱(图7-b)。在弯管中,饲料颗粒由于离心力作用紧贴外侧管道壁面前进,在弯管末端颗粒被抛向水平管道上表面,此时与壁面产生冲击力作用,颗粒回弹到水平管道下表面,此过程中可以看到颗粒速度有所下降,产生了能量损失(图7-c)。对管道结构改进时,要减小弯管曲率,以降低弯管与直管连接处颗粒碰撞的能量损失。在图7-d中,颗粒运动速度达到最大,同时出现悬浮状态,颗粒位置布满整个管道截面,颗粒间空隙率也达到最大。整个过程中并未见到颗粒堆积现象发生,管道内排料顺畅。

      图  7  饲料颗粒不同运动段形态分布

      Figure 7.  Morphological distribution of feed particles at different motion stages

    • 饲料颗粒沿着管道前进所需的能量是由气流压降来实现的。在仿真过程中,根据能量守恒定律,一定质量的流体在入口截面所具有的能量等于颗粒运动的能量、沿程压力损失、局部损失、出口截面能量之和。仿真采用简化的管道模型,长度较短,并且没有不同段的辅助衔接装置。因此,忽略气体流动过程的局部压力损失。将饲料颗粒从供料器下落至管道内的位置作为起点,运动方向为正向,每隔695 mm做一管道总压云图(图8)。随着饲料颗粒运动距离增加,管道内总压下降明显。饲料颗粒在管道内刚开始发生运动时,气流总压为4 600 Pa。颗粒发生水平距离2 780 mm时速度达到26 m·s–1,管道总压3 614 Pa,压降986 Pa。这个过程饲料颗粒运动速度从0增加到26 m·s–1,压降能量大部分用于转化为颗粒运动的动能及颗粒群提升的重力势能。

      图  8  颗粒运动不同距离处管道截面总压云

      Figure 8.  Total pressure cloud of pipe cross section at different distance of particle movement

    • 物料所受耦合力有曳力、马格努斯旋转提升力、赛夫曼剪切提升力、压力梯度力等,在仿真中后3项作用力微小,因此计算过程中将其忽略。饲料颗粒所受曳力公式为:

      $ {F_D} = \frac{1}{8}{C_D}{\text {π}} d_p^2\overline \rho \mu _p^2 $

      式中$ {C_D}$为阻力系数;$ {\mu _p}$为球形颗粒在静止流体中的运动速度(m·s–1);$ {d_p}$为球形颗粒直径(m);$ \overline \rho $为流体密度(kg·m–3)。

      将非球形饲料颗粒等效为球形,$ {d_p}$=7 mm,Re=2.02×106对应的$ {C_D}$近似为0.2,$ {\mu _p}$=60 m·s–1,代入公式(12)可得$ {F_D}$=1.79×10–2 N,则单个饲料颗粒运动1 m气流作功为1.79×10–2 J。

    • 在此工况下,空气为不可压牛顿流,因此,由雷诺数计算公式$ {\mathop{\rm Re}\nolimits} = \displaystyle\frac{{\rho vd}}{\mu }$可得Re=2.02×106。所以,管道内空载时压降计算公式为:

      $ \Delta p = \frac{{2{f_{sw}}{\rho _w}{v_w}{L_y}}}{D} $

      式中$ \Delta p$为直管压降(Mpa);$ {f_{sw}}$为流体摩阻系数;$ {\rho _w}$为流体密度 (kg·m–3);$ {v_w}$为流体流经管道的速度(m·s–1);$ {L_y}$为直管长度(m);D为直管直径(m)。

      此工况下Re > 2 100为湍流运动,摩阻系数公式为:

      $ {f_{sw}} = \frac{1}{4}{\left( {\frac{1}{{0.87\ln \frac{D}{{2{e_p}}} + 1.74}}} \right)^2} $

      式中$ {e_p}$为直管的表面粗糙度(m)。

      由以上公式代入数据计算可得1 m管道内$ \Delta p$=24 Pa。

    • 1) 利用CFD-DEM耦合法对饲料颗粒在气力输送设备管道中初始阶段的运动过程进行了分析,得到颗粒从初始状态到运动稳定阶段的颗粒位置分布情况,以及不同阶段颗粒的平均动能大小。整个初始阶段过程中并未见到颗粒堆积现象发生。

      2) 饲料颗粒运动速度从0 m·s–1增加到26 m·s–1时,颗粒饲料管道入口到颗粒运动稳定状态时刻的压降为986 Pa,压降能量大部分用于转化为颗粒运动的动能及颗粒群提升的重力势能。

      3) 根据仿真分析结果,对设备结构改进提出了意见。对供料器壳体做进一步优化,使其可以对颗粒运动有良好的导向性。减小弯管曲率,以降低弯管与直管连接处颗粒碰撞能量损失。

      4) 经现场测试,投饵机水平管道出口风速为24.03 m·s–1,与仿真模拟结果比较接近,说明采用CFD-EDM两相流耦合技术进行饲料颗粒气力输送运动特性的模拟仿真是可行的。

      5) 后续研究中,将着重讨论饲料尺寸、管道尺寸、管道弯曲程度、风速、风压等工艺参数对颗粒饲料运动特性的影响规律,获得管道尺寸、输送参数和饵料特性的最优化方案。

参考文献 (22)

目录

    /

    返回文章
    返回