2. 反应器类型与放大对象建模
放大计算的第一步不是套公式,而是把“被放大的对象”描述清楚。同一个体积数字,在搅拌槽、气升式反应器和鼓泡塔中对应的混合机制、剪切来源和传质路径完全不同。如果建模阶段把反应器抽象成一个只有体积和通气量的黑箱,后续无论用单位体积功率还是kLa关联式,都会因为几何与操作变量定义不一致而失真。本节的任务是:识别三类典型生物反应器的结构特征,明确放大计算中必须固定的几何参数与操作变量,并据此整理出一份可直接用于计算的输入输出清单。
三类反应器的结构特征与放大敏感点
机械搅拌槽由罐体、搅拌轴、叶轮、挡板、分布器和换热构件组成。它的核心特征是能量由叶轮集中输入,因此局部剪切、循环流量和功率消耗都强烈依赖叶轮直径D、转速N、叶轮型式与层数。放大时最敏感的量是单位体积功率P/V与叶端线速度πND,两者往往不能同时保持。搅拌槽的几何相似通常以D/T(叶轮直径与罐径之比)、H/T(液位高度与罐径之比)、挡板数及宽度来约束。
气升式反应器没有机械搅拌,依靠上升区与下降区的表观气速差形成液体循环。其结构关键是导流筒(或隔板)把床层分为上升区和下降区,气液分离区决定循环回流的脱气程度。放大敏感点是循环时间、气含率与下降区是否出现气窜。几何参数包括导流筒直径与罐径之比、导流筒高度、上升区与下降区截面积比、底部间隙和顶部浸没深度。
鼓泡塔结构最简单:无内构件或仅有分布器与挡板,气体从底部经分布器进入,靠气泡浮力驱动液体宏观循环。其放大难点在于高径比增大后轴向返混加剧、气泡聚并行为随直径变化,导致kLa与气含率的经验关联式外推风险高。关键几何量为塔径T、静液高H、分布器开孔数与孔径、高径比H/T。
几何参数与操作变量的定义必须“先约定、后计算”
放大计算中最常见的错误,是不同来源的关联式使用了不同的参数定义。例如“通气量”可以指体积流量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 / m | 0.30 | 小试罐内径 |
| 装液高 | H / m | 0.60 | H/T=2.0 | |
| 叶轮直径 | D / m | 0.10 | D/T=0.33 | |
| 挡板数 | n / 个 | 4 | 宽度T/10 | |
| 操作输入 | 转速 | N / r·min⁻¹ | 300 | 小试条件 |
| 通气量 | Q / vvm | 1.0 | 需换算为m³/s | |
| 温度 | T / ℃ | 30 | 恒温控制 | |
| 装液体积 | V / L | 42 | 由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⁻¹ | 待算 | 关联式外推需谨慎 |
从清单到计算:一个可复现的示例
以搅拌槽为例,假设小试罐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。三类反应器共用一张清单框架,但字段按类型启用,这样在后续章节做具体放大计算时,可以直接引用同一套定义,避免参数口径漂移。
完成本节建模后,读者应能对任意一台候选反应器回答三个问题:它属于哪一类结构、它的放大敏感量是什么、它的几何与操作参数是否已经按统一口径记录完毕。只有这三个问题都有明确答案,后续的功率、传质与混合计算才有可靠起点。