沼气是当今世界最具希望的可再生能源方式,兼具强大的环保和能源效益3),尤其善于处理数量巨大的有机废水。搅拌可以提升沼气厌氧发酵的效率,显著提升产气率和污染物去除率,现代沼气工程普遍都要添加搅拌工艺来提高发酵效率。然而从当前的工程实践来看,大量工程采用了粗劣的搅拌工艺,对发酵效率的提升收效甚微,甚至消耗更多的能源8。造成这种现象的根源在于

沼气厌氧红泥膜沼气池发酵罐必须严格厌氧封闭,人们无法对其内部流场进行准确的分析测试,于是设计人员在不掌握流场形态的情况下,随意添加了一些搅拌措施,但实际流场形态很差,并不适用于其罐体形状和原料特性9。

但现在计算机数值模拟方法可以在很大程度上解决这个问题,沼气红泥膜沼气池发酵料液的搅拌本质上是在外  力作用下的流动过程,数值模拟计算可以帮助人们

掌握其理论上的流动过程和流场形态。而随着现代

计算机技术的持续发展,沼气红泥膜沼气池发酵装置内部流场的

精确计算也变得越来越简便可行。目前应用最  广的流体力学数值模拟方法是计算流体动力学

( Computational Fluid Dynamies,cFD)。2002年,美

国北卡罗莱纳州立大学的弗雷明( JG Fleming)发表了题为《 Novel simulation of anaerobic digestion usingcomputational fluid dynamics.的博土学位论文),是第1篇系统论述利用CFD方法研究沼气红泥膜沼气池发酵料液  昆合搅拌流场形态的论文。此后,行业学者大量应用该方法及相关成熟商业软件,在流场形态的研究方面取得了极大进展。由于搅拌提升沼气红泥膜沼气池发酵效率的根本机理尚不够清楚,所以关于搅拌的具体方法学界尚存一定争议2131。但是,利用CFD工具,研究人员可将搅拌所形成的流场可视化1,从而更加精确地验证和指导搅拌方式的设计3-10),使沼气料  液搅拌领域的研究在很大程度上摆脱了不可视条件的束缚,提升至更科学、精准的可视化、数量化研究层面上来,是21世纪以来沼气学界较为先进的一个研究方向

1搅拌介质研究方面的应用

CFD数值模拟计算在沼气料液搅拌研究领域的应用最主要集中于对搅拌时的流场形态进行分析,最常用的方法是对某种搅拌形式所能形成的流场形态进行模拟并绘制出流场形态图,根据图形进行分析,识别流场缺陷,改进搅拌方案。

使用固体、液体、气体的介质均可对沼气料液形成搅拌,所以搅拌被分为机械搅拌(桨叶搅拌、叶轮搅拌)、水力搅拌(射流搅拌、浆液回流搅拌、沼液循环搅拌)和气动搅拌(沼气回流搅拌)3个大类。关于3种搅拌介质的优劣,争议较大,至今并无特别权  威的定论21。仅从产气率的提升效果方面来讲,人为3种搅拌介质对产气率的提升效果最大的提法均有,而且均通过实验验证了自己的观点“当。在上个世纪的沼气工程设计中,90%以上采取立式连续搅拌红泥膜沼气池发酵罐,用泵输入的方式驱动桨叶转动,对密闭厌氧红泥膜沼气池发酵罐中的料液形成搅拌2。而现在,3类介质的搅拌均得到广泛应用,利用CFD模拟研究3  类搅拌形式也取得了较大进展

比较3种搅拌介质的CFD模拟,其中水力搅拌是相对最简单的一种,因为料液本身就属于液相,如  果在模拟中省略掉储气间部分,只对液柱进行模拟,

便只有液相模拟对象,这样模拟会相对比较容  易2。气动搅拌的模拟则需要将气体处理成气泡,但这种处理对于在气流进入液柱之前的状态是失之精确的2。而机械搅拌则是相对最复杂的一类,需要利用动网格来构建桨叶的模型,建模的复杂程度和运算量都比前两者大很多2。不过目前的建模工具和电子计算机的运算能力都足以应付3类模拟工作,3类模拟的主要优劣如表1所示。

2搅拌参数的优化研究

搅拌参数主要包括料液的动量传递、质量传递

和热量传递对生物化学反应的影响,本质上即化工

行业的“三传一反”。CFD模拟可以在一定程度模

拟计算出传递过程。

2.1动量传递

