基于Fluent流场模拟的半干式连续厌氧发酵反应器的设计及中试试验

李 靖，张 重，赵国明，袁存亮，李爱芹，丁 伟

(吉林省农业机械研究院, 吉林 长春 130022)

随着畜禽养殖业的飞速发展，畜禽粪便的产生量也十分巨大，据统计，我国畜禽粪污每年所产生的COD的排放量是我国工业总COD排放量的5倍以上[1-2]，如何合理地将畜禽粪便无害化处理已经成为当下一种主要的研究方向。目前厌氧发酵处理工艺已被广泛应用于畜禽粪便的无害化处理，并且还可在处理畜禽粪便的同时产生沼气作为清洁能源而被广泛利用。常见的畜禽粪便的厌氧发酵处理工艺可分为湿式与干式两种[3]，近几年又将发酵物料TS在15%～25%的称为半干式厌氧发酵技术。基于环保与节能基础上，我国畜禽舍现多为干清粪，其物料TS浓度在10%～25%，水分含量少[4]，故而半干式厌氧发酵技术更适合于畜禽干清粪便的处理。

根据文献显示，厌氧发酵技术虽然存在很多优点，但是反应器的设计与推广上仍存在许多问题[5]，对于畜禽干清粪便来说，一是反应物料浓度过高，易出现中间产物例如乙酸等物质的积累，从而造成反应物料的酸化，进而阻止厌氧反应的进行；二是反应基质的不均匀性[6]，使得厌氧发酵反应体系不易控制，从而造成连续运行的不稳定性；三是搅拌阻力大，畜禽干清粪便厌氧发酵时的物料浓度高，搅拌阻力大，物料的搅拌混合困难[7]。

国内外学者近年来针对畜禽粪便厌氧发酵反应器进行了深入研究，赵兰兰[8]等在传统CSTR反应器的基础上对外形进行卵型结构设计，提高厌氧发酵混合效果进而提升厌氧发酵产气效率，但整体物料浓度基本维持在8%～10%，对于干清粪便物料浓度来说仍需进行稀释，进而增加反应器负荷。杨开宇[9]等在畜禽粪污厌氧发酵产甲烷反应器选型中也认为在畜禽粪便湿法发酵中，应选择CSTR式厌氧反应器进行，但随着技术的不断发展，考虑节能减排的因素，传统的湿法发酵技术在处理干清粪便中已出现很大的局限性，而国外一些研究学者则主要针对干式厌氧反应器进行了大量研究，其中最常见的3种典型连续式单相干发酵反应器为Dranco、Kompogas、Valorga。其中比利时的Dranco[10]干法反应器多处理经分类后的有机废弃物，其最大特点就是可将发酵剩余物回流至反应器内二次发酵，但由于未设置机械搅拌装置，难以解决干清粪便厌氧发酵时的发料问题。此外，BEKON[11]工艺与BIOFerm[12]工艺都是采用高温干式厌氧发酵技术，在发酵过程中并未设置搅拌装置，并且无好氧堆肥阶段，但针对干清粪便的处理来说，高温厌氧发酵会加速干清粪便的水解酸化进程，无搅拌装置，易造成畜禽粪便的酸化及发料，从而抑制厌氧发酵的进行。

尽管国内外许多学者对厌氧发酵反应器都开展了大量的研究，但针对干清粪便处理的半干式厌氧发酵反应器的研究尚存在空白。因此本文设计了1种适用于处理干清粪便的300 m3的半干式连续厌氧反应器，并对其进行结构参数的计算设计及强度分析，模拟仿真物料在罐体内的运动情况，进而为半干式厌氧发酵方式提供高效有利的厌氧发酵反应器。

1 半干式连续厌氧反应器总体设计与工作原理

经过设计的半干式厌氧反应器其结构如图1所示，厌氧反应器主要由发酵罐体、增温换热盘管、搅拌装置、进出料口及出料绞龙、搅拌电机、保温层组成。当厌氧反应器正常工作时，畜禽干清粪便由进料口匀速进入厌氧发酵罐体中，当装料浓度达到总发酵罐体体积的80%时停止进料，并且增温换热盘管进行工作，在太阳能辅助增温条件下对盘管内的热水进行循环加热，使其内部物料温度维持在(35±2)℃范围内，之后每隔4 h搅拌电机开始工作，采用星型卧式布置的搅拌轴在工作时能将物料实现沿搅拌轴周向运动与沿罐体周向的运动，从而使物料与底物充分接触，进而提高厌氧发酵产气效率，当完成厌氧发酵进程后，物料将在底部绞龙与出料口的共同作用下完成出料。

1.进料口；2.发酵罐体；3.出料口；4.搅拌电机；5.出料绞龙；6.搅拌装置；7.换热盘管；8.罐体顶盖

2 半干式厌氧反应器关键部件设计及有限元分析

