向闯1,2,李莉莉1,2,王琨1,2,赵庆良1,2,姜珺秋1,2,魏亮亮1,2,汤云榕1,2
(1.城市水资源与水环境国家重点实验室(哈尔滨工业大学),哈尔滨 150090; 2.哈尔滨工业大学 环境学院,哈尔滨 150090)
摘要:
为研究厨余垃圾固态发酵反应器中机械搅拌装置的搅拌流场特性及优选机械搅拌设计参数,采用计算流体力学技术(CFD)方法模拟5种机械搅拌桨型的混合流场,并对搅拌桨的桨径比、转速等参数通过搅拌功率、混合时间、混合能和死区百分比等指标进行定量分析.模拟结果表明:二折叶搅拌桨、六折叶涡轮搅拌桨死区百分比大,对物料扰动范围小,不能有效进行混合;双螺带搅拌桨死区百分比最小,对物料扰动范围最大,搅拌功率也最大.对于双螺带搅拌桨桨径,物料死区随桨径比的增加而减小,当桨径比增大至0.75后不再减小,而搅拌功率随桨径比增大而一直增大.对于双螺带搅拌桨,随着转速增加,搅拌功率呈直线式增加,混合时间以抛物线式降低,混合能逐渐增加.综合考虑混合效果和搅拌功耗,转速为20 r/min、桨径比为0.75的双螺带搅拌桨是本研究反应器中的最佳搅拌设计.
关键词: 固态发酵 非牛顿流体 双螺带搅拌桨 死区 数值模拟 计算流体力学
DOI:10.11918/201911179
分类号:S216.4
文献标识码:A
基金项目:国家重点研发计划项目(2018YFC1900902)
Numerical analysis of flow field and optimization of agitation parameters in solid state fermentation reactor
XIANG Chuang1,2,LI Lili1,2,WANG Kun1,2,ZHAO Qingliang1,2,JIANG Junqiu1,2,WEI Liangliang1,2,TANG Yunrong1,2
(1.State Key Lab of Urban Water Resource and Environment(Harbin Institute of Technology), Harbin 150090, China;2.School of Environment, Harbin Institute of Technology, Harbin 150090, China)
Abstract:
To investigate the characteristics of the flow field of high solids in solid state fermentation reactor with mechanical agitators and optimize the design parameters of mechanical agitators, computational fluid dynamics (CFD) method was used to simulate the flow field of five types of mechanical agitators, and parameters such as ratio of impeller diameter and stirring speed were analyzed quantitatively through the indexes of stirring power, mixing time, mixing energy, and the percentage of dead zone. Simulation results show that the high solids had a small range of disturbance and a large percentage of dead zone, indicating that it cannot be effectively mixed through double blade and turbine agitator. The high solids had the largest range of disturbance, the smallest percentage of dead zone, and the biggest stirring power through double helical agitator. For the impeller diameter of the double helical agitator, the dead zone decreased with the increase of the ratio of impeller diameter, and when the ratio of impeller diameter increased to 0.75, it no longer decreased, while the stirring power increased with the increase of the ratio of impeller diameter all along. With regard to the stirring speed of the double helical agitator, the stirring power increased linearly with the increase of the stirring speed, the mixing time decreased along a parabolic curve, and the mixing energy gradually increased. Considering a balance of mixing performance and mixing power, the double helical agitator with a stirring speed of 20 r/min and the ratio of impeller diameter of 0.75 was the best agitator design in this reactor.
Key words: solid state fermentation non-Newtonian fluid double helical agitator dead zone numerical simulation CFD
向闯, 李莉莉, 王琨, 赵庆良, 姜珺秋, 魏亮亮, 汤云榕. 固态发酵反应器流场数值分析及搅拌参数优选[J]. 哈尔滨工业大学学报, 2021, 53(5): 65-71. DOI: 10.11918/201911179.
![复制到剪切板](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/copy.gif)
XIANG Chuang, LI Lili, WANG Kun, ZHAO Qingliang, JIANG Junqiu, WEI Liangliang, TANG Yunrong. Numerical analysis of flow field and optimization of agitation parameters in solid state fermentation reactor[J]. Journal of Harbin Institute of Technology, 2021, 53(5): 65-71. DOI: 10.11918/201911179.
![复制到剪切板](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/copy.gif)
基金项目 国家重点研发计划项目(2018YFC1900902) 作者简介 向闯(1994—),男,硕士研究生;
赵庆良(1962—),男,教授,博士生导师;
魏亮亮(1985—),男,副教授,博士生导师 通信作者 王琨, wang02kun@163.com 文章历史 收稿日期: 2019-11-20
Abstract Full text Figures/Tables PDF
固态发酵反应器流场数值分析及搅拌参数优选
向闯1,2, 李莉莉1,2, 王琨1,2
![通信作者](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/REcor.gif)
![wang02kun@163.com](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/REemail.gif)
1. 城市水资源与水环境国家重点实验室(哈尔滨工业大学), 哈尔滨 150090;
2. 哈尔滨工业大学 环境学院, 哈尔滨 150090
收稿日期: 2019-11-20
基金项目: 国家重点研发计划项目(2018YFC1900902)
作者简介: 向闯(1994—),男,硕士研究生; 赵庆良(1962—),男,教授,博士生导师; 魏亮亮(1985—),男,副教授,博士生导师
通信作者: 王琨, wang02kun@163.com
摘要: 为研究厨余垃圾固态发酵反应器中机械搅拌装置的搅拌流场特性及优选机械搅拌设计参数,采用计算流体力学技术(CFD)方法模拟5种机械搅拌桨型的混合流场,并对搅拌桨的桨径比、转速等参数通过搅拌功率、混合时间、混合能和死区百分比等指标进行定量分析.模拟结果表明:二折叶搅拌桨、六折叶涡轮搅拌桨死区百分比大,对物料扰动范围小,不能有效进行混合;双螺带搅拌桨死区百分比最小,对物料扰动范围最大,搅拌功率也最大.对于双螺带搅拌桨桨径,物料死区随桨径比的增加而减小,当桨径比增大至0.75后不再减小,而搅拌功率随桨径比增大而一直增大.对于双螺带搅拌桨,随着转速增加,搅拌功率呈直线式增加,混合时间以抛物线式降低,混合能逐渐增加.综合考虑混合效果和搅拌功耗,转速为20 r/min、桨径比为0.75的双螺带搅拌桨是本研究反应器中的最佳搅拌设计.
关键词: 固态发酵 非牛顿流体 双螺带搅拌桨 死区 数值模拟 计算流体力学
Numerical analysis of flow field and optimization of agitation parameters in solid state fermentation reactor
XIANG Chuang1,2, LI Lili1,2, WANG Kun1,2
![Corresponding author](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/REcor.gif)
![wang02kun@163.com](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/REemail.gif)
1. State Key Lab of Urban Water Resource and Environment(Harbin Institute of Technology), Harbin 150090, China;
2. School of Environment, Harbin Institute of Technology, Harbin 150090, China
Abstract: To investigate the characteristics of the flow field of high solids in solid state fermentation reactor with mechanical agitators and optimize the design parameters of mechanical agitators, computational fluid dynamics (CFD) method was used to simulate the flow field of five types of mechanical agitators, and parameters such as ratio of impeller diameter and stirring speed were analyzed quantitatively through the indexes of stirring power, mixing time, mixing energy, and the percentage of dead zone. Simulation results show that the high solids had a small range of disturbance and a large percentage of dead zone, indicating that it cannot be effectively mixed through double blade and turbine agitator. The high solids had the largest range of disturbance, the smallest percentage of dead zone, and the biggest stirring power through double helical agitator. For the impeller diameter of the double helical agitator, the dead zone decreased with the increase of the ratio of impeller diameter, and when the ratio of impeller diameter increased to 0.75, it no longer decreased, while the stirring power increased with the increase of the ratio of impeller diameter all along. With regard to the stirring speed of the double helical agitator, the stirring power increased linearly with the increase of the stirring speed, the mixing time decreased along a parabolic curve, and the mixing energy gradually increased. Considering a balance of mixing performance and mixing power, the double helical agitator with a stirring speed of 20 r/min and the ratio of impeller diameter of 0.75 was the best agitator design in this reactor.
Keywords: solid state fermentation non-Newtonian fluid double helical agitator dead zone numerical simulation CFD
厌氧发酵产沼气是最有前景的垃圾处理与资源化方法之一.厌氧发酵根据含固率的不同分为湿式和干式工艺.其中物料含固率在20%~30%的称为干发酵,也叫固态发酵.与传统的湿式发酵相比,干发酵可用于处理秸秆、生活垃圾和畜禽粪便等固体的有机废物,具有原料利用范围广、容积负荷大、污水产生量低、更加节能降耗、节省基建面积等优势[1-2].但是干发酵物料含固率高,属于高黏度假塑性流体,混合极为困难,大大降低了酶、微生物、底物之间的传质效率[3-5].搅拌有助于优化传热传质,但针对不同流体的具体搅拌方式与搅拌强度上仍有许多争议[6].Wu等通过模拟证明对于低黏度非牛顿流体,气动搅拌混合效率较高[7-8];Karim等通过实验证明低强度机械搅拌是对含固率为15%的牛粪最有效的搅拌方式[9].为了评价搅拌效果和优化搅拌装置设计,反应器流场研究通常是必须的,但传统实验如激光多普勒测速、粒子示踪法等限制条件较多且成本高昂[10-11].计算液体力学(CFD)方法基于有限容积法能对流场精确模拟计算,且成本低廉,越来越受到重视而广泛用于搅拌装置设计的研究中.Wu等[12]建立了畜禽粪便非牛顿流体的数学模型,并模拟对比了6种搅拌桨型的高含固流场特性.Yu等[13]比较了A-310和螺带搅拌桨对TS < 5%和TS=10%物料的搅拌特性,证明了螺带搅拌桨搅拌高含固物料的优越性.曹亚东等[14]考察了卧式干发酵反应器中搅拌转速、桨叶直径和搅拌桨层数等因素对流场和混合时间的影响.但是针对高径比较大的垂直流高含固反应器流场的研究仍极为有限.
针对厨余垃圾垂直流固态发酵反应器,物料含固率为25%,采用CFD方法模拟比选不同搅拌桨型的搅拌流场,并针对适用于高黏度流体的桨型[15-16],通过改变搅拌桨桨径比和转速,分析流场特性、搅拌功率、混合时间、混合能和死区百分比等因素,确定了最佳设计参数,为实际生产中搅拌系统的设计提供了借鉴.
1 试验 1.1 物理模型本研究反应器借鉴了成熟的已商业化运行的厨余垃圾固态发酵反应器[17],预先设计了如图 1所示的实验室规模反应器,该反应器高H为600 mm,直径D为200 mm,下部卸料锥角为45°.反应器所搭配的搅拌桨为3层二折叶搅拌桨,具有一定的径向和轴向混合作用,有效防止物料下沉,搅拌桨直径DJ为100 mm,3层桨桨间距离为200 mm.基于反应器所构建的物理模型如图 2所示.为了研究不同类型搅拌桨的混合效果,针对高含固非牛顿流体特征与反应器高径比较大的特点,根据行业标准(HG-T 3796.1~12—2005)另外构建了六折叶涡轮搅拌桨、螺杆搅拌桨、单螺带与双螺带搅拌桨与二折叶搅拌桨进行对比分析,搅拌桨设计参数如表 1所示.
Fig. 1
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-1.jpg)
Fig. 2
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-2.jpg)
表 1
![点击查看表格内容](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/table-icon.gif)
二折叶式 100 4 2 2 2
六折叶涡轮式 100 6 2 6 2
螺杆式 100 20 2 1 8
单螺带式 190 6 2 1 2
双螺带式 190 6 2 2 2
表 1 5种搅拌桨设计参数 Tab. 1 Design parameters of five agitators
1.2 网格划分搅拌桨结构较为复杂,手动划分结构化网格较困难且费时费力,因此,采用非结构化四面体网格对反应器进行网格划分.非结构化网格内部没有规则的拓扑关系,内部毗邻单元均不相同,对不同模型适应性强,通过网格独立性测试,以3层双折叶搅拌桨反应器为例,网格划分如图 3所示,网格划分数量为477 783,节点数量为89 348.
Fig. 3
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-3.jpg)
1.3 模拟方法对于搅拌反应器问题,Fluent提供了多种计算方法,包括多重参考系(MRF)、滑移网格与动网格模型.其中,MRF的基本思路是将搅拌桨所在的旋转区域使用运动参考系进行模拟,反应器无搅拌桨的区域成为静止区域,使用静止参考系模拟.旋转区域和静止区域通过interface进行质量动量交换.MRF计算速度快,搅拌运动模拟准确性高,因此,本文使用MRF方法.
高含固厨余垃圾属于假塑性非牛顿流体,黏度随着剪切速率的增加而减小.其本构方程为
$\eta = \kappa \cdot {\gamma ^{n - 1}}.$ (1)
式中:η为流体表观黏度,Pa·s;κ为稠度系数;γ为剪切速率,s-1;n为流变指数.参照文献将物料流变性质设为κ=35.4,n=0.29 [18].
反应器内雷诺数决定了模拟过程湍流模型的选择.当雷诺数小于1 000时,流体处于层流区,应选择层流模型,反之应选择标准k-ε、RNG k-ε湍流方程等湍流模型.搅拌反应器内雷诺数可通过式(2)进行计算:
$Re = \frac{{\rho N{d^2}}}{\eta }.$ (2)
式中:ρ为流体密度,kg/m3;N为搅拌桨转速,r/min;d为搅拌桨直径,mm.
使用商业软件Fluent17.0,基于压力求解器,采用稳态MRF模型、层流模型和标准k-ε、RNG k-ε湍流模型,并使用收敛性好、计算速度快的Simple算法对流场进行求解,当各项指标残差达到1×10-4时认为计算收敛,求解完成.
Fluent中使用的连续性方程、动量方程以及湍流方程如下:
连续性方程
$\frac{{\partial \rho }}{{\partial t}} + \frac{\partial }{{\partial {x_i}}}\left( {\rho {u_i}} \right) = 0.$ (3)
动量方程
$\frac{\partial }{{\partial t}}\left( {\rho {u_i}} \right) + \frac{\partial }{{\partial {x_j}}}\left( {\rho {u_i}{u_j}} \right) = {F_i} - \frac{{\partial \rho }}{{\partial {x_i}}} + \frac{{\partial {\tau _{ij}}}}{{\partial {x_{ij}}}} + \rho {g_i}.$ (4)
标准k-ε湍流方程
$\begin{array}{l}\frac{\partial }{{\partial t}}(\rho \kappa ) + \frac{\partial }{{\partial {x_i}}}\left( {\rho \kappa {u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left( {{\sigma _\kappa }{\mu _{{\rm{eff}}}}\frac{{\partial \kappa }}{{\partial {x_j}}}} \right) + \\{G_\kappa } + {G_{\rm{b}}} - \rho \varepsilon - {Y_{\rm{M}}} + {S_\kappa }, \end{array}$ (5)
$\begin{array}{l}\frac{\partial }{{\partial t}}(\rho \varepsilon ) + \frac{\partial }{{\partial {x_i}}}\left( {\rho \varepsilon {u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left( {{\sigma _\xi }{\mu _{{\rm{eff}}}}\frac{{\partial \xi }}{{\partial {x_j}}}} \right) + \\{C_{1\xi }}\frac{\xi }{\kappa }\left( {{G_\kappa } + {C_{3\xi }}{G_{\rm{b}}}} \right) + {C_{2\xi }}\rho \frac{{{\xi ^2}}}{\kappa } - {R_\xi } - S.\end{array}$ (6)
式中:ui为湍动黏度,可表示为k和ε的函数;Gк为平均速度梯度引起的湍流动能产生项; Gb为浮力引起的湍流动能产生项; YM为可压缩湍流中的波动扩张对整体耗散率的贡献; Sκ与S为用户定义的原项;模型常量C1ξ=1.44,C2ξ=1.92,σκ=1.0,σξ=1.3.
1.4 分析指标 1.4.1 搅拌功率搅拌功率是指在单位时间内搅拌反应器消耗的能量,该能量提供给反应器内物料进行剪切和破坏.影响搅拌功率的因素主要包括搅拌装置的结构参数、转速、固态发酵反应器的构型和反应器物料的流变特性等.计算搅拌对搅拌桨设计和电动机选择意义突出.搅拌功率的计算如式(7)所示,其中扭矩M可以在Fluent模拟计算后得到[19].
$P = M\omega = M \times \frac{{2\pi n}}{{60}}.$ (7)
式中:P为搅拌功率,W;M为扭矩,N·m;ω为角速度,rad/s;n为转速,r/min.
1.4.2 混合时间混合时间是物料达到规定混合程度所需要的时间,反映了物料混合的快慢[20].均匀度M(T)小于5%时消耗的时间即为混合时间,M(T)的定义为
$M(T) = \frac{{\left| {{c_i} - {c_{{\rm{average }}}}} \right|}}{{{c_{{\rm{average }}}}}}.$ (8)
式中:T为时间,ci为在某一计算域内的示踪剂浓度,caverage为反应器内示踪剂的平均浓度.
在Fluent中可通过Patch定义进料球体,将球体内物料初始浓度百分比设为1,球体外初始浓度百分比为0,通过欧拉模型模拟球体内物料进入反应器后的混合过程.同时在反应器内设定监测点,当监测点物料浓度稳定达到规定浓度即认为混合完成,此时所消耗时间即为混合时间.
1.4.3 混合能混合能是评价混合效率的综合性指标,综合考虑了搅拌功率和混合时间的影响.混合能越大,通常意味着混合效率越低.混合能计算如下:
$Q = P \times T.$ (9)
式中Q为混合能,J.
1.4.4 死区百分比反应器内各部位物料速度具有不均一性,根据王令闪等[21]的研究将反应器内流体速率低于流体最大速率0.01的区域定义为死区.将死区百分比定义为死区物料体积占反应器物料总体积的比例,该区域属于搅拌动能的死角,死区物料基本不进行动量和能量交换,容易结壳黏附在反应器壁,造成反应器有效区域大大减少甚至堵塞等问题,因此,必须通过搅拌器设计尽可能减少死区体积.
2 结果与讨论 2.1 搅拌桨选型 2.1.1 宏观流场针对大高径比反应器中的高黏度非牛顿流体,选取既产生径向流又能产生轴向流的混流型搅拌器:二折叶搅拌桨、六折叶涡轮搅拌桨、螺杆搅拌桨、单螺带搅拌桨、双螺带搅拌桨进行分析比较.在搅拌转速均为20 r/min的条件下,Fluent求解生成的轴向截面速度场和矢量速度场分别如图 4,5所示.显然,二折叶搅拌桨与六折叶涡轮式搅拌桨由于其较小的桨叶面积,仅仅在桨叶周围的流体受到扰动,在剪切力的作用下黏度减小,能够实现桨叶局部物料的循环与充分混合;对于远离桨叶的流体,未受到强烈的扰动,流体依然保持着高黏度,难以实现运动,整体来看,二折叶搅拌桨与六折叶涡轮式搅拌桨对高黏度非牛顿流体的混合效果非常不理想.对于螺杆搅拌反应器,流体在搅拌桨的轴流作用下开始整体向上运动,随后由于桨叶的挤压作用沿反应器壁向下运动,完成反应器内的整体循环过程,但是在远离搅拌桨的上部和下部,流体未受到扰动,处于静止状态.对于单螺带搅拌反应器,流体在桨的带动下形成了局部涡流,而未能形成整体循环.当增加一条螺带后,在桨局部形成涡流的基础上,流体整体自下而上运动,随后沿着反应器壁向下,不断进行循环.从宏观流场运动来看,得益于与物料较大的接触面积,螺杆搅拌桨与双螺带搅拌桨能够带动物料在反应器内整体实现循环,因而混合效果较好.
Fig. 4
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-4.jpg)
Fig. 5
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-5.jpg)
2.1.2 死区百分比各反应器在轴向截面的死区分布如图 6所示,红色区域代表速度大于0.01um的有效混合区,红色区域外的部分即代表未能充分混合的死区.通过模拟计算,各搅拌桨的死区百分比如表 2所示.由图 6和表 2可以看出,二折叶搅拌桨与六折叶涡轮式搅拌反应器内大部分区域属于死区,死区百分比分别达到了80.9%和71.4%,因此不能使物料得到充分混合.螺杆搅拌反应器内,在反应器上下部以及器壁存在少量死区,死区百分比为24.5%.而在螺带搅拌反应器中仅有反应器上下部存在少量死区,其中单螺带搅拌桨的死区百分比为23.9%,双螺带搅拌桨死区百分比最小,为19.9%.混合效果明显优于其他搅拌装置.
Fig. 6
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-6.jpg)
表 2
![点击查看表格内容](http://hit.alljournals.cn/html/hitxb_cn/2021/5/images/table-icon.gif)
a 80.9 5.381
b 71.4 16.225
c 24.5 31.133
d 23.9 72.213
e 19.9 119.771
表 2 不同搅拌桨型死区百分比、搅拌功率模拟计算 Tab. 2 Simulation calculation of the percentage of dead zone and stirring power
2.1.3 搅拌功率通过Flunet模拟计算得到的各搅拌桨消耗的功率见表 2.不难发现,由于二折叶搅拌桨与六折叶涡轮式搅拌反应器结构相对简单,桨叶叶片数量少,与流体接触面积小,搅拌时阻力较小,消耗功率也较小,单位体积搅拌功率分别为5.381和16.225 W/m3.而螺杆与螺带搅拌桨结构复杂,与物料接触面大,搅拌时阻力大,消耗功率也较大,其中双螺带搅拌桨消耗功率最大,达119.771 W/m3.但是由于反应器间歇运行,为了实现充分混合,在可接受范围内较高的搅拌功率是能接受的.
2.2 搅拌桨叶直径为了研究桨叶直径对混合物料的影响,在搅拌速度恒定为20 r/min的情况下,定义桨径比为桨叶直径与反应器直径的比,预设桨径比分别为0.55,0.65,0.75,0.85和0.95,则桨叶直径对应110,130,150,170和190 mm.通过Fluent求解,得到如图 7和8所示流场速度云图和矢量图.从轴向看,反应器远离搅拌桨的上部和下部均未得到扰动,因而不能进行混合.从径向看,桨叶周围的物料因剧烈搅拌速度最大,而且在桨叶周围形成强大的涡流,与桨叶距离越远的物料混合效果越差,物料在整体上沿着搅拌桨壁形成了由下而上到达反应器顶部、再沿反应器壁向下到达反应器底部的循环.桨径比为0.55时,在反应器壁的物料没有得到较好的混合,容易造成物料在器壁结壳的现象.但增大桨叶直径后,对物料的扰动范围也增大,混合效果越好.当桨径比为0.95时,反应器壁物料也能得到较好的混合.
Fig. 7
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-7.jpg)
Fig. 8
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-8.jpg)
反应器内死区百分比情况如图 9所示.桨径比为0.55和0.65时,反应器壁仍存在少量死区,此处物料完全得不到扰动,未参与混合.当桨径比增大到0.75以后,反应器壁几乎不存在死区,混合效果优越.反应器搅拌功率如图 10所示.随着桨径比增大,与物料接触面也越大,搅拌功率呈直线式递增.综上,当桨径比为0.75时,反应器壁物料能够得到扰动而进行混合,且此时搅拌功率相比桨径比为0.85和0.95时更小,能耗更少.因此,将桨径比定为0.75,此时搅拌桨直径为150 mm.
Fig. 9
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-9.jpg)
Fig. 10
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-10.jpg)
2.3 搅拌转速为了确定所选双螺带搅拌桨的最佳转速,在桨叶直径为150 mm的情况下,将转速分别设为5,10,20,50和100 r/min.通过Fluent模拟流场,得到如图 11所示轴向截面速度云图和图 12的轴向截面死区分布图.不难发现,由于搅拌桨结构不变,增加搅拌转速只是增加了搅拌强度,而不会改变物料流动结构.在死区方面,即使是在低转速条件下,除了上、下部分的死区,仅能避免存在一定的死区,随着搅拌转速增加,壁面部分死区逐渐减小;当转速增加到20 r/min以后,壁面部分已消除死区.死区是搅拌死角,即使转速达到极端情况,也只能使有效混合区内动能更加充沛,而无法减小死区空间,这一部分要改善流场特性必须通过优化反应器和搅拌装置的设计.图 12中反应器上、下部的死区,可在今后考虑增加搅拌桨高度进行优化.
Fig. 11
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-11.jpg)
Fig. 12
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-12.jpg)
另外,在不同转速下模拟得到的搅拌曲线如图 13所示,可以看出,搅拌功率随着转速增加呈直线式递增,而混合时间随着搅拌转速以抛物线形式降低.在搅拌转速为0~20 r/min时混合时间迅速降低,在20 r/min以后降低幅度很小.而混合能随着转速的增加而增加.由于混合能是搅拌功率与混合时间的乘积,转速为5~50 r/min时,混合能以抛物线式增加.不难预测,随着混合时间趋于平稳,混合能在100 r/min以后将直线式递增.虽然搅拌转速越高,物料的混合效果越好,能有效促进传质,但针对高黏度非牛顿流体,高转速意味着巨大的搅拌阻力和功耗,实际生产中没有电机能够支撑这一搅拌转速.同时,高转速带来的高剪切速率会破环微生物絮凝体结构,高转速下充分混合反而造成了反应器沼气产量的下降.为了降低系统的成本,Rivard等[22-23]尝试寻找满足有效消化率的混合要求,将搅拌器转速设为1和25 r/min的实验结果表明,两者在发酵性能上没有显著差异,证明在1 r/min条件下,消化池内的微生物和亚硝酸盐也会进行迁移.因此,低强度间歇混合是确保反应器成功运行的关键.综合考虑本文混合效果与搅拌功耗的研究,20 r/min为目前最佳转速,但针对本文反应器找到合适的间歇搅拌时间,仍有待今后继续研究混合条件与产气量之间的关系.
Fig. 13
![点击查看原图](http://hit.alljournals.cn/html/hitxb_cn/2021/5/PIC/hebgydxxb-53-5-65-13.jpg)
3 结论1) 针对TS=25%的高含固非牛顿流体,在5种搅拌桨型中,从宏观流场看,双螺带搅拌桨带动物料整体自下而上运动,随后沿着反应器壁向下进行循环;从死区看,双螺带搅拌桨死区百分比最小,仅为19.9;从搅拌功率看,双螺带搅拌桨搅拌功率最大为119.771 W/m3.为了确保物料的充分混合,综合来看,双螺带搅拌桨最符合要求.
2) 对于双螺带搅拌桨,增大反应器桨径比,物料死区百分比逐渐减小,混合效果变好.同时搅拌功率也随反应器桨径比呈直线式上升.当桨径比为0.75即直径为150 mm时,死区百分比最小,既保证了物料的有效混合,同时搅拌功率较小.
3) 对于双螺带搅拌桨,随着搅拌转速的增加,物料流动特征基本未发生变化,扰动范围并未增加,只是混合区内物料运动速度增加.而搅拌功率随着搅拌转速增加呈直线式增长,混合时间呈抛物线式下降,混合能逐渐增加.综合考虑混合效果与搅拌功耗,搅拌转速20 r/min为最佳转速.
参考文献
[1] 吴健, 华银锋, 张海涛, 等. 高含固餐厨垃圾半干法厌氧发酵系统运行分析[J]. 科学技术与工程, 2018, 18(35): 94.
WU Jian, HUA Yinfeng, ZHANG Haitao, et al. Operation analysis of semi-dry anaerobic digestion system of high solid content food waste[J]. Science Technology and Engineering, 2018, 18(35): 94. DOI:10.3969/j.issn.1671-1815.2018.35.014
[2] 叶小梅, 常志州. 有机固体废物干法厌氧发酵技术研究综述[J]. 生态与农村环境学报, 2008, 24(2): 76.
YE Xiaomei, CHANG Zhizhou. State of arts and perspective of dry anaerobic digestion of organic solid waste[J]. Journal of Ecology and Rural Environment, 2008, 24(2): 76. DOI:10.3969/j.issn.1673-4831.2008.02.017
[3] LEBRANCHU A, DELAUNAY S, MARCHAL P, et al. Impact of shear stress and impeller design on the production of biogas in anaerobic digesters[J]. Bioresource Technology, 2017, 245(Part A): 1139. DOI:10.1016/j.biortech.2017.07.113
[4] MAO Liwei, ZHANG Jingxin, DAI Yanjun, et al. Effects of mixing time on methane production from anaerobic co-digestion of food waste and chicken manure: Experimental studies and CFD analysis[J]. Bioresource Technology, 2019, 294: 122177. DOI:10.1016/j.biortech.2019.122177
[5] 戴晓虎, 盖鑫, 董滨. 高含固厌氧消化污泥流变特性[J]. 环境工程学报, 2014, 8(9): 3912.
DAI Xiaohu, GAI Xin, DONG Bin. Rheological characteristics of sludge in high-solids anaerobic digestion[J]. Chinese Journal of Environmental Engineering, 2014, 8(9): 3912.
[6] 黄如一, 赵鑫, 李江, 等. 利用CFD方法优化沼气发酵罐内流场形态的研究综述[J]. 中国沼气, 2018, 36(4): 23.
HUANG Ruyi, ZHAO Xin, LI Jiang, et al. A review on optimization of flow pattern in biogas anaerobic digester by CFD method[J]. China Biogas, 2018, 36(4): 23. DOI:10.3969/j.issn.1000-1166.2018.04.005
[7] WU Binxin. Computational fluid dynamics investigation of turbulence models for non-Newtonian fluid flow in anaerobic digesters[J]. Environmental Science & Technology, 2010, 44(23): 8989. DOI:10.1021/es1010016
[8] WU Binxin. CFD simulation of gas and non-Newtonian fluid two-phase flow in anaerobic digesters[J]. Water Research, 2010, 44(13): 3861. DOI:10.1016/j.watres.2010.04.043
[9] KARIM K, HOFFMANN R, KLASSON T, et al. Anaerobic digestion of animal waste: Waste strength versus impact of mixing[J]. Bioresource Technology, 2005, 96(16): 1771. DOI:10.1016/j.biortech.2005.01.020
[10] DAPELO D, BRIDGEMAN J. Assessment of mixing quality in full-scale, biogas-mixed anaerobic digestion using CFD[J]. Bioresource Technology, 2018, 265: 480. DOI:10.1016/j.biortech.2018.06.036
[11] DAPELO D, ALBERINI F, BRIDGEMAN J. Euler-Lagrange CFD modelling of unconfined gas mixing in anaerobic digestion[J]. Water Research, 2015, 85: 497. DOI:10.1016/j.watres.2015.08.042
[12] WU Binxin. CFD simulation of mixing for high-solids anaerobic digestion[J]. Biotechnology and Bioengineering, 2012, 109(8): 2116. DOI:10.1002/bit.24482
[13] YU Liang, MA Jingwei, CHEN Shulin. Numerical simulation of mechanical mixing in high solid anaerobic digester[J]. Bioresource Technology, 2011, 102(2): 1012. DOI:10.1016/j.biortech.2010.09.079
[14] 曹亚东, 李欣欣, 黄振峰, 等. 干式发酵搅拌器流场数值分析及结构参数优选[J]. 化学反应工程与工艺, 2018, 34(1): 11.
CAO Yadong, LI Xinxin, HUANG Zhenfeng, et al. Flow field numerical analysis of dry-fermentation stirrer and its structure parameter optimization[J]. Chemical Reaction Engineering and Technology, 2018, 34(1): 11. DOI:10.11730/j.issn.1001-7631.2018.01.0011.09
[15] 朱桂华, 柳颖娇, 徐洪威, 等. 高黏度污泥螺带螺杆搅拌混合特性的数值模拟[J]. 安全与环境学报, 2018, 18(3): 1121.
ZHU Guihua, LIU Yingjiao, XU Hongwei, et al. Numerical simulation for the mixture features of the ribbon-screw impeller with high viscosity sludge[J]. Journal of Safety and Environment, 2018, 18(3): 1121. DOI:10.13637/j.issn.1009-6094.2018.03.054
[16] 董敏, 夏晨亮, 李想. 双螺带及六斜叶组合桨搅拌槽内部流场及混合时间数值模拟[J]. 排灌机械工程学报, 2019, 37(1): 43.
DONG Min, XIA Chenliang, LI Xiang. Numerical simulation of internal flow field and mixing time in stirred tank with combined paddles[J]. Journal of Drainage and Irrigation Machinery Engineering, 2019, 37(1): 43. DOI:10.3969/j.issn.1674-8530.17.0160
[17] 胡鑫, 冯晶, 赵立欣, 等. 干法厌氧发酵反应器及过程控制技术研究进展[J]. 中国沼气, 2018, 36(2): 68.
HU Xin, FENG Jing, ZHAO Lixin, et al. Research progress of dry anaerobic fermentation reactors and its process control technique[J]. China Biogas, 2018, 36(2): 68. DOI:10.3969/j.issn.1000-1166.2018.02.012
[18] 刘刈. 高浓度物料沼气发酵过程混合搅拌及其影响因素的研究[D]. 北京: 中国农业科学院, 2009
LIU Yi. The study of mixing and strring for high solids biogas fermentation and its influence factors[D]. Beijing: Chinese Academy of Agricultural Sciences, 2009
[19] WU Binxin, KENNEDY S, ESHTIAGHI N, et al. CFD modeling of active volume creation in a non-Newtonian fluid agitated by submerged recirculating jets[J]. Chemical Engineering & Technology, 2018, 41(7): 1441. DOI:10.1002/ceat.201700323
[20] ZHAI Xiaodong, KARIYAMA I D, WU Binxin. Investigation of the effect of intermittent minimal mixing intensity on methane production during anaerobic digestion of dairy manure[J]. Computers and Electronics in Agriculture, 2018, 155: 121. DOI:10.1016/j.compag.2018.10.002
[21] 王令闪, 苏红军, 徐世艾. 高黏体系中最大叶片式搅拌桨直径的CFD优化[J]. 化学工程, 2011, 39(7): 9.
WANG Lingshan, SU Hongjun, XU Shiai. Optimization of maxblend impeller diameter in high viscous fluid by CFD[J]. Chemical Engineering, 2011, 39(7): 9. DOI:10.3969/j.issn.1005-9954.2011.07.003
[22] RIVARD C J, HIMMEL M E, VINZANT T B, et al. Development of a novel laboratory scale high solids reactor for anaerobic digestion of processed municipal solid wastes for the production of methane[J]. Applied Biochemistry and Biotechnology, 1989, 20(1): 461. DOI:10.1007/BF02936503
[23] RIVARD C J, HIMMEL M E, VINZANT T B, et al. Anaerobic digestion of processed municipal solid waste using a novel high solids reactor: Maximum solids levels and mixing requirements[J]. Biotechnology Letters, 1990, 12(3): 235. DOI:10.1007/bf01026806