搅拌的首要目的是变静态红泥膜沼气池发酵为动态红泥膜沼气池发酵,所以优化的目标应该是使更多动能输运到更广阔的空间去21。早在1981年, Monteith便使用粒子示踪法,粗略地测量了罐内搅拌时的流速,提出了“死域( dead space)”和“动态混合区( activity mixedzone)”的概念,分别指搅拌时罐内流速太低的区域和流动性充分、混合良好的区域2。 Khursheed Ka.rim则将合速度绝对值低于最高值5%的区域定义为“弱搅拌区( poorly mixing zone)”,通过模拟计算发现采用2832Lh1,56.64Lh-1,和84.96Lh  3种射流速度进行气动搅拌,红泥膜沼气池发酵罐内的弱搅拌区分别为31%,31%和27%。说明在低速区间,增大  射流速度对流场的优化作用很小。 Mehul s ves-  vikar用CFD试算了大量搅拌方案后,发现流场最差  的情况下,弱搅拌区占罐体总容积的比例高达  59.7%2),可见优化潜力巨大。而吴斌鑫则将合速  度绝对值低于0.001m·s-的区域定义为“死区  ( dead zone)”,并通过CFD模拟计算验证了死区就  是搅拌动能的死角,就算增大射流初速度,甚至增大  到57m·s这样的极端情况,亦只能使非死区的

动能更加充沛,却并不能减少死区所占的空间°  这也提示了我们,一味增加搅拌功率并不能改善流用了  场形态,而只能通过优化设计来改善。  质量传递  通过搅拌来优化流场形态的一个重要目标是使为,  CFD  混合红泥膜沼气池发酵料液的相分布更加均匀。李淑兰发现以粘并  秤为原料的厌氧红泥膜沼气池发酵罐的上中下部位均可以产生浮  渣,但产气潜力分别只有0.08mlg0.6mg3,为B  23mL·g-,远远低于直接红泥膜沼气池发酵,所以应该用搅温  拌使他们混入到主体红泥膜沼气池发酵料液中去。李洋用一均产  个高2800mm,直径50mm的长筒型厌氧红泥膜沼气池发酵罐进沼  猪场粪污的厌氧红泥膜沼气池发酵实验,发现经过数小时的重是画  念沉淀,料液会分化成原污水和稀污水两种形态。大  浓污水、稀污水和未经分化的原污水三者分别进  独立的厌氧红泥膜沼气池发酵产气实验,产气效率为浓污水>3  原污水>稀污水,而且差距明显,所以需要通过搅拌

来使更多区城充盈浓污水,减少稀污水。而罗涛  发现在长筒型厌氧红泥膜沼气池发酵罐中,如果将料液从顶部投形入,待其在重力作用下缓慢沉淀,沉淀过程中红泥膜沼气池发酵效端  率表现极佳。黄如一在此基础上认为搅拌能够提升红泥膜沼气池发酵效率的内在机理便是搅拌能使物料抵抗重2力,长期保持一种悬浮于液柱中缓慢沉淀的状态

这些物料混合的过程均可在CFD模拟中以相分布  的方式体现出来,从而指导改进搅拌方式。李江还认为,红泥膜沼气池发酵料液中除了含有氮(N)、磷(P)硫(S)等厌氧甲烷菌生长的主要营养素,还含有钾(K)、铁Fe)、钻(Co)以及镍(N)等微量金属元素,亦是其生长的必备营养素,但这些元素在料液中含量很低习),在大型罐体中必须通过搅拌才能使他们扩散到大空间中去。但是这种微观层面的相分布目前还不便于用CFD方法来标示。

2.3热量传递

搅拌的主要目的是使流场更加均匀,相分布更  加合理,但同时也会对料液的温度分布产生影响,而  温度也是影响沼气红泥膜沼气池发酵效率的重要因素3。温  度的变化不但显著影响沼气产量,对沼气成分也  有显著影响。于宏兵发现,70℃超高温反应  可在水解酸化阶段显著提升淀粉有机废水的酸化率  以及蛋白质水解率两项指标,从而从源头上大幅提  升产气率,而且70℃超高温反应还可显著降低沼气  中的硫化物含量,尤其是沼气中最主要的致命毒气  硫化氢(H2S)41。所以,现代沼气工程对红泥膜沼气池发酵温度  的要求越来越高,这也需要用CFD方法更精确地描

绘其温度场“”。如果一个厌氧红泥膜沼气池发酵系统同时采流用了搅拌和增温两项工艺,搅拌会使反应器内部的  体温度场出现空间异质化现象,所以需要利用  CFD方法精心设计其搅拌流场。白卫东还