2.1 发酵罐体结构参数设计及有限元分析

2.1.1 发酵罐体结构参数设计

通过对比国内外畜禽粪便厌氧发酵反应器及物料基本特性，最终确定反应器的发酵罐体设计成圆柱形结构，外形设计为矮胖型，既能节约构造材料，还可以解决垂直方向上传质不均的问题。

本次研究以总体积300 m3为例设计反应器，考虑厌氧发酵罐体顶部需要留有一定的空间用于排气与防止溢流，故而填料量为总体积的80%，所以罐体有效容积为240 m3，如图2所示。本厌氧发酵罐采用卧式搅拌装置星型结构布置，通过搅拌使罐内料液围绕圆柱罐体的中心轴形成转动，则发酵罐体内传质相对不充分区域沿垂直方向的投影为图3中的阴影部分。

图2 发酵罐体结构示意图(透视)

图3 厌氧发酵罐二维简图

需满足传质相对不充分区域投影面积占总面积比≤10%，选取占比为8%，则最终可确定，发酵罐体的直径D=11.04 m，料液高度H料=2.51 m。在保持厌氧发酵罐体直径不变的情况下最终确定厌氧发酵罐高度H=3.13 m，满足发酵罐体要求。厌氧发酵罐体材料选择Q235钢材，根据圆柱型罐体的相关最小壁厚计算公式同时考虑到钢板的负偏差与腐蚀余量，最终选择厌氧发酵罐罐体材料最小壁厚为8 mm的Q235钢板。

2.1.2 发酵罐体的有限元分析

厌氧反应中，发酵罐体为厌氧发酵设备中的核心功能部件，利用三维设计软件SolidWorks中创建了发酵罐体的三维图，由于在实际工作过程中罐体内部设置正负压保护器，因此发酵料液只会对发酵罐体底部与侧壁产生压力，故而在分析时简化发酵罐体顶盖，采用无顶盖时三维图进行力学分析，赋予发酵罐体模型材料参数。根据厌氧发酵罐内在实际工作情况下装料80%时，罐体底部与侧壁焊缝处受到物料产生的重力载荷与压力载荷，在此种情况下，罐体底部受到压力为固定约束，侧壁焊缝处受到压力为均布载荷，方向垂直于侧壁。

由图4和图5的应力图与总变形量云图中可以看到，厌氧发酵罐体最大应力在罐体的底部与侧壁的焊接处，在发酵罐体正常工作下会产生垂直侧壁向外的拉力，因此受到最大应力为0.6007 MPa，最大应力小于材料Q235的屈服强度185 MPa，符合设计要求。

图4 厌氧发酵罐体有限元应力云图

图5 厌氧发酵罐体有限元总变形量云图

2.2 搅拌装置的结构参数设计及有限元分析

2.2.1 搅拌装置的结构参数设计

实现物料的搅拌需要输入动力，即为搅拌功率。由于搅拌功率受叶翅形状，物料粘度等因素的影响，故而，在进行搅拌装置的设计时需计算搅拌功率。根据公式(1)可计算出搅拌功率[13]，其中，搅拌功率受功率准数、物料密度、搅拌速度与搅拌直径等因素有关，故而需计算功率准数。

P=NPρN3d5

(1)

根据设计要求，本次设计所选择的叶翅结构形状近似于MIG式搅拌叶翅，在计算功率准数时近似于多层平直桨叶，利用永田进治的搅拌功率计算式[14]，得:

(2)

式中：μ为物料粘度；b为桨叶宽度；d为搅拌直径；D为搅拌槽内径。

以干清牛粪为例，利用NDJ-5S数显粘度计测得物料粘度为19668 mpa·s，雷诺数Re=60.86，属于较高粘度，所设计的桨叶类型可近似为4层平直桨叶，其搅拌介质粘度范围通常为<50000 mPa·s[15]。由于本论文所设计反应器搅拌类型为卧式搅拌，故而，搅拌槽内径可近似选用料液高度，即D=2.51 m，根据规定，搅拌直径与搅拌槽内径的比值d/D在0.5～0.98，为增加搅拌解除面积，提高搅拌效率，故而将比值选为0.98，即搅拌直径d=2.4598 m。经过公式(2)的计算可得功率准数Np=5.79[15]。由于，干清粪便的粘度过高，故而设计转速条件时选取10 r·min-1。经过公式(1)的计算，单个搅拌轴功率为2.87 kW。

为实现搅拌轴的受力平衡，本论文在设计时将叶翅以螺旋形式分布，并且设有螺带，可以使物料沿轴向运动。故而，螺带在搅拌时消耗功率需进行计算，可将其近似于螺杆式搅拌器，其中，功率计算公式同(1)，但功率准数的计算略有不同，根据前苏联OB B BacH ∏ C等人对螺杆式搅拌器的功率准数的研究提出计算式公式为:

