2. 反应器类型与放大对象建模

放大计算的第一步不是套公式,而是把“被放大的对象”描述清楚。同一个体积数字,在搅拌槽、气升式反应器和鼓泡塔中对应的混合机制、剪切来源和传质路径完全不同。如果建模阶段把反应器抽象成一个只有体积和通气量的黑箱,后续无论用单位体积功率还是kLa关联式,都会因为几何与操作变量定义不一致而失真。本节的任务是:识别三类典型生物反应器的结构特征,明确放大计算中必须固定的几何参数与操作变量,并据此整理出一份可直接用于计算的输入输出清单。

三类反应器的结构特征与放大敏感点

机械搅拌槽由罐体、搅拌轴、叶轮、挡板、分布器和换热构件组成。它的核心特征是能量由叶轮集中输入,因此局部剪切、循环流量和功率消耗都强烈依赖叶轮直径D、转速N、叶轮型式与层数。放大时最敏感的量是单位体积功率P/V与叶端线速度πND,两者往往不能同时保持。搅拌槽的几何相似通常以D/T(叶轮直径与罐径之比)、H/T(液位高度与罐径之比)、挡板数及宽度来约束。

气升式反应器没有机械搅拌,依靠上升区与下降区的表观气速差形成液体循环。其结构关键是导流筒(或隔板)把床层分为上升区和下降区,气液分离区决定循环回流的脱气程度。放大敏感点是循环时间、气含率与下降区是否出现气窜。几何参数包括导流筒直径与罐径之比、导流筒高度、上升区与下降区截面积比、底部间隙和顶部浸没深度。

鼓泡塔结构最简单:无内构件或仅有分布器与挡板,气体从底部经分布器进入,靠气泡浮力驱动液体宏观循环。其放大难点在于高径比增大后轴向返混加剧、气泡聚并行为随直径变化,导致kLa与气含率的经验关联式外推风险高。关键几何量为塔径T、静液高H、分布器开孔数与孔径、高径比H/T。

三类典型生物反应器的结构与放大敏感点 机械搅拌槽 H 叶轮 D 气体分布器 能量集中输入 · 剪切强 敏感量:P/V、πND 气升式反应器 导流筒 上升区 下降区 气液分离区 气差驱动循环 · 无机械剪切 敏感量:循环时间、气含率 鼓泡塔 分布器开孔 结构简单 · 轴向返混大 敏感量:H/T、kLa外推

几何参数与操作变量的定义必须“先约定、后计算”

放大计算中最常见的错误,是不同来源的关联式使用了不同的参数定义。例如“通气量”可以指体积流量Q(m³/s),也可以指空塔表观气速vs=Q/A(m/s),还可以指vvm(每分钟通入体积与液体积之比)。三者数值差异巨大,代入关联式前必须统一。类似地,P/V中的V应明确是装液体积还是罐体公称体积;kLa的单位是s⁻¹还是h⁻¹;叶端线速度用πND还是πN(D+2δ)。建议在建模阶段就建立一张参数定义表,把符号、单位、参照基准固定下来,全文与图示保持一致。

几何参数方面,搅拌槽至少需要:罐径T、装液高H、叶轮直径D、叶轮层数与间距、挡板数与宽度、分布器型式。气升式需要:导流筒直径与高度、上升/下降区截面积、底部间隙、顶部浸没深度。鼓泡塔需要:塔径T、静液高H、分布器开孔数与孔径、高径比。操作变量方面,三类反应器共有的是通气流量与表观气速、温度、pH、进料速率;搅拌槽额外需要转速N与功率P;气升式与鼓泡塔则以表观气速和循环特性为主要操作量。

放大计算的输入输出清单

把上述内容整理成清单,是本节最直接的交付物。清单的作用不是罗列名词,而是保证每次计算时输入完整、输出可追溯。下表给出一个可填写的示例框架,其中数值为教学示例假设,实际项目需按真实设备与工艺填写。

类别 参数 符号/单位 示例值(假设) 说明
几何输入罐径T / m0.30小试罐内径
装液高H / m0.60H/T=2.0
叶轮直径D / m0.10D/T=0.33
挡板数n / 个4宽度T/10
操作输入转速N / r·min⁻¹300小试条件
通气量Q / vvm1.0需换算为m³/s
温度T / ℃30恒温控制
装液体积V / L42由T、H计算
中间量单位体积功率P/V / W·m⁻³待算由功率准数关联
叶端线速度πND / m·s⁻¹待算剪切判据
表观气速vs / m·s⁻¹待算Q除以截面积
输出放大后转速N₂ / r·min⁻¹待算按P/V或πND准则
放大后通气量Q₂ / m³·s⁻¹待算按vs或vvm准则
kLa预测值kLa / s⁻¹待算关联式外推需谨慎
清单中的“中间量”一栏非常关键。放大准则往往不是直接给出目标转速,而是先要求P/V或πND相等,再反解N₂。把中间量显式列出,可以避免在多个关联式之间跳步时丢失单位与基准。

从清单到计算:一个可复现的示例

以搅拌槽为例,假设小试罐T=0.30 m、D=0.10 m、装液体积V=42 L、N=300 r/min。若采用Rushton涡轮,功率准数Np约为5.5(湍流区),则P=Np·ρ·N³·D⁵。取ρ=1000 kg/m³,N=5 s⁻¹,计算得P≈5.5×1000×125×1×10⁻⁵≈6.9 W,P/V≈164 W/m³。叶端线速度πND≈1.57 m/s。若放大到V₂=2000 L并保持P/V不变,则需先由几何相似确定D₂与V₂的关系,再反解N₂;此时πND₂通常会显著上升,必须回头检查剪切是否超过菌体耐受上限。这个“算完再校验”的环节,正是输入输出清单要预留输出项的原因。

对于气升式与鼓泡塔,清单中的几何项要换成导流筒或分布器参数,操作项以表观气速为主,输出项则重点关注循环时间与kLa。三类反应器共用一张清单框架,但字段按类型启用,这样在后续章节做具体放大计算时,可以直接引用同一套定义,避免参数口径漂移。

不要用“几何相似”一句话代替几何参数记录。放大中真正被保持的往往只是部分相似(如H/T、D/T),而分布器开孔率、导流筒间隙等细节在小试与中试间经常被改动。清单中应记录这些“非相似项”,否则放大偏差无法归因。

完成本节建模后,读者应能对任意一台候选反应器回答三个问题:它属于哪一类结构、它的放大敏感量是什么、它的几何与操作参数是否已经按统一口径记录完毕。只有这三个问题都有明确答案,后续的功率、传质与混合计算才有可靠起点。