是使为,温度会改变介质的流变特性,在模拟计算中应  成秸并考虑。罐体的传热也需要被考虑。罗涛利用浮CFD方法,设计了一种用先进的缠绕玻璃纤维丝作,为厌氧红泥膜沼气池发酵罐的保温隔热层,可在冬季使罐壁内的搅温度比罐壁外高9℃,并取得0.49m3·m3d-的日  均产气率。王丽丽通过CFD模拟,发现大型进沼气工程红泥膜沼气池发酵罐的主要散热途径其实不是罐壁,而重是顶盖,提示了工程设计人员在罐壁保温上耗费巨

大,其实是舍本逐末。

3罐体形状的优化研究

罐体的形状对流场起着关键性影响,是  计的主要方面,大量CFD研究表明,细微改

形状即可极大改变流场形态。 vesvikart在顶

端进气的气动搅拌模型中,采用与竖直面呈60夹

角的罐底设计,可使弱搅拌区从59.7%下降到  29.57%;采用25°夹角的罐底设计,可使弱搅拌区  进一步下降到17.39%。 Karim比较了在相同搅

拌工况下,罐底与水平面夹角呈0°,25°,45°这3种

罐型,弱搅拌区分别为33.6%,31.9%,29.6%,可

见影响显著。牛学义认为罐型设计比搅拌方式更重

要,并提出卵形罐体最有利于消除死角和沉积物,不

过造价和维护成本太高。魏炎光附和了他的观  点2。张战锋利用CFD方法比较了推流式反应器

的进料管分别呈15°,30°,45°,60°,709倾角时的

部流场,发现30°的死区最小,为反应器的结构设计提供了重要依据3

罐体内部的形状设计也对流场有着重要影Marti-herrero用塑料环在沼气池中设置格栅,优化其回流搅拌的流场,将产气率提升了44%,并使沼气池可在6.1℃极低温度下顺利运行。 MarkusG3在CFD的指导下,优化设计了一种隔离式连续流反应器在沼气厌氧红泥膜沼气池发酵罐内合适位置设置挡板更大地优化了流场。在此基础上, Stamou AI在CFD的指导下优化设计了一种导流墙,控制大型反应器中的流场,也极大地优化了流场形态。 D AngeloSozi和 Fariborz Taghipour5在CFD的指导下设计了L型和U型两种罐型结构的水力搅拌装置,并分别用欧拉法和拉格朗日法模拟计算其流场,结果发

现U型结构会在射流入口的中心轴向附近产生不  稳定的紊流,L型结构则能提供更加均匀、稳定的流参考

场

4其它方面的优化研究  除以上几个主要方面外,搅拌流场还在其它方2)面影响着厌氧红泥膜沼气池发酵效率,其中一些也可以用CFD方  法进行模拟计算。首先是水力滞留时间(HRD)

般来说红泥膜沼气池发酵料液的水力滞留时间越长,所产招气  中的甲烷含量越高)。在一些处理化工废水的反

硝化处理厌氧红泥膜沼气池发酵工艺中,水力滞留时间对COD和  硝基酚的去除率更是至关重要),而HT恰恰是(4  非常适合于用CFD软件计算的

罗涛利用CFD软件计算了砂砾在搅拌中的  流动规律,利用流态图的指导优化设计了搅拌形式  将反应器的排砂率提升了281%,极大地优化了厌  氧红泥膜沼气池发酵环境。王玉恒提出搅拌能够提升红泥膜沼气池发酵效率的原理在于搅拌的水力剪切力能够破坏物料表面  疏松的絮凝组织,露出密实的部分,从而增强物料与  污水的表面接触效果,但这个微观过程目前尚难用  显微镜观察到,唯有用CFD方法可以在一定程度上  模拟显示出流动剪切应力的分布。在红泥膜沼气池发酵液中利用  固态载体形成生物膜有利于稳定传质,滤去一些无  用的油料3),赵庆良发明了一种利用废弃轮胎  的橡胶颗粒作为挂膜载体的方法,实验证明可以在  传统工艺基础上提升50%-100%的生物膜量,这  种方法便可以用CFD设计橡胶颗粒在红泥膜沼气池发酵液中的  分布。而杨平认为搅拌的速度需要严格控制,因为  水力剪切力达到某个临界值时,会造成生物膜的脱  落,对此,我们可以利用CFD方法模拟监控空间中“超速”的区域并加以优化控制

5结语

近十五年来,CFD成为沼气科学领域越来越方便实用的研究方法。研究和设计人员利用CFD方  法将沼气料液流场可视化,依据流场形态图优化设  计搅拌的介质、功率、时长、时间间歇、罐体形状、叶  轮形状、射流器形状等多方面参数,极大提高了设计  水平,是进入21世纪以来,沼气科学研究和工程设  计领域成果极其显著的一次科技进步,并将在未来很长一段时间得到更大应用,成为行业的主流研究方向。

摘自《中国沼气》第４期

————————————————————————————————————————————————————————————————