(3)

式中：Re为雷诺数；h为螺带高度；dj为螺带直径。

本次所设计的螺带条数为Z=1，dj=1.36 m，由于反应器中心需设计一个内中柱作为星型卧式搅拌轴的支撑底座，故而螺带高度h=4.91 m。经过计算可得，螺带的功率准数Np=4.15，故而单条螺带的搅拌功率P=0.11 kW。故而，经过计算单个搅拌设备的总功率P=2.98 kW，取整为3 kW。通常电机的传动效率为0.91，则单个搅拌设备的3.3 kW，根据常规电机功率规格表，取电机功率为4 kW的电机，搅拌轴及搅拌叶翅示意图如图6和图7所示。

图6 搅拌轴结构示意图

图7 搅拌叶翅结构示意图

通过最终计算，可以求得抗扭截面系数9.285×105 mm3，对于空心轴内外径之比通常取0.5～0.7，可求得最小外径D=116.98 mm，对最小外径进行取整得D=184 mm，则内径d=128 mm。

2.2.2 搅拌轴与搅拌叶翅的有限元分析

利用SolidWorks软件，创建搅拌轴及搅拌叶翅的三维模型，如图8所示。

图8 搅拌装置三维模型

将图8所示的搅拌装置的三维模型导入SolidWorks Simulation模块内，赋予搅拌轴与搅拌叶翅的材料参数。约束和压力添加完毕后，利用SolidWorks Simulation的Direct Sparse求解器进行求解。由图9和图10可以看出，搅拌轴与搅拌叶翅的最大应力集中在搅拌轴与电机连接端，最大应力为34.76 MPa，小于Q235钢材的屈服强度185 MPa，符合设计要求。

图10 搅拌装置有限元总变形量云图

3 半干式厌氧反应器流场数值模拟

3.1 workbench中模型的建立与网格划分

模型形状与实际形状基本相同，计算域选取搅拌区域与反应器其他区域，网格划分如图11和图12所示，546216个网格数，108150个节点数。

图11 整体网格划分图

图12 搅拌装置区域网格划分图

3.2 边界条件的设置及初始参数确定

本次选取求解器湍流模型，将流体的计算区域与边界进行初始化，初始速度为0，残差的收敛曲线设置为0.001，进行迭代计算，直至结果收敛。对搅拌装置的转速、物料的粘度等进行设置，不考虑反应的进行，只考虑混合流动情况。设置转速为10 r·min-1，转动方向沿搅拌轴顺时针，设置负方向重力加速度g=9.8 m·s-2，物料粘度设定为19.668 Pa·s，密度1187 kg·m-3。

3.3 结果及分析

当转速为10 r·min-1时、粘度为19.668 Pa·s时，反应器内物料的速度矢量图及运动轨迹图如图13和图14所示。沿搅拌轴轴向方向可以看出，搅拌轴附近的物料速度基本一致，可稳定在0.0104 m·s-1，而在搅拌叶片外径处速度达到最大，可达0.0174 m·s-1。由于物料混合的多样性与不同物料具有的复杂流变特性，在进行搅拌时，由于物料的粘度较高，在粘滞力的影响下会出现部分的层流状态。此外，从图中可以看出，物料以角速度旋转方向做沿搅拌轴的周向运动，并且速度大小会沿着搅拌直径方向呈线性增大，而在模拟过程中考虑重力场的影响，故而在底部出现转速大的现象；此外在搅拌轴作用之外的其他区域，由于物料存在的粘滞力，可使物料随着惯性出现沿罐体的周向运动，经过模拟后的流场速度分布基本与设计目标相符合，能够实现物料的全混合无死角搅拌。

图13 反应器物料运动速度矢量图

图14 反应器物料运动轨迹图

4 中试试验验证

4.1 工程概述

为验证该反应器针对畜禽粪便连续运行实际情况，吉林省农业机械研究院于乾安县牧羊厂建立300 m3沼气示范工程，该工程依托于吉林省科技发展计划项目“畜禽粪便高浓度沼气发酵及沼渣基质化利用技术装备研究”，总工程占地1300 m2。工程运行处理畜禽干清粪便产生天然气可供示范园区生产生活，沼渣渣液可固态、液态有机肥还田使用。实际工程如图15所示。

图15 300 m3沼气示范工程

4.2 实际工程运行情况

本次工程运行原料以干清牛粪作为主要发酵原料，物料浓度可达16.4%，按接种比例30%进行投料，装料容积为总体积的80%，装料完毕后，开启循环换热系统与搅拌装置，进行10 min低转速搅拌使发酵罐体内部物料均匀匀质，之后循环换热系统开启，当物料温度高于或低于(35±2)℃时，自动开启或停止循环换热以维持厌氧发酵罐体温度；保持每4 h 1次，转速为10 r·min-1，时长为2 min的间隙搅拌强度，当厌氧发酵反应进入产气高峰时，开始连续进出料，并每天记录产气情况与设备的运行检查。

本次半干式厌氧发酵反应器连续运行试验产气结果如图16所示，本次干清牛粪连续厌氧发酵试验，截取进入连续运行后的15天连续日产量，在厌氧发酵反应进入第7天开始进行连续进出料，在连续运行期间，平均日产气量可达410.6 m3，并且平均容积产气率可达1.37 m3·m-3d-1，整个连续厌氧发酵期内，干清牛粪能够正常高效完整厌氧发酵，并且维持较高的产气能力，可充分证明本次论文所研究的半干式厌氧反应器的工作性能，并值得推广应用。

图16 半干式连续厌氧发酵日产气量变化

5 讨论

5.1 对前人研究结果有所突破、创新的阐释

(1)针对畜禽粪便厌氧发酵，目前我国常见的厌氧反应器为CSTR型，本次论文设计的半干式厌氧发酵反应器与CSTR型反应器相比结构与外型有很大差异，首先，常规的CSTR型厌氧反应器罐体采用细高型立式圆柱结构[16]，而本次论文所设计的半干式厌氧发酵反应器采用矮胖型立式圆柱罐体结构，在合理范围内较小的高径比能减少物料沿垂直方向上的分层，并且在大型沼气工程中还可节省板材厚度，减少开支；其次，常规CSTR型反应器的机械搅拌为单一搅拌的小直径搅拌叶翅方式，需利用低浓度物料自身的水力特点在机械搅拌参与时实现物料自下而上的运动[17]，但对于畜禽干清粪便来说属高浓度物料，其粘度过高，CSTR型反应器的搅拌无法实现预期效果，而本文设计的半干式厌氧发酵反应器采用星型布置的卧式机械搅拌，不受物料粘度影响，能实现均匀传热传质，搅拌效果较佳。因此，本文设计的半干式厌氧反应器可针对干清粪便进行高效处理。

(2)本文所设计的半干式厌氧发酵反应器在中试试验中运行较为稳定，因半干式厌氧发酵反应器可满足高浓度物料的厌氧发酵，并且厌氧发酵环境适宜，接种物比例合适，因此，反应周期较短，效率较高，较传统的CSTR型厌氧发酵反应器的水利滞留期30 d相比，半干式厌氧发酵反应器的水利滞留期可缩短至20 d，并且在连续厌氧发酵过程中，畜禽粪便的降解率可达56.9%，连续厌氧发酵平均容积产气率也可达1.37 m3·m-3d-1，较CSTR型厌氧发酵反应器的平均容积产气率0.8 m3·m-3d-1有较高的提升[18]。

5.2 本研究需要说明的问题

本次论文所设计的半干式厌氧发酵反应器能高效处理浓度在12%～25%的农业废弃物，不仅仅只适用于干清粪便的厌氧发酵，对于秸秆类原料，合理的搅拌方式可以避免物料的结壳及分层，进而达到高效运行效果。此外，由于设备可进行半干式厌氧发酵，发酵浓度较高，能避免经过厌氧发酵后产生大量沼液造成沼液难以消纳的问题，因此本设备同时具有很好的推广应用价值。在今后的研究中将会继续加大对反应器运行时间的监测，从而确保反应器的稳定运行，并且还会增加对玉米秸秆连续厌氧发酵的研究，用以验证反应器的处理能力。

6 结论

本文设计了1种针对畜禽干清粪便连续厌氧发酵的反应器，为畜禽干清粪便的无害化处理与资源化利用提供了1种新的处理设备及厌氧发酵工艺，并基于300 m3沼气工程进行了实际运行试验，最终得到以下结论。

(1)利用三维建模软件SolidWorks对半干式厌氧反应器进行总体结构设计，基于发酵物料的理化性质、半干式厌氧发酵工艺及基本工作原理，对总体及关键部件进行设计计算机有限元分析，均符合实际生产需要与设计目标，值得推广应用。

(2)通过流体数值仿真分析对半干式厌氧反应器进行流体力学仿真分析，进而模拟实际运行结果，在转速为10 r·min-1的条件下，物料能够实现在罐体内的全混合式搅拌，使物料沿搅拌轴周向旋转的同时，实现物料沿罐体的周向运动，使得物料与底物的接触更加充分，完全可避免在实际厌氧反应进程中发酵原料发料、结壳、分层等问题。

(3)通过300 m3沼气工程实际运行情况来看，半干式厌氧反应器整机性能良好，能维持稳定的厌氧发酵环境，并且对于干清粪便的处理效果十分显著，在连续运行期间，平均日产气量能稳定在410.6 m3左右，平均容积产气率可达1.37 m3·m-3d-1，工程运行稳定，值得推广